lm è il nome della funzione del modello lineare (cioè di regressione multipla) nel pacchetto delle statistiche R. Per i modelli lineari in generale usa invece il tag `linear-model`.
Per la mia comprensione, sono interessato a replicare manualmente il calcolo degli errori standard dei coefficienti stimati poiché, ad esempio, viene fornito l'output della lm()funzione R, ma non sono stato in grado di fissarlo. Qual è la formula / implementazione utilizzata?
In particolare, voglio sapere se c'è una differenza tra lm(y ~ x1 + x2)e glm(y ~ x1 + x2, family=gaussian). Penso che questo caso particolare di glm sia uguale a lm. Ho sbagliato?
Voglio presumere che la temperatura della superficie del mare del Mar Baltico sia la stessa anno dopo anno, e quindi descriverlo con un modello funzione / lineare. L'idea che ho avuto è stata quella di inserire solo l'anno come un numero decimale (o num_months / 12) e capire quale dovrebbe …
Qual è la formula esatta usata in R lm() per il rettangolo R rettificato? Come posso interpretarlo? Formule rettangolari rettificate Sembra che esistano diverse formule per calcolare il rettangolo R rettificato. Formula di Wherry: 1 - ( 1 - R2) ( n - 1 )( n - v )1-(1-R2)(n-1)(n-v)1-(1-R^2)\frac{(n-1)}{(n-v)} Formula …
Possiamo usare lm()per prevedere un valore, ma in alcuni casi abbiamo ancora bisogno dell'equazione della formula del risultato. Ad esempio, aggiungere l'equazione ai grafici.
Come prequel di una domanda sui modelli lineari misti in R e di condividere come riferimento per gli appassionati di statistiche per principianti / intermedi, ho deciso di pubblicare come "domande e risposte" indipendenti i passaggi coinvolti nel calcolo "manuale" del coefficienti e valori previsti di una semplice regressione lineare. …
Ho una trama di valori residui di un modello lineare in funzione dei valori adattati in cui l'eteroscedasticità è molto chiara. Tuttavia non sono sicuro di come procedere ora perché, per quanto ho capito, questa eteroscedasticità rende il mio modello lineare non valido. (È giusto?) Usa un robusto raccordo lineare …
Vorrei fare un test W di Shapiro Wilk e un test Kolmogorov-Smirnov sui residui di un modello lineare per verificare la normalità. Mi stavo solo chiedendo quali residui dovrebbero essere usati per questo - i residui grezzi, i residui di Pearson, i residui studentizzati o i residui standardizzati? Per un …
sfondo Sto cercando di capire il primo esempio in un corso sui modelli di adattamento (quindi questo può sembrare ridicolmente semplice). Ho fatto i calcoli a mano e corrispondono all'esempio, ma quando li ripeto in R, i coefficienti del modello sono disattivati. Ho pensato che la differenza potrebbe essere dovuta …
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. …
Potete darmi un esempio dell'uso degli stimatori sandwich per eseguire inferenze di regressione robuste? Riesco a vedere l'esempio in ?sandwich, ma non capisco bene come possiamo passare da lm(a ~ b, data)( codice r ) a una stima e un valore p risultanti da un modello di regressione usando la …
Sto semplicemente cercando di ricalcolare con dnorm () la verosimiglianza log fornita dalla funzione logLik da un modello lm (in R). Funziona (quasi perfettamente) per un elevato numero di dati (ad es. N = 1000): > n <- 1000 > x <- 1:n > set.seed(1) > y <- 10 + …
Ho un GLMM del modulo: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Quando uso drop1(model, test="Chi"), ottengo risultati diversi rispetto a quelli che utilizzo Anova(model, type="III")dal pacchetto auto o summary(model). Questi ultimi due danno le stesse risposte. Usando un mucchio di dati fabbricati, …
Un mio collega mi ha inviato questo problema apparentemente facendo il giro su Internet: If $3 = 18, 4 = 32, 5 = 50, 6 = 72, 7 = 98$, Then, $10 =$ ? La risposta sembra essere 200. 3*6 4*8 5*10 6*12 7*14 8*16 9*18 10*20=200 Quando eseguo una …
I seguenti innesti sono presi da questo articolo . Sono un novizio di bootstrap e sto cercando di implementare il bootstrap bootstrap parametrico, semiparametrico e non parametrico per il modello misto lineare con R bootpacchetto. Codice R Ecco il mio Rcodice: library(SASmixed) library(lme4) library(boot) fm1Cult <- lmer(drywt ~ Inoc + …
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.