Please also provide the following information in addition to your question:
- Operating System: OSX 10.11.6
- brms Version: 2.4.0
I am using a cumulative() family to model membership for 3 discrete states (state) based on a single continuous variable (pheno).
When plotting the marginal effects I am confused by the x-axis. The probabilities do not match up with my continuous data. My “pheno” variable has a range:
> range(data$pheno)  0.4676923 0.9285714
yet the highest probability for inclusion in state 1 is from 0.0 - 0.2, which is out of the range of my variable.
I am not sure if this is an error in the plotting of the x-axis scale or a broad misinterpretation of the plot on my part. Does something need to be re-scaled? I know from looking at the data that there is a trend. Here is a plot of the raw data:
Below is the model, and the data are attached.
Thanks for your time,
o_test <- brm( state ~ pheno, data = data, family = cumulative(), iter = 5000)
ordinal_test.csv (3.1 KB)