Sto iniziando a dilettarsi con l'uso di glmnetcon LASSO Regressione dove il mio risultato di interesse è dicotomica. Di seguito ho creato un piccolo frame di dati finti: age <- c(4, 8, 7, 12, 6, 9, 10, 14, 7) gender <- c(1, 0, 1, 1, 1, 0, 1, 0, 0) …
Ho un background statistico di livello introduttivo (suppongo di conoscere statistiche matematiche e probabilità a livello universitario (ad es. Wackerly et al., Probabilità di Ross) e avere una certa conoscenza della teoria delle misure). Di recente ho iniziato un lavoro facendo progettazione sperimentale e rapporti statistici nelle statistiche sull'istruzione e …
Supponiamo che due gruppi, composti da e n 2, ciascuno classifichi un insieme di 25 elementi dal più al meno importante. Quali sono i modi migliori per confrontare queste classifiche?n1n1n_1n2n2n_2 Chiaramente, è possibile eseguire 25 test di Mann-Whitney U, ma ciò comporterebbe 25 risultati di test da interpretare, che possono …
Sto recensendo un articolo di un giornale accademico e gli autori hanno scritto quanto segue come giustificazione per non riportare alcuna statistica inferenziale (ho deidentificato la natura dei due gruppi): In totale, 25 dei 2,349 (1,1%) rispondenti presentate X . Ci asteniamo in modo appropriato dalla presentazione di analisi che …
Ho un set di dati molto grande e mancano circa il 5% di valori casuali. Queste variabili sono correlate tra loro. Il seguente set di dati R è solo un esempio di giocattolo con dati correlati fittizi. set.seed(123) # matrix of X variable xmat <- matrix(sample(-1:1, 2000000, replace = TRUE), …
Il mgcvpacchetto per Rha due funzioni per adattare le interazioni del prodotto tensore: te()e ti(). Comprendo la divisione di base del lavoro tra i due (adattamento di un'interazione non lineare rispetto alla scomposizione di questa interazione in effetti principali e un'interazione). Quello che non capisco è perché te(x1, x2)e ti(x1) …
Quindi, prima di tutto, ho fatto alcune ricerche su questo forum e so che sono state poste domande molto simili, ma di solito non hanno ricevuto una risposta corretta o, a volte, la risposta non è abbastanza dettagliata da poter essere compresa. Quindi questa volta la mia domanda è: ho …
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 …
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, …
Sto cercando di adattare un modello a tempo discreto in R, ma non sono sicuro di come farlo. Ho letto che puoi organizzare la variabile dipendente in diverse righe, una per ogni osservazione temporale e utilizzare la glmfunzione con un collegamento logit o cloglog. In questo senso, ho tre colonne: …
Supponiamo di avere un campione di due popolazioni: Ae B. Supponiamo che queste popolazioni siano composte da individui e scegliamo di descriverli in termini di caratteristiche. Alcune di queste funzioni sono categoriche (ad esempio, guidano per funzionare?) E alcune sono numeriche (ad esempio, la loro altezza). Chiamiamo queste funzionalità: . …
Supponiamo di avere un campione di frequenze di 4 possibili eventi: Event1 - 5 E2 - 1 E3 - 0 E4 - 12 e ho le probabilità attese dei miei eventi: p1 - 0.2 p2 - 0.1 p3 - 0.1 p4 - 0.6 Con la somma delle frequenze osservate dei …
Come posso verificare se la media (ad esempio, la pressione sanguigna) di un sottogruppo (ad esempio, quelli che sono morti) differisce da tutto il gruppo (ad esempio, tutti coloro che hanno avuto la malattia, compresi quelli che sono morti)? Chiaramente, il primo è un sottogruppo del secondo. Quale test di …
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.