Is there an easy way to get confidence intervals from "glmm" in Jim Lindsey's library(repeated)? Consider the following slight modification of an example from the help page: > df <- data.frame(r=rbinom(10,10,0.5), n=rep(10,10), x=c(rep(0,5), + rep(1,5)), nest=1:10) > fit <- glmm(cbind(r,n-r)~x, family=binomial, nest=nest, data=df) > summary(fit) Error in print.summary.glmm(structure(list(call = glmm(cbind(r, n - r) ~ : couldn't find function "print.summary.glm" > confint(fit) Waiting for profiling to be done... Error in glmm(cbind(r, n - r) ~ x, family = binomial, nest = nest, data = df, : unused argument(s) (method ...) I'm running R 1.9.0pat under Windows 2000. Thanks, spencer graves
Spencer Graves <spencer.graves at pdf.com> writes:> Is there an easy way to get confidence intervals from "glmm" in > Jim Lindsey's library(repeated)? Consider the following slight > modification of an example from the help page: > df <- > data.frame(r=rbinom(10,10,0.5), n=rep(10,10), x=c(rep(0,5), > > + rep(1,5)), nest=1:10) > > fit <- glmm(cbind(r,n-r)~x, family=binomial, nest=nest, data=df) > > summary(fit) > Error in print.summary.glmm(structure(list(call = glmm(cbind(r, n - r) ~ : > couldn't find function "print.summary.glm"It's a namespace problem. The S3 method, which is being called directly, is now in a namespace.> getAnywhere("print.summary.glm")A single object matching 'print.summary.glm' was found It was found in the following places registered S3 method for print from namespace stats namespace:stats Check the form of the call in question to see if the object being printed does have class summary.glm, in which case the generic function print should be used instead of print.summary.glm.