s-news
[Top] [All Lists]

nlme, extract stand deviation value

To: s-news@lists.biostat.wustl.edu
Subject: nlme, extract stand deviation value
From: "Li, Qishan" <qishan@email.arizona.edu>
Date: Thu, 26 Apr 2001 11:29:46 -0700
Dear s-newsers:
 
I should be grateful for suggestions on how to do the following:

I am using Splus (nlme)  and I want to extract a parameter and use it later.
The parameter value is 4.077946, which is the StdDev of Asym.
                  > My_buffer<-  StdDev of Asym , the value is  4.077946

I do not know which function I should use. I tried
                      >random.effects(fm2)
but, it does not work.

Many thanks.

Qishan Li
---------------------------------------------
>fm2 <- nlme(weight ~ SSlogis(Time, Asym, xmid, scal), data = Soybean,
            fixed = Asym + xmid + scal ~ 1, random = pdDiag(Asym + xmid
~ 1),
            start = c(18, 52, 7.5))

> fm2
Nonlinear mixed-effects model fit by maximum likelihood
  Model: weight ~ SSlogis(Time, Asym, xmid, scal)
  Data: Soybean
  Log-likelihood: -755.7323
  Fixed: Asym + xmid + scal ~ 1
     Asym     xmid    scal
 18.96426 55.14188 8.55292

Random effects:
 Formula: list(Asym ~ 1, xmid ~ 1)
 Level: Plot
 Structure: Diagonal
            Asym    xmid Residual
StdDev: 4.077946 2.52487   1.1924

Number of Observations: 412
Number of Groups: 48
>





<Prev in Thread] Current Thread [Next in Thread>