lme4 e nlme sono pacchetti R utilizzati per il montaggio di modelli di effetti misti lineari, lineari generalizzati e non lineari. Per domande generali sui modelli misti utilizzare il tag [modello misto].
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 …
Il montaggio di una regressione logistica con lme4 termina con Error in mer_finalize(ans) : Downdated X'X is not positive definite. Una probabile causa di questo errore è apparentemente la carenza di rango. Cos'è la carenza di rango e come devo affrontarla?
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 …
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 …
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 …
Derivare numericamente gli MLE di GLMM è difficile e, in pratica, lo so, non dovremmo usare l'ottimizzazione della forza bruta (ad esempio, usando optimin modo semplice). Ma per il mio scopo educativo, voglio provarlo per assicurarmi di capire correttamente il modello (vedi il codice qui sotto). Ho scoperto che ottengo …
EffectsIl pacchetto fornisce un modo molto veloce e conveniente per tracciare risultati di modelli a effetto misto lineare ottenuti attraverso il lme4pacchetto . La effectfunzione calcola gli intervalli di confidenza (IC) molto rapidamente, ma quanto sono affidabili questi intervalli di confidenza? Per esempio: library(lme4) library(effects) library(ggplot) data(Pastes) fm1 <- lmer(strength …
Di recente mi sono reso conto che un modello misto con solo soggetto come fattore casuale e gli altri fattori come fattori fissi equivalgono a un ANOVA quando si imposta la struttura correlativa del modello misto su simmetria composta. Pertanto, vorrei sapere cosa significa simmetria composta nel contesto di un …
Ho i dati raccolti da un esperimento organizzato come segue: Due siti, ciascuno con 30 alberi. 15 sono trattati, 15 controllano ogni sito. Da ogni albero, campioniamo tre pezzi dello stelo e tre pezzi delle radici, quindi 6 campioni di livello 1 per albero che è rappresentato da uno dei …
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.