[R] test of significance for nlme coefficients
Lindsey Root
lroot at umiami.ir.miami.edu
Fri Mar 4 19:49:54 CET 2005
Hello,
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
More information about the R-help
mailing list