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].
Sto sperimentando l'algoritmo della macchina per aumentare il gradiente tramite il caretpacchetto in R. Utilizzando un piccolo set di dati di ammissione al college, ho eseguito il seguente codice: library(caret) ### Load admissions dataset. ### mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") ### Create yes/no levels for admission. ### mydata$admit_factor[mydata$admit==0] <- "no" mydata$admit_factor[mydata$admit==1] <- …
La domanda del CV Sto cercando di fornire una (a) rappresentazione matematica dettagliata e concisa di un modello di effetti misti. Sto usando il lme4pacchetto in R. Qual è la rappresentazione matematica corretta per il mio modello? I dati, la domanda scientifica e il codice R. Il mio set di …
I test di permutazione (chiamati anche test di randomizzazione, test di ri-randomizzazione o test esatto) sono molto utili e sono utili quando l'assunzione della distribuzione normale richiesta da per esempio t-testnon è soddisfatta e quando la trasformazione dei valori per classifica del test non parametrici come Mann-Whitney-U-testquesto porterebbero alla perdita …
Vorrei abbinare le uscite di lmer (davvero glmer) con un esempio binomiale di giocattoli. Ho letto le vignette e credo di capire cosa sta succedendo. Ma a quanto pare non lo faccio. Dopo essermi bloccato, ho corretto la "verità" in termini di effetti casuali e ho seguito solo la stima …
Ho un set di dati in cui la variabile che vorrei usare come effetto casuale ha una sola osservazione per alcuni livelli. Sulla base delle risposte alle domande precedenti, ho raccolto che, in linea di principio, questo può andare bene. Posso adattare un modello misto a soggetti che hanno solo …
Attualmente sto esaminando alcuni lavori e ho riscontrato quanto segue, il che mi sembra sbagliato. Due modelli misti sono montati (in R) usando lmer. I modelli non sono nidificati e vengono confrontati mediante test del rapporto di verosimiglianza. In breve, ecco un esempio riproducibile di ciò che ho: set.seed(105) Resp …
Ho un set di dati di misure ripetute sbilanciate da analizzare e ho letto che il modo in cui la maggior parte dei pacchetti statistici gestisce questo con ANOVA (ovvero la somma dei quadrati di tipo III) è sbagliato. Pertanto, vorrei utilizzare un modello di effetti misti per analizzare questi …
Nello specificare un modello di effetti misti incrociati, sto cercando di includere interazioni. Tuttavia, viene visualizzato il seguente messaggio di errore: Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) …
Ho una domanda sul mio uso di un modello misto / lmer. Il modello base è questo: lmer(DV ~ group * condition + (1|pptid), data= df) Il gruppo e la condizione sono entrambi fattori: il gruppo ha due livelli (gruppo A, gruppo B) e la condizione ha tre livelli (condizione1, …
Ho eseguito un esperimento in cui ho cresciuto diverse famiglie provenienti da due diverse popolazioni di origine. Ad ogni famiglia è stato assegnato uno di due trattamenti. Dopo l'esperimento ho misurato diversi tratti su ciascun individuo. Per verificare l'effetto del trattamento o della fonte, nonché la loro interazione, ho usato …
Supponiamo che io abbia NNN partecipanti, ognuno dei quali dà una risposta YYY 20 volte, 10 in una condizione e 10 in un'altra. Adatto un modello lineare di effetti misti confrontando YYY in ogni condizione. Ecco un esempio riproducibile che simula questa situazione usando il lme4pacchetto in R: library(lme4) fml …
Supponiamo di avere una certa certa variabile di risposta yijyijy_{ij} che è stato misurato da jjj esima sibling iii esima famiglia. Inoltre, alcuni dati comportamentali xijxijx_{ij} sono stati raccolti contemporaneamente da ciascun soggetto. Sto cercando di analizzare la situazione con il seguente modello lineare a effetti misti: yij=α0+α1xij+δ1ixij+εijyij=α0+α1xij+δ1ixij+εijy_{ij} = \alpha_0 …
Qualcuno può dirmi come fare in modo che R valuti il punto di interruzione in un modello lineare a tratti (come parametro fisso o casuale), quando devo anche stimare altri effetti casuali? Di seguito ho incluso un esempio di giocattolo che si adatta a una regressione di un bastone da …
Douglas Bates afferma che i seguenti modelli sono equivalenti "se la matrice varianza-covarianza per gli effetti casuali con valori vettoriali ha una forma speciale, chiamata simmetria composta" ( diapositiva 91 in questa presentazione ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor + …
Sto eseguendo un modello di regressione logistica nel modulo: lmer(response~1+(1|site), family=binomial, REML = FALSE) Normalmente calcolerei l'ICC dall'intercettazione e dalle varianze residue, ma il riepilogo del modello non include la varianza residua. Come posso calcolarlo?
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.