理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

残差分析

残差分析

Rのサンプル ChickWeight を使用
y<-ChickWeight$weight[1:10]
x<-ChickWeight$Time[1:10]
summary(lm(y~x))
plot(x,y,xlab = "生後日数",ylab = "体重")
abline(30.327 ,7.030)         #回帰直線の挿入

f:id:yoshida931:20170802161419p:plain

残差 ei をプロットすることで、回帰モデルからのズレを視覚化してみます

標準化残差Sie

f:id:yoshida931:20170802154446p:plain

ei<-y-(30.327+7.03*x)
eis<-ei/sqrt( (sum(ei^2) )/(10-2) )  #標準化残差
plot(eis,xlab="生後日数",ylab="残差", xaxt="n")
name<-c(0:9)*2
axis(side=1,1:10,labels=name)
abline(h=0, lty=2, col=2)

f:id:yoshida931:20170802161440p:plain

こうやって見たら、回帰直線からの「ズレ」がより理解できます