Betrachten Sie den folgenden R-Code: example <- function(n) { X <- 1:n Y <- rep(1,n) return(lm(Y~X)) } #(2.13.0, i386-pc-mingw32) summary(example(7)) #R^2 = .1963 summary(example(62)) #R^2 = .4529 summary(example(4540)) #R^2 = .7832 summary(example(104))) #R^2 = 0 #I did a search for n...
Ich habe viele Probleme mit einem Datensatz, auf den ich SEM anwenden möchte. Wir nehmen die Existenz von 5 latenten Faktoren A, B, C, D, E mit Indikatoren bzw. A1 bis A5 (geordnete Faktoren), B1 bis B3 (quantitativ), C1, D1, E1 (alle letzten drei geordneten Faktoren mit nur 2 Ebenen für E1. Wir...
Ich modelliere einige Daten, bei denen ich denke, dass ich zwei gekreuzte zufällige Effekte habe. Der Datensatz ist jedoch nicht ausgewogen, und ich bin mir nicht sicher, was getan werden muss, um dies zu berücksichtigen. Meine Daten sind eine Reihe von Ereignissen. Ein Ereignis tritt auf, wenn...
In den meisten Situationen beschäftigen wir uns nur mit einer Ergebnis- / Antwortvariablen wie . In einigen Szenarien, insbesondere in den klinischen Daten, können die Ergebnisvariablen jedoch hochdimensional / multivariat sein. Beispiel: , wobei die , und enthält und diese Ergebnisse alle...
In der Ökologie verwenden wir häufig die logistische Wachstumsgleichung: N.t= K.N.0er tK.+ N.0er t - 1Nt=KN0ertK+N0ert−1 N_t = \frac{ K N_0 e^{rt} }{K + N_0 e^{rt-1}} oder N.t= K.N.0N.0+ ( K.- N.0) e- r tNt=KN0N0+(K−N0)e−rt N_t = \frac{ K N_0}{N_0 + (K -N_0)e^{-rt}} wobei die Tragfähigkeit ist...
Ich habe einige Zehntausende von Beobachtungen, die in einer Zeitreihe liegen, aber nach Orten gruppiert sind. Zum Beispiel: location date observationA observationB --------------------------------------- A 1-2010 22 12 A 2-2010 26 15 A 3-2010 45 16 A 4-2010 46 27 B 1-2010 167 48 B 2-2010 134 56 B...
Dies ist nur ein Beispiel, auf das ich mehrmals gestoßen bin, daher habe ich keine Beispieldaten. Ausführen eines linearen Regressionsmodells in R: a.lm = lm(Y ~ x1 + x2) x1ist eine stetige Variable. x2ist kategorisch und hat drei Werte, z. B. "Niedrig", "Mittel" und "Hoch". Die von R gegebene...
Angenommen, ich habe 'Experten', von denen ich eine vorherige Verteilung auf eine Variable auslösen möchte . Ich möchte sie mit echtem Geld motivieren . Die Idee ist, die Priors zu entlocken, Realisierungen der Zufallsvariablen und dann eine vorgegebene "Geldbörse" unter den Experten aufzuteilen,...
Angenommen, wir haben den folgenden Datensatz: Men Women Dieting 10 30 Non-dieting 5 60 Wenn ich den Fisher-Exakt-Test in R durchführe, was bedeutet dann alternative = greater(oder weniger)? Beispielsweise: mat = matrix(c(10,5,30,60), 2,2) fisher.test(mat, alternative="greater") Ich bekomme das...
Ich habe eine Drei-Wege-ANOVA mit wiederholten Messungen durchgeführt. Welche Post-hoc-Analysen sind gültig? Dies ist ein vollständig ausgewogenes Design (2x2x2), bei dem einer der Faktoren eine wiederholte Messung innerhalb der Probanden aufweist. Ich kenne multivariate Ansätze für ANOVA mit...
Es gibt viele Methoden zur Durchführung der Regularisierung - beispielsweise die Regularisierung , und . Laut Friedman Hastie & Tibsharani hängt der beste Regularisierer vom Problem ab: nämlich der Art der tatsächlichen Zielfunktion, der jeweils verwendeten Basis, dem Signal-Rausch-Verhältnis...
Wenn ich meine Daten mit etwas wie lm(y~a*b)in der R-Syntax anpasse, wobei aes sich um eine binäre Variable und beine numerische Variable handelt, ist der a:bInteraktionsterm die Differenz zwischen der Steigung von y~bat a= 0 und at a= 1. Nehmen wir nun an, die Beziehung zwischen yund bist...
Die Statistik des Wahrscheinlichkeitsverhältnisses (auch bekannt als Abweichung) und der Test auf mangelnde Anpassung (oder Anpassungsgüte) sind für ein logistisches Regressionsmodell (Anpassung unter Verwendung der Funktion) in R ziemlich einfach zu erhalten . Dies kann jedoch sein Es ist leicht,...
Ist es für eine lineare Regression mit mehreren Gruppen (natürliche Gruppen, die a priori definiert wurden) akzeptabel, zwei verschiedene Modelle mit demselben Datensatz auszuführen, um die folgenden zwei Fragen zu beantworten? Hat jede Gruppe eine Steigung ungleich Null und einen Achsenabschnitt...
Ich versuche herauszufinden, ob einfache Wahrscheinlichkeiten für mein Problem funktionieren oder ob es besser ist, komplexere Methoden wie die logistische Regression zu verwenden (und etwas darüber zu lernen). Die Antwortvariable in diesem Problem ist eine binäre Antwort (0, 1). Ich habe eine...
Ich versuche, die Arbeit eines Kollegen zu replizieren, und verschiebe die Analyse von Stata nach R. Die von ihr verwendeten Modelle rufen die Option "Cluster" innerhalb der Funktion nbreg auf, um die Standardfehler zu gruppieren. Unter http://repec.org/usug2007/crse.pdf finden Sie eine ziemlich...
Ich versuche, die Vorlaufverzögerung zwischen Zeitreihen zweier Aktienkurse zu analysieren. In der regelmäßigen Zeitreihenanalyse können wir Cross Correlaton, VECM (Granger Causality) durchführen. Wie geht man jedoch in unregelmäßig verteilten Zeitreihen damit um? Die Hypothese ist, dass eines der...
Ich möchte zwei verschiedene Zeitvariablen modellieren, von denen einige in meinen Daten stark kollinear sind (Alter + Kohorte = Periode). Dabei hatte ich einige Probleme mit lmerund und Interaktionen von poly(), aber es ist wahrscheinlich nicht darauf beschränkt lmer, dass ich mit nlmeIIRC die...
Das Bild unten stammt aus diesem Artikel in Psychological Science . Ein Kollege wies auf zwei ungewöhnliche Dinge hin: Gemäß der Beschriftung zeigen die Fehlerbalken "± 2,04 Standardfehler, das 95% -Konfidenzintervall". Ich habe bisher nur ± 1,96 SE für den 95% CI gesehen, und ich kann nichts über...
Ich mache eine Einweg-ANOVA (pro Art) mit benutzerdefinierten Kontrasten. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 wo ich Intensität 0,5 gegen 5, 5 gegen 12,5 und so weiter vergleiche. Dies sind die Daten, an denen ich arbeite mit den folgenden Ergebnissen...