Hey,
I was wondering about the output of the conditional_effects() function, especially about the interpretation of the standard error.
Until now I was assuming, that the standard error was the same as the standard deviation so that I could get the upper and lower bounds of a 95% (0.025 - 0.975) interval with mean ± 2se.
However, when I tried this for one of my models, I had a small error in there:

(0.01663656 + 2*0.002212549) - 0.02146316 should then give me 0 but results in -0.000401502 instead. I made sure to use robust=FALSE to get the mean and not the median so my guess is that my interpretation of the standard error is off.

Hmm, well it’s very close to zero, so maybe a rounding error, which doesn’t make much of a difference? Let’s see what @paul.buerkner has to say about this.

@paul.buerkner Thanks for the reply. I want to present the conditional effect for one of my population-level effects with mean, median, sd and 95% interval and I thought conditional_effects(fit, effects=c("effect)$effect would give me that information.
I was assuming, that the se__ column of the resulting table was the same as the standard deviation, but when I tried to use the mean+2sd formula the result was not the same as the upper bound which confuses me.

The approach mean +2se is a frequentist way to obtain CI bounds assuming normality of the estimate. We don’t need such approximation in Bayesian statistics when we have posterior draws from which we can simply compute CI bounds based on the quantiltes of the draws.

Is my understanding correct then, that if the effect’s distribution were of perfect gaussian shape, the mean +2se approach would still work as expected?
Thank you for the quick answer!