## Illustrating the construction of a "fixed" local odds ratios structure
## using the arthritis dataset. Here, we assume a uniform local odds ratios
## structure equal to 2 for each time pair.
## Create the uniform local odds ratios structure.
lorterm <- matrixLOR(matrix(2, 4, 4))
## Create the LORterm argument.
lorterm <- c(lorterm)
lorterm <- matrix(c(lorterm), 3, 25, TRUE)
## Fit the marginal model.
data(arthritis)
fitmod <- ordLORgee(y ~ factor(trt) + factor(time) + factor(baseline),
data = arthritis, id = id, repeated = time, LORstr = "fixed",
LORterm = lorterm)
fitmod
Run the code above in your browser using DataLab