diff --git a/R/format.data/Canada.R b/R/format.data/Canada.R index 909a51823397c4b57e6f672f7f2c62939ac2a4ba..c0c621956ef2e6602bb8f0b8d0492501745e35a6 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 69d301709306636bb43a7a9bd156589bab677933..4deefae321707a88d3496f64637f1926b15ec5e4 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