for (i in 1:100){

total[i] = inv_logit(b_0 + b_titer*(loglower_trans_1[n]+(i*diff_logtiter[n])/100));

y_step[i] = bernoulli_pmf(infection[n] | total[i]);

}

I would like bernoulli_pmf, but seems there is no this function, I only could get bernoulli_lpmf…but I don’t need log at this step