I modelli misti (aka multilivello o gerarchici) sono modelli lineari che includono sia effetti fissi che casuali. Sono utilizzati per modellare dati longitudinali o nidificati.
Sto lavorando su un set di dati per valutare l'impatto dell'essiccazione sulle attività microbiche dei sedimenti. L'obiettivo è determinare se l'impatto dell'essiccazione varia tra i tipi di sedimenti e / o la profondità all'interno del sedimento. Il design sperimentale è il seguente: Il primo fattore Sedimento corrisponde a tre tipi …
Di recente mi è stato detto che non era possibile incorporare covariate variabili nel tempo in modelli misti longitudinali senza introdurre un ritardo per queste covariate. Puoi confermare / negare questo? Hai qualche riferimento su questa situazione? Propongo una situazione semplice per chiarire. Supponiamo di avere ripetute misure (diciamo più …
Sto cercando di riprodurre diversi test di interazione tra con entrambi lme lmersu misure ripetute (2x2x2). Il motivo per cui voglio confrontare entrambi i metodi è perché il GLM di SPSS per misure ripetute produce esattamente gli stessi risultati lmdell'approccio presentato qui, quindi alla fine voglio confrontare SPSS vs R-lmer. …
Sto provando a calcolare manualmente le previsioni di effetti casuali da un modello misto lineare, e usando la notazione fornita da Wood in Generalized Additive Models: un'introduzione con R (pag. 294 / pag. 307 di pdf), mi sto confondendo su ciò che ogni parametro rappresenta. Di seguito è riportato un …
Come controparte di questo post , ho lavorato sulla simulazione di dati con variabili continue, prestandomi a intercettazioni e pendenze correlate. Sebbene sul sito e al di fuori del sito siano presenti ottimi post su questo argomento , ho avuto difficoltà a trovare un esempio dall'inizio alla fine con dati …
Nella lmerfunzione all'interno lme4di Rc'è un invito a costruire una matrice modello di effetti casuali, , come spiegato qui , pagine 7 - 9.ZZZ Il calcolo di comporta i prodotti KhatriRao e / o Kronecker di due matrici, J_i e X_i . ZZZJiJiJ_iXiXiX_i La matrice è un boccone: " Matrice …
Sto eseguendo test post-hoc su un modello lineare a effetti misti in R( lme4pacchetto). Sto usando il multcomppacchetto ( glht()funzione) per eseguire i test post-hoc. Il mio disegno sperimentale è misure ripetute, con un effetto blocco casuale. I modelli sono specificati come: mymod <- lmer(variable ~ treatment * time + …
Date tre variabili, ye x, che sono continue positive e z, che è categorico, ho due modelli candidati dati da: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) e fit.fe <- lm( y ~ 1 + x ) Spero di confrontare …
Prima di porre questa domanda, ho cercato nel nostro sito e ho trovato molte domande simili (come qui , qui e qui ). Ma ritengo che quelle domande correlate non abbiano avuto una risposta o una discussione, quindi vorrei sollevare di nuovo questa domanda. Penso che dovrebbe esserci un grande …
Sto analizzando i risultati di un esperimento sul tempo di reazione in R. Ho eseguito una misurazione ripetuta ANOVA (1 fattore entro soggetto con 2 livelli e 1 fattore tra soggetto con 2 livelli). Ho eseguito un modello misto lineare simile e volevo riassumere i risultati più bassi sotto forma …
Nel modello misto, supponiamo che gli effetti casuali (parametri) siano variabili casuali che seguono le normali distribuzioni. Sembra molto simile al metodo bayesiano, in cui tutti i parametri sono considerati casuali. Quindi il modello a effetti casuali è una specie di caso speciale del metodo bayesiano?
Mi è sempre stato insegnato che gli effetti casuali influenzano solo la varianza (errore) e che gli effetti fissi influenzano solo la media. Ma ho trovato un esempio in cui gli effetti casuali influenzano anche la media - la stima del coefficiente: require(nlme) set.seed(128) n <- 100 k <- 5 …
Esistono diversi thread su R-sig-ME su come ottenere intervalli di confidenza per le previsioni usando lme4e nlmein R. Ad esempio qui e qui nel 2010, tra cui alcuni commenti di Dougals Bates, uno degli autori di entrambi i pacchetti. Esito a citarlo alla lettera, per paura che vengano portati fuori …
Sto lavorando su un set di dati. Dopo aver usato alcune tecniche di identificazione del modello, sono uscito con un modello ARIMA (0,2,1). Ho usato la detectIOfunzione nel pacchetto TSAin R per rilevare un valore anomalo innovativo (IO) alla 48a osservazione del mio set di dati originale. Come posso incorporare …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.