Measurement error for mixed models
Thanks for the code and suggestions. Trying to run the analysis in metafor was on the TODO list. Thank you for saving me the time on figuring out the syntax. If I may ask for clarification, what do you mean by "the variance component is larger than 0" ? Best wishes Krzysztof ??????? Original Message ???????
On Thursday, September 12, 2019 9:04 PM, Viechtbauer, Wolfgang (SP) <wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
Thanks for the feedback.
For what it's worth, you can fit the model with those 0 variances with metafor:
df_Dummy<-as.data.frame(matrix(rnorm(20),ncol=2,nrow=10))
df_Dummy$id<-as.factor(1:nrow(df_Dummy))
df_Dummy$merror<-rexp(10)
colnames(df_Dummy)[1:2]<-c("x","y")
df_Dummy$merror[1:4]<-0
library(metafor)
rma(y ~ x, merror, data=df_Dummy, method="ML")
This will work as long as the variance component is larger than 0.
Best,
Wolfgang