I can think of two ways to get confidence intervals on intraclass
correlations (ICCs) and more accurate intervals for variance
components: (1) modifying 'simulate.lme' to store the estimated
variance components as well as logLik and (2) using 'lmer' and
'mcmcsamp' in library(lme4).
On Sun, 2006-10-29 at 11:06 -0800, Spencer Graves wrote:
I can think of two ways to get confidence intervals on intraclass
correlations (ICCs) and more accurate intervals for variance
components: (1) modifying 'simulate.lme' to store the estimated
variance components as well as
see in line
Rick Bilonick wrote:
On Sun, 2006-10-29 at 11:06 -0800, Spencer Graves wrote:
I can think of two ways to get confidence intervals on intraclass
correlations (ICCs) and more accurate intervals for variance
snip
The way I've done the bootstrapping
I'm using the lme function in nmle to estimate the variance components
of a fully nested two-level model:
Y_ijk = mu + a_i + b_j(i) + e_k(j(i))
lme computes estimates of the variances for a, b, and e, call them v_a,
v_b, and v_e, and I can use the intervals function to get confidence
intervals.