Displaying 1 result from an estimated 1 matches for "fm3r".
Did you mean:
fm3
2008 Apr 13
2
prediction intervals from a mixed-effects models?
...he 'interval' argument to the 'predict' function was ignored.
# It works works for an 'lm' object, but not an 'lme' object.
One way to do this might be via mcmcsamp of the corresponding
'lmer' model:
library(lme4)
set.seed(3)
samp3r <- mcmcsamp(fm3r, n=10000)
samp3r[1:2,]
Then use library(coda) to check convergence and write a function
to simulate a single observation from each set of simulated parameters
and compute quantile(..., c(.025, .975)) for each prediction level
desired.
However, before I coded that, I thought I woul...