Newer
Older
#################
################
### PROCESS THE DATA FOR NVS
rm(list = ls())
source("process.fun.R")
library(reshape)
data.tree <- read.csv("../../output/formatted/NVS/tree.csv", stringsAsFactors = FALSE)
data.tree <- subset(data.tree,subset=!is.na(data.tree[["D"]]))
## change sp
data.traits <- read.csv("../../output/formatted/NVS/traits.csv")
data.traits$sp <- paste("sp",data.traits$sp,sep=".")
data.tree$sp <- paste("sp",data.tree$sp,sep=".")
# vector of ecoregion name
ecoregion.unique <- unique(data.tree[["ecocode"]])
#### lapply function
system.time(lapply(ecoregion.unique, FUN = fun.data.per.ecoregion, data.tot = data.tree,
weight.full.plot = NA, site.name = "NVS", data.TRAITS = data.traits))