powered by
data(ustemp)
y
x1
x2
data(ustemp) lmod<-lm(y~x1+x2, data=ustemp) x1<-ustemp[,"x1"] x2<-ustemp[,"x2"] res<-resid(lmod) hctest(x1, res)
Run the code above in your browser using DataLab