Statistiche e Big Data r

5

Esempio di forte coefficiente di correlazione con un alto valore p

Mi chiedevo, è possibile avere un coefficiente di correlazione molto forte (diciamo .9 o superiore), con un valore di p elevato (diciamo .25 o superiore)? Ecco un esempio di un basso coefficiente di correlazione, con un valore p elevato: set.seed(10) y <- rnorm(100) x <- rnorm(100)+.1*y cor.test(x,y) cor = 0,03908927, …

21 r hypothesis-testing correlation

1

Regressione logistica per serie storiche

Vorrei utilizzare un modello di regressione logistica binaria nel contesto dei dati di streaming (serie temporali multidimensionali) al fine di prevedere il valore della variabile dipendente dei dati (ovvero riga) appena arrivati, date le osservazioni passate. Per quanto ne so, la regressione logistica viene tradizionalmente utilizzata per l'analisi post mortem, …

21 r time-series logistic

2

Come usare i pesi nella funzione lm in R?

Bloccato . Questa domanda e le sue risposte sono bloccate perché la domanda è fuori tema ma ha un significato storico. Al momento non accetta nuove risposte o interazioni. Qualcuno potrebbe offrire alcuni suggerimenti su come utilizzare l' weightsargomento nella lmfunzione di R ? Ad esempio, stavi cercando di adattare …

21 r regression

1

Generazione di variabili casuali binomiali correlate

Mi chiedevo se fosse possibile generare variabili binomiali casuali correlate seguendo un approccio di trasformazione lineare? Di seguito, ho provato qualcosa di semplice in R e produce una certa correlazione. Ma mi chiedevo se esiste un modo di principio per farlo? X1 = rbinom(1e4, 6, .5) ; X2 = rbinom(1e4, …

21 r correlation binomial random-generation bernoulli-distribution

2

Come funziona il metodo di trasformazione inversa?

Come funziona il metodo di inversione? Supponiamo di avere un campione casuale con densità over e quindi con cdf su . Quindi con il metodo di inversione ottengo la distribuzione di X come F_X ^ {- 1} (u) = u ^ \ theta . f ( x ; θ ) …

21 r distributions inference random-generation inverse-cdf

1

Perché il quasi-Poisson nella GLM non è trattato come un caso speciale di binomio negativo?

Sto cercando di adattare modelli lineari generalizzati ad alcune serie di dati di conteggio che potrebbero essere o meno sovradispersi. Le due distribuzioni canoniche che si applicano qui sono Poisson e Negative Binomial (Negbin), con EV e varianzaμμ\mu Va rP= μVun'rP=μVar_P = \mu Va rNB= μ + μ2θVun'rNB=μ+μ2θVar_{NB} = \mu …

21 r generalized-linear-model negative-binomial poisson-regression quasi-likelihood

4

Come testare se la mia distribuzione è multimodale?

Quando tracciamo un istogramma dei miei dati, ha due picchi: Significa una potenziale distribuzione multimodale? Ho eseguito in dip.testin R ( library(diptest)) e l'output è: D = 0.0275, p-value = 0.7913 Posso concludere che i miei dati hanno una distribuzione multimodale? DATI 10346 13698 13894 19854 28066 26620 27066 16658 …

21 r hypothesis-testing distributions self-study histogram

3

Come calcolare la bontà di adattamento in glm (R)

Ho il seguente risultato dall'esecuzione della funzione glm. Come posso interpretare i seguenti valori: Deviazione nulla Devianza residua AIC Hanno qualcosa a che fare con la bontà di adattarsi? Posso calcolare la bontà della misura di adattamento da questi risultati come R-quadrato o qualsiasi altra misura? Call: glm(formula = tmpData$Y …

21 r regression generalized-linear-model

1

Ruolo del parametro n.minobsinnode di GBM in R [chiuso]

È improbabile che questa domanda aiuti eventuali futuri visitatori; è rilevante solo per una piccola area geografica, un momento specifico nel tempo o una situazione straordinariamente stretta che non è generalmente applicabile al pubblico mondiale di Internet. Per assistenza nel rendere questa domanda più ampiamente applicabile, visitare il centro assistenza …

21 r gbm

3

Perché nls () mi dà errori di "matrice gradiente singolare alle stime iniziali dei parametri"?

Ho alcuni dati di base sulla riduzione delle emissioni e sul costo per auto: q24 <- read.table(text = "reductions cost.per.car 50 45 55 55 60 62 65 70 70 80 75 90 80 100 85 200 90 375 95 600 ",header = TRUE, sep = "") So che questa è …

21 r self-study exponential starting-values

4

Come proiettare un nuovo vettore nello spazio PCA?

Dopo aver eseguito l'analisi dei componenti principali (PCA), voglio proiettare un nuovo vettore nello spazio PCA (ovvero trovare le sue coordinate nel sistema di coordinate PCA). Ho calcolato PCA in linguaggio R utilizzando prcomp. Ora dovrei essere in grado di moltiplicare il mio vettore per la matrice di rotazione PCA. …

21 r pca r variance heteroscedasticity misspecification distributions time-series data-visualization modeling histogram kolmogorov-smirnov negative-binomial likelihood-ratio econometrics panel-data categorical-data scales survey distributions pdf histogram correlation algorithms r gpu parallel-computing approximation mean median references sample-size normality-assumption central-limit-theorem rule-of-thumb confidence-interval estimation mixed-model psychometrics random-effects-model hypothesis-testing sample-size dataset large-data regression standard-deviation variance approximation hypothesis-testing variance central-limit-theorem kernel-trick kernel-smoothing error sampling hypothesis-testing normality-assumption philosophical confidence-interval modeling model-selection experiment-design hypothesis-testing statistical-significance power asymptotics information-retrieval anova multiple-comparisons ancova classification clustering factor-analysis psychometrics r sampling expectation-maximization markov-process r data-visualization correlation regression statistical-significance degrees-of-freedom experiment-design r regression curve-fitting change-point loess machine-learning classification self-study monte-carlo markov-process references mathematical-statistics data-visualization python cart boosting regression classification robust cart survey binomial psychometrics likert psychology asymptotics multinomial

2

Avviso "Il modello non è riuscito a convergere" in lmer ()

Con il seguente set di dati, volevo vedere se la risposta (effetto) cambia in relazione a siti, stagione, durata e loro interazioni. Alcuni forum online sulle statistiche mi hanno suggerito di continuare con i modelli lineari a effetti misti, ma il problema è che, poiché i replicati sono randomizzati all'interno …

21 r mixed-model lme4-nlme

4

Come creare una matrice di covarianza arbitraria

Ad esempio, in R, la MASS::mvrnorm()funzione è utile per generare dati per dimostrare varie cose nelle statistiche. Prende un Sigmaargomento obbligatorio che è una matrice simmetrica che specifica la matrice di covarianza delle variabili. Come potrei creare una matrice simmetrica con voci arbitrarie?n×nn×nn\times n

21 r random-generation covariance-matrix

1

lme () e lmer () che danno risultati contrastanti

Ho lavorato con alcuni dati che hanno alcuni problemi con misurazioni ripetute. In tal modo ho notato un comportamento molto diverso tra lme()e lmer()utilizzando i miei dati di test e voglio sapere perché. Il set di dati falsi che ho creato ha misurazioni di altezza e peso per 10 soggetti, …

20 r mixed-model lme4-nlme

3

Trovare un modo per simulare numeri casuali per questa distribuzione

Sto cercando di scrivere un programma in R che simula numeri pseudo casuali da una distribuzione con la funzione di distribuzione cumulativa: F(x)=1−exp(−ax−bp+1xp+1),x≥0F(x)=1−exp⁡(−ax−bp+1xp+1),x≥0F(x)= 1-\exp \left(-ax-\frac{b}{p+1}x^{p+1}\right), \quad x \geq 0 dove a,b>0,p∈(0,1)a,b>0,p∈(0,1)a,b>0, p \in (0,1) Ho provato il campionamento della trasformata inversa ma l'inverso non sembra risolvibile dal punto di vista …

20 r random-generation

Domande taggate «r»