####################################################################################
# Read data
data(Smoking)
if (FALSE) {
dat <- Smoking
####################################################################################
# Table TXBR
# matrix producing 4x2x3x6 table TXBR
atTXBR <- MarginalMatrix(var = c("X", "B", "R1", "R2", "R3", "R4"),
marg = list(c("X", "B", "R1"), c("X", "B", "R2"), c("X", "B", "R3"), c("X", "B", "R4")),
dim = c(2, 3, 5, 5, 5, 5))
bt <- ConstraintMatrix(var = c("T", "X", "B", "R"), suffconfigs = list(c("T", "X", "B"), c("R")),
dim = c(4, 2, 3, 5))
model = list(bt, "log", atTXBR)
fit = MarginalModelFit(dat = dat, model = model, MaxStepSize = .3, MaxSteps = 100,
ShowProgress = 5)
}
Run the code above in your browser using DataCamp Workspace