R Script for the analysis of relationships between PhoGV and Potato tuber moth distributions

Bioclimatic predictor raster stacks are loaded, glm models for densities of insests and virus prevalence are fitted, and correlations between these predicted variables at the sampling points are performed. Finaly anova is performed to evaluate the hierarchical effects of site, storage bags and tubercule on virus infection. Libraries MASS, raster, mapplots, rgdal and matrixStats are used. (R file)