Commit 2224624c authored by Daniel Laughlin's avatar Daniel Laughlin
Browse files

Minor tweaks to New Zealand scripts

No related merge requests found
Showing with 5 additions and 3 deletions
+5 -3
......@@ -27,7 +27,7 @@ data.trait$SLA.mean <- 1/data.trait$lma.gm2; data.trait$SLA.mean <- data.trait$S
data.trait$SLA.sd <- NA
data.trait$Wood.density.mean <- data.trait$wood; data.trait$wood <- NULL
data.trait$Wood.density.sd <- NA
data.trait$Max.height.mean <- log10(data.trait$height.m); data.trait$height.m <- NULL
data.trait$Max.height.mean <- data.trait$height.m; data.trait$height.m <- NULL
data.trait$Max.height.sd <- NA
## read traits from TRY
......
......@@ -67,8 +67,8 @@ rm(data.sp, data.sp2)
# library(RColorBrewer)
# library(rworldmap)
# newmap <- getMap(resolution = 'coarse')
# ## # different resolutions available
# plot(newmap,xlim=c(166,177),ylim=c(-34,-47))
# # ## # different resolutions available
# plot(newmap,xlim=c(166,177),ylim=c(-47,-34))
# ecoreg <- factor(data.nz$ecocode); levels(ecoreg) <- mycols <- brewer.pal(length(levels(ecoreg)), "Set1")
# points(data.nz[['Lon']],data.nz[['Lat']],pty='.',cex=.2,col = as.character(ecoreg))
# legend('bottomleft', col = mycols, legend = names(table(data.nz$ecocode)), pch = rep(19,length(levels(ecoreg))))
......@@ -87,4 +87,6 @@ vec.abio.var.names <- c("MAT", "MAP")
vec.basic.var <- c("obs.id","tree.id", "sp", "sp.name","cluster","plot", "ecocode", "D", "G", "dead",
"year", "htot", "Lon", "Lat", "perc.dead","weights","census")
data.tree <- subset(data.nz, select = c(vec.basic.var, vec.abio.var.names))
write.csv(data.tree,file="output/formatted/NVS/tree.csv",row.names = FALSE)
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment