資源描述：
《《r語言課件》第六課回歸分析》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、回歸分析殘差.吐1??宀■■■X,A'YZ???——X>AY??????■■X?■/■殘差殘差(c)(d)非線性情況有異常值點的情況Box-Cox變換：在做回歸分析時，殘差不是像上圖仗)所示的正常情況，不滿足齊性，前面介紹加權(quán)方法，這里介紹數(shù)據(jù)變換方法——Box-Cox變換。在出現(xiàn)異方差情況下，Box-Cox變換可使回歸方程的殘差滿足齊性要求，它對Y做如下變換：Y入-1Y(x)_r入>IInY,入為待定參數(shù)，boxcox()函數(shù)可以繪出不同參數(shù)下對數(shù)似然函數(shù)的目標(biāo)值，可以通過圖形來選擇boxcox
2、(object,lambda=seq(-2,2,1/10),plotit=TRUE,interp,eps=1/50,xlab=expression(lambda),ylab=tflog-Likelihoodft,???)參數(shù)object為lm生成的對象,lambda為參數(shù)默認(rèn)值為-2到2,間隔值0.1,plotit為邏輯變量，表示是否畫出圖形，默認(rèn)值為TRUE,interp為邏輯變量，表示是否使用三次樣條插值，eps為控制精度，默認(rèn)值為0?02,???為附加參數(shù)。例某公司為了研究產(chǎn)品的營銷策略，對
3、產(chǎn)品的銷售情況進(jìn)行了調(diào)查。設(shè)Y為某地區(qū)該產(chǎn)品的家庭人均購買量(單位：元)，X為家庭人均收入(單位：元)，下表給出了53個家庭的數(shù)據(jù),試通過這些數(shù)據(jù)建立Y與X的關(guān)系式。X:{679292101249358211569972189109720781818170074720301643414354127674543554087415431029710143483717481381142812551777370231611304637707248087907834061242658174646811144
4、1317873560149522211526}Y:{0.790.440.560.792.703.644.739.505.346.855.845.213.254.433.160.500.171.880.771.390.561.565.280.644.000.314.204.883.487.582.634.990.598.194.790.511.744.103.940.963.290.443.242.145.710.641.900.518.3314.945.113.853.93}解：(程序名exam0
5、607.R)##%%輸入數(shù)據(jù)，作回歸方程X<-scan()6792921012493582115699721891097207818181700747203016434143541276745435540874154310297101434837174813811428125517773702316113046377072480879078340612426581746468111441317873560149522211526Y<-scan()0.790.440.560.792.703.644.
6、739.505.346.855.845.213.254.433.160.500.171.880.771.390.561.565.280.644.000.314.204.883.487.582.634.990.598.194.790.511.744.103.940.963.290.443.242.145.710.641.900.518.3314.945.113.853.93lm?solv-lm(Y?X);summary(lm.sol)##%%加載MASS程序包library(MASS)##%%作圖，
7、共4張op<-par(mfrow=c(2,2),mar=.4+c(4,4,14)?oma=c(0,0,2,0))##%%第1張，殘差與預(yù)測散點圖plot(fitted(lnieSol),resid(lm>sol),cex=1.2,pch=21,col二''red'；bjg=norangeH,xlab=nFittedValue”，ylab=nResidualsn)###fitted(Im.sol)是擬合值(由Y二0.831304+0.003683X求出)，###resid(lm.sol)是殘差##%
8、%第2張，確定參數(shù)lambdaboxcox(lm>sol,lambda=seq(0,1,by=O>l))##%%Box-Cox變換后，作回歸分析Iambda<-0.55;YIam<-(YAlambda-l)/lambdalm?lamv-lm(Ylam~X);summary(lm.lam)##%%第3張，變換后殘差與預(yù)測散點圖plot(fltted(lm.lam)9resid(lm.lam),cex=1.2,pch=21,col=MredH,bg=HorangeH,xlab=HFit
