Commit 822f6f82 authored by Georges Kunstler's avatar Georges Kunstler
Browse files

further progress on improving jags

parent 3d41364d
......@@ -37,15 +37,15 @@ SUMABSBBN[i] <- (vec_sumTnTfBn_abs_biomes)*sumTnTfBn_abs[i]
### species random param
for (n in 1:N_species)
{
vec_logD_species[n] ~ dnorm(0,tau_logD_species)T(-7, 7)
vec_sumBn_species[n] ~ dnorm(0,tau_sumBn_species)T(-7, 7)
intercept_species[n] ~ dnorm(0,tau_inter_species)T(-7, 7)
vec_logD_species[n] ~ dnorm(0,tau_logD_species)
vec_sumBn_species[n] ~ dnorm(0,tau_sumBn_species)
intercept_species[n] ~ dnorm(0,tau_inter_species)
}
## plot effect
for (j in 1:N_plot)
{
intercept_plot[j] ~ dnorm(0,tau_inter_plot)T(-7, 7)
intercept_plot[j] ~ dnorm(0,tau_inter_plot)
}
## tree effect
## for (j in 1:N_tree)
......@@ -56,7 +56,7 @@ for (j in 1:N_plot)
## set effect
for (j in 1:N_set)
{
intercept_set[j] ~ dnorm(0,tau_inter_set)T(-7, 7)
intercept_set[j] ~ dnorm(0,tau_inter_set)
}
## biomes
......
......@@ -37,26 +37,26 @@ SUMABSBBN[i] <- (vec_sumTnTfBn_abs_biomes[biomes_id[i]])*sumTnTfBn_abs[i]
### species random param
for (n in 1:N_species)
{
vec_logD_species[n] ~ dnorm(0,tau_logD_species)T(-7, 7)
vec_sumBn_species[n] ~ dnorm(0,tau_sumBn_species)T(-7, 7)
intercept_species[n] ~ dnorm(0,tau_inter_species)T(-7, 7)
vec_logD_species[n] ~ dnorm(0,tau_logD_species)
vec_sumBn_species[n] ~ dnorm(0,tau_sumBn_species)
intercept_species[n] ~ dnorm(0,tau_inter_species)
}
## plot effect
for (j in 1:N_plot)
{
intercept_plot[j] ~ dnorm(0,tau_inter_plot)T(-7, 7)
intercept_plot[j] ~ dnorm(0,tau_inter_plot)
}
## tree effect
## for (j in 1:N_tree)
## {
## intercept_tree[j] ~ dnorm(0,tau_inter_tree)T(-7, 7)
## intercept_tree[j] ~ dnorm(0,tau_inter_tree)
## }
## set effect
for (j in 1:N_set)
{
intercept_set[j] ~ dnorm(0,tau_inter_set)T(-7, 7)
intercept_set[j] ~ dnorm(0,tau_inter_set)
}
## biomes
......
......@@ -11,6 +11,7 @@ run.multiple.model.for.set.one.trait(model.files.jags.Tf.1, run.jags,
var.sample = 'ecocode.id',
iter = 200, warmup = 20, chains = 3,
init.TF = TRUE)
resb <- run.jags(model.files.jags.Tf.1b, trait = 'Wood.density', type.filling = 'species',
init.TF = FALSE, cat.TF = FALSE,
sample.size = 10000, iter = 2000,
......
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