To account for clustering, multivariable 2-level multinomial regression models were constructed to assess for factors associated with colonization of major MRSA clones, using SAS Proc Glimmix.