首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >R中的PCA与prcomp

R中的PCA与prcomp
EN

Stack Overflow用户
提问于 2021-12-09 08:47:31
回答 1查看 242关注 0票数 1

我试图用PCA计算一个指数。我使用‘mtcar’作为我的数据集。这是我的密码:

代码语言:javascript
复制
 pca_mt<- prcomp(mtcars[,c(1:7,10,11)], center = TRUE,scale. = TRUE)

我如何获得每个变量的权重,以求它们之和以计算索引?如果有人能解释一下pca_mt$rotationpca_mt$x是什么,我很感激。

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2021-12-09 09:20:57

也许这有助于您更好地理解pca_mt$rotationpca_mt$x是什么:

代码语言:javascript
复制
data <- mtcars[,c(1:7,10,11)]
pca_mt<- prcomp(data, center = F,scale. = F)
sum(data[1,] * pca_mt$rotation[,1]) == pca_mt$x[1,1]
[1] TRUE

请注意,这只是因为center=Fscale.=F才起作用。

为了进一步阅读,我建议https://builtin.com/data-science/step-step-explanation-principal-component-analysis首先了解PCA是什么。

因此,要了解PCA组件可以解释多少方差,可以使用

代码语言:javascript
复制
summary(pca_mt)
Importance of components:
                            PC1      PC2      PC3     PC4     PC5     PC6    PC7    PC8    PC9
Standard deviation     310.1169 40.88309 15.83774 2.12821 0.99854 0.74271 0.4412 0.2913 0.2168
Proportion of Variance   0.9803  0.01704  0.00256 0.00005 0.00001 0.00001 0.0000 0.0000 0.0000
Cumulative Proportion    0.9803  0.99738  0.99993 0.99998 0.99999 1.00000 1.0000 1.0000 1.0000

第一组分解释了数据的98,03%,下一组(98,03%-99,738%)等。因此,使用第一个组件,您可以解释数据中几乎所有的差异。

票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/70286925

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档