The data was analyzed using
generalized linear models and
generalized estimating equations, which are specifically used to address the
multilevel design of data in which schools with participating schoolchildren were randomized (rather than individual participants).