我没有从R(版本2.15.2)中的嵌套for循环中收到正确的结果。
GCP <- matrix(nrow=143,ncol=4)
for (i in 1:13) {
for (j in 1:11) {
GCP[i*j,1] = 410*(j-1)
GCP[i*j,2] = 400*(i-1)
# print(GCP[i*j,1])
}
}
当我取消注释打印(GCP [i * j,1])时,我收到了我期望的内容:
[1] 0
[1] 410
[1] 820
[1] 1230
[1] 1640
[1] 2050
[1] 2460
[1] 2870
[1] 3280
[1] 3690
[1] 4100
[1] 0
[1] 410
[1] 820
[1] 1230
[1] 1640
[1] 2050
[1] 2460
[1] 2870
[1] 3280
[1] 3690
[1] 4100
[1] 0
...
但是当我在循环结束后输入GCP [,1]时:
[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 410 820
[16] 410 NA 410 NA 410 820 410 NA 410 1640 410 820 1230 NA 820
[31] NA 1230 820 NA 1640 820 NA NA 820 1230 NA 2050 NA 1230 1640
[46] NA NA 1230 2460 1640 NA 1230 NA 2050 1640 2460 NA NA NA 1640
[61] NA NA 2460 2870 1640 2050 NA NA NA 2460 NA 2050 NA NA NA
[76] NA 2460 2050 NA 2870 3280 NA NA 2460 NA NA NA 2870 NA 3280
[91] 2460 NA NA NA NA 2870 NA NA 3280 3690 NA NA NA 2870 NA
[106] NA NA 3280 NA 3690 NA NA NA NA NA NA 3280 NA NA 3690
[121] 4100 NA NA NA NA NA NA NA NA 3690 NA 4100 NA NA NA
[136] NA NA NA NA NA NA NA 4100
我该如何更正我的代码?
答案 0 :(得分:2)
问题是,当你这样做时:
GCP <- matrix(nrow = 143, ncol = 4)
您实际上正在使用(143,4)
对维度NA
进行初始化。然后,您只为row = i * j
构建值。而你只是打印这些价值观。也许这段代码可以帮助你清理一下:
GCP <- matrix(nrow = 143, ncol = 4)
for (i in 1:13) {
for (j in 1:11) {
GCP[i*j,1] = 410*(j-1)
GCP[i*j,2] = 400*(i-1)
cat("i =", i, "j =", j, "i*j =", i*j, "GCP[i*j] = ", GCP[i*j,1], "\n")
}
}
您会看到循环遍历i
和j
的所有值。但是,您的索引仅填充特定值= i*j
。因此,其余值仍为NA
。由于您只打印已设置的值,因此在打印时永远不会获得NA。
除非你解释你想要做什么,否则我无法提供解决方案。首先,您要设置一个ncol=4
的矩阵,但是您只填充前两列,并且正如@redmode指出的那样,您将重写一些值。
编辑:也许那么,这就是你要找的东西?
GCP <- matrix(nrow = 143, ncol = 4)
for (i in 1:13) {
for (j in 1:11) {
GCP[11*(i-1) + j, 1] = 410*(j-1)
GCP[11*(i-1) + j, 2] = 400*(i-1)
}
}
这将通过1到143之间的值。