data/hepatitis.R

In pmlr: Penalized Multinomial Logistic Regression

"hepatitis" <-
structure(list(
city = structure(c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 1, 1, 1, 2, 2, 2, 3, 3, 1, 1, 1, 1, 2, 3, 3, 3), class = "factor", .Label = c("Toronto", "Hamilton", "Winnipeg")),
group = structure(c(1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 1, 2, 2), class = "factor", .Label = c("no-withhold", "withhold")), 
time = structure(c(1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 1, 2, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2), class = "factor", .Label = c("pre", "post")), 
HCV = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0),
nonABC = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1),
counts = c(144, 960, 160, 965, 172, 629, 160, 626, 73, 275, 80, 305, 2, 1, 2, 2, 2, 1, 1, 2, 3, 10, 1, 6, 3, 1, 1, 1)),
.Names = c("city", "group", "time", "HCV", "nonABC", "counts"),
row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28"),
class = "data.frame")