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.
Ci sono molte discussioni in corso su questo forum sul modo corretto di specificare vari modelli gerarchici usando lmer. Ho pensato che sarebbe bello avere tutte le informazioni in un unico posto. Un paio di domande per iniziare: Come specificare più livelli, in cui un gruppo è nidificato nell'altro: è …
Ecco come ho compreso gli effetti casuali nidificati o incrociati: Gli effetti casuali nidificati si verificano quando un fattore di livello inferiore appare solo all'interno di un determinato livello di un fattore di livello superiore. Ad esempio, gli alunni delle classi in un determinato momento. In lme4ho pensato che rappresentiamo …
Se rielaboriamo un bagliore, potremmo ricevere un avviso che ci dice che il modello sta trovando difficoltà a convergere ... ad es >Warning message: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, : Model failed to converge with max|grad| = 0.00389462 (tol = 0.001) un altro modo per verificare la convergenza …
Ho adattare alcuni modelli effetti misti (in particolare modelli longitudinali) utilizzando lme4in Rma vorrei padroneggiare davvero i modelli e il codice che va con loro. Tuttavia, prima di immergermi con entrambi i piedi (e acquistare alcuni libri) voglio essere sicuro di imparare la biblioteca giusta. Ho usato lme4fino ad ora …
Ho letto nell'abstract di questo documento che: "La procedura di massima verosimiglianza (ML) di Hartley aud Rao viene modificata adattando una trasformazione di Patterson e Thompson che suddivide la verosimiglianza in due parti, una delle quali priva degli effetti fissi. Massimizzare questa parte produce ciò che viene chiamato verosimiglianza massima …
Considera i seguenti tre fenomeni. Paradosso di Stein: dati alcuni dalla distribuzione normale multivariata in , la media campionaria non è un ottimo stimatore della media vera. Si può ottenere una stima con errore quadratico medio inferiore se si riducono tutte le coordinate della media campionaria verso zero [o verso …
Da un po 'di tempo uso felicemente modelli di effetti misti con dati longitudinali. Vorrei poter inserire le relazioni AR in lmer (penso di aver ragione nel non poterlo fare?) Ma non penso che sia disperatamente importante, quindi non mi preoccupo troppo. Ho appena incontrato equazioni di stima generalizzate (GEE) …
Uso lme4 in R per adattarmi al modello misto lmer(value~status+(1|experiment))) dove il valore è continuo, lo stato e l'esperimento sono fattori e ottengo Linear mixed model fit by REML Formula: value ~ status + (1 | experiment) AIC BIC logLik deviance REMLdev 29.1 46.98 -9.548 5.911 19.1 Random effects: Groups …
Pensavo che il "modello a effetti casuali" in econometria corrispondesse a un "modello misto con intercettazione casuale" al di fuori dell'econometria, ma ora non ne sono sicuro. Vero? L'econometria usa termini come "effetti fissi" e "effetti casuali" in modo leggermente diverso dalla letteratura sui modelli misti, e questo provoca una …
Di recente ho misurato il modo in cui il significato di una nuova parola viene acquisito su esposizioni ripetute (pratica: dal giorno 1 al giorno 10) misurando ERP (EEG) quando la parola veniva vista in contesti diversi. Ho anche controllato le proprietà del contesto, ad esempio, la sua utilità per …
Mi chiedo se fa differenza nell'interpretazione se solo le variabili dipendenti, dipendenti e indipendenti, o solo le variabili indipendenti, vengono trasformate in log. Considera il caso di log(DV) = Intercept + B1*IV + Error Posso interpretare il IV come l'aumento percentuale, ma come cambia quando lo faccio log(DV) = Intercept …
EDIT 2: Inizialmente pensavo di dover eseguire un ANOVA a due fattori con misure ripetute su un fattore, ma ora penso che un modello lineare a effetto misto funzionerà meglio per i miei dati. Penso di sapere quasi cosa deve succedere, ma sono ancora confuso da alcuni punti. Gli esperimenti …
Ho calcolato AIC e AICc per confrontare due modelli misti lineari generali; Gli AIC sono positivi con il modello 1 con un AIC inferiore rispetto al modello 2. Tuttavia, i valori per AICc sono entrambi negativi (il modello 1 è ancora <modello 2). È valido utilizzare e confrontare valori AICc …
Voglio ottenere un intervallo di previsione attorno a una previsione da un modello lmer (). Ho trovato alcune discussioni su questo: http://rstudio-pubs-static.s3.amazonaws.com/24365_2803ab8299934e888a60e7b16113f619.html http://glmm.wikidot.com/faq ma sembrano non tenere conto dell'incertezza degli effetti casuali. Ecco un esempio specifico. Sto correndo pesce d'oro. Ho dei dati sulle ultime 100 gare. Voglio prevedere il …
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.