| diagplot {PKtools} | R Documentation |
diagplot creates plots of observed versus predicted values and residuals (ordinary and standardized) versus predicted values for both the population (marginal) and individual (conditional) levels.
diagplot(x,...)
x |
object of class, NONMEM, PKNLME, WinBUGS |
... |
additional arguments to be passed to lower level functions |
Plots of observed versus predicted values and residuals (ordinary and standardized) versus predicted values for both the population (marginal) and individual (conditional) levels.
M.S. Blanchard <sblanchard@coh.org>
residplot, obvsprplot, identify
if (.Platform$OS.type =="windows"){
curwd=getwd()
if (file.exists("C:/nmv/run")) {
setwd("C:/nmv/run")
library(nlme)
library(PKtools)
data(Theoph)
Theoph<-Theoph[Theoph$Time!=0,]
id<-as.numeric(as.character(Theoph$Subject))
dose<-Theoph$Dose
time<-Theoph$Time
conc<-round(sqrt(Theoph$conc),4)
Theo<-data.frame(cbind(id,dose,time,conc))
names(Theo)<-c("id","dose","time","conc")
wt.v<-Theoph$Wt
nameData<-list(covnames=c("wt"),
yvarlab="Sqrt(Theop. Conc.) Sqrt(mg/L)",
xvarlab="Time since dose (hrs)",
reparams=c("Ka", "V", "Cl"),
params=c("Ka", "V", "Cl"),
tparams=c("log(Ka)", "log(V)", "log(Cl)"),
varnames=c("D[1,1]","D[1,2]","D[2,2]","D[1,3]","D[2,3]","D[3,3]")
)
data<-list(pkvar=Theo, cov=wt.v)
NM<-RunNM(inputStructure="control.model5", data=data, nameData=nameData)
diagplot(NM)
setwd(curwd)
}
else{
"You do not have NONMEM."
}
}