| 查看: 3413 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助/交流】虫虫们，有没有做PCA主成分分析的啊？ 已有12人参与
|
|||
| 各位虫虫们，有没有做PCA主成分分析的啊，一般是用什么软件呢，有可以免费下载的网址吗？ |
回复此楼» 猜你喜欢
生物质超级电容器
已经有0人回复
-
ALK5 高选择性小分子 LY-3200882：可放大合成路线与质控要点解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有183人回复
-
AM630（CAS 164178-33-0）：CB2 选择性配体的实验使用与供应要点
已经有0人回复
-
ADC标杆分子T-DM1全解析：从结构设计到实验室应用的技术要点
已经有0人回复
-
HILIC磁珠研究蛋白组学怎么做？
已经有0人回复
-
His标签蛋白纯化怎么才能快？
已经有0人回复
-
WRN变构抑制剂HRO761：从D1/D2界面锁定到MSI合成致死的工具化合物深度解析
已经有0人回复
-
KRAS G12V(ON) 选择性抑制剂 RMC-5127：从三元复合物机制到科研应用
已经有0人回复
-
SHP2 变构抑制剂 Sitneprotafib：从作用机制到科研试剂选型要点
已经有0人回复
-
PRMT5·MTA 抑制剂 MRTX1719 的科研应用
已经有0人回复
快乐的米粒
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1838.7
- 帖子: 164
- 在线: 48.9小时
- 虫号: 1956333
- 注册: 2012-08-26
- 性别: MM
- 专业: 女性生殖系统肿瘤
送鲜花一朵 
12楼2012-09-13 17:22:11
2楼2010-08-25 09:18:13
reasonspare
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.894
- 金币: 2137.1
- 红花: 11
- 帖子: 1909
- 在线: 72.5小时
- 虫号: 180928
- 注册: 2006-02-10
- 专业: All in one
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
silicare(金币+2):专家级~ 2010-08-25 12:28:28
小木虫(金币+0.5):给个红包，谢谢回帖交流
silicare(金币+2):专家级~ 2010-08-25 12:28:28
|
最经典的是用SAS得princomp过程。另来SAS 得的主因子分析 Factor 因不错。 此外SPSS.也不错。其他的诸如 Metlab，也有现成的命令 princomp - Principal component analysis on data Syntax [COEFF,SCORE] = princomp(X) [COEFF,SCORE,latent] = princomp(X) [COEFF,SCORE,latent,tsquare] = princomp(X) [...] = princomp(X,'econ') Description COEFF = princomp(X) performs principal components analysis on the n-by-p data matrix X, and returns the principal component coefficients, also known as loadings. Rows of X correspond to observations, columns to variables. COEFF is a p-by-p matrix, each column containing coefficients for one principal component. The columns are in order of decreasing component variance. princomp centers X by subtracting off column means, but does not rescale the columns of X. To perform principal components analysis with standardized variables, that is, based on correlations, use princomp(zscore(X)). To perform principal components analysis directly on a covariance or correlation matrix, use pcacov. [COEFF,SCORE] = princomp(X) returns SCORE, the principal component scores; that is, the representation of X in the principal component space. Rows of SCORE correspond to observations, columns to components. [COEFF,SCORE,latent] = princomp(X) returns latent, a vector containing the eigenvalues of the covariance matrix of X. [COEFF,SCORE,latent,tsquare] = princomp(X) returns tsquare, which contains Hotelling's T2 statistic for each data point. The scores are the data formed by transforming the original data into the space of the principal components. The values of the vector latent are the variance of the columns of SCORE. Hotelling's T2 is a measure of the multivariate distance of each observation from the center of the data set. When n <= p, SCORE(:,n:p) and latent(n:p) are necessarily zero, and the columns of COEFF(:,n:p) define directions that are orthogonal to X. [...] = princomp(X,'econ') returns only the elements of latent that are not necessarily zero, and the corresponding columns of COEFF and SCORE, that is, when n <= p, only the first n-1. This can be significantly faster when p is much larger than n. Examples Compute principal components for the ingredients data in the Hald data set, and the variance accounted for by each component. load hald; [pc,score,latent,tsquare] = princomp(ingredients); pc,latent pc = 0.0678 -0.6460 0.5673 -0.5062 0.6785 -0.0200 -0.5440 -0.4933 -0.0290 0.7553 0.4036 -0.5156 -0.7309 -0.1085 -0.4684 -0.4844 latent = 517.7969 67.4964 12.4054 0.2372 The following command and plot show that two components account for 98% of the variance: cumsum(latent)./sum(latent) ans = 0.86597 0.97886 0.9996 1 biplot(pc(:,1:2),'Scores',score(:,1:2),'VarLabels',... {'X1' 'X2' 'X3' 'X4'})  |
» 本帖已获得的红花（最新10朵）
3楼2010-08-25 09:59:53
4楼2010-08-25 10:54:10