Commit 6a316e07 authored by Kunstler Georges's avatar Kunstler Georges
Browse files

new figures for presentations

parent 5f18d2cc
......@@ -1399,12 +1399,14 @@ plot.param.mean.and.biomes.fixed <- function(list.res,
'Compet response x trait',
'Compet x trait dissimilarity'),
param.print = 1:5,
traits.print = 1:3,
col.names = fun.col.param() ,
data.type = "simple",
col.vec,
pch.vec,
names.bio,
intra.TF = FALSE,
legend.TF = TRUE,
...){
if(!intra.TF) x.line <- -0.91
if(intra.TF) x.line = -0.91
......@@ -1418,7 +1420,7 @@ traits_letters <- c('a', 'b', 'c')
names(traits_letters) <- c('Wood.density', 'SLA', 'Max.height')
##################################
## model fixed biomes
for (i in traits){
for (i in traits[traits.print]){
list.temp <- list.res[[paste0(data.type, "_", i ,
"_", models[2])]]
list.temp.m <- list.res[[paste(data.type, "_", i ,
......@@ -1485,7 +1487,7 @@ names(traits_letters) <- c('Wood.density', 'SLA', 'Max.height')
mtext(traits.names[i], side=3, cex =2.25, line = 1)
}
fun.legend( biomes.names, biomes.c, col.vec, pch.vec)
if (legend.TF) fun.legend( biomes.names, biomes.c, col.vec, pch.vec)
}
......
......@@ -84,7 +84,7 @@ plot.param.mean.and.biomes.fixed(list.all.results.set , data.type = "simple",
dev.off()
pdf('figs/figres12_TP_intra.pdf', height = 14, width = 18)
pdf('figs/figres12_TP_intrab.pdf', height = 14, width = 18)
plot.param.mean.and.biomes.fixed(list.all.results.intra , data.type = "intra",
models = c('lmer.LOGLIN.ER.AD.Tf.MAT.MAP.intra.r.set.species',
'lmer.LOGLIN.ER.AD.Tf.MAT.MAP.intra.r.set.fixed.biomes.species'),
......@@ -107,12 +107,12 @@ plot.param.mean.and.biomes.fixed(list.all.results.intra , data.type = "intra",
intra.TF = TRUE)
dev.off()
biomes.white <- rep('white', length.out = length(biomes.white))
biomes.white <- rep('white', length.out = length(fun.col.pch.biomes()$col.vec))
names(biomes.white) <- names(fun.col.pch.biomes()$col.vec)
biomes.names.vec.w <- rep('', length.out = length(fun.biomes.names()))
names(biomes.names.vec.w) <- names(fun.biomes.names())
pdf('figs/figres12_TP_intra_m.pdf', height = 14, width = 16)
pdf('figs/figres12_TP_intra_mb.pdf', height = 14, width = 18)
plot.param.mean.and.biomes.fixed(list.all.results.intra , data.type = "intra",
models = c('lmer.LOGLIN.ER.AD.Tf.MAT.MAP.intra.r.set.species',
'lmer.LOGLIN.ER.AD.Tf.MAT.MAP.intra.r.set.fixed.biomes.species'),
......@@ -132,10 +132,35 @@ plot.param.mean.and.biomes.fixed(list.all.results.intra , data.type = "intra",
col.vec = biomes.white,
pch.vec = fun.col.pch.biomes()$pch.vec,
names.bio = names.biomes ,
xlim = c(-0.30, 0.35),
xlim = c(-0.30, 0.365),
intra.TF = TRUE)
dev.off()
pdf('figs/figres12_TP_intra_ma.pdf', height = 14, width = 18)
plot.param.mean.and.biomes.fixed(list.all.results.intra , data.type = "intra",
models = c('lmer.LOGLIN.ER.AD.Tf.MAT.MAP.intra.r.set.species',
'lmer.LOGLIN.ER.AD.Tf.MAT.MAP.intra.r.set.fixed.biomes.species'),
traits = c('Wood.density' , 'SLA', 'Max.height'),
param.vec = c("sumBn.intra", "sumBn.inter",
"sumTnTfBn.abs", "sumTfBn","sumTnBn",
"Tf"),
biomes.names = biomes.names.vec.w,
param.print = 1:6,
traits.print = 1,
param.names = c(expression('Intra'(alpha['0 intra'])),
expression('Inter'(alpha['0 inter'])),
expression('Dissim '(alpha['d'])),
expression('Tolerance '(alpha['t'])),
expression('Effect '(alpha['e'])),
expression("Trait link "(m[1])),
expression("Size "(gamma %*% log('D'))) ) ,
col.vec = biomes.white,
pch.vec = fun.col.pch.biomes()$pch.vec,
names.bio = names.biomes ,
xlim = c(-0.30, 0.365),
intra.TF = TRUE, legend.TF = FALSE)
dev.off()
pdf('figs/figres12_ecocode_TP_intra.pdf', height = 14, width = 16)
......
Markdown is supported
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