基于R语言进行Monolix建模后的Boostrap验证

Monolix建模后进行Boostrap验证可以基于GUI界面进行，如果会R语言的话，非常简单。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
library(Rsmlx)
library(RJSONIO)
library(lixoftConnectors)

# 检查是否正常运行
initializeLixoftConnectors(software = "monolix")
loadProject("./***.mlxtran") # 载入模型
runScenario()
getEstimatedPopulationParameters()
# 一旦本代码验证过后，上述代码可以删除


# Final boosting
initializeLixoftConnectors(software = "monolix")

project <- "./project_file_path/***.mlxtran" # 修改对应的文件路径
boostdat <- bootmlx(project, nboot = 1000, parametric = T, settings = list(plot =T, level = 0.95))

print(boostdat)

dat <- read.table("./project_file_path/boostrap/parametric/PopulationParameters.txt",header = T, sep = ",") # 修改对应的文件路径

head(dat)
colMeans(dat)
median(dat$v_pop)