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
>>
>>
>>
>
>
>
> The way I've done the bootstrapping
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 "l
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).
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. M