Hi,

maybe try 'intervals()':



> fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc),
+             data = Loblolly,
+             fixed = Asym + R0 + lrc ~ 1,
+             random = Asym ~ 1,
+             start = c(Asym = 103, R0 = -8.5, lrc = -3.3))


> intervals(fm1)
Approximate 95% confidence intervals

 Fixed effects:
         lower       est.      upper
Asym 96.625880 101.449600 106.273320
R0   -9.250358  -8.627331  -8.004303
lrc  -3.300904  -3.233751  -3.166598
attr(,"label")
[1] "Fixed effects:"

 Random Effects:
  Level: Seed 
            lower     est.   upper
sd(Asym) 2.479449 3.650642 5.37506

 Within-group standard error:
    lower      est.     upper 
0.6090031 0.7188625 0.8485398 
> 


HTH,

Luca





-----Original Message-----
From: SHAILLY MEHROTRA <shaillymehrotra@gmail.com>
To: r-sig-mixed-models@r-project.org
Date: Tue, 7 Feb 2012 09:40:58 -0500
Subject: [R-sig-ME] nlme -95% CI of parameter estimates

Dear R users,

I am using nlme to fit my model and estimate fixed and random effects. I
want to calculate  confidence intervals of my parameter estimates. Is there
some generic function to calculate 95%CI of the estimates in nlme?

 I get standard errors of the parameter estimates using summary command. I
want to know how to output it as a data set so that I can use the standard
error to calculate CI.

Best Regards

Shailly Mehrotra

Center for Translational Medicine

School of Pharmacy

University of Maryland

Baltimore, MD, USA

   [[alternative HTML version deleted]]

_______________________________________________
R-sig-mixed-models@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
	[[alternative HTML version deleted]]


