2019-12-20

data<-read.csv("C:/Users/ZWPC/Desktop/pls.csv",header=F)
pls <- plsr(V1~.,data=data,30,validation="CV")#全波段最小二乘法
selectNcomp(pls,plot=T)
cor(data$V1,pls$fitted.values[,,10])^2
summary(pls)#查看
plot(data$V1,pls$fitted.values[,,10])
write.csv(data,"C:/Users/ZWPC/Desktop/pls1.csv",row.names=F)
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容