From d3f3e86b2637451e42f54dccf0866cdf4e99fc29 Mon Sep 17 00:00:00 2001
From: markvanderwel <mark.vanderwel@ufl.edu>
Date: Mon, 16 Sep 2013 22:27:13 -0400
Subject: [PATCH] Remove all trees <10 cm dbh and Canadian plots <0.02 ha
---
R/format.data/Canada.R | 6 ++++++
R/format.data/US.R | 6 ++++++
2 files changed, 12 insertions(+)
diff --git a/R/format.data/Canada.R b/R/format.data/Canada.R
index 909a518..c0c6219 100644
--- a/R/format.data/Canada.R
+++ b/R/format.data/Canada.R
@@ -10,6 +10,11 @@ dir.create("output/formatted/Canada", recursive=TRUE,showWarnings=FALSE)
######################### READ DATA read individuals tree data
data.canada <- read.csv("data/raw/Canada/Canada_Data2George_20130910.csv",
header = TRUE, stringsAsFactors = FALSE)
+
+#MCV remove plots smaller than 0.02 ha (plot size with average of <20 trees 10cm dbh)
+#MCV and remove trees <10 cm dbh
+data.canada <- subset(data.canada,subset=(data.canada$SubPlotSize>=0.02 & data.canada$InitDBH>=10) )
+
sum(is.na(data.canada$species))
data.canada$sp = data.canada$species_FIACode
data.canada$sp[data.canada.sp==8349] <- 762 #MCV replace incorrect sp code for Prunus serotina
@@ -80,3 +85,4 @@ vec.basic.var <- c("obs.id","tree.id", "sp", "sp.name","cluster","plot", "ecocod
data.tree <- subset(data.canada, select = c(vec.basic.var, vec.abio.var.names))
write.csv(data.tree,file="output/formatted/Canada/tree.csv",row.names = FALSE)
+
diff --git a/R/format.data/US.R b/R/format.data/US.R
index 69d3017..4deefae 100644
--- a/R/format.data/US.R
+++ b/R/format.data/US.R
@@ -27,6 +27,12 @@ write.csv(data.max.height, file = "output/formatted/US/max.height.csv") ## I wa
#### READ DATA US
data.us <- read.csv("data/raw/US/FIA51_trees_w_supp.csv",header=TRUE,stringsAsFactors =FALSE)
+
+#MCV remove trees <10 cm dbh
+data.us <- subset(data.us,subset=(data.us$InitDbh>=10) )
+
+
+
############## MASSAGE TRAIT DATA HEIGHT DATA FOR TREE MISSING BRING US DATA FOR HEIGHT OVER
############## WHEN WE ANALYZE THAT DATASET LATER ON
############## FORMAT INDIVIDUAL TREE DATA
--
GitLab