样本矩阵的均值协方差

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：协方差 means 协方差矩阵 mean LME 矩阵 样本

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

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

eigen-decomposition

子衿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

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

谢啦。

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

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

子衿1219 发表于 2012-12-16 15:24
我还想问下，在R中，如果使用了一个文件包，必须加library（）吗？那这个命令加在哪里呢？

在你用这个包之前都可以。一般在程序的最开始把要用的所有包都library一下