To account for clustering at the household level , we analysed all injuries from falls at home per person-year with a negative binomial generalised linear model with generalised estimating equations .