powered by
data(windmill)
Y
x
data(windmill) g1<-lm(Y~x,data=windmill) res<- resid(g1) x<- windmill[,2] hctest(x,res)
Run the code above in your browser using DataLab