###################### #Allelic richness ###################### AR=c(5.89,5.83,5.76,5.85,5.87,5.7,5.74,5.69,5.45,5.66,4.91,5.28,5.41,5.5) Latitude=c(34.49,36.13,36.22,36.34,36.41,37.4,41.18,41.11,41.04,40.36,45.35,47.1,46.55,43.42) res=lm(AR~Latitude) plot(Latitude, AR, xlab="Latitude", ylab="Allelic Richness”, pch=16) res abline(res) summary(res) ############################### #Private Allelic Richness ############################### PAR=c(0.27,0.22,0.19,0.21,0.22,0.2,0.23,0.18,0.14,0.13,0.08,0.11,0.13,0.12) plot(Latitude, PAR, xlab="Latitude", ylab=“Private Allelic Richness”, pch=16) res=lm(PAR~Lat) res abline(res) summary(res)