From 3af17152add6d2faf6ef9e1e01bb88844fbd3b33 Mon Sep 17 00:00:00 2001 From: Douglas Ezra Morrison Date: Sun, 4 Jan 2026 23:03:57 -0800 Subject: [PATCH] Increment version number to 0.0.0.9045 --- DESCRIPTION | 2 +- inst/extdata/model.jags | 12 +++++++----- 2 files changed, 8 insertions(+), 6 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index 251af531..a813d1d1 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,6 +1,6 @@ Package: serodynamics Title: What the Package Does (One Line, Title Case) -Version: 0.0.0.9044 +Version: 0.0.0.9045 Authors@R: c( person("Peter", "Teunis", , "p.teunis@emory.edu", role = c("aut", "cph"), comment = "Author of the method and original code."), diff --git a/inst/extdata/model.jags b/inst/extdata/model.jags index d38909e2..81ed58e2 100644 --- a/inst/extdata/model.jags +++ b/inst/extdata/model.jags @@ -55,11 +55,13 @@ model { } # these are random effects - y0[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,1]) - y1[subj,cur_antigen_iso] <- y0[subj,cur_antigen_iso] + exp(par[subj,cur_antigen_iso,2]) # par[,,2] must be log(y1-y0) - t1[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,3]) - alpha[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,4]) # `nu` in the paper - shape[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,5]) + 1 # `r` in the paper + y0[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,1]) + # par[,,2] is log(y1-y0): + change_y[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,2]) + y1[subj,cur_antigen_iso] <- y0[subj,cur_antigen_iso] + change_y[subj,cur_antigen_iso] + t1[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,3]) + alpha[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,4]) # `nu` in the paper + shape[subj,cur_antigen_iso] <- exp(par[subj,cur_antigen_iso,5]) + 1 # `r` in the paper # `n_params` is the number of model parameters; y0, y1, t1, alpha (aka nu), and r # this is the prior distribution