# save all output in a list
source("R/analysis/lmer.output-fun.R")
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.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,
##                               model.files.lmer.Tf.intra.2),
##                    traits = c("SLA", "Wood.density", "Max.height"),
##                    data.type = 'intra')