24小时热门版块排行榜    

查看: 3293  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

一笑35

铁虫 (初入文坛)

[交流] 【求助/交流】虫虫们,有没有做PCA主成分分析的啊? 已有12人参与

各位虫虫们,有没有做PCA主成分分析的啊,一般是用什么软件呢,有可以免费下载的网址吗?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

努力学习,充实自己!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

reasonspare

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
silicare(金币+2):专家级~ 2010-08-25 12:28:28
引用回帖:
Originally posted by 一笑35 at 2010-08-25 09:12:36:
各位虫虫们,有没有做PCA主成分分析的啊,一般是用什么软件呢,有可以免费下载的网址吗?

最经典的是用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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

yanruoke

银虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请问分析对象具体是?
戒嗔怒以养肝气,省言语以养神气,多读书以养质气,顺时令以养元气,不拘节以养大气,观天变以养灵气,莫强求规于运气。
2楼2010-08-25 09:18:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一笑35

铁虫 (初入文坛)

我是用来做DGGE的PCA主成分分析的,可以从较官方网上免费下载吗?
努力学习,充实自己!
4楼2010-08-25 10:54:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sqhnsd

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by 一笑35 at 2010-08-25 09:12:36:
各位虫虫们,有没有做PCA主成分分析的啊,一般是用什么软件呢,有可以免费下载的网址吗?

大名鼎鼎的R软件业可以做,不过需要自己去学习,关键R是免费的
5楼2010-08-25 10:55:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 7/350 2026-08-25 03:24 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见