推 bugle:你所謂的 c1'c1 不就是 sum(x[,1]^2) ? 10/09 21:31
→ bcs:謝了，我知道了^^" 10/09 22:05
x<-matrix(1:15,nrow=5,3)
ci:ith column vector of x
想把x的 c1'c1+c2'c2+c3'c3
理論上應該等於
tcrossprod(x) = x'x
若一步步算
tcrossprod(x[,1])+
tcrossprod(x[,2])+
tcrossprod(x[,3])
但是若matrix有n個column時，是否可以利用apply的方式
得到上式。
apply(x,2,function(x) sum( tcrossprod(x) )
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.86.139