added more covariates,
removed some sites with missing data,
updated model statement