网状Meta分析在JAGS软件中的实现.pdfVIP

数理医药学杂志 2014年第 27卷第 4期 or[c，k]~--exp(d[k]--d[c]) 2 NMA 的实现步骤 lor[c，k]一(d[k]一d[c]) JAGS软件虽是基于BUGS语言而设定的，但其操作需要 } 全程编程，且在模块、数据及初始值的格式及输入方面与 } BUGS软件有所不同。因此，为保证正常运行操作，需要对模 处理完毕后 ，将其保存为：“．bug”或 “．jags”格式。本例以 块、数据及初始值的格式在其加载之前进行相关处理。 “ networkmodel”为名、以 “．jags”格式保存在桌面的Rwork文 2．1 模块加载 件夹下。至此，就完成了该软件的模块处理部分。 与BUGS软件 相比，JAGS软件 自身没有排列函数。其 在加载模块之前 ，为了方便操作，可以改变其原始工作 目 模块加载的代码命令为： 录为桌面的Rwork文件夹，具体命令如下： model cd“C：＼Users＼Administrator＼Desktop＼Rwork” ( 至此，所有的操作将默认在桌面的Rwork文件夹中，并可 for(iin1：ns) 以通过下述命令加载模块 ： { modelinnetworkmode1．jags w[i，1]一0 上述命令输入后 ，按回车键没有出现任何提示时，则说明 delta-[i，1]一O 模块已经加载完成。 mu[i]～dnorm(o，0．0001) 2．2 数据加载 for(kin1：na[i-]) 加载之前，需要进行数据的排列。JAGS软件中数据的排 { 列格式与R软件数据排列格式一致 ，本例如下： r[i，k]~dbin(p[i，k]，nil，k-t) t 一 structure(c(1，1，1，1，1，1，1，1，1，1，1，1，1， logit(p[i，k1)~--mu[i-l-t-delta[i，k] 1，1，1，l，1，1，l，1，1，1，1，1，1，1，1，1，2，2，2，3， rhat[i，k]--p[i，k]*nil，k] 3，5，5，5，6，6，6，6，7，7，7，7，7，7，7，7，7，7，7，7， dev[i，k]一2*(r[i，k]*(1og(r[i，k])一log(rhat[i， 7，7，9，9，1O，10，11，12，12，12，12，2，2，2，4，4，4，4， k-I))+(nl-i，k]一r[i，k])*(1og(n-[i，k]一r[i，k])--log(nil，k] 4，5，5，5，5，5，5，5，7，7，7，7，8，1O，10，l1，12，12， --rhat-[i。k-I)))