myplotnls <- function(x){
par(mfrow=c(2,2))
plot(fitted(x), resid(x), xlab="Fitted Values", ylab="Residual", main="Residuals vs Fi
tted
Values")
abline(h=0, lty=2)
qqnorm(resid(x))
plot(fitted(x), sqrt(abs(resid(x))), xlab="Fitted Values", ylab="Sqrt of abs(Residuals
)", main =
"Sqrt of abs(Residuals) vs Fitted Values")
plot(fitted(x)+resid(x), fitted(x), xlab="Response (Y)", ylab="Fitted values", main =
"Fitted
Values vs Response (Y)")
abline(0,1, lty=2)
}