Re: [R] lme and confidence intervals

From: Gregory Warnes <gregory.warnes_at_mac.com>
Date: Thu, 10 Apr 2008 10:09:33 -0400

FWIW, the ci() function in the gmodels package supports generating confidence intervals for the fixed effects of lme objects.

-G

On Apr 8, 2008, at 1:34PM , Dieter Menne wrote:
> Cristian Carranza <cristiancarranza_1 <at> hotmail.com> writes:
>> After fitting a mixed effects model to repeated measurements data
>> set, and
> after several unsuccessful
>> atempts to make a simple plot of the confidence interval for the
>> fitted model,
> I gave up and now I am asking
>> for help in this useful list.
>>
>> Could anyone be so kind to give me some code lines in order to
>> make a plot of
> the fitted equation and the
>> correspondent confidence interval?
>
> The fitted equation is simple: use predict in its variations. For the
> "confidence interval", most people on the list will ask back what
> it means. But
> I know my colleagues, they or their reviewers insist on having
> error bar in the
> graphics, and don't care about the interpretation.
>
> Dimitri has given an approach that could be used to produce some
> reviewer-satisfaction plots limited borborygmy.
>
> http://finzi.psych.upenn.edu/R/Rhelp02a/archive/118174.html
>
> Note that his method is meant to estimate well-defined confidence
> intervals for
> the coefficients, but you can use the profiled fits and compute the
> "confidence
> intervals" at the data points. Alternatively, you could use mcmc to
> get a range
> of curves for averaging.
>
> Dieter
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 10 Apr 2008 - 14:16:35 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 10 Apr 2008 - 14:30:27 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive