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期
   77   78   79   80   81   82   83   84   85   86   87