From 2224624c89892bd469a3c3baa0a1aa5be882cd6f Mon Sep 17 00:00:00 2001 From: Daniel Laughlin <d.laughlin@waikato.ac.nz> Date: Tue, 17 Sep 2013 16:56:12 +1200 Subject: [PATCH] Minor tweaks to New Zealand scripts --- R/find.trait/NVS.R | 2 +- R/format.data/NVS.R | 6 ++++-- 2 files changed, 5 insertions(+), 3 deletions(-) diff --git a/R/find.trait/NVS.R b/R/find.trait/NVS.R index 7336f4b..eebacb1 100644 --- a/R/find.trait/NVS.R +++ b/R/find.trait/NVS.R @@ -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 diff --git a/R/format.data/NVS.R b/R/format.data/NVS.R index 7f09974..e9d47dd 100644 --- a/R/format.data/NVS.R +++ b/R/format.data/NVS.R @@ -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) -- GitLab