search for: samp3r

Displaying 1 result from an estimated 1 matches for "samp3r".

Did you mean: samp3
2008 Apr 13
2
prediction intervals from a mixed-effects models?
...7 34.45574 # NOTE: The '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 th...