Project examining community dynamics on Portal krat exclosure plots. Based off original project which focused on control plots only: https://github.com/emchristensen/Extreme-events-LDA This work expands on the original project by asking whether the rodent community present on the experimentally altered plots (altered by removing krats) experienced changes at the same time as the rodent community on the control plots.
- compare_changepoint_distributions.R code for analyzing overlap of changepoint distributions from control plots and krat exclosures
- main_results_figures.R code for making Fig 1 and 2 in ms (multipart LDA and changepoint results)
- run-hg.R code for running analyses on hipergator
- whole-analysis.R main script for analyzing rodent community change using LDA. Pulls data from PortalData repo, manipulates it for input to LDA model, runs LDA models and model selection, runs changepoint model and model selection
- paper_dat.csv Portal rodent data that was used in the original Extreme-events-LDA project
- baileys-population.R plot population time series of C. baileyi on control and krat exclosure plots (Fig 4 in ms)
- species-populations.R plot popultion time series of all species around 1990 changepoint (supplemental fig in ms)
- total_abundance_plots.r plot total abundance time series on control and krat exclosure plots (not used in ms, but was used in Extreme-events-LDA project)
Contains scripts used by whole-analysis.R
- changepoint_histogram_plot.R function for plotting changepoint histograms (used by main_results_figures.R)
- convert_summary_to_dates.R function to convert output of changepoint model (in time steps) to dates for interpretation and plotting
- eval_changepoint_model.R functions for evaluating candidate changepoint models and doing model selection
- get-data.R function "get_rodent_lda_data" for extracting data from PortalData repo, using portalr package, and preparing it for input to LDA model
- lda_plot_function.R functions for visualizing output of LDA model: barplots of species composition of topics (from previous work in Extreme-events-LDA)
- plots_from_ldats.R functions for visualizing output of LDA model: timeseries (from LDATS package)
- run_changepoint_model.R function for running changepoint model, potentially with multiple numbers of changepoints
- select_LDA.R function "run_rodent_LDA" to run a bunch of LDA models and do model selection. Returns only the best LDA model.
Contains .Rdata results files from runs done on hipergator
Contains code and data from Extreme-events-LDA repo (may or may not be completely up to date).
Rmarkdown files and pdfs of results for facilitating discussion