We used this location category as a predictor in a generalized linear model with the binomial distribution using the glm function and compared it to a null model with a likelihood ratio test using lrtest from the lmtest (Zeileis and Hothorn, 2002) package to test for significant differences in the frequency of first and last seasonal presence depending upon station location.