ecocode output for TP only 5

R/analysis/lmer.output-fun.R

@@ -1586,7 +1586,8 @@ fun.plot.wd.sl.param <-  function(list.res,
                                   dir.root = '.',
                                   data.type = 'simple',
                                   big.m = 1.0,
-                                  small.m = 0.42
+                                  small.m = 0.42,
+                                  MAT.MAP.TF = FALSE
                                   ){
 
 
@@ -1596,7 +1597,8 @@ fun.plot.wd.sl.param <-  function(list.res,
                                 model = model,
                                 dir.root = dir.root,
                                 list.res = list.res,
-                                data.type = data.type)
+                                data.type = data.type,
+                                MAT.MAP.TF = MAT.MAP.TF)
 require(dplyr)
 # Layout
  m <- matrix(c(1:4), 2, 2)

scripts/analysis/lmer.all.no.log.output.R

@@ -8,7 +8,7 @@ source("R/analysis/lmer.run.R")
 
 format.all.output.lmer(file.name = "NA.koppen.results.nolog.all.rds",
                    list.file.name = 'list.lmer.out.all.NA.simple.ecocode.koppen.rds',
-                   models = c(model.files.lmer.Tf.4,# model.files.lmer.Tf.2),
+                   models = c(model.files.lmer.Tf.4),# model.files.lmer.Tf.2),
                    traits = c("SLA", "Wood.density", "Max.height"))
 
 ## format.all.output.lmer(file.name = "NA.wwf.results.nolog.all.rds",

scripts/analysis/results.fig.R

@@ -186,6 +186,12 @@ pdf('figs/figres4.pdf', width = 8, height = 7)
 fun.plot.wd.sl.param(list.all.results.set)
 dev.off()
 
+pdf('figs/figres4.TP.pdf', width = 8, height = 7)
+fun.plot.wd.sl.param(list.all.results.set,
+                     model = 'lmer.LOGLIN.ER.AD.Tf.MAT.MAP.r.set.species',
+                     MAT.MAP.TF = TRUE)
+dev.off()
+
 
 pdf('figs/figres4b.pdf', width = 8, height = 11)
 fun.plot.all.param(list.all.results.set)