devel/make_paramcorr_tables.R

In atredennick/community_synchrony: Calculates syunchrony and stability metrics from time series data

##  Script to create tables for parameter 
##  correlations from vital rate regressions

rm(list=ls())

library("xtable")

in_dir <- "../results/param_covariance/"
tmp_files <- list.files(in_dir)
all_files <- tmp_files[grep(".RDS", tmp_files)]

if(grep("surv_grow_paramcorrs", list.files("../results/")) != 0){
  file.remove("../results/surv_grow_paramcorrs.txt")
}
newpages <- seq(3,length(all_files), by=3)
for(i in 1:length(all_files)){
  tmp_corrs <- readRDS(paste0(in_dir,all_files[i]))
  vital <- strsplit(all_files[i], "_")[[1]][1]
  site <- strsplit(all_files[i], "_")[[1]][2]
  species <- strsplit(all_files[i], "_")[[1]][3]
  cap <- paste0("Parameter correlations for site = ", site, ", vital rate = ",vital, ", and species = ",species,".")
  cat(print(xtable(tmp_corrs, digits=-2, caption = cap), caption.placement = "top", comment=FALSE), file = "../results/surv_grow_paramcorrs.txt", append = TRUE, sep="\n")
  cat("\n\n", file = "../results/surv_grow_paramcorrs.txt", append = TRUE)
  if(i %in% newpages){cat("\\newpage{} \n\n", file = "../results/surv_grow_paramcorrs.txt", append = TRUE)}
}