data(hivsurv)
fit1 <- gtReg(type = "sp",
formula = groupres ~ AGE + EDUC.,
data = hivsurv, groupn = gnum,
sens = 0.9, spec = 0.9,
method = "Xie")
summary(fit1)
# 5x6 and 4x5 array
set.seed(9128)
sa2a <- gtSim(type = "array", par = c(-7, 0.1),
size1 = c(5, 4), size2 = c(6, 5),
sens = 0.95, spec = 0.95)
sa2 <- sa2a$dframe
# \donttest{
fit2 <- gtReg(type = "array",
formula = cbind(col.resp, row.resp) ~ x,
data = sa2, coln = coln, rown = rown,
arrayn = arrayn, sens = 0.95, spec = 0.95,
linkf = "logit", n.gibbs = 1000, tol = 0.005)
summary(fit2)# }
Run the code above in your browser using DataLab