Ich habe PCA mit 25 Variablen ausgeführt und die 7 besten PCs mit ausgewählt prcomp. prc <- prcomp(pollutions, center=T, scale=T, retx=T) Ich habe dann Varimax-Rotation an diesen Komponenten durchgeführt. varimax7 <- varimax(prc$rotation[,1:7]) Und jetzt möchte ich die PCA-gedrehten Daten...
Bearbeiten Ich habe ein Dokument gefunden , das genau das beschreibt, was ich brauche. Der einzige Unterschied besteht darin, dass das Papier die Monatsmittelwerte in Tagesdaten interpoliert, während die Monatsmittelwerte beibehalten werden. Ich habe Probleme, den Ansatz zu implementieren R....
Ich verwende R Lavaan-Paket , um ein Strukturgleichungsmodell zu schätzen. Angenommen, das Modell besteht aus 1 endogenen Manifestvariablen mit 1 latenten und 2 manifest erklärenden Variablen: group = {0,1} attitude1 = latent,scale age = respondent's age Das gewünschte Lavamodell ist dann...
Ich möchte ausgewählte Modelle mit First, Lasso und elastischem Netz vergleichen. Die folgende Abbildung zeigt die Koeffizientenpfade mit allen drei Methoden: Grat (Abb. A, alpha = 0), Lasso (Abb. B; alpha = 1) und elastisches Netz (Abb. C; alpha = 0,5). Die optimale Lösung hängt vom ausgewählten...
Ich habe eine Menge von Werten und y, die theoretisch exponentiell zusammenhängen:xxxyyy y=axby=axby = ax^b Eine Möglichkeit, die Koeffizienten zu erhalten, besteht darin, natürliche Logarithmen auf beiden Seiten anzuwenden und ein lineares Modell zu erstellen: > fit <- lm(log(y)~log(x)) >...
Bei annähernd normal verteilten Daten sind Boxplots eine großartige Möglichkeit, den Median und die Verbreitung der Daten sowie das Vorhandensein von Ausreißern schnell zu visualisieren. Bei stärker schwanzförmigen Verteilungen werden jedoch viele Punkte als Ausreißer angezeigt, da Ausreißer als...
Ich habe Ergebnisse aus dem gleichen Test auf zwei unabhängige Proben angewendet: x <- c(17, 12, 13, 16, 9, 19, 21, 12, 18, 17) y <- c(10, 6, 15, 9, 8, 11, 8, 16, 13, 7, 5, 14) Und ich möchte einen Wilcoxon-Rang-Summen-Test berechnen. Wenn ich die Statistik von Hand berechne , erhalte...
Gibt es Literatur, die die Wahl der Minibatch-Größe bei der Durchführung eines stochastischen Gradientenabfalls untersucht? Nach meiner Erfahrung scheint es sich um eine empirische Entscheidung zu handeln, die normalerweise durch Kreuzvalidierung oder unter Verwendung unterschiedlicher Faustregeln...
Ich muss die Grundlinien-Gefährdungsfunktion in einem zeitabhängigen Cox-Modell schätzenλ0( t )λ0(t)\lambda_0(t) λ ( t ) = λ0( t ) exp( Z.( t )'β)λ(t)=λ0(t)exp(Z.(t)'β)\lambda(t) = \lambda_0(t) \exp(Z(t)'\beta) Während ich den Überlebenskurs belegte, erinnere ich mich, dass die direkte Ableitung...
Ich habe ein Glm.nb von durchgeführt glm1<-glm.nb(x~factor(group)) wobei group eine kategoriale und x eine metrische Variable ist. Wenn ich versuche, eine Zusammenfassung der Ergebnisse zu erhalten, werden geringfügig unterschiedliche Ergebnisse angezeigt, je nachdem, ob ich summary()oder...
Ist es möglich, eine Leistungsanalyse für den Kruskal-Wallis- und den Mann-Whitney-U-Test durchzuführen? Wenn ja, gibt es R-Pakete / -Funktionen, die dies
Ich habe einige Erklärungen zu den Eigenschaften von linearen und nichtlinearen Modellen gelesen, bin mir aber manchmal nicht sicher, ob es sich bei dem vorliegenden Modell um ein lineares oder ein nichtlineares Modell handelt. Ist beispielsweise das folgende Modell linear oder
Ich arbeite an einem Logistikmodell und habe einige Schwierigkeiten, die Ergebnisse zu bewerten. Mein Modell ist ein Binomial Logit. Meine erklärenden Variablen sind: eine kategoriale Variable mit 15 Ebenen, eine dichotome Variable und 2 stetige Variablen. Mein N ist groß> 8000. Ich versuche,...
Mir ist klar, dass dies wahrscheinlich eine sehr einfache Frage ist, aber nach der Suche kann ich die gesuchte Antwort nicht finden. Ich habe ein Problem, bei dem ich die Variablen standardisieren muss, die die (Ridge-Regression) ausführen, um die Ridge-Schätzungen der Betas zu berechnen. Ich...
Ich habe einige Daten, die ich glatt benutze loess. Ich möchte die Wendepunkte der geglätteten Linie finden. Ist das möglich? Ich bin mir sicher, dass jemand eine ausgefallene Methode entwickelt hat, um dieses Problem zu lösen. Ich meine, schließlich ist es R! Ich kann die Glättungsfunktion, die...
Ich möchte die Imputation verwenden, um fehlende Werte in meinem Datensatz unter bestimmten Bedingungen zu ersetzen. Zum Beispiel möchte ich, dass die unterstellte Variable x1größer oder gleich der Summe meiner beiden anderen Variablen ist, sagen wir x2und x3. Ich möchte x3auch von entweder 0oder...
In R (2.15.2) habe ich einmal eine ARIMA (3,1,3) auf eine Zeitreihe und einmal eine ARMA (3,3) auf die einmal differenzierten Zeitreihen gepasst. Die angepassten Parameter unterscheiden sich, was ich der Anpassungsmethode in ARIMA zuschrieb. Auch das Anpassen einer ARIMA (3,0,3) an dieselben Daten...
Ich habe eine multiple Regression durchgeführt, bei der das Modell als Ganzes signifikant ist und ungefähr 13% der Varianz erklärt. Ich muss jedoch den Betrag der Varianz finden, der von jedem signifikanten Prädiktor erklärt wird. Wie kann ich das mit R machen? Hier sind einige Beispieldaten und...
Angenommen, X0, X1, … , XnX0,X1,…,Xn X_{0},X_{1},\ldots,X_{n} sind Zufallsvariablen, die der Poisson-Verteilung mit dem Mittelwert folgen . Wie kann ich nachweisen, dass es keinen unvoreingenommenen Schätzer für die Menge ?λλ \lambda 1λ1λ \dfrac{1}{\lambda}
Ich benutze die R "Multcomp" -Bibliothek ( http://cran.r-project.org/web/packages/multcomp/ ), um Dunnetts Test zu berechnen. Ich benutze das folgende Skript: Group <- factor(c("A","A","B","B","B","C","C","C","D","D","D","E","E","F","F","F")) Value <-