Data were analyzed using a mixed-effects multilevel Poisson regression model ( primary outcome ) and a Cox proportional hazards model ( secondary outcome ) , adjusting for facility-level clustering and resident-level variables .