knitr::opts_chunk$set(echo = TRUE)
This Markdown file provides the documentation for updating the Nutrient Biotic Index (NBI) (See Smith et al. 2007) to include:
The methods used rely on the use of existing nutrient bins, weighted average equation, and tolerance value assignment ranges as published in Smith et al. 2007.
References:
Smith, A.J., R.W. Bode, and G.S. Kleppel. 2007. A nutrient biotic index (NBI) for use with benthic macroinvertebrate communities. Ecological Indicators 7:371-386
Import in the raw data files from the working directory. These include a species data table, water chemistry data table, and a site data table.
species<-read.csv("Species_2002_2016.csv") chemistry<-read.csv("Chemistry_2002_2016.csv") chembiojoin<-read.csv("RIBSchem_BioSiteJoin_Final_2019-01-09.csv")
Create two new tables for species and chemistry data that add the column "RIBS_ID" to the species table and "BASIN_LOC_RIV" to the chemistry table. By providing this set of joins and adding the additional columns from matching records in each table, we have the necessary data to compare nutrient values to macroinvertebrate species.
joinedspecies <- merge(x = species, y = chembiojoin, by = "BASIN_LOC_RIV", all = FALSE) joinedchemistry <- merge(x = chemistry, y = chembiojoin, by = "RIBS_ID", all = FALSE) summary(joinedspecies) summary(joinedchemistry)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.