Skip to content

Commit

Permalink
jags model.R
Browse files Browse the repository at this point in the history
  • Loading branch information
XuchenF authored Jun 5, 2024
1 parent baa0ea1 commit 395e658
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion jags model.R
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@ nyr<-length(ct)
N<-length(SSB)
jags.data<- c('ct','btj', "SSB","N",'nyr', 'prior.k', 'startbio', 'q.priorj',"Q.priorj","init.Q",'init.q','init.h','init.k','pen.bk','pen.F','b.yrs','bk.beta','CV.C','CV.cpue','nbks')

jags.save.params <- c('h','k','q',"Q",'P','ct.jags','cpuem', "cpuemm",'proc.logB','B','F')#if run the Schsefer function, use r to replace h.
jags.save.params <- c('h','k','q',"Q",'P','ct.jags','cpuem', "cpuemm",'proc.logB','B','F')#if run the Schsefer function, use r to replace h. And x is added in BHDPF2.

j.inits <- function(){list("h"=rlnorm(1, mean=log(init.h), sd=0.17), "k"=rlnorm(1, mean=log(init.k), sd=0.2),"q"=rlnorm(1,mean=log(init.q),sd=0.2),"Q"=rlnorm(1,mean=log(init.Q),sd=0.2),"itau2"=1000,"isigma2"=1000)}
#if run the Schsefer function, use r and init.r to replace h and init.h.
Expand Down

0 comments on commit 395e658

Please sign in to comment.