Remove all trees <10 cm dbh and Canadian plots <0.02 ha

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)
 
+

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