基础模型的建立

基础模型的建立应尽量简洁,其中可以包括基本的药代动力学参数以及个体间差异(IIV)以及残留误差

当这些参数符合模型诊断图后,基本上就可以认为达到稳定的基础结构模型。(base structural model)

一般来说,基础模型不应该包括协变量。但对于那些显著影响药代动力学参数(parameter)的协变量,在基础模型中就应该引入(比如儿童群体的体重)。

那么,该如何确定你的基础模型建立的是合适的呢?

一、通过结果表格

Stderr:

CV%:

CI置信区间:不要经过0


一、模型诊断图

几个概念:

DV:一个数据集中因变量(dependent)的数值。

PRED: 模型根据DV模拟的数值。

RES:DV与PRED的绝对差值。

WRES:加权后的DV与PRED的相对差值,一般可以认为加权后为均值为0方差为1的正态分布。一个好的模型WRES应该在-3到+3之间。

1. PRED vs. DV

图形点应该在截距为0,斜率为1的单位线的两侧对称随机分布。

像这个图就显示点在线两边的随机分布,但是线的斜率不为1,且横纵坐标不是一样的。当把坐标改为一样时就会发现预测值并不能全部覆盖真实值。预测值只覆盖了0到200,而真实值到了接近400。如下图。

2. WRES vs. PRED

当误差模型合理时,点应该在WRES=0的线上匀称分布。

拟合良好的的误差图
当出现U字型时,说明模型需要增加一个房室 
表明数据最好用对数残留误差模型

3. RES vs. PRED

RES是预测值与实际值的实际差异,数据点应该在零位线的两侧均匀分布。由比例型误差模型得到的数据往往允许数据较大的点的误差大一些,这样算出来的相对误差才能一致化的最小,因此呈现出扇形。

用比例型误差模型的图

4. WRES vs. TIME

可以了解这些残留误差是否随时间而变化。

模拟良好的图

二、随机效应的估计

选择个体间随机效应的原则是一开始先选择这个数据集可以描述的较小的参数的随机效应,然后再根据情况逐一加入其它的随机效应。即由少到多逐一添加。

比如说:如果一个数据集在谷浓度的点相对于峰浓度来得多,那么就应该先加入的是清除率的个体间差异η。

残留误差与随机效应可能会相互影响。有时候错误的残差模型的估计会导致ETA没有办法被合理估计,但是一旦残差模型设置合理后,这个ETA就可以成功估计了,以下是一些随机效应的诊断图。

几个概念:

IPRED:个体预测值,计算时包括了η。

IRES:个体预测值与实际值的差值。

IWRES:加权后的IRES。

1. IPRED vs. DV

由于考虑到了个体间差异η,因此图形拟合的会比 PRED vs. DV 好。图形差不多,就不放这里了

2. IWRES vs. IPRED and |IWRES| vs. IPRED

下面是一张|IWRES| vs. IPRED的图,由于是绝对值,因此纵坐标大于0,模拟良好的图应该是那条模拟的线基本呈一条水平光滑的直线,且大部分点落于小于3的位置上。

模拟良好的图

而在模型拟合中,加和型误差模型能够反映差值的绝对值。如果数据集的跨度范围过大,对于数据小的点来说,由于数值小,因此误差也会小,而对于数值大的点来说,由于本身的点就大,因此数值大的点误差也会大,导致加和型残留误差模型不足以解决这么大的差异,就会呈现出下图的模样。此时最好用比例型误差模型,即每个误差的绝对值再除以本身的值。这样就能形成相对误差。

需要比例误差模型却用了加和模型

错误的比例型误差模型应用时,由于小的数值的分母小,造成算出来的相对误差会增大,在图上显示出小的误差相对大的误差来得大,此时可能用加和型的误差模型会更好。

需要用加和模型却用了比例模型


2. 推测参数估计的准确度

当协方差步骤完成后,THETA,OMEGA,SIGMA的标准误都可以得到。这个标准误除以估计值得到的是percent standard error of the mean(%SEM)或者是 percent relative standard error(%RSE)。

一般来说,对于固定效应参数的RSE应当在30%以下,而随机效应的参数应该在40~50%之间。

随机效应OMEGA or SIGMA也有其他的计算方法,比如OMEGA的RSE的计算通常是?完成。

一般来说,所有参数的相关性应该要小于0.9 ,这样的话每个估计值才能是独立的。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,313评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,369评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,916评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,333评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,425评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,481评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,491评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,268评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,719评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,004评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,179评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,832评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,510评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,153评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,402评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,045评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,071评论 2 352