Hello,

I have specified the following model (a varying intercepts model for a reaction time (see figure for the sample data distribution) as a function of condition = [baseline, treatment1, treatment2]) with the bambi package:

b_model = Model(data)

results = b_model.fit(āreaction_time ~ 1 + conditionā, random = [ā1|subject_idā],

samples = 3000, tune = 2000)

As the title says: is there a possibility to model the response variable (in my case a reaction time) as lognormal (or any other right skewed option)? I have seen the possibility to specify a family, but canāt understand how to do it.

Thanks