diag（）を返すこの関数の説明はR

Question

でこれはなぜですか？ y <- min(dim(x))？

このコードを試してみたところ、c <- vector()となっていました（c <- x[i,i]がループしていました）。ループ内でc[i]を使用するとエラーが表示されます）。

x <- matrix(1:25, 5,5) 
> diag(x) 
> n <- function(x){ 
    y <- min(dim(x)) 
    c <- vector() 
    for(i in 1:y){ 
     c[i] <- x[i,i] 
    } 
    c 
} 
> n(x) 

Answer 1

y <- min(dim(x))は、マトリックスのより小さい寸法を決定するために使用される。したがって、yを使用すると、forループの反復回数は、行列のこの次元に等しくなるように固定されます。次に、[i,i]を使用して対角要素にアクセスします。

この例では、行列のサイズは5x5なので、両方の次元が同じです。したがって、yが5の場合、iを1から5までループします。今度は、[i,i]をマトリックスインデックスに使用しているので、[1,1]の後に[2,2]が得られます。本質的には、[i,i]は、行番号と列番号の両方が同じ値である場合に要素を選択します。

c <- vector()を削除すると、cは1つの変数として扱われ、したがってc[i]はエラーをスローします。また、cの値は上書きされ続けるので、ループを終了すると、最後の値が得られます。この場合、この値は25です。