NONMEN 控制文件组成

每一条命令行或者模块以特定符号“$”开头，可缩写。

控制文件主要由 $PROBLEM、 $DATA、 $INPUT 、 $PRED 、 $THETA、 $OMEGA、 $SIGMA、 $ESTIMATION 、 $TABLE、等模块组成。

一、$PROBLEM模块

描述文件的信息，作为控制文件的标题或者注释

$PROB 3cp_PK

二、$DATA模块

指定所需数据文件及读取路径，后面和附加选项 IGNORE：

$DATA infusion_3cp.csv IGNORE=C

IGNORE还可以指定特定数据：

$DATA infusion_3cp.csv IGNORE= (AGE.LT.12)

FORTRAN 常用的判断逻辑

==或.eq. 判断是否相等
/=或.ne. 判断是否不等
>或.gt. 判断是否大于
>=或 .ge. 判断是否大于等于
<或.lt. 判断是否小于
<=或.le. 判断是否小于等于

IGNORE还可以增加条件选项，如：

$DATA infusion_3cp.csv IGNORE=C  IGNORE= (AGE.LT.12)

$DATA infusion_3cp.csv IGNORE=C  ACCEPT= (WK.GT.40,AGE.GT.12)

三、$INPUT模块

指定数据文件的数据结构。模块定于的变量顺序需和数据文件的每列

$INPUT C ID TIME TAD NMT CONC=DV AMT RATE MDV DOSE GROUP sex=DROP age weight height bmi

CONC=DV
表示用CONC取代保留变量DV 即为CONC为药代动力学浓度的观测值
sex =DROP
表示忽略改变量

四、$SUBROUTINES模块

$SUBROUTINES用于确定调用哪种子程序。

ADVAN指定房室模型种类

TRANS指定基本药动学参数

ADVAN1-4、10-12是内置的药动学模型，分别指一室静脉、一室口服、二室静脉、二室口服、米氏模型、三室静脉、三室口服。选择了以上的ADVAN子程序后，还需选择相应的TRANS子程序。
如：

$SUBROUTINES ADVAN11 TRANS4

ADVAN和TRANS所需的参数表：

五、$PK模块

ADVAN和TRANS所需的必须参数和附件参数，$PK模块指定这些参数的典型值和体间变异。

定义参数：

TVCL = THETA(1)

定义体间变异和类型：

CL = TVCL +ETA(1)            ;加和型
CL = TVCL * (1+ETA(1))       ;比例型
CL = TVCL * EXP(ETA(1))      ;指数型

以TV开头的变量名表示参数的典型值，如TVCL表示清除率的典型值，CL表示特定的个体值。

如：

$PKTVCL=THETA(1)CL=TVCL*EXP(ETA(1)) TVV1=THETA(2)V1=TVV1*EXP(ETA(2))TVV2=THETA(3)V2=TVV2*EXP(ETA(3))TVQ2=THETA(4)Q2=TVQ2*EXP(ETA(4))TVV3=THETA(5)V3=TVV3*EXP(ETA(5))TVQ3=THETA(6)Q3=TVQ3*EXP(ETA(6))

附加参数定义

   TVALAG1=THETA(7)ALAG1=TVALAG1*EXP(ETA(7))

IF= (DOSE.GT.3)  F1=1
IF= (DOSE.LT.3)  F1=THETA(8)

定义换算系数

定义 Sn 可使给药量和药物浓度单位纲保持一致。 n表示观测事件发生的隔室, Sn 表示分布容积和无单位标量值得乘积：
Sn=V×usvS_n = V ×usv Sn=V×usv
Cn=amountSn=amountV×usvC_n = \frac{amount}{S_n} = \frac{amount}{V ×usv} Cn=Snamount=V×usvamount
如：给药剂量为mg，分布体积单位为L,浓度观测单位为ng/ml则
Cn=mgL×usv=ngmlC_n = \frac{mg}{L×usv} = \frac{ng}{ml} Cn=L×usvmg=mlng
此时：usvusvusv =1000

   S1=V1/1000S2=V2/1000S3=V3/1000

IF语句应用

选择结构语句

if (logical expression) then      statement(s)
elseother_statement(s)
end if

嵌套选择结构语句

if (logical expression) then      statement(s)
elseif  (logical expression2)  then
statement2(s)
elseif  (logical expression3)  then
statement3(s)
······
other_statement(s)
end if

示例：

IF (RATE.EQ.1) THENF1=1
ELSEF1=THETA(1)
ENDIF

六、$ERROR模块

描述预测值F与实测值Y的差异，称为残差变异。

简单残差变异模型：加和型，比例型，指数型，结合型

$ERRORIPRED=F               ;individual predictionY=F*+EPS(1)           ;加和型Y=F*(1+EPS(1))        ;比例型Y=F*EXP(EPS(1))       ;指数型Y=F*(1+ERR(1))+ERR(2) ;结合型

简单残差变异模型:考虑了不同分析方法

七、$PRED模块

$PRED模块用于不需要调用PREDPP子程序的模型，位于 $DATA和 $INPUT 之后，代替了 $SUBROUTINES 、 $PK和 $ERROR模块。参数（包括个体间变异）、表达式、残差模型完全自己定义，变量可任意命名。

$PREDA=THETA(1)*EXP(ETA(1))B=THETA(2)EFF=A*CMAX+BY=EFF*(1+EPS(1))
$THETA(1 FIXED)      ;KA(0,5)          ;CL(0,60)         ;V
$OMEGA0.09           ;BSV of CL0.09           ;BSV of V
$SIGMA0.09           ;RV

八、 $THETA、 $OMEGA、 $SIGMA模块

用于设定参数的初值，上下限，以及固定参数。

$THETA(下限，初值，上限)
$THETA(0，45，100)
$THETA(45 FIX)     ;固定$OMEGA
(0.09 FIX)   ;OMEGA1
0.09    ;OMEGA2
0.09    ;OMEGA3$SIGMA
0.09   ;add
0.02   ;eps

九、 $ESTIMATION、 $COVARIANCE模块

$ESTIMATION

METHOD指定估计方法，MAXEVAL指定最大迭代次数，PRINT指定每多少次运算输出一次结果，$SIGDIGITS模块指定有效数字位数。

$EST METHOD=1 INTER NOABORT PRINT=5 MAX=9999 SIG=3

NONMEN 估算方法：

一阶估算法（FO）
一阶条件估算法（FOCE）
含个体间变异-残差变异交互作用的一阶条件估算法（FOCE-I）
拉普拉斯法（LAPL）
最大期望算法（EM），包括：
蒙特卡洛重要抽样法（IMP）
基于后验估计的重要抽样法（IMPMAP）
随机渐进最大期望法（SAEM）
马尔科夫链蒙特卡洛贝叶斯法（MCMC）
迭代两步法（ITS）

$COVARIANCE模块对估计结果进行协方差分析。MATRIX指定估计R矩阵还是S矩阵，PRINT=E指定输出方差-协方差矩阵的特征值。

$COVARIANCE MATRIX=R PRINT=E CONDITIONAL

十、 $TABLE、 $SCATTERPOL模块

$TABLE 以指定的文本形式输出模型参数结果

NOAPPEND指定不输出观测值（DV）、群体预测值（PRED）、残差（RES）、加权残差（WRES），默认输出
NOPRINT指定不在.lst文件里输出数据列表；
FILE指定结果输出的文件名；
ONEHEADER指定列名只出现在第一行。

$SCATTERPOL模块输出绘制散点图的数据

NONMEN 控制文件相关推荐

oracle重建实例_记一次误删Oracle控制文件并恢复过程
概述当你在数据库运行时误删除了控制文件怎么办?很不幸有一次我就有这个情况,虽然是测试环境,这里因为我有事先把控制文件分别备份,所以恢复还是比较简单的.下面简单记录下怎么恢复. 问题控制文件版本不一 ...
ORACLE11g 没有控制文件如何通过rman备份恢复数据的详细实战过程
1.副总裁需要裸恢复的严峻现实集团总部的信息部负责人给我打电话说为了找一年前的记录,所以需要对一年前2015年5月1日的数据进行恢复.而2016年初因为进行迁移,所以有些文件可能丢失,手上只有rma ...
oracle某个表丢失,丢失一个控制文件并恢复数据库
只丢失或损坏一个控制文件的情况下来恢复数据库,相对来说简单一点.一般来说,控制文件都需要形成一个多路径冗余策略,来提高数据库的安全性.这样的话只需将完好的控制文件复制一个副本放到丢失或者损坏了的控 ...
oracle00214,oracle添加控制文件,ORA-00214: 错误
今天学习oracle添加控制文件,控制文件是一个二进制文件,记录了物理数据库的当前状态,每一个数据文件只属于一个数据库,为了防止文件损坏丢失等情况,一个数据库不止一个控制文件.一个数据库的控制文件中的 ...
重建控制文件具体解释
http://note.youdao.com/share/?id=d3991566cad33723db8b92183dcd9161&type=note 完整图文最新版创造"重建 ...
oracle控制文件发生坏块,控制文件坏块处理方法
1 说明控制文件毫无疑问是一个数据库正常工作的必不可少的一个文件.非常重要,所以控制文件出现坏块需要第一时间进行修复.通常会以ORA-00227来表示控制文件出现坏块.https://www.cnd ...
oracle log block size,案例:Oracle无法启动报错ORA-00218: block size 0 重建控制文件
天萃荷净 oracle数据库无法启动报错ORA-00218故障,通过重建控制文件解决文件损坏问题遇到一个案例在数据库启动的时候报ORA-00218错误,而这个故障的引起原因是因为重建控制文件的时候, ...
oracle 控制文件在哪里设置_从Oracle到PostgreSQL：最全控制文件
原文: 从Oracle到PostgreSQL:最全控制文件(上) https://www.enmotech.com/web/detail/1/770/1.html 从Oracle到PostgreSQL ...
Oracle控制文件的重建
场景:存储介质被破坏,造成控制文件丢失,需更换硬盘,并重新建立控制文件. (1)修改初始化参数control_files,改变控制文件位置 SQL> ALTER SYSTEM SET CO ...

NONMEN 控制文件