Plausible values in over argument
Hi, it seems that the results are incorrect (at least for S.E.) when running Rrepest with plausible values in the "over" argument. I have some examples working with PISA data.
For example, with PISA data: Rrepest::Rrepest(data =PISA_2022%>% drop_na(PV1READ_quart) , svy ="PISA2015", est = est("mean","PV@READ"), by = c("CNT"), over="PV@READ_quart", average =c(grp("OECDa", "CNT",OECD2022)) )
or with PISA for school data Rrepest::Rrepest(data =f_Std_Stq_AZE , svy ="PBTS", est = est("freq",c("ST129Q01PA_rec","ST129Q02PA_rec","ST129Q03PA_rec","ST129Q04PA_rec","ST129Q05PA_rec","ST129Q06PA_rec","ST129Q07PA_rec","ST129Q08PA_rec")), by = c("SchoolId") , over="PV@SCIE_quart", n.pvs = 10)
Happy to discuss this
Vanessa