贝叶斯网络元分析 (gemtc)

Question

我正在尝试使用 R 包 gemtc 进行 NMA。我有三个比较： A v D: MD: 119 (17.6) B v D：MD：73（23.98） C v D：MD：92（21.94）

我正在关注一个似乎可以正常工作的示例代码：

library(gemtc)
library(rjags)
library(dmetar)

data(TherapyFormatsGeMTC)

head(TherapyFormatsGeMTC$data)

TherapyFormatsGeMTC$treat.codes
network <- mtc.network(data.re  = TherapyFormatsGeMTC$data,
                       treatments = TherapyFormatsGeMTC$treat.codes)


model <- mtc.model(network,
                   likelihood = "normal",
                   link = "identity",
                   linearModel = "random",
                   n.chain = 4)

但是当我尝试使用我的代码时：

library(gemtc)
library(rjags)
library(dmetar)

treatments <- read.table(textConnection('
  id  description
  A   "Treatment A"
  B   "Treatment B"
  C   "Treatment C"
  D   "Treatment D"'), header=TRUE)
data <- read.table(textConnection('
  study                 diff         std.err            treatment  
  "Study A"             119          17.60              A          
  "Study A"             NA           NA                 D          
  "Study B"             73           23.98              B          
  "Study B"             NA           NA                 D          
  "Study C"            92           21.94              C          
  "Study C"            NA           NA                 D'), header=TRUE)

data$diff<-as.numeric(data$diff)

network1 <- mtc.network(data, description="Example", treatments=treatments)

model <- mtc.model(network1,
                   likelihood = "normal",
                   link = "identity",
                   linearModel = "random",
                   n.chain = 4)

我收到以下错误：

Error in validate.data.normal.identity(list(study = c(1L, 1L, 2L, 2L,  : 
  all(data.ab[["std.err"]] > 0) is not TRUE

我无法弄清楚我做错了什么

Answer 1

也许您可以尝试使用data.re 参数而不是默认的data.ab：

network1 <- mtc.network(data.re=data, description="Example", treatments=treatments)