Usa questo tag per qualsiasi domanda * sull'argomento * che (a) coinvolga `R` come parte critica della domanda o risposta prevista, e (b) non è * solo * su come usare` R`.
Sto cercando di utilizzare lme4::glmer()per adattare un modello misto generalizzato binomiale (GLMM) con una variabile dipendente che non è binaria, ma una variabile continua tra zero e uno. Si può pensare a questa variabile come a una probabilità; in effetti è la probabilità riportata da soggetti umani (in un esperimento …
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 …
Dopo aver letto questo post sul blog sui modelli strutturali delle serie temporali bayesiane, ho voluto esaminare l'implementazione nel contesto di un problema per il quale avevo precedentemente utilizzato ARIMA. Ho alcuni dati con alcuni componenti stagionali noti (ma rumorosi) - ci sono sicuramente componenti annuali, mensili e settimanali a …
Ho un esempio funzionante (in R), che sto cercando di capire ulteriormente. Sto usando Limma per creare un modello lineare e sto cercando di capire cosa sta succedendo passo dopo passo nei calcoli del cambio di piega. Sto principalmente cercando di capire cosa succede per calcolare i coefficienti. Da quello …
Vorrei sapere se la foresta casuale di Breiman (foresta casuale nel pacchetto R randomForest) utilizza come criterio di suddivisione (criterio per la selezione degli attributi) guadagno di informazioni o indice Gini? Ho provato a scoprirlo su http://www.stat.berkeley.edu/~breiman/RandomForests/cc_home.htm e nella documentazione per il pacchetto randomForest in R. Ma l'unica cosa che …
Ho due variabili: X e Y e ho bisogno di rendere il cluster massimo (e ottimale) = 5. Facciamo un diagramma ideale delle variabili come il seguente: Vorrei creare 5 cluster di questo. Qualcosa come questo: Quindi penso che questo sia un modello misto con 5 cluster. Ogni cluster ha …
Ho provato a calcolare l'AIC di una regressione lineare in R ma senza usare la AICfunzione, in questo modo: lm_mtcars <- lm(mpg ~ drat, mtcars) nrow(mtcars)*(log((sum(lm_mtcars$residuals^2)/nrow(mtcars))))+(length(lm_mtcars$coefficients)*2) [1] 97.98786 Tuttavia, AICfornisce un valore diverso: AIC(lm_mtcars) [1] 190.7999 Qualcuno potrebbe dirmi cosa sto facendo di sbagliato?
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 …
Sono nuovo nel glmnetpacchetto e non sono ancora sicuro di come interpretare i risultati. Qualcuno potrebbe aiutarmi a leggere il seguente diagramma di traccia? Il grafico è stato ottenuto eseguendo quanto segue: library(glmnet) return <- matrix(ret.ff.zoo[which(index(ret.ff.zoo)==beta.df$date[2]), ]) data <- matrix(unlist(beta.df[which(beta.df$date==beta.df$date[2]), ][ ,-1]), ncol=num.factors) model <- cv.glmnet(data, return, standardize=TRUE) op <- …
Ho alcuni dati limitati tra 0 e 1. Ho usato il betaregpacchetto in R per adattare un modello di regressione con i dati limitati come variabile dipendente. La mia domanda è: come interpretare i coefficienti dalla regressione?
Ho stimato la matrice di covarianza del campione di un campione e ho ottenuto una matrice simmetrica. Con , vorrei creare -variate rn distribuita normale, ma quindi ho bisogno la decomposizione di Cholesky di . Cosa devo fare se non è definito positivo?C n C CCCCCCCnnnCCCCCC
Sto utilizzando una serie temporale giornaliera di dati sulle vendite che contiene circa 2 anni di punti dati giornalieri. Sulla base di alcuni tutorial / esempi online ho cercato di identificare la stagionalità nei dati. Sembra che vi sia una periodicità / stagionalità settimanale, mensile e probabilmente annuale. Ad esempio, …
Il pacchetto di stampa R ggplot2 ha una fantastica funzione chiamata stat_smooth per tracciare una linea (o curva) di regressione con la banda di confidenza associata. Tuttavia, faccio fatica a capire esattamente come viene generata questa banda di confidenza, per ogni volta della linea di regressione (o "metodo"). Come posso …
Ho un database di eventi (cioè una variabile di date) e covariate associate. Gli eventi sono generati dal processo di Poisson non stazionario con il parametro che è una funzione sconosciuta (ma possibilmente lineare) di alcune covariate. Penso che il pacchetto NHPoisson esista proprio per questo scopo; ma dopo 15 …
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.