Statistiken und Big Data r

1

Hervorheben signifikanter Ergebnisse aus nicht parametrischen Mehrfachvergleichen auf Boxplots

Ich habe Box-Plots von 13 Gruppen, die ich in einem Plot zeige. Die Gruppen haben unausgeglichene Bevölkerungsgruppen und sind nicht normal verteilt. Ich möchte zeigen, welche Paare statistisch ähnlich sind (dh kruskal.test p-Wert <0,05 haben), indem ich a, b, c usw. auf übereinstimmende Kästchen setze. Hier ist ein Pseudocode, um …

9 r statistical-significance boxplot

1

Ist es möglich, ein lmer-Modell ohne feste Effekte anzugeben?

Wie spezifiziere ich in R ein früheres Modell ohne globalen festen Effekt? Zum Beispiel, wenn ich so etwas sage lmer(y ~ (1 | group) + (0 + x | group), data = my_df) das angepasste Modell wird sein yi j= a + αich+ βichxi jyij=a+αi+βixijy_{ij} = a + \alpha_{i} + …

9 r mixed-model lme4-nlme

3

Schätzen zufälliger Effekte und Anwenden einer benutzerdefinierten Korrelations- / Kovarianzstruktur mit dem Paket R lme4 oder nlme

Ich habe folgende Art von Daten. Ich habe 10 Personen bewertet, die jeweils 10 Mal wiederholt wurden. Ich habe eine 10x10 Beziehungsmatrix (Beziehung zwischen allen Kombinationen der Individuen). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld = rnorm(10, 5, 0.5)) Bei diesem Gen …

9 r mixed-model

1

Statistische Projektverzeichnisstruktur mit mehreren Sprachen (z. B. R und Splus)?

Aufbauend auf der Post Wie man ein statistisches Analyseprojekt und das ProjectTemplatePaket in R ... F: Wie erstellen Sie Ihre statistische Projektverzeichnisstruktur, wenn mehrere Sprachen stark vertreten sind (z. B. R AND Splus)? Die meisten Diskussionen zu diesem Thema beschränkten sich auf Projekte, die hauptsächlich eine Sprache verwenden. Ich beschäftige …

9 r project-management splus

1

Kontinuitätskorrektur für den Chi-Quadrat-Test von Pearson und McNemar

Ich denke, das ist eine ziemlich grundlegende Frage, aber mir ist nur klar, dass ich den Begriff Kontinuitätskorrektur nicht ganz verstehe. Ich benutze R und habe die gleiche Syntax correct=TRUEfür beide chisq.testund gefunden mcnemar.test. Beziehen sie sich auf unterschiedliche Kontinuitätskorrekturmethoden? Ich habe gehört, dass die Kontinuität des Yate für den …

9 r chi-squared yates-correction

2

Anpassung des exponentiellen Zerfalls mit negativen y-Werten

Ich versuche, eine exponentielle Abklingfunktion an y-Werte anzupassen, die bei hohen x-Werten negativ werden, kann meine nlsFunktion jedoch nicht richtig konfigurieren . Ziel Ich interessiere mich für die Steigung der Abklingfunktion (λλ\lambdanach einigen Quellen ). Wie ich diese Steigung erhalte, ist nicht wichtig, aber das Modell sollte so gut wie …

9 r model exponential

2

ANOVA zum Vergleichen von Modellen

Ich suche auf dieser Seite nach einem Workshop zu GAM in R: http://qcbs.ca/wiki/r_workshop8 Am Ende des Abschnitts 2. Multiple smooth termszeigen sie ein Beispiel, in anovadem drei verschiedene Modelle verglichen werden, um das am besten passende Modell zu ermitteln. Die Ausgabe ist Analysis of Deviance Table Model 1: y ~ …

9 r anova gam mgcv

1

So wählen Sie den Typ der GAM-Parameter

Ich habe begonnen, mit GAM in R zu arbeiten, und ich habe Simon Woods ausgezeichnetes Buch zum Thema ( "Generalisierte additive Modelle Eine Einführung mit R" ) erworben. Anhand eines seiner Beispiele betrachte ich Folgendes: library(mgcv) data(trees) ct1<-gam(log(Volume) ~ Height + s(Girth), data=trees) Ich habe zwei allgemeine Fragen zu diesem …

9 r regression machine-learning gam mgcv

1

Anpassen der negativen Binomialverteilung an Daten mit großer Anzahl

Ich habe ~ 1 Million Datenpunkte. Hier ist der Link zur Datei data.txt. Jeder von ihnen kann einen Wert zwischen 0 und 145 annehmen. Es handelt sich um einen diskreten Datensatz. Unten ist das Histogramm des Datensatzes. Auf der x-Achse ist die Zählung (0-145) und auf der y-Achse ist die …

9 r statistical-significance goodness-of-fit negative-binomial

1

Simulieren Sie aus einer Null-aufgeblasenen Poisson-Verteilung

Ich versuche aus beobachteten Daten zu simulieren, dass ich zu einem Poisson-Regressionsmodell ohne Inflation passt. Ich passe die Daten in R an, indem ich sie zeroinfl()aus dem Paket verwende pscl, aber ich habe Probleme herauszufinden, wie ich die ZIP-Verteilung aus den Koeffizientenschätzungen ableiten kann. Ich weiß, wie ich die vorhergesagten …

9 r simulation poisson-regression zero-inflation

1

Zufällige Waldregression - Residuen korrelieren mit der Reaktion

Ich versuche, die Random Forest-Regression zu verwenden. Ich habe eine Antwortvariable: y = rnorm(10000, mean=0, sd=3) Und ein paar Prädiktorvariablen (die nur die Antwort mit zusätzlichem Rauschen sind): x = data.frame(v1=y + rnorm(10000, mean=0, sd=3), v2=y + rnorm(10000, mean=0, sd=3), v3=y + rnorm(10000, mean=0, sd=3)) Ich baue den zufälligen Wald: …

9 r regression machine-learning random-forest ensemble

2

Ridge-Regression in R mit p-Werten und Anpassungsgüte

Gratregression in RI haben entdeckt linearRidgeim ridgePaket - das zu einem Modell passt, Koeffizienten und p-Werte angibt, aber nichts, um die allgemeine Anpassungsgüte zu messen lm.ridgeim MASSPaket - das Koeffizienten und GCV, aber keine p-Werte für Parameter meldet Wie kann ich all diese Dinge (Anpassungsgüte, Koeffizienten und p-Werte) aus derselben …

9 r p-value r-squared ridge-regression

2

R: Augmented Dickey Fuller (ADF) -Test

Ich habe ein Problem mit den Dickey-Fuller-p-Werten und der Teststatistik für den Einheitswurzeltest in R. Ich habe versucht, Funktionen zu verwenden: urca::ur.df() fUnitRoots::adfTest() tseries::adf.test() Alle zeigten unterschiedliche Ergebnisse für die gleichen Testeinstellungen (Verzögerung, Typ) im Vergleich zur gretl-Ausgabe. Zum Beispiel: set.seed(1) x <- rnorm(50, 0, 3) schwert.param <- trunc(12 * …

9 r time-series unit-root augmented-dickey-fuller

1

Wie kann man die Anzahl der Personen in einem Bild einer Menschenmenge messen?

Hintergrund : Israel (und der Nahe Osten im Allgemeinen) ist voller Proteste. Ich bin neugierig, wenn ich ein Bild bekomme, um zu schätzen, wie viele Personen sich auf diesem Bild befinden (oft ein Bild einer großen Menge). Welche Modellierung kann eine Lösung für dieses Problem bieten? (und natürlich kann es …

8 r image-processing pattern-recognition

2

Ursachen für bimodale Verteilungen beim Bootstrapping eines Metaanalysemodells

Ich helfe einem Kollegen, ein Metaanalyse-Modell mit gemischten Effekten mithilfe des von @Wolfgang verfassten Metafor R-Paket-Frameworks zu booten. Interessanterweise und besorgniserregend erhalte ich für einen der Koeffizienten des Modells beim Bootstrapping eine bimodale Verteilung (siehe unten rechts in der Abbildung unten). Ich denke, eine der Hauptursachen könnte die Tatsache sein, …

8 r mixed-model bootstrap meta-analysis bimodal

Als «r» getaggte Fragen