我已获得数据集Usedata
的协方差矩阵。我还得到了特征值的对角矩阵和R的输出,并将其保存为csv。
我无法在csv文件的顶部看到每个变量名称,而是只显示“X1”,“X2”等。我想在每列的顶部看到变量名,这样我就可以看到哪个变量具有最大的特征值。
我的代码:
Usedata <- structure(list(X1 = c(1, 0, 0, 0.244012404, 0, 0, 6, 0, 0, 0,
0, 0), X2 = c(52.72564729, 2, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0),
X3 = c(0, 0, 3, 0, 0.142522511, 0, 0, 0, 8, 0, 0, 0), X4 = c(0,
0.341103073, 0, 4, 0, 0, 0, 0, 0, 9, 0, 0), X5 = c(0, 0,
0, 0, 5, 0.091644475, 0, 0, 0, 0, 10, 0)), .Names = c("X1",
"X2", "X3", "X4", "X5"), class = "data.frame",
row.names = c(NA, -12L))
smallcov <- cov(Usedata)
lam <- eigen(smallcov)$values
LamM <-diag(lam)
diagresult <- data.frame(LamM)
write.csv(diagresult, file = "myoutput.csv")