一、数据集描述
某厂生产的一种电器年销售量Y与竞争对手的价格X1及本厂的价格X2有关,数据如下:
10个城市某种电器的年销售量和竞争对手的价格(单位:元) | |||||
X1 | X2 | Y | X1 | X2 | Y |
120 | 100 | 102 | 140 | 110 | 100 |
190 | 90 | 120 | 130 | 150 | 77 |
155 | 210 | 46 | 175 | 150 | 93 |
125 | 250 | 26 | 145 | 270 | 69 |
180 | 300 | 65 | 180 | 250 | 85 |
二、建立Y与X1与X2的回归关系
x1<-c(120,190,155,125,180,140,130,175,145,150)
x2<-c(100,90,210,250,300,110,150,150,270,250)
y<-c(102,120,46,26,65,100,77,93,69,85)
df<-data.frame(y,x1,x2)
lm.reg<-lm(y~x1+x2,data=df)
summary(lm.reg)
运行得到:
二、计算本例的残差和标准残差,并绘制残差散点图
y.res<-residuals(lm.reg)
print(y.res)
y.rst<-rstandard(lm.reg)
print(y.rst)
y.fit<-predict(lm.reg)
op<-par(mfrow=c(1,2))
plot(y.res~y.fit)
plot(y.rst~y.fit)
par(op)
运行得到:
标准化残差的绝对值都小于2,可见本例样本数据并无异常点。