Hi all,

I’m new to Bayesian statistics using brms and stan, hence apologies if my question seems very basic.

I’m trying to fit a regression model in brms with informative priors. The prior knowledge is based on a previously conducted pilot study, for which I produced posterior distributions of the intercept, group level effects and population level effects. I am then using this information of the posteriors from my pilot study as priors for my actual study (including of course only the recent data collected).

My pilot model was:

fitdur <- brm(formula = Duration ~ DSI + rankdiff + agediff + sexcode +(1|Initiator), data=entryduration, family=lognormal(), warmup= 1000, iter = 2000, chains = 4, control = list(adapt_delta = 0.99, stepsize=0.01, max_treedepth=15))

I am, however, struggling to understand how to set the right distribution of priors for my final updated model; I see that, by using prior_summary, the default priors of my pilot model were all student-t-priors. But since I used a model with the family “lognormal”, I wonder whether I should change the default priors of my final, updated model (including the new data collected and the prior of the pilot study) to lognormal priors for all parameters (i.e., intercept, b’s, sd, sigma)? Or, should I leave the default prior distributions as they are and just insert the posterior means and sd of my pilot model? If so, should I transform these values (means and sd), considering they stem from a lognormal model?

As a more general question - should you inform your model by including the distribution of your previous posteriors as your new prior distributions?

Any advice would be of great help, thank you!

Best, Raphaela