Re: [R] Nonlinear mixed effects model

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Sun 15 Oct 2006 - 02:11:07 GMT

      Have you read Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer), especially the discussion of 'corStruct' classes in ch. 5 and the nonlinear mixed-effects models in Part II? Bates has kept a small army of graduate students busy for over 20 years developing good ways to do things like this, and this book is probably the best-known product of that work.

      Hope this helps. 
      Spencer Graves

Cam wrote:
> In nonlinear mixed effects models, SAS doesn't allow for free manipulation of
> the covariance matrix (you can only specify a "type", and our "type" doesn't
> exist). Can R accomplish this? For example:
>
> Parameters:
> B1= Beta 1i
> B2= Beta 2i
> G1= Gamma i
>
>
> y = B1 -(B1 - B2) exp { - G1 time} + e
>
> the covariance matrix for
> (B1 [( covB1? covB1B2 covB1G1
> B2 ~ covB2B1 covB2? covB2G1
> G1) covG1B1 covG1B2 covG1? )]
>
> **If we want to specify covG1B1 and make everything else 0's, for example,
> what would the code look like?
>
>
>
> Thanks for your time.
>



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 and provide commented, minimal, self-contained, reproducible code. Received on Sun Oct 15 12:17:04 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sun 15 Oct 2006 - 02:30:10 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.