• Add plot functions for Hierarchical Gaussian Filter models: plot_hgf_ibrb, plot_hgf_ibrb_single.

• Added a Hierarchical Gaussian Filter model for binary inputs and binary responses: hgf_ibrb for hierarchical Bayesian analysis and hgf_ibrb_single for individual Bayesian analysis.
• Added new article: Hierarchical Bayesian Analysis on Hierarchical Gaussian Filter

Fixed a pkgdown error.

• Added a drift diffusion model and two reinforcement learning-drift diffision models for the probabilistic selection task: pstRT_ddm, pstRT_rlddm1, and pstRT_rlddm6.
• Added multiple models for the banditNarm task: banditNarm_2par_lapse, banditNarm_4par, banditNarm_delta, banditNarm_kalman_filter, banditNarm_lapse, banditNarm_lapse_decay, and banditNarm_singleA_lapse.
• Fixed bart_ewmv to avoid dividing by zero.

• Fix the symlink error in the Python version due to releasing with poetry
• Fix minor errors in both R and Python

• Added the cumulative model for the Cambridge gambling task: cgt_cm.
• Added two new models for aversive learning tasks: alt_delta and alt_gamma.
• Added exponential-weight mean-variance model for BART task: bart_ewmv.
• Added simple Q learning model for the probabilistic selection task: prl_Q.
• Added signal detection theory model for 2-alternative forced choice task: task2AFC_sdt.

• Fix an error on using data.frame objects as data (#112).

• Minor fix on R and Python codes (R, #111).

Major changes

• Now, hBayesDM has both R and Python version, with same models included!
  You can run hBayesDM with a language you prefer!
• Models in hBayesDM are now specified as YAML files. Using the YAML files,
  R and Python codes are generated automatically. If you want to contribute
  hBayesDM by adding a model, what you have to do is just to write a Stan file
  and to specify its information! You can find how to do in the hBayesDM wiki
  (https://github.com/CCS-Lab/hBayesDM/wiki).
• Model functions try to use parameter estimates using variational Bayesian
  methods as its initial values for MCMC sampling by default (#96). If VB
  estimation fails, then it uses random values instead.
• The data argument for model functions can handle a data.frame object (#2, #98).
• choiceRT_lba and choiceRT_lba_single are temporarily removed since their codes
  are not suitable to the new package structure. We plan to re-add the models
  in future versions.
• The Cumulative Model for Cambridge Gambling Task is added (cgt_cm; #108).

Minor changes

• The tau parameter in all models for the risk aversion task is modified to
  be bounded to [0, 30] (#77, #78).
• bart_4par is fixed to compute subject-wise log-likelihood (#82).
• extract_ic is fixed for its wrong rep function usage (#94, #100).
• The drift rate (delta parameter) in choiceRT_ddm and choiceRT_ddm_single is
  unbounded and now it is estimated between [-Inf, Inf] (#95, #107).
• Fix a preprocessing error in choiceRT_ddm and choiceRT_ddm_single (#95, #109).
• Fix igt_orl for a wrong Matt trick operation (#110).

• Add three new models for the bandit4arm task: bandit4arm_2par_lapse,
  bandit4arm_lapse_decay and bandit4arm_singleA_lapse.
• Fix various (minor) errors.