From: Joerg Schaber (Joerg.Schaber@uv.es)
Date: Fri 14 Nov 2003 - 02:49:34 EST
I naively thought when I can give estimates of the random effects I
should also be able to calculate confidence levels of these estimates
(that's what statistics is about, isn't it?)
For example, similar to the fixed case, I can calculate a
variance-covariance matrix (C) for the random effects (e.g. following
Hemmerle and Hartley,TECHNOMETRICS 15 (4): 819-831 1973) and using the
t-value for the given confidence level and degrees of freedom (t), I can
estimate confidence intervals for random effect i (r[i]) as something
like r[i] +- t*sqrt(C[i][i]).
What does the statistician say?
Douglas Bates wrote:
>Joerg Schaber <Joerg.Schaber@uv.es> writes:
>>I have a linear mixed-effects model object and want to extract the 95%
>>confidence intervals for the fixed and random effects, respectively. I
>>found the function intervals() for confidence intervals for the fixed
>>effects but no corresponding function for the random effects. Does it
>>exist or do I have to calculate the confidence intervals for the
>>random effects myself?
>You have to calculate them yourself, partly because it is not clear
>what such an interval should be. Technically, the random effects are
>not parameters and defining a "confidence interval" on a random
>variable that is part of the model is, at the very least, awkward.
-- ---------------------------------------------------------- Jörg Schaber Instituto Cavanilles de Biodiversidad y Biologia Evolutiva Universidad de Valencia Tel.: ++34 96 354 3666 A.C. 22085 Fax.: ++34 96 354 3670 46071 Valencia, España email : firstname.lastname@example.org
______________________________________________ Remail@example.com mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
This archive was generated by hypermail 2.1.3 : Thu 01 Jan 2004 - 09:20:41 EST