MarcLavielle / Rsmlx

Other
4 stars 2 forks source link

error in pkbuild with stat = T #12

Open Vincent813 opened 1 year ago

Vincent813 commented 1 year ago

I have tested the warfarin dataset and got an error :

Any idea to make it work?

my code

library(tidyverse) # 1.3.2 + R-4.2.3 library(lixoftConnectors) # 2023.1 library(Rsmlx) # 2023.1.1 + Monolix2023R1

dir_data = "./../00_DATA/"

warfarinPK <- list( dataFile = paste0(dir_data, "warfarinPK.csv"), # from the "Rsmlx_demos" headerTypes = c("id", "time", "observation", "amount", "contcov", "contcov", "catcov"), administration = "oral")

out.pk.all <- pkbuild(data = warfarinPK, param = "clearance", criterion = "BIC", new.dir = "./warfarinPK", linearization = T, level = 2, stat = TRUE, settings.stat = list(covToTransform = c("wt", "age"), criterion ="BIC", model=c("covariate", "residualError"), linearization=T))

my code

(truncated) output with error message

... [1] "./warfarinPK/pk_TlagkaV1V2QCl.mlxtran" Estimation of the population parameters... Estimation of the log-likelihood... [INFO] The linearization methods is more relevant around the conditional mode, we encourage you to make this individual mode task in your scenario OFV AIC BIC BICc 646.3027 674.3027 694.8230 711.1722 [INFO] Results have been successfully loaded [INFO] Results have been successfully loaded [INFO] Results have been successfully loaded [INFO] Results have been successfully loaded [INFO] Results have been successfully loaded


Building:


Covariate model: sex age wt Tlag 0 0 0 Tk0 0 0 0 V 0 0 0 Cl 0 0 0

Residual error model: y "combined2" Sampling of the conditional distribution using the initial model ... Error in if (!is.na(ll["standardError"])) llr["s.e."] <- ll["standardError"] : argument is of length zero