location/treatment/species1
BC/yes/20
BC/no/11
BC/yes/880
DW/no/3
DW/yes/230
DW/yes/0
YZ/yes/390
YZ/no/0
YZ/yes/540
(etc...)

sp1calls=glm(species1~treatment,family=poisson)

sp1callswithblock=glm(sp1~treatment+block,family=poisson)

lmersp1loc=lmer(sp1~treatment + (1|location/treatment), data = mydataframe)