> setwd("C:/Users/W.H.Y/Desktop")
> library(gemtc)
载入需要的程辑包:coda
> library(rjags)
Linked to JAGS 4.3.1
Loaded modules: basemod,bugs
> data <- read.csv("test1.csv",sep=",",header = TRUE)
> treatments <- read.csv("treatments1.csv",sep=",",header = TRUE)
> network <- mtc.network(data, description="Example",treatments=treatments)
> plot(network)
> model <- mtc.model(network,type = "consistency",factor = 2.5,n.chain = 4,linearModel = "random")
>results <- mtc.run(model,sampler = NA,n.adapt = 5000,n.iter = 20000,thin = 1)
forest(relative.effect(results,"E"))
参照链接https://zhuanlan.zhihu.com/p/467276292 ,把里面的参数和它设的一样试一下。尤其是network <- mtc.network和mtc.run这两个函数。如果这两个试了还不行,那就对换个和他一样的,如果再不行的话可以留言。
这不好搞