setwd("C:/Users/IBM/Desktop/研一课程/2.2回归分析/回归作业") #设定当前的工作目录
shuju=read.table("shuju.txt",header=T)
shuju #读取数据
#主成分回归,用函数pcr
library(pls) #pls包提供pcr函数,也提供偏最小二乘回归函数plsr,而函数mvr是二者的集成
attach(shuju) #pcr可以在scale=T的设置下对自变量进行标准化,但注意对因变量并未标准化
set.seed(1)
shuju.pcr<-pcr(y~x1+x2+x3+x4+x5+x6+x7,4,scale=T,validation="CV")#主成分回归
summary(shuju.pcr) #这里采用的最优量化指标为均方根误差RMSEP,要得到MSE,需要对RMSEP做平方
coef(shuju.pcr)
#还原到原始变量的模型系数
cx1<-coef(shuju.pcr)[1]/sd(shuju$x1)
cx2<-coef(shuju.pcr)[2]/sd(shuju$x2)
cx3<-coef(shuju.pcr)[3]/sd(shuju$x3)