MarcLavielle / mlxR

mlxR
Other
19 stars 5 forks source link

SimulX - Compute continuous output (AUC) from a Monolix Project #15

Open NicolasAzzopardi opened 6 years ago

NicolasAzzopardi commented 6 years ago

Dear Marc,

I used to compute AUC over time in the OUTPUT block of Monolix 4.x (table={AUC}). In (The excellent) Monolix2018R1, it is not possible anymore outside observation times...

So I tried to do it with SimulX as explained here (http://simulx.webpopix.org/userguide/monolix1/):

In "~/Lixoft/MonolixSuite2018R1/resources/demos/monolix/1.creating_and_using_models/1.1.libraries_of_models" folder I ran in R:

library(mlxR)
project.file <- 'theophylline_project.mlxtran'
out1  <- list(name = 'Cc', time = seq(0, 25, by=0.1))
sim.param <- "mode"
res5  <- simulx(project = project.file, output = out1, parameter = sim.param)

Works just fine!!

I tried to use the same parameter "mode" for my Monolix project but it felt with the not so obvious error:

Error in `$<-.data.frame`(`*tmp*`, id, value = 23L) : 
  replacement has 1 row, data has 0

After some tests, I figure out that the cause of this error was the presence of "table={AUC}" in the monolix structural model file...

It may be interesting to have a more explicit error message.

Kind regards, Nicolas

MarcLavielle commented 6 years ago

Hi Nicolas,

I was not able to reproduce this bug... :-) or :-( ? can you send me the project please?

best, Marc

2018-04-10 15:31 GMT+02:00 Nicolas Azzopardi notifications@github.com:

Dear Marc,

I used to compute AUC over time in the OUTPUT block of Monolix 4.x (table={AUC}). In (The excellent) Monolix2018R1, it is not possible anymore outside observation times...

So I tried to do it with SimulX as explained here ( http://simulx.webpopix.org/userguide/monolix1/):

In "~/Lixoft/MonolixSuite2018R1/resources/demos/monolix/1. creating_and_using_models/1.1.libraries_of_models" folder I ran in R:

if (isAwesome){ library(mlxR)project.file <- 'theophylline_project.mlxtran'out1 <- list(name = 'Cc', time = seq(0, 25, by=0.1))sim.param <- "mode"res5 <- simulx(project = project.file, output = out1, parameter = sim.param) }

Works just fine!!

I tried to use the same parameter "mode" for my Monolix project but it felt with the not so obvious error:

if (isAwesome){Error in $<-.data.frame(*tmp*, id, value = 23L) : replacement has 1 row, data has 0 }

After some tests, I figure out that the cause of this error was the presence of "table={AUC}" in the monolix model file...

It may be interesting to have a more explicit error message.

Kind regards, Nicolas

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/MarcLavielle/mlxR/issues/15, or mute the thread https://github.com/notifications/unsubscribe-auth/AJ3yugzyIrRDuObitsWVNLlZIcTDZ9irks5tnLQ5gaJpZM4TOQx0 .

-- https://soutienvenezuela.org/en/ https://soutienvenezuela.org/en/

Marc Lavielle

CMAP, Ecole Polytechnique91128 Palaiseau Cedex

Francehttps://soutienvenezuela.org/en/ Tel: (33) 1 69 33 46 00 (33) 6 37 31 93 05

Marc.Lavielle@inria.fr http://www.cmap.polytechnique.fr/~lavielle

"""""""""""""""""""""""""""""""""""""""""""""""""""""""""