# NOT RUN {
## --- Continuing the Example from '?sim.mp' and '?gtreg.mp':
# 5*6 and 4*5 matrix
set.seed(9128)
sa1a<-sim.mp(par=c(-7,0.1), n.row=c(5,4), n.col=c(6,5),
sens=0.95, spec=0.95)
sa1<-sa1a$dframe
# }
# NOT RUN {
fit1mp <- gtreg.mp(formula = cbind(col.resp, row.resp) ~ x, data = sa1,
coln = coln, rown = rown, arrayn = arrayn,
sens = 0.95, spec = 0.95, linkf = "logit", n.gibbs = 1000, tol = 0.005)
summary(fit1mp)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab