Tableaux de résultats

Un ou plusieurs tableaux de résultats peuvent être créés selon ce qui est spécifié dans le block TABLE.

Ces tableaux sont être créés sous la forme de fichier(s) ASCII.

Script R pour ouvrir une table créée par nonmem...

# chemin complet
ch <-"c:/nmvi/run/theo/base.df"

# nom du fichier à ouvrir
filename<-"base.fit"

cd<-paste(ch,filename,sep="/")

data<-read.table(cd,header=TRUE,skip=1)

# pour visualiser les noms des colonnes
names(data)

... et pour faire un graphique des PRED en fonction des DV

# plot de base
plot(data$DV,data$PRED, xlab="Observed concentrations (mg/L)", ylab="Predicted concentrations (mg/L)", xlim=c(0,10000), ylim=c(0,10000))
options(na.action=na.exclude)

# ajout de la régression linéaire (après avoir éliminé les valeurs à t0)
data2 <- data[data$TIME>0,]
lm.obs <- lm(data2$PRED~data2$DV)
# lty peut prendre plusieurs valeurs: 0=blank, 1=solid, 2=dashed, 3=dotted
abline(lm.obs, lty=3, col="red")
# ajout de la droite théorique
abline(a=0, b=1, lty=1, col="blue")