sum(mites$pa) / nrow(mites)