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`.
È possibile eseguire un'analisi della potenza per un test Kolmogorov Smirnov su 2 lati in R? Sto testando se due distribuzioni empiriche differiscono usando ks.test () e sto cercando di aggiungere un'analisi di potenza. Non sono stato in grado di trovare analisi di potenza integrate per i test KS in …
Sto usando la formula del test ab bayesiano per calcolare i risultati del test AB usando la metodologia bayesiana. Pr ( pB> pUN) = ∑i = 0αB- 1B ( αUN+ i , βB+ βUN)( βB+ i ) B ( 1 + i , βB) B ( αUN, βUN)Pr(pB>pA)=∑i=0αB−1B(αA+i,βB+βA)(βB+i)B(1+i,βB)B(αA,βA) \Pr(p_B > …
Sto analizzando alcuni dati comportamentali usando lme4in R, principalmente seguendo gli eccellenti tutorial di Bodo Winter , ma non capisco se sto gestendo correttamente le interazioni. Peggio ancora, nessun altro coinvolto in questa ricerca utilizza modelli misti, quindi sono un po 'alla deriva quando si tratta di assicurarsi che le …
Sto eseguendo l'analisi delle serie storiche usando R. Devo scomporre i miei dati in componenti di tendenza, stagionali e casuali. Ho dati settimanali per 3 anni. Ho trovato due funzioni in R - stl()e decompose(). Ho letto che stl()non va bene per la decomposizione moltiplicativa. Qualcuno può dirmi in quale …
Dalla Quick-R di Robert Kabacoff che ho # Bootstrap 95% CI for regression coefficients library(boot) # function to obtain regression weights bs <- function(formula, data, indices) { d <- data[indices,] # allows boot to select sample fit <- lm(formula, data=d) return(coef(fit)) } # bootstrapping with 1000 replications results <- boot(data=mtcars, …
Confronto di due campioni di proporzioni, stima della dimensione del campione: R vs Stata Ho ottenuto risultati diversi per le dimensioni del campione, come segue: In R power.prop.test(p1 = 0.70, p2 = 0.85, power = 0.90, sig.level = 0.05) Risultato: (quindi 161) per ciascun gruppo.n = 160.7777n=160.7777n = 160.7777 In …
Non ho trovato alcuna documentazione sull'applicazione di Random Forests a MNIST, CIFAR, STL-10, ecc. Quindi ho pensato di provarli con il MNIST invariante alla permutazione . In R , ho provato: randomForest(train$x, factor(train$y), test$x, factor(test$y), ntree=500) Questo ha funzionato per 2 ore e ha avuto un errore di prova del …
Nel codice seguente eseguo una regressione logistica su dati raggruppati usando glm e "a mano" usando mle2. Perché la funzione logLik in R mi dà una probabilità di log logLik (fit.glm) = - 2.336 che è diverso da quello logLik (fit.ml) = - 5.514 che ottengo a mano? library(bbmle) #successes …
Sto cercando di usare SMOTE per correggere lo squilibrio nel mio problema di classificazione multi-classe. Sebbene SMOTE funzioni perfettamente sul set di dati dell'iride secondo il documento della guida di SMOTE, non funziona su un set di dati simile. Ecco come appaiono i miei dati. Nota che ha tre classi …
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: …
Sto riscontrando problemi con la definizione dei termini di errore per una misura ripetuta a due vie ANOVA in R. I miei dati consistono in stime della densità del legno per tre posizioni radiali (interna, media ed esterna) lungo un nucleo estratto da un albero. Ci sono un totale di …
Mi sento sopraffatto dopo aver tentato di approfondire la letteratura su come eseguire la mia analisi di modelli misti seguendola con l'utilizzo di AIC per selezionare il modello oi modelli migliori. Non credo che i miei dati siano così complicati, ma sto cercando conferma che ciò che ho fatto sia …
Cosa succede quando non hai idea della distribuzione dei parametri? Quale approccio dovremmo usare? Il più delle volte ci proponiamo di sottovalutare se una determinata variabile ha qualche influenza sulla presenza / assenza di una determinata specie e la variabile è accettata o meno in base all'importanza della variabile. Ciò …
Le mie domande riguardano i GAM nel pacchetto mgcv R. A causa delle dimensioni ridotte del campione, desidero determinare l'errore di previsione utilizzando la convalida incrociata senza esclusione di dati. È ragionevole? Esiste un pacchetto o un codice come posso farlo? La errorest()funzione nel pacchetto ipred non funziona. Un semplice …
Nel seguente esempio ho una cornice di dati che consiste in una serie temporale di misurazioni della temperatura dell'acqua registrate a 5 profondità nell'oceano in cui ciascun valore in Tempcorrisponde alla data in DateTimee alla profondità in Depth. set.seed(1) Temp <- rnorm(43800,sd=20) AirT <- rnorm(8760,sd=20) Depth <- c(1:5) DateTime = …
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.