[R-sig-ME] lme4 categorical outcome probability values at each time point

[Igor Yakovenko]

Hello Everyone,

 

Does anyone know how to graph the results of a lme4 longitudinal multilevel
model as a bar graph of Time (3 points) on x-axis with probability of
getting the 1 value on the outcome after taking into account the predictors
as the y-axis? Alternatively, if there is a way to derive the mean
probability of the outcome being 1 (vs. 0) during an intervention at time 1,
2 and 3 divided by Condition (two levels), I can graph those values
separately. The model is based on data from a two-group RCT evaluating an
intervention over 3 time points. Here is a sample model with a dichotomous
outcome and fixed effects of time (continuous), condition (2 levels) and
whether an additional workshop was completed (2 levels).

 

workshop<-glmer(CurrentAsst~Time*Condition + Workshop_completion_reordered +
(1|id), data = vse, family = binomial, nAGQ = 25, control =
glmerControl(optimizer = "bobyqa"))

 

Thank you,

Igor Yakovenko


	[[alternative HTML version deleted]]