powered by
data(ozone)
O3
vh
wind
humidity
temp
ibh
dpg
ibt
vis
doy
data(ozone) alm <- lm(O3 ~ vis+doy+ibt+humidity+temp,data=ozone) res <- resid(alm) hctest(ozone[,10], res)
Run the code above in your browser using DataLab