Displaying 2 results from an estimated 2 matches for "phenomodel".
2003 May 12
1
plot.ranef.lme (PR#2986)
library(nlme)
data(Phenobarb)
na.include <- function(x)x
phe1 <- nlme(conc~phenoModel(Subject, time, dose, lCl, lV),
data = Phenobarb,
fixed = lCl+lV~1,
random= pdDiag(lCl+lV~1),
start = c(-5,0),
na.action = na.include,
naPattern = ~!is.na(conc))
phe.ranef <- ranef(phe1,augFrame=TRUE)
plot(phe.ranef, form=lCl~Wt+ApgarInd)
[Error in max(length(x0),...
2004 Sep 18
0
Pharmacokinetics using R
...d for PK analysis. However, I have to use R
for some purposes anyways. Here is a really strange problem I am running
into and I hope to get some assistance from the group. The following code
works fine in S-Plus (Population estimates and MAP estimates look reasonable
).
cs.fit<-nlme ( CONC~phenoModel(CID,TIME, DOSE, lCl, lV),
fixed = lCl+lV ~1, random = pdDiag(lCl+lV ~1),
data = cs1grp, start = c(-5, 0), na.action=na.include, naPattern = ~
!is.na(CONC)
)
But if I run the same problem in R, I get the following error: "Error in
chol((value + t(value))/2) : the lea...