I have a model like the following:

m <- brm(y | se (ySE) ~ 1 + me(x, sdx = xSE), d, family=gaussian, save_mevars = T,

prior=prior, sample_prior = T, etc)

I am a bit at a loss as to how to test the hypothesis that the relation between x and y is positive.

hypothesis(m, “x > 0”)

understandably gives a “Error: Some parameters cannot be found in the model”.

however, even using the predictor from the model summary fails with the same error:

hypothesis(m, “mexsdxEQxSE > 0”)

How should I specify the command?

The `class`

of these parameters is `bsp`

not `b`

. You see that when running `parnames(m)`

. I understand that’s possibly confusing, but I can’t change it until brms 3.0 and even then I am undecided whether to allign all these effects to have class `b`

.

Thanks! It makes sense. Maybe just an example with sbp in the help file would suffice, instead of a change of class.

Hi Paul,

Is it now possible to test hypotheses about effects with measurement error in brms using hypothesis()?

Best wishes,

Ladislas

Yes it is. What code did you try that didnt work?

I have tried `hypothesis(model, "bsp_meesT1seT1 = 0")`

where `esT1`

is the name of the predictor with known measurement error.

It throws out the following error message: `Error: Some parameters cannot be found in the model: 'b_bsp_meesT1seT1'`

I have also tried replacing `bsp_meesT1seT1`

with `esT1`

or `meesT1`

without much more success…

Aaah my bad, I should have read the doc… my apologies for digging up this old thread. Thank you!