powered by
lmLikSI(lm.mod, level)
set.seed(10) x <- rnorm(50, 0 , 5) y <- sapply(1:length(x), function(z) 3+.5*x[z]+rnorm(1, 0, 5)) lm.obj <- lm(y~x) lmLikSI(lm.obj, 8)
Run the code above in your browser using DataLab