powered by
hmctest
hmctest(formula, T, data=list())
"htest"
lm
x <- c(1:30) err <- c(rnorm(10,0,1), rnorm(20,0,10)) y <- x^2 + err formular <- y ~ x hmc <- hmctest(formular, 10)
Run the code above in your browser using DataLab