一、示例文件
本示例使用的示例文件来自于Phoenix8.0的示例文件,文件路径为C:\Program Files (x86)\Certara\Phoenix\application\Examples\NLME\Theo_Model.phxproj
图表 1
二、了解数据
1.双击“Theo_Model.phxproj”项目文件,在Phoenix中将该示例文件打开。
图表 2
2.该示例中“Partial Diagonal Model tlag”对象是该数据所对应的最终模型,在左侧对象浏览器中鼠标左键单击该对象,之后在右侧内容展示窗中的“Results(结果)”子选项卡下,选择“Pop DV. IPRED vs IVAR”图表,可查看个体预测曲线与观测点对时间的散点图。
图表 3
三、绘制VPC图(不对数据分组)
1.选中项目文件中的最终模型“Partial Diagonal Model tlag”,复制并粘贴该模型,并将名称更改为“VPC1”
2.导航到“VPC1”对象的“Run Options”选项卡下,将运行模式切换为“Sim./Pred/Check”,
#replicates: (重复次数)设定为10,然后执行该对象。
图表 4
3.执行完成,在其“Results(结果)”子选项卡下导航至其图表结果“Pop PredCheck ObsQ_SimQ”,可在右侧查看到其VPC图形
图表 5
为什么图形看起来会这么糟糕呢?
因为没有对图中的数据分组(或者说分段),接下来我们重新绘制一下VPC图形,在绘制的同时进行分组。
四、绘制VPC图(对数据分组)
1.首先复制操作对象VPC1,将新的操作对象命名为VPC2;
2.导航至其Run Options标签页,勾选“Keep?”,固定随机种子;
3.接着来到“Binning(分组)”子标签页,下拉“Binning Option(分组选项)”下拉框,选择“Explicit Centers(指定分组的中心)”,然后在出现的对话框中输入数值序列“0,0.25,0.5,1,2,4,5,7,9,12,24”(输入时不需要输入引号,逗号需要为英文半角逗号);
图表 6
“Binning Option(分组选项)”提供了四种分组方式:None不分组,K-means应用K-means算法对数据分组,Explicit centers用户指每组的中心对数据进行分组,Explicit boundaries用户指定组与组的边界值进行分组,之后执行操作对象;
4.执行完成,在其“Results(结果)”子选项卡下导航至其图表结果“Pop PredCheck ObsQ_SimQ”,可在右侧查看到其VPC图形;
图表 7
这样我们就绘制得到了一个典型的VPC图。
五、绘制VPC图(同时绘制预测值的PI的置信区间)
1.复制VPC2对象,并命名为VPC3
2.导航的VPC3的Run Options标签页
3.选择“Predictive Check Options”的子选项卡“Qunantiles(分位数)”,勾选“Qunantile%”选项,之后执行操作对象;
图表 8
4.执行完成,在其“Results(结果)”子选项卡下导航至其图表结果“Pop PredCheck ObsQ_SimQCI”,可在右侧查看到其VPC图形;
图表 9
六、图形绘制细节调整
如果希望修改图表中的绘制风格,可双击拟修改的图表将其在新的窗口中打开,然后就可以通过图表下部的菜单来修改图片。
图表 10
图表 11
七、小节
到此我们就绘制得到了们平常所见到的几种VPC图,通过Phoenix Model,我们可以方便快速的绘制出我们想要的VPC图,以满足报告要求;在下次案例分享,将通过手工绘制的方式,一步一步的绘制出VPC图,使大家可以更加容易理解VPC图的绘制原理。