par(mfrow=c(2,1))
plot(fit1$resid)



qqnorm(fit1$resid)

abline(a=0,b=1)