library("gemtc")
library("rjags")
library("igraph")
library("coda")
setwd("D:\Rlanguage\wangzhuangtu")
data <- read.csv("data.csv", sep=",", header=TRUE)
treatments <- read.csv("treatments.csv", sep=",", header=TRUE)
network <- mtc.network(data=data,description = "Network data",treatments = treatments)
summary(network)
plot(network,use.description=FALSE,
vertex.label.cex=1.5,
vertex.size=20,
vertex.label.color="black",
vertex.label.dist=2,
vertex.label.degree=-pi/3,
vertex.shape="circle",
vertex.color="royalblue",
dynamic.edge.width = TRUE,
edge.color="black",
vertex.label.font=2)
model <- mtc.model(network,type="consistency",n.chain=4,likelihood="binom",link="logit",linearModel="fixed")
最后一个model命令一直报错如何改