[R] test of significance for nlme coefficients

From: Lindsey Root <lroot_at_umiami.ir.miami.edu>
Date: Sat 05 Mar 2005 - 05:49:54 EST


I am using the nlme package to fit an exponential decay model to a longitudinal data set:

> ##Model a variable called RASCH16A as an exponential function ## of an
initial value "A" and a decay rate "K"
> fm2 <- nlme(model = RASCH16A ~ A* (exp(-K*DAYS)),
+ ###A has nonzero fixed and random effects, K has nonzero fixed and random effects

+ fixed = A + K ~ 1,
+ data = rasch,
+ random = A + K ~ 1 | SUBJECT,
+ ##Give a start value for the fixed effects, expressed using the
concatenate function
+ start = c(60, 0.007))

However, 20 of the 161 estimates decay rates are NOT decay rates, but are very small positive growth rates (all less than a standard deviation (.03) from zero). I would like to run a significance test to see how many of these positive growth rates are significantly different from zero. However, the nlme package only offers confidence intervals for the model parameters, not the person-specific coefficients. I would welcome any and all suggestions.

Thanks if you can help,
Lindsey Root

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Mar 05 05:59:27 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:30:40 EST