rpartSoweit ich weiß, hilft das Argument cp für die Funktion dabei, den Baum auf die gleiche Weise wie die Argumente minsplit oder minbucket vorab zu bereinigen. Was ich nicht verstehe, ist, wie CP-Werte berechnet werden. Zum Beispiel df<-data.frame(x=c(1,2,3,3,3,4), y=as.factor(c(TRUE, TRUE,...
library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1...
Ich möchte Inzidenzraten berechnen, die entlang der Gefährdungsquoten dargestellt werden, um sowohl relative als auch absolute Risikomessgrößen darzustellen. Ich habe in anderen Studien gesehen, dass solche Inzidenzraten unter Verwendung von Poisson-Modellen mit Nachbeobachtungszeit im Modell als...
Wie berechnen Sie den entsprechenden Prior, wenn Sie den Messfehler eines Instruments haben? Dieser Absatz stammt aus Cressies Buch "Statistik für räumlich-zeitliche Daten": Es ist häufig der Fall, dass einige vorherige Informationen bezüglich der Messfehlervarianz verfügbar sind, so dass ein...
Betrachten Sie den folgenden R-Code: > data <- data.frame( a=c(NA,2,3,4,5,6),b=c(2.2,NA,6.1,8.3,10.2,12.13),c=c(4.2,7.9,NA,16.1,19.9,23)) > data a b c 1 NA 2.20 4.2 2 2 NA 7.9 3 3 6.10 NA 4 4 8.30 16.1 5 5 10.20 19.9 6 6 12.13 23.0 Wie Sie sehen, habe ich die Daten so konstruiert,...
Ich habe eine Verteilung von Proben mit einer kleinen Anzahl von Werten in jeder (weniger als ). Ich habe den Median für jede Stichprobe berechnet, den ich mit einem Modell vergleichen und die Differenz zwischen dem Modell und dem Median jeder Stichprobe ermitteln möchte. Um ein konsistentes...
Ich habe einige Probleme, gleichwertige Ergebnisse zwischen einem aovModell zwischen wiederholten Messungen und einem lmergemischten Modell zu erhalten. Meine Daten und mein Skript sehen wie folgt
Ich möchte die Faltung berechnen / bewerten G( x ) = ∫D.f( x - t ) ϕ ( t ) dt ,g(x)=∫Df(x−t)ϕ(t)dt,g(x)=\int_D f(x-t) \phi(t) dt, wobei eine Dichte und ist eine glatte Funktion mit kompaktem Träger . Die Faltung ist nicht in geschlossener Form verfügbar und ich muss sie numerisch integrieren....
Ich habe die Funktion 'polr' im MASS-Paket verwendet, um eine ordinale logistische Regression für eine ordinale kategoriale Antwortvariable mit 15 kontinuierlichen erklärenden Variablen auszuführen. Ich habe den Code (siehe unten) verwendet, um zu überprüfen, ob mein Modell die Annahme der...
Die nicht kontextbezogene Kurzversion Sei eine Zufallsvariable mit CDF yyyF.( ⋅ ) ≡ { θθ + ( 1 - θ ) × CDFlog-normal( ⋅ ; μ , σ) y = 0 y> 0F(⋅)≡{θ y = 0 θ+(1−θ)×CDFlog-normal(⋅;μ,σ) y > 0 F(\cdot) \equiv \cases{\theta & y = 0 \\ \theta + (1-\theta) \times...
In einem logistischen verallgemeinerten linearen gemischten Modell (Familie = Binomial) weiß ich nicht, wie die Varianz der zufälligen Effekte zu interpretieren ist: Random effects: Groups Name Variance Std.Dev. HOSPITAL (Intercept) 0.4295 0.6554 Number of obs: 2275, groups: HOSPITAL, 14 Wie...
Ich habe drei makroökonomische Variablen (IKS - Verbraucherstimmung, ER - Beschäftigungsquote, DGO - Bestellung langlebiger Güter) und habe Granger-Kausaltests in R durchgeführt. Ich weiß nicht wirklich, wie ich die Ergebnisse eines Granger-Tests interpretieren soll. Könnte mir jemand helfen, die...
Ich versuche, den durchschnittlichen Behandlungseffekt anhand von Beobachtungsdaten mithilfe der Neigungsbewertung (speziell IPTW) abzuschätzen. Ich denke, ich berechne die ATE korrekt, aber ich weiß nicht, wie ich das Konfidenzintervall der ATE unter Berücksichtigung der inversen...
Ich habe Daten aus dem folgende experimentellen Design: meine Beobachtungen zählen die Anzahl der Erfolge sind ( K) aus Anzahl der Versuche von entsprechenden ( N), gemessen für zwei Gruppen von jeweils umfassten IIndividuen, von TBehandlungen, wobei in jeder solchen Faktorkombination gibt es...
Ich habe den folgenden einfachen Datensatz mit zwei kontinuierlichen Variablen; dh: d = data.frame(x=runif(100,0,100),y = runif(100,0,100)) plot(d$x,d$y) abline(lm(y~x,d), col="red") cor(d$x,d$y) # = 0.2135273 Ich muss die Daten so umordnen, dass die Korrelation zwischen Variablen ~ 0,6 beträgt....
Ich habe versucht, die Bayes'schen Statistiken neu zu lernen (jedes Mal, wenn ich dachte, ich hätte sie endlich bekommen, taucht etwas anderes auf, das ich vorher nicht in Betracht gezogen habe ...), aber es war (für mich) nicht klar, wie der Datengenerierungsprozess ablief im Bayesianischen Rahmen...
Bei der Auswahl einer geeigneten Anzahl von Knoten für ein GAM sollte möglicherweise die Anzahl der Daten und Inkremente auf der x-Achse berücksichtigt werden. Was ist, wenn wir 100 Inkremente auf der x-Achse mit 1000 Datenpunkten bei jedem Inkrement haben? Die Info hier sagt: Wenn sie nicht...
Warum wird bei der Schätzung von Parametern mit MAP geschrieben, dass wir den "Modus" schätzen? Ich dachte, es wäre der Mittelwert der posterioren
Ich suche nach einem Paket, das mir hilft, einige quadratische Optimierungsprobleme zu lösen, und ich sehe, dass es mindestens ein halbes Dutzend verschiedene Pakete gibt. Laut dieser Seite: QP (Quadratic Programming, 90C20): cplexAPI , Kernlab , limSolve , LowRankQP , Quadprog , Rcplex , Rmosek...
Meine Rohdaten bestehen aus einer 60-Tage-Zeitreihe mit einem Abwärtstrend. Die Daten sind wöchentlich, daher wird die Frequenz auf 7 eingestellt. Ich habe die Differenz der Daten berechnet, die so aussieht Wenn ich ACF- und PACF-Diagramme über den Unterschied durchführe, erhalte ich scheinbar...