0

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])
 }
}

print(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

コードを修正するにはどうすればよいですか?

4

1 に答える 1

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いますが、最初の 2 つの列のみを埋めており、@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 までの値を通過します。

于 2013-02-03T11:51:11.377 に答える