［R］行列のトレース（固有和）を求める

Rには行列のトレースを直接計算する関数は標準では搭載されていない。

sum関数とdiag関数を組み合わせると簡単に求めることができる。

> d <- 1:9
> mx <- matrix(d, 3, 3)
> mx
     [,1] [,2] [,3]
[1,]    1    4    7
[2,]    2    5    8
[3,]    3    6    9
> sum(diag(mx))
[1] 15
> d <- 1:16
> mx <- matrix(d, 4, 4)
> mx
     [,1] [,2] [,3] [,4]
[1,]    1    5    9   13
[2,]    2    6   10   14
[3,]    3    7   11   15
[4,]    4    8   12   16
> sum(diag(mx))
[1] 34