fit_jabba: Fit JABBA model

Description Usage Arguments Value Examples

View source: R/fit_jabba.R

Description

Fits JABBA model in JAGS and produce output object as list()

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
fit_jabba(
  jbinput,
  ni = 30000,
  nt = 5,
  nb = 5000,
  nc = 2,
  init.values = FALSE,
  init.K = NULL,
  init.r = NULL,
  init.q = NULL,
  peels = NULL,
  save.all = FALSE,
  save.jabba = FALSE,
  save.csvs = FALSE,
  save.prjkobe = FALSE,
  output.dir = getwd(),
  silent = FALSE
)

Arguments

jbinput

List of input variables as output by build_jabba() MCMC settings

ni

number of iterations

nt

thinning interval of saved iterations

nb

burn-in

nc

number of mcmc chains Initial values

init.values

= FALSE,

init.K

= NULL,

init.r

= NULL,

init.q

= NULL,# vector

peels

= NULL, # retro peel option

save.all

= FALSE

save.jabba

= FALSE

output.dir

path to save plot. default is getwd()

silent

option mute messages

Value

A result list containing estimates of model input, settings and results

Examples

1
2
data(bet)
jbinput <- build_jabba()

Henning-Winker/JABBApkg documentation built on Oct. 29, 2020, 9:11 a.m.