I am working on a dataset of longitudinal (baseline and control) data about some values in blood in two treatment groups
I performed a propensity score matching with the matchthem methodology in R (first imputing the data with mice
function, and then matching with the matchthem
package). After that the package suggest using svyglm
for computing the tests on a dichotomic variable as outcome. However I am trying to compute a linear mixed model, adjusted by id as random variable and time and treatment in this matched dataset, but I can't find the way of doing it.
Any ideas?