Introduction
This vignette will walk you through the different functions for cleaning quadrat data within the ‘quadcleanR’ package.
library(quadcleanR)
data("corals")
Cleaning Functions
There are 5 cleaning functions within quadcleanR.
1. change_names()
Using a new data frame of labels, change column names in one function. Helpful if column names are shorthands or contain spaces and characters that are not supported in column names in R.
data("coral_labelset")
corals_change_names <- change_names(corals, coral_labelset, from = "short_name", to = "full_name")
2. change_values()
Using two vectors, change the values in one column to a new set of values. Helpful if you need to change many values at once, like updating changes to site names or taxonomy.
corals_change_values <- change_values(corals_change_names, "Field.Season", c("KI2015a", "KI2016a"), c("KI2015", "KI2016"))
3. keep_rm()
Using a character, or part of character select rows or columns of the data frame to either keep or remove. A more customizable way to subset your data as you can keep or remove based on partial matches, or cells containing select characters.
corals_keep_rm <- keep_rm(corals_change_values, "DEEP" , select = "row", keep = FALSE, drop_levels = TRUE, exact = FALSE, "Site")
4. rm_chr()
Parts of characters can be removed based on a vector of removal characters. When working with images, this can be helpful to remove extra characters from image IDs, or anywhere else where you want to remove specific characters from your data.
Adding Functions
There are 2 adding functions within quadcleanR.
1. add_data()
Using key identifying columns, add additional columns to an existing data frame. Helpful for adding environmental or taxonomic data to your quadrat data
2. categorize()
Using a column within the data frame, categorize rows in a binary of yes or no, or customize with a set of category names. Especially useful for turning ID tags into useful categories for analysis such as morphology, taxonomy etc.
corals_categorize <- categorize(corals_add_data, "Field.Season", values = c("2015", "2016", "2017"), name = "ElNino", binary = FALSE, exact = FALSE, categories = c("Before", "During", "After"))
Processing Functions
There are 3 processing functions within quadcleanR.
1. cover_calc()
Convert the number of observations for each species or non-species to proportion or percent cover within each row based on the total number of observations in each row. Useful for quadrats with varying numbers of observations to calculate each row’s percent cover all at once.
corals_cover_calc <- cover_calc(corals_categorize, spp = colnames(corals_categorize[,c(7:15)]), prop = TRUE, total = TRUE)
2. crop_area()
Using the location of annotated points within quadrats and the size of the quadrat, crop quadrat data to a smaller area, while maintaining the spatial relationships between points. Useful for making different sized quadrat data comparable.
This function is covered in depth with the Why to Crop Quadrats by Area vignette.
3. usable_obs()
Sum columns containing unusable observations and remove rows that contain more than the specified cutoff number of unusable points. Helpful if there are annotations that were unidentifiable and you want to remove them from the total usable observations, and you can remove quadrats with too many unusable observations.
corals_usable_obs <- usable_obs(corals_cover_calc, c("Unclear", "Shadow"), max = TRUE, cutoff = 0.9, above_cutoff = FALSE, rm_unusable = TRUE)
Assessing Functions
There are 2 assessing functions within quadcleanR.
1. sample_size()
Specify which columns to use to produce a table with sample sizes. Helpful to visualize number of samples in your data.
corals_sample_size <- sample_size(corals_usable_obs, dim_1 = "Site", dim_2 = "Field.Season", count = "Quadrat")
2. visualize_app()
Using an interactive shiny app, visualize and explore cleaned quadrat data.
A good combination to examine this shiny with is : - y-axis: Sarcophyton - x-axis: Field.Season - color: ElNino (treat as discrete) - facet: HD_cat - group by: Field.Season, ElNino, Site and HD_Cat - view as a violin plot
corals_usable_obs$ElNino <- factor(corals_usable_obs$ElNino, levels = c("Before", "During", "After"))
corals_usable_obs$HD_Cat <- factor(corals_usable_obs$HD_Cat, levels = c("Very High", "High", "Medium", "Low", "Very Low"))
visualize_app(data = corals_usable_obs, xaxis = colnames(corals_usable_obs[,c(1:6)]), yaxis = colnames(corals_usable_obs[,c(7:13)]))