Dear colleagues,
I'm analysing a randomized controlled trial in Nepal. I built a linear multilevel model with a score as the outcome, time as a dummy variable (i.time) and the interaction term i.intervention#i.time and id as second level in Stata:
xtmixed outcome i.time i.intervention#i.time || id:, var
According to Twisk 2018, excluding the treatment variable intervention should adjust for baseline score ( 10.1016/j.conctc.2018.03.008 ).
I then used the fitted predicted values
predict predri, fitted level(id)
to make a graph of the predicted scores over time. Predicted scores at baseline are different between control and intervention group, even though if they were adjusted for baseline they should be the same.
Am I doing something wrong?
Many thanks for your help!
Sebastián