An official website of the United States government.

This is not the current EPA website. To navigate to the current EPA website, please go to www.epa.gov. This website is historical material reflecting the EPA website as it existed on January 19, 2021. This website is no longer updated and links to external websites and some internal pages may not work. More information »

CADDIS Volume 4

R Script: Compare Taxa Names in Tolerance Value and Assessment Data

# Compare taxa names in tolerance value and assessment data.

names.tv <- toupper(names(site.species)[-1])
names.assess <- toupper(names(site.species.or)[-1])

names.all <- c(names.tv, names.assess)
names.match <- names.all[duplicated(names.all)]

names.miss.tv <- character(0)
names.miss.assess <- character(0)

for (i in 1:length(names.tv)) {
     if (is.na(match(names.tv[i], names.match))) {
          names.miss.tv <- c(names.miss.tv, names.tv[i])
     }
}

for (i in 1:length(names.assess)) {
     if (is.na(match(names.assess[i], names.match))) {
          names.miss.assess <- c(names.miss.assess, names.assess[i])
     }
}

print("Taxa in both databases")
print(sort(names.match))

print("Taxa in tolerance value db but not observed in assessment db")
print(sort(names.miss.tv))

print("Taxa in assessment db but not observed in tolerance value db")
print(sort(names.miss.assess))