程序三：投影寻踪回归分析

##################################################################################################################ͶӰѰ×ٻعé·ÖÎöProjectionPursuitRegression###rm(list=ls(all=TRUE))#ɾ³ıRÈí¼şÖĞËùÓжÔÏósetwd('d:/qsardata')#ÉèÖù¤×÷Ŀ¼£¬½«¾­¹ıÔ¤´¦ÀíµÄÊı¾İ·ÅÔÚ¸ÃĿ¼Ï£¬ºóĞø¹¤×÷Ëù²úµÄÎļş¶¼½«±£´æÔÚ¸ÃĿ¼ÏÂgetwd()##########################################################train.data-read.csv(file=file.choose(),header=T)#¶ÁȡѵÁ·¼¯Êı¾İtraindata-train.data[,-1]#½«ÑµÁ·¼¯Êı¾İµÄ×Ô±äÁ¿ÓëÒò±äÁ¿ÌáÈ¡µ½ÎļştraindataÖĞdim(traindata);colnames(traindata)trainactivity-traindata$activitytest.data-read.csv(file=file.choose(),header=T)#¶ÁÈ¡²âÊÔ¼¯Êı¾İtestdata-test.data[,-1]#½«²âÊÔÊı¾İµÄ×Ô±äÁ¿ÓëÒò±äÁ¿ÌáÈ¡µ½ÎļştestdataÖĞtestactivity-testdata$activitydim(testdata);colnames(testdata)##############################################################################################################################################################################require(graphics)ppr.model-ppr(activity~.,data=traindata,nterms=4,max.terms=10)#½¨Á¢Í¶Ó°Ñ°×ٻعéÄ£ĞÍppr.modelsummary(ppr.model)#²é¿´Í¶Ó°Ñ°×ٻعéÄ£Ğͽá¹ûppr.model$xnamesppr.model$dfPredicted.Value-ppr.model$fitted.values;Predicted.Value#¸ù¾İͶӰѰ×ÙÄ£Ğ͵õ½µÄÔ¤²âÖµPredicted.Residual-ppr.model$residuals;Predicted.Residual#²é¿´Ô¤²âÖµµÄ²Ğ²î£¨ÎÄÏ×Öµ-Ô¤²âÖµ£©trainactivity.and.Predicted-cbind(trainactivity,Predicted.Value,Predicted.Residual)trainactivity.and.Predictedjpeg(filename=ͶӰѰ×ٻعéÄ£ĞÍѵÁ·¼¯½á¹û.jpeg,units=px,width=800,height=600,restoreConsole=TRUE,quality=75)#Êä³öͼĞÎÃüÁîplot(Predicted.Value,trainactivity,xlab='ͶӰѰ×ٻعéÄ£ĞÍѵÁ·¼¯Ô¤²âÖµ',ylab='ʵÑéÖµ',col=red,type=p)#ÎÄÏ×ÖµÓëÔ¤²âÖµ×÷ͼabline(lm(Predicted.Value~trainactivity),col=blue)#Ìí¼ÓÉ¢µãͼµÄÇ÷ÊÆÏßdev.off()#¹Ø±ÕͼĞÎÊä³ör-cor(trainactivity,Predicted.Value);r;r^2#ÇóÎÄÏ×ÖµÓëÔ¤²âÖµµÄÏà¹ØϵÊı#Ïà¹ØϵÊıµÄƽ·½£¬¾ö¶¨ÏµÊı£¬ÔÚºóÃæµÄlm.solÏßĞԻعé·ÖÎöÖĞÓĞlm.sol-lm(trainactivity~Predicted.Value)#¶ÔÎÄÏ×ÖµºÍÔ¤²âÖµ½øĞлعé·ÖÎölm.solsummary(lm.sol)#²é¿´ÎÄÏ×ÖµºÍÔ¤²âÖµ»Ø¹é·ÖÎö½á¹ûpar(mfrow=c(2,2))plot(lm.sol)#ÎÄÏ×ÖµºÍÔ¤²âÖµ»Ø¹é·ÖÎöÕï¶Ïopar-par()write.table(trainactivity.and.Predicted,file='ͶӰѰ×ٻعéÄ£ĞÍѵÁ·¼¯½á¹û.csv',sep=',',col.names=TRUE,row.names=FALSE)#½«ÎÄÏ×Öµ¡¢Ô¤²âÖµÒÔ¼°²Ğ²î½á¹ûÁí´æΪcsvÎļş##############################################################################################################################################################################test.Predicted.Value-predict(ppr.model,testdata)#¸ù¾İͶӰѰ×ÙÄ£Ğ͵õ½µÄÔ¤²âÖµtestactivity.and.Predicted-cbind(testactivity,test.Predicted.Value);testactivity.and.Predictedwrite.table(testactivity.and.Predicted,file='ͶӰѰ×ٻعéÄ£ĞͲâÊÔ¼¯½á¹û.csv',sep=',',col.names=TRUE,row.names=FALSE)#½«ÎÄÏ×Öµ¡¢Ô¤²âÖµÒÔ¼°²Ğ²î½á¹ûÁí´æΪcsvÎļşr-cor(testactivity,test.Predicted.Value);r;r^2#ÇóÎÄÏ×ÖµÓëÔ¤²âÖµµÄÏà¹ØϵÊılm.sol.2-lm(testactivity~test.Predicted.Value)#¶ÔÎÄÏ×ÖµºÍÔ¤²âÖµ½øĞлعé·ÖÎölm.sol.2summary(lm.sol.2)#²é¿´ÎÄÏ×ÖµºÍÔ¤²âÖµ»Ø¹é·ÖÎö½á¹ûpar(mfrow=c(2,2))plot(lm.sol.2)#ÎÄÏ×ÖµºÍÔ¤²âÖµ»Ø¹é·ÖÎöÕï¶Ïopar-par()##############################################################################################################################################################################jpeg(filename=ͶӰѰ×ٻعéÄ£ĞͲâÊÔ¼¯½á¹û.jpeg,units=px,width=800,height=600,restoreConsole=TRUE,quality=75)#Êä³öͼĞÎÃüÁîplot(test.Predicted.Value,testactivity,xlab='ͶӰѰ×ٻعéÄ£ĞͲâÊÔ¼¯Ô¤²âÖµ',ylab='ʵÑéÖµ',col=red,type=p)#ÎÄÏ×ÖµÓëÔ¤²âÖµ×÷ͼabline(lm(test.Predicted.Value~testactivity),col=blue)#Ìí¼ÓÉ¢µãͼµÄÇ÷ÊÆÏßdev.off()#¹Ø±ÕͼĞÎÊä³ö