method.R

<br> ### function for calculating growth synchrony ### fun_syn = function(chron, interval){ syn = data.frame() for(i in 1900:(2014 - interval + 1)){ res = filter(chron, year %in% c(i : (i + interval - 1))) res = res[, colSums(is.na(res)) == 0] %&gt;% select(-year) %&gt;% cor() n = ncol(res) res = res[lower.tri(res)] res = data.frame(year = i + interval - 1, mean = mean(abs(res)), sd = std.error(abs(res)), n = n) syn = rbind(syn, res)} return(syn)} <br> ###### within-population growth syn calculation ###### wps = list() res1 = data.frame() res2 = data.frame()