pengchy 还是有点问题: 关于score的计算,根据《实用现代统计分析方法与SPSS应用——张文璋》一书的第250页,某主成分的得分计算公式应为 <br /> x=scale(x); <br /> cor=cor(x); <br /> eig=eigen(cor); <br /> scores <- x %*% eig$vec<br /> 也可能是我理解的有问题,多讨论。
smalltail 回复 第1楼 的 dengyishuo:我还发现用你的code算出来的其实和princomp的意义是一样的,不同的是数值的大小,可能princomp()的结果是你的code输出结果绝对值的几倍...成比例的,但是解释的意义我想是一样的..
dengyishuo 回复 第16楼 的 smalltail:主贴中的code只是为了使得princomp()的输出结果SPSS和sas的输出结果相一致。如果说主成分分析的原理和意义的话,两者的确是一致的。你说的成比例应该等于sqrt(value)。
zxcv_bnm R做主成分,用eigen或者princomp都应该注意,最大问题是,与SPSS在载荷或者特征向量上,相差一个负号哦(这才是实质)。很多地方没发觉,包括《统计建模与R书》,我一个博客,可以解决大家的疑惑。这里我有点怀疑R里面的特征向量实现的算法,其他人没发觉? 再提醒,相差一个符号,会导致综合评价排序成相反。