NONMEM Users Network Archive

Hosted by Cognigen

RE: R/S+ nlme ODE

From: Denney, William S. <william_denney>
Date: Mon, 23 Nov 2009 13:46:54 -0500

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: owner-nmusers
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 infor=
mation of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jerse=
y, USA 08889), and/or its affiliates Direct contact information for affilia=
tes is available at http://www.merck.com/contact/contacts.html) that may be=
 confidential, proprietary copyrighted and/or legally privileged. It is int=
ended 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 e=
rror, please notify us immediately by reply e-mail and then delete it from =
your system.
Received on Mon Nov 23 2009 - 13:46:54 EST

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.