-
Notifications
You must be signed in to change notification settings - Fork 0
/
lvb_jags2_multi_t0_free
32 lines (29 loc) · 1.11 KB
/
lvb_jags2_multi_t0_free
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
model {
for(i in 1:n) {
for(j in 1:4) {
L[i,j] ~ dnorm(mu[i,j], tau[j])
ypp[i,j] ~ dnorm(mu[i,j], tau[j])
}
mu[i,1] <- L_inf[1]*(1-exp(-k[1]*(a[i]-t0[1]))) # VBM - VB
mu[i,2] <- L_inf[2]*(1-exp(-k[2]*(a[i]-t0[2])))^P[2] # GGM - generalized VB
mu[i,3] <- L_inf[3]/(1+exp(-k[3]*(a[i]-t0[3]))) # LGM - logistic growth model
mu[i,4] <- L_inf[4]*exp((-1/k[4])*exp(-k[4]*(a[i]-t0[4]))) # PGM - Gompertz
}
for(j in 1:4) {
tau[j] <- pow(sig[j], -2)
sig[j] ~ dunif(0, 1000)
L_inf[j] ~ dnorm(400, .0001)
k[j] ~ dnorm(0,.1)T(0,)
t0[j] ~ dnorm(0,.01)T(,6)#<-0#
P[j] ~ dlnorm(0,1)
}
for(ifit in 1:nfit) {
mufit[ifit,1] <- L_inf[1]*(1-exp(-k[1]*(agefit[ifit]-t0[1]))) # VBM - VB
mufit[ifit,2] <- L_inf[2]*(1-exp(-k[2]*(agefit[ifit]-t0[2])))^P[2] # GGM - generalized VB
mufit[ifit,3] <- L_inf[3]/(1+exp(-k[3]*(agefit[ifit]-t0[3]))) # LGM - logistic growth model
mufit[ifit,4] <- L_inf[4]*exp((-1/k[4])*exp(-k[4]*(agefit[ifit]-t0[4]))) # PGM - Gompertz
for(j in 1:4) {
ppfit[ifit,j] ~ dnorm(mufit[ifit,j], tau[j])
}
}
}