Page 82 - 《中国药房》2024年1期
P. 82
表2 协变量筛选过程 80
步骤 协变量筛选 OFV ΔOFV P 备注
1 无 3 276.31 基础模型
60
正向纳入
2 CL-weight 3 233.09 -43.22 <0.01
3 CL-weight,V 2-ALB 3 205.18 -27.91 <0.01 ( μg/mL ) 40
4 CL-weight,V 2-ALB,Q-ALB 3 183.62 -21.56 <0.01 DV/
5 CL-weight,Q-weight,V 2-ALB,Q-ALB 3 174.84 -8.78 <0.01
6 CL-weight,Q-weight,V 1-ALB,V 2-ALB,Q-ALB 3 165.44 -9.40 <0.01 20
7 V 2-weight,CL-weight,Q-weight,V 1-ALB,V 2-ALB,Q-ALB 3 156.89 -8.55 <0.01 全量模型
逆向剔除
8 V 2-weight,Q-weight,V 1-ALB,V 2-ALB,Q-ALB 3 158.21 1.32 >0.001 0
9 V 2-weight,Q-weight,V 2-ALB,Q-ALB 3 167.10 8.89 >0.001 最终模型 0 20 40 60 80
10 Q-weight,V 2-ALB,Q-ALB 3 203.61 25.68 <0.001 PRED/(μg/mL)
11 V 2-weight,V 2-ALB,Q-ALB 3 202.72 24.79 <0.001 A. PRED-DV散点图
12 V 2-weight,Q-weight,Q-ALB 3 203.28 25.35 <0.001 80
13 V 2-weight,Q-weight,V 2-ALB 3 186.30 19.21 <0.001
表3 MPA最终模型的PPK参数及自举结果 60
最终模型 Bootstrap
参数 相对偏差/%
估计值 RSE/% 中位值 95%置信区间 ( μg/mL )
θ Ka/(/h) 1.29 8.95 1.31 1.06~1.66 2.19 40
θ V1/(L) 5.79 18.80 5.79 3.42~8.97 -0.02 DV/
θ V2/(L) 430.93 14.59 426.45 221.97~641.05 -1.04
θ CL/(L/h) 4.06 16.04 4.04 2.45~6.47 -0.46 20
θ Q/(L/h) 15.40 7.11 15.24 12.18~17.8 -1.01
1.68 13.09 1.73 1.20~2.39 3.09
θ V2,weight
1.14 14.04 1.15 0.85~1.60 0.55 0
θ Q,weight
0 20 40 60 80
-1.64 -22.87 -1.68 -3.17~-0.51 2.18
θ V2,ALB
IPRED/(μg/mL)
-1.45 -18.91 -1.48 -2.12~-0.38 2.17
θ Q,ALB B. IPRED-DV散点图
个体间变异 4
2
ω Ka 0.20 32.51 0.16 0.16~0.17 -18.24
2 1.07 38.49 1.05 0.98~1.05 -2.19
ω V1
2 y=2
ω CL 0.26 23.02 0.23 0.22~0.24 -9.45 2
2
ω V2 0.15 36.59 0.13 0.11~0.13 -11.82
2 0.13 21.26 0.12 0.12~0.13 -4.35
ω Q
残差变异 CWRES 0 y=0
σ 2 0.58 3.86 0.58 0.56~0.57 -0.98
注:RSE为相对标准误差;θ Ka、θ V1、θ V2、θ CL、θ Q分别为各药动学参数群 y=-2
体典型值;θ V2,weight为体重对V 2的固定效应参数,θ Q,weight为体重对Q的固定 -2
效应参数;θ V2,ALB为ALB对V 2的固定效应参数,θ Q,ALB为ALB对Q的固定
效应参数;ω Ka、ω V1、ω CL、ω V2、ω Q为各药动学参数个体间变异方差;σ 2
2
2
2
2
2
-4
为残差变异方差。相对偏差(%)=(自举验证参数中位值-最终模型
0 5 10 15 20 25
参数估算值)/最终模型参数估算值×100%。 PRED/(μg/mL)
2.2.3 模型评价结果 C. CWRES-PRED散点图
4
最终模型的拟合优度图(图 1)显示 MPA 最终模型
具有良好的拟合优度。其中,图1A和图1B分别为最终 y=2
2
模型群体预测值(predicted value,PRED)与观测值(de‐
pendent variable,DV)的 散 点 图 及 个 体 预 测 值(indi‐ y=0
vidual predicted value,IPRED)与 DV 的散点图,图中数 CWRES 0
据点基本在y=x直线两侧对称分布,说明模型拟合度良
y=-2
好;图 1C 和图 1D 分别展示了最终模型条件加权残差 -2
(conditional weighted residuals-predicted values,CWRES)
与 PRED 的散点图及 CWRES 与霉酚酸酯给药后时间 -4
(time after dose,TAD)的散点图,图中大多数数据点在 0 3 6 9 12
TAD/h
CWRES=0两侧对称随机分布,大部分数据在-2~2之 D. CWRES-TAD散点图
间,无明显的趋势性分布,提示模型预测性能良好。 图1 最终模型的拟合优度图
· 72 · China Pharmacy 2024 Vol. 35 No. 1 中国药房 2024年第35卷第1期