Ecology in R

Lesson 2 Overview

•Occurrence Data

•Packages: rgbif, mapview, scrubr, sp, dplyr

•Data mining species occurrence data from GBIF

•Creating spatial points data objects

•Exploratory visualization techniques

•Species Density Maps

•Packages: sp, raster, usdm, mapview, rgbif, scrubr, GISTools, maps, ggplot2, RColorBrewer, ggspatial

•Spatial data cleaning techniques

•Querying country polygon outlines

•Clustering techniques

•Basic species density and spatial distance analyses

•Multiple data visualization techniques

A walk-through of Lesson 2, exercise 1

Density and distance analysis using species occurrence data queried from GBIF, various mapping techniques, and simple statistical analyses.

Lesson 3 Overview

•Environmental Data

•Packages: raster, mapview, ggplot2, dplyr, rgbif, maptools, scrubr

•Querying environmental data (elevation, temperature) from various online sources

•Forming datasets with species occurrence records and extracted environmental data

•Raster manipulation techniques

•Field Guide Maps

•Packages: rnaturalearth, rnaturalearthdata, ggspatial, sf, dismo, jsonlite, mapdata, raster, mapview, ggplot2, dplyr, rgbif, maptools, rgdal

•Use environmental data to project current climate suitability models

•Raster to polygon manipulation for range map generation

•Species Distribution Models (SDM)

•Packages: dismo, jsonlite, mapdata, raster, mapview, ggplot2, dplyr, rgbif, maptools, rgdal

•Using climate data to create current and future projections of climate suitability models

Mining environmental data using the raster package, basic raster manipulation techniques, extracting data from species occurrences, and building species data frames.

Basic techniques to create species distribution maps  based on climate suitability

Correction: In the video I say the highest AUC value indicates the best fit model. I accidentally misspoke in this regard. The best fit model is indicated by the LOWEST AUC value. I must have been still thinking of AIC at the time, in which the best fit model is the highest value. Sorry for the confusion everyone!

A basic introduction to MaxEnt modelling in R using current and future climate data from WorldClim and species occurence data from GBIF.

Note: The models generated in this exercise and the last exercise are simply to show you the basic functionality of MaxEnt in R. SDMs should incorporate multiple parameters and as much data as possible to make robust models when used for scientific publication.

Lesson 4 Overview

•Movement Dynamics

•Packages: adehabitatHR, sp, ggplot2, dplyr, lubridate, mapview

•Analysis of time series, tracking data

•Tabulating movement summaries and saving to .csv

•Home Ranges

•Packages: move, adehabitatHR, caTools, spatialEco, reshape2, tibble, sp, ggplot2, dplyr, lubridate, mapview, cowplot, ggspatial

•Testing several utilization density (UD) estimators

•Visual and mathematical methods of examining Type I (overfitting) and Type II (under-fitting) errors in UD estimators

A tutorial on basic exploratory analyses on movement data

An in-depth lecture on movement home range models and the limitations of Utilization Density (UD) estimators regarding Type I (over-fitting) and Type II (under-fitting) errors.

Lesson 5 Overview

•Trapping Analysis

•Packages: rgbif, mapview, sp, dplyr, ggplot2, Viridis, scatterpie, raster, reshape2, tibble, Vegan, BiodiversityR, ggspatial

•Query trapping study datasets from GBIF

•Data wrangling techniques to transpose data columns to long-form

•Visualizations of species richness

•Map richness based pie-charts for each study site

•Analysis of species richness, density, and abundance

•Species richness indices

•Rarefaction curves, and rank abundance plots

A simple walk-through on querying datasets from GBIF, exploring trapping data, and using data visualizations and statistical indexes to analyze the data.

Note: When discussing the best-fit model of the radlattice() function call on line 274, I incorrectly state the NULL model would be the best-fit model in this case, therefore disqualifying the others. The best-fit model is in fact the Mandelbrot test with an AIC of 350.17.

Sorry for my mistake and any confusion it may have caused!

Lesson 6 Overview

•Trait Analysis

•Packages: ape, taxise, rentrez, phytools, Select, treeio, ggtree, data.tree, tidytree, ggplot2, dplyr, traits, stringr, cowplot

•Create phylogenetic trees using a species list and taxonomic data from NCBI

•Examine various branch distance simulation options

•Combine trait data with tree data

•Subset and visualize trait data in trees for comparative analysis

•Sequence Analysis

•Packages: sequin, adegenet, ape, ggtree, DECIPHER, viridis, ggplot2

•Query raw Internal Transcribed Spacer (ITS) sequences from NCBI

•Read fasta data into R

•Form alignments and distance matrices

•Visualize tree data using baseR

•Modify and customize tree plots in ggtree

An end-to-end tutorial on gathering sequence data from NCBI in fasta format, importing/exporting, aligning, and building customized trees using ggtree.

Basic guide to tree building and including trait data in R for comparative analysis