powered by
data(tensile)
Y
x
data(tensile) tmod1 <- lm(Y~x+I(x^2), tensile) x<-tensile[,2] res <- resid(tmod1) hctest(x, res)
Run the code above in your browser using DataLab