Hi Leonid, This is how to do it for fixed effects in R with a regular nlme model (I've not used nlmeODE). I don't have an example with generating new random variables.
library(MASS) ## make a model named "model" here coef <- model$coefficients$fixed cov.matrix <- model$varFix new.coef <- as.data.frame(mvrnorm(n=10000, coef, cov.matrix, empirical=TRUE)) Thanks, Bill -----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Leonid Gibiansky Sent: Monday, November 23, 2009 11:50 AM To: nmusers Subject: [NMusers] R/S+ nlme ODE Dear All, This is off-topic but I tried to ask this question in R-help list and it was not successful: 132 email/day and no reply to this question: I have a population PK model in nlmeODE (R or S+ versions) and I need a posthoc (in Nonmem terminology) estimate: estimate individual parameters of new subjects given population parameters and variance-covariance matrix. Is this possible? If yes, could you send me a sample code? Thanks! Leonid -------------------------------------- Leonid Gibiansky, Ph.D. President, QuantPharm LLC web: www.quantpharm.com e-mail: LGibiansky at quantpharm.com tel: (301) 767 5566 Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.
