各路大神，帮我看看这个程序计算残差部分哪里错了，怎么出不来结果

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

#### 主成分分析
D=scale(DA[1:50])
R=cov(D)
eigen(R)$values
pr=princomp(D)
summary(pr, loadings=F)
# 预测样本主成分, 并作主成分分析
pre<-predict(pr)
DA$z1<-pre[,1]
DA$z2<-pre[,2]
DA$z3<-pre[,3]
lm.z<-lm(DA$Y~DA$z1+DA$z2+DA$z3)
summary(lm.z)

#### 作变换, 得到原坐标下的关系表达式
beta<-coef(lm.z); A<-loadings(pr)
x.bar<-colMeans(DA[,1:50]); x.sd<-sapply(DA[,1:50], sd)
coef<-(beta[2]*A[,1]+ beta[3]*A[,2]+ beta[4]*A[,3])/x.sd
beta0 <- beta[1]- sum(x.bar * coef)
c(beta0, coef)
#计算残差
x<-as.matrix(DA[,1:50])
coef=as.vector(coef)
y<-DA[,51]
r<- y-beta0-x%*%coef
Res<-t(r)
Time<-1:78
plot(Time,Res,type="l")
lines(c(-15,515),c(0,0))
Res%*%r

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：princomp Loadings Summary loading values 表达式 程序 样本

What is DA?