Hi,

I'm trying to extract a histogram over the results from a bootstrap. However
I keep receiving the error message "Error in hist.default(boot.lrtest$ll,
breaks = "scott") : 'x' must be numeric".
The bootstrap I'm running looks like:

> boot.test <- function(data, indeces, maxit=20) {
+ y1 <- fit1+e1[indeces]
+ mod1 <- glm(y1 ~ X1-1, maxit=maxit)
+ y2 <- fit2+e2[indeces]
+ mod2 <- glm(y2~1, maxit=maxit)
+ ll <- 2*(logLik(mod1)-logLik(mod2))
+ ll
+ }
> boot.lrtest <- boot(data=M1, statistic=boot.test, R=2000, maxit=100);

> hist(boot.lrtest$ll, breaks="scott")  # This results in the error message
stated above

So my question is: what am I doing wrong?

//Anders

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to