样本矩阵的均值协方差

13046

收藏 2012-12-14

如题，如何求样本矩阵的均值协方差呢?
自变量含有X1，X2，X3，X4，X5，X6观察值比方说50个
就组成50*6的样本矩阵
估计样本矩阵的一阶矩，二阶矩
是不是一阶矩就是列平均，用colMeans()？
那二阶怎么做呢？应该好像是个协方差矩阵的形式

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

子衿1219

2012-12-14 16:12:16

补充问下矩阵A的1/2次方如何求啊？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

qoiqpwqr

2012-12-14 23:32:16

子衿1219 发表于 2012-12-14 16:12
补充问下矩阵A的1/2次方如何求啊？

eigen-decomposition

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

qoiqpwqr

2012-12-14 23:35:00

子衿1219 发表于 2012-12-14 16:12
补充问下矩阵A的1/2次方如何求啊？

> A <- matrix(c(2,0.6,0.6,2), 2, 2)
> b <- eigen(A)
> B <- b$vectors %*% diag(sqrt(b$values)) %*% t(b$vectors)
> B %*% B
   [,1] [,2]
[1,]  2.0  0.6
[2,]  0.6  2.0

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

子衿1219

2012-12-16 15:23:00

qoiqpwqr 发表于 2012-12-14 23:35
> A  b  B  B %*% B
   [,1] [,2]
[1,]  2.0  0.6

谢啦。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

子衿1219

2012-12-16 15:24:40

qoiqpwqr 发表于 2012-12-14 23:35
> A  b  B  B %*% B
   [,1] [,2]
[1,]  2.0  0.6

我还想问下，在R中，如果使用了一个文件包，必须加library（）吗？那这个命令加在哪里呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群