《r语言课件》第六课回归分析

《《r语言课件》第六课回归分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、回归分析残差.吐1••宀■■■X,A'YZ•••——X>AY••••••■■X•■/■残差残差(c)(d)非线性情况有异常值点的情况Box-Cox变换：在做回归分析时，残差不是像上图仗)所示的正常情况，不满足齐性，前面介绍加权方法，这里介绍数据变换方法——Box-Cox变换。在出现异方差情况下，Box-Cox变换可使回归方程的残差满足齐性要求，它对Y做如下变换：Y入-1Y(x)_r入>IInY,入为待定参数，boxcox()函数可以绘出不同参数下对数似然函数的目标值，可以通过图形来选择boxcox

2、(object,lambda=seq(-2,2,1/10),plotit=TRUE,interp,eps=1/50,xlab=expression(lambda),ylab=tflog-Likelihoodft,•••)参数object为lm生成的对象,lambda为参数默认值为-2到2,间隔值0.1,plotit为逻辑变量，表示是否画出图形，默认值为TRUE,interp为逻辑变量，表示是否使用三次样条插值，eps为控制精度，默认值为0・02,・・・为附加参数。例某公司为了研究产品的营销策略，对

3、产品的销售情况进行了调查。设Y为某地区该产品的家庭人均购买量(单位：元)，X为家庭人均收入(单位：元)，下表给出了53个家庭的数据,试通过这些数据建立Y与X的关系式。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)##%%输入数据，作回归方程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张，残差与预测散点图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张，确定参数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张，变换后残差与预测散点图plot(fltted(lm.lam)9resid(lm.lam),cex=1.2,pch=21,col=MredH,bg=HorangeH,xlab=HFit