FR

Table

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 du fichier "base.fit"
ch <-"c:/nmvi/run/theo/base.nm7/base.fit"
cd<-paste(ch,filename,sep="/")

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

# visualiser les noms des colonnes et verifier:
names(data)
str(data)

# faire un graphique des PRED en fonction des DV

# plot de base DV en fc de PRED, bornes 0 et 10000:
plot(data$DV,data$PRED, xlab="Observed concentrations (ng/mL)", ylab="Predicted concentrations (ng/mL)", 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")

img