La libreria languageR fornisce un metodo (pvals.fnc) per eseguire test di significatività MCMC degli effetti fissi in un modello di regressione ad effetti misti adattandosi usando lmer. Tuttavia, pvals.fnc genera un errore quando il modello lmer include pendenze casuali.
C'è un modo per fare un test di ipotesi MCMC di tali modelli?
Se é cosi, come? (Per essere accettata, una risposta dovrebbe avere un esempio funzionante in R) In caso contrario, c'è un motivo concettuale / di calcolo perché non c'è modo?
Questa domanda potrebbe essere correlata a questa, ma non ho capito bene il contenuto abbastanza per esserne certo.
Modifica 1 : una dimostrazione del concetto che mostra che pvals.fnc () fa ancora "qualcosa" con i modelli lme4, ma che non fa nulla con i modelli a pendenza casuale.
library(lme4)
library(languageR)
#the example from pvals.fnc
data(primingHeid)
# remove extreme outliers
primingHeid = primingHeid[primingHeid$RT < 7.1,]
# fit mixed-effects model
primingHeid.lmer = lmer(RT ~ RTtoPrime * ResponseToPrime + Condition + (1|Subject) + (1|Word), data = primingHeid)
mcmc = pvals.fnc(primingHeid.lmer, nsim=10000, withMCMC=TRUE)
#Subjects are in both conditions...
table(primingHeid$Subject,primingHeid$Condition)
#So I can fit a model that has a random slope of condition by participant
primingHeid.lmer.rs = lmer(RT ~ RTtoPrime * ResponseToPrime + Condition + (1+Condition|Subject) + (1|Word), data = primingHeid)
#However pvals.fnc fails here...
mcmc.rs = pvals.fnc(primingHeid.lmer.rs)
Dice: Errore in pvals.fnc (primingHeid.lmer.rs): il campionamento MCMC non è ancora implementato in lme4_0.999375 per i modelli con parametri di correlazione casuali
Domanda aggiuntiva: pvals.fnc funziona come previsto per il modello di intercettazione casuale? Gli output dovrebbero essere attendibili?