24小时热门版块排行榜

>论坛更新日志 (1727)
>虫友互识 (165)
>考博 (31)
>导师招生 (26)
>硕博家园 (19)
>考研 (18)
>论文道贺祈福 (17)
>基金申请 (17)
>博后之家 (16)
>文献求助 (15)
>论文投稿 (14)
>教师之家 (13)
>找工作 (13)
>休闲灌水 (12)
>公派出国 (11)
>医学 (4)

返回列表

jiaxin11

铜虫 (正式写手)

应助: 5 (幼儿园)
金币: 11.3
散金: 35
红花: 2
帖子: 433
在线: 100.7小时
虫号: 2106229
注册: 2012-11-04
性别: MM
专业: 信息理论与信息系统

[求助] 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)

回复此楼