| 查看: 840 | 回复: 1 | ||
[求助]
principle analysis components
|
|
各位， 你们好！ 因为我有很多变量，想做Multi-linear regression， 但是变量之间存在着multicollinearity, 所以之后 我用principle analysis components 分析出了主要的components, 然后用Multi-linear regression analyses 得到了的以components 为因子的线性模型。然后我用R算出了原来那些变量的系数，我不知道合不合理。 请各位大神帮我看看。非常感谢 ！ pcal=princomp(sub_agri211, scores=TRUE,cor=TRUE) summary(pcal) loadings(pcal) plot(pcal) pre=predict(pcal) sub_agri21$Z1=pre[,1] sub_agri21$Z2=pre[,2] lm_sub_agri21=lm(sta_day~Z1+Z2,data=sub_agri21) summary(lm_sub_agri21) Anova(lm_sub_agri21) lm_sub_agri211=update(lm_sub_agri21,.~.-Z1) summary(lm_sub_agri211) Anova(lm_sub_agri211) beta=coef(lm_sub_agri211) A=loadings(pcal) x.bar=pcal$center x.bar x.sd=pcal$scale x.sd coef=beta[2]*A[,2]/x.sd beta0=beta[1]-sum(x.bar*coef) c(beta0,coef) |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
churchery
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2475.8
- 散金: 250
- 红花: 2
- 帖子: 225
- 在线: 202.2小时
- 虫号: 2866067
- 注册: 2013-12-12
- 性别: GG
- 专业: 控制理论与方法
2楼2016-06-06 20:39:39
20