- Have you released a vignette yet on measurement error using the brms package? I am trying to do a multilevel measurement error model but within a hierarchical structure.
The hierarchy is over time (year_cat5) and the measurement error is on the variable “Q” which we are using to explain the investment rate of firms (inv_rf).
- In, principle, in this formulation, for a hierarchical measurement error model in brms correct?
Links: mu = identity; sigma = identity; nu = identity
Formula: inv_rf ~ me(Q_manual, Q_sd) + cashflow_rf + (1 + me(Q_manual, Q_sd) + cashflow_rf | year_cat5) + capital_output + fic_advanced_nousa
Data: data_STAN_lender_sam (Number of observations: 584)
Samples: 3 chains, each with iter = 4000; warmup = 2000; thin = 1;
total post-warmup samples = 6000
I ask because when adding the hierarchical structure the impact of the measurement error vanishes basically, Ie when running a ‘fixed effects’ measurement error it makes a big difference; but the hierarchical measurement error makes little difference.
- One issue is that I have a large dataset (360,000 data points. But it doesn’t seem to work with more 10,000 points or so. Understandable given that every datapoint will now be treated as an unknown parameter. What do you recommend to try work around this?
- Operating System:
- brms Version: