The script (habitat_density_model.r) in this code repository filters data, sets up model input, and runs an N-mixture model with two components of detection (availability and perceptibility) to estimate bird density from point-count data. This analysis was developed to evaluate covariates of population density for three species of boreal birds: Lesser Yellowlegs, Olive-sided Flycatcher, and Rusty Blackbird.