process output new numm model with MAT MAP

R/analysis/model.lmer/model.lmer.LOGLIN.MAT.MAP.r.set.fixed.biomes.species.R

+load.model <- function () {
+    list(name="lmer.LOGLIN.MAT.MAP.r.set.fixed.biomes.species",
+         var.BLUP = 'set.id',
+         lmer.formula.tree.id=formula("logG~1+(1|set.id)+(1|species.id)+(1|plot.id)+MAT+MAP+biomes.id+logD+sumBn+sumBn:biomes.id +(logD-1|species.id) +(sumBn-1|species.id)+(sumBn-1|set.id)"))
+}
+

R/analysis/model.lmer/model.lmer.LOGLIN.MAT.MAP.r.set.species.R

+load.model <- function () {
+    list(name="lmer.LOGLIN.MAT.MAP.r.set.species",
+         var.BLUP = 'set.id',
+         lmer.formula.tree.id=formula("logG~1+(1|set.id)+(1|species.id)+(1|plot.id)+MAT+MAP+logD+sumBn +(logD-1|species.id) +(sumBn-1|species.id)+(sumBn-1|set.id)"))
+}
+
+
+

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

@@ -6,17 +6,17 @@ source("R/analysis/lmer.run.R")
 
 ## ## save list of all output for  NA
 
-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),
-                   traits = c("SLA", "Wood.density", "Max.height"))
-
-## format.all.output.lmer(file.name = "NA.wwf.results.nolog.all.rds",
-##                    list.file.name = 'list.lmer.out.all.NA.simple.set.rds',
-##                    models = c(model.files.lmer.Tf.0, model.files.lmer.Tf.1,
-##                               model.files.lmer.Tf.2, model.files.lmer.Tf.3),
+## 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),
 ##                    traits = c("SLA", "Wood.density", "Max.height"))
 
+format.all.output.lmer(file.name = "NA.wwf.results.nolog.all.rds",
+                   list.file.name = 'list.lmer.out.all.NA.simple.set.rds',
+                   models = c(model.files.lmer.Tf.0, model.files.lmer.Tf.0b,
+                              model.files.lmer.Tf.1, model.files.lmer.Tf.3),
+                   traits = c("SLA", "Wood.density", "Max.height"))
+
 ## format.all.output.lmer(file.name = "NA.wwf.results.nolog.all.rds",
 ##                    list.file.name = 'list.lmer.out.all.NA.intra.set.rds',
 ##                    models = c(model.files.lmer.Tf.intra.1,