data(schizo)
model.10.3.1 <-
olmm(imps79o ~ tx + sqrt(week) + re(1|id),
data = schizo, family = cumulative())
summary(model.10.3.1)
(newCoef <- coef(model.10.3.1))
newCoef[1:3] <- c(-7, -4, -2) # thresholds
newCoef["tx"] <- 1.5
newCoef["sqrt(week)"] <- 1.6
newCoef[6] <- 2
model.10.3.1.mod <- olmm_overwritePar(object = model.10.3.1, coefficients = newCoef)
summary(model.10.3.1.mod)
Run the code above in your browser using DataLab