Hi,

I am trying to analyze data from a matched-pair cluster randomized trial.
The group-level outcomes were analyzed using a random effect meta-analysis
but I now want to investigate individual-level outcomes. I have used glmer
to take into account the school clustering but I am unsure about the correct
way to take into account the matching...

glmer(Outcome ~ Trt + (1|SCHOOL), data=Data, family = binomial)

Is there an option to include the matching as a strata?

(There is a total of 18 schools that were matched on geographic area and
size. In each pair (9), one school was randomized to an intervention. There
is a total of 1159 observations (but cluster size varies from 20 to 180).
The primary outcome is binary.)

Thanks,
Mathieu MG

	[[alternative HTML version deleted]]


