Ich habe einige EEG-Datensätze, die ich gegen zwei Klassen teste. Ich kann eine anständige Fehlerrate von LDA erhalten (die klassenbedingten Verteilungen sind nicht Gaußsch, haben aber ähnliche Schwänze und eine ausreichend gute Trennung), und deshalb möchte ich den ROC des LDA-Prädiktors gegen...
Dies ist eine Anfängerfrage zu einer Übung in Jim Alberts „Bayesian Computation with R“. Beachten Sie, dass dies zwar Hausaufgaben sein können, in meinem Fall jedoch nicht, da ich Bayes'sche Methoden in R lerne, weil ich denke, dass ich sie in meinen zukünftigen Analysen verwenden könnte. Obwohl...
Ich schreibe diese Frage unter Bezugnahme auf ein Beispiel auf Seite 138-142 des folgenden Dokuments: ftp://ftp.software.ibm.com/software/analytics/spss/documentation/amos/20.0/de/Manuals/IBM_SPSS_Amos_User_Guide.pdf . Hier sind illustrative Abbildungen und eine Tabelle: Ich verstehe, dass die...
Wenn ich diesen Code ausführe: require(nlme) a <- matrix(c(1,3,5,7,4,5,6,4,7,8,9)) b <- matrix(c(3,5,6,2,4,6,7,8,7,8,9)) res <- lm(a ~ b) print(summary(res)) res_gls <- gls(a ~ b) print(summary(res_gls)) Ich bekomme die gleichen Koeffizienten und die gleiche statistische Signifikanz für...
Ich möchte Leistungsanalysen für Hypothesentests der (Nicht-) Proportionsgleichheit durchführen, bei denen die Proportionen sehr klein sind. Ich möchte dies tun, ohne normale (oder Poisson-) Näherungen der Binomialverteilung zu verwenden. Es gibt verschiedene allgemeine Arten von Machtfragen, die...
Bitte erlauben Sie mir, eine grundlegende Frage zu stellen. Ich verstehe die Mechanik von Naive Bayes für diskrete Variablen und kann die Berechnungen "von Hand" wiederholen. (Code von HouseVotes84 ganz unten). Ich habe jedoch Schwierigkeiten zu sehen, wie die Mechanik für kontinuierliche...
Ich habe einen Parameter der zwischen . Nehmen wir an, ich kann ein Experiment durchführen und , wobei ein Standard-Gaußscher ist. Was ich brauche, ist eine Schätzung von θ, die 1) unvoreingenommen 2) fast sicher begrenzt ist. Voraussetzung (2) ist für mich
Ich versuche, eine Cox-Regression für einen Beispieldatensatz von 2.000.000 Zeilen wie folgt mit nur R auszuführen. Dies ist eine direkte Übersetzung eines PHREG in SAS. Die Stichprobe ist repräsentativ für die Struktur des Originaldatensatzes. ## library(survival) ### Replace 100000 by...
Ich arbeite an einem Modell, das auf einer hässlichen parametrisierten Funktion beruht, die als Kalibrierungsfunktion für einen Teil des Modells fungiert. Bei Verwendung einer Bayes'schen Einstellung muss ich nicht informative Prioritäten für die Parameter erhalten, die meine Funktion beschreiben....
Ich habe 2 Studien, die die Reaktion des Patienten auf dasselbe Medikament untersuchen. In Studie 1 wurden 10.000 über dem Hintergrund exprimierte Gene gefunden, von denen 500 differentiell exprimiert und als Signatur der Arzneimittelantwort bezeichnet werden. In Studie 2 wurden 1.000 Gene...
Ich benötige einen einzelnen unterstellten Datensatz (z. B. um einen Ländergruppen-Dummy aus den unterstellten Pro-Kopf-Einkommensdaten des Landes zu erstellen). R bietet ein Paketpaket zum Erstellen mehrerer unterstellter Daten (z. B. Amelia) und zum Kombinieren von Ergebnissen aus mehreren...
Ich möchte Stichprobengewichte in mein Quantil-Regressionsmodell aufnehmen, bin mir aber nicht sicher, wie ich das machen soll. Ich habe bereits mein Gewicht definiert. Hierbei handelt es sich um replizierte Gewichte, die bereits im Umfragedatensatz angegeben sind (berechnet im Umfragepaket):...
Ich selbst würde immer einen geometrischen Mittelwert verwenden, um einen logarithmischen Median zu schätzen. In der Industrie führt die Verwendung des Stichprobenmedians jedoch manchmal zu besseren Ergebnissen. Die Frage ist also, gibt es einen Grenzbereich / Punkt, ab dem der Stichprobenmedian...
Kennt jemand eine Methode zum Vergleichen zweier Variablen mit einem Chi-Quadrat-Test, wenn die Variablen aus verschiedenen Umfragen mit unterschiedlichen svydesign()Aussagen stammen? Ich möchte einen Unterschied in einer variablen Verteilung über zwei Wellen einer Umfrage testen, aber die...
Ich habe eine Clusterbildung von Koordinatenpunkten (Längen- und Breitengrad) durchgeführt und überraschende, nachteilige Ergebnisse aus Clustering-Kriterien für die optimale Anzahl von Clustern gefunden. Die Kriterien sind dem clusterCrit()Paket entnommen . Die Punkte, die ich auf einem Plot zu...
Ich habe einen täglichen Wetterdatensatz, der wenig überraschend einen sehr starken saisonalen Effekt hat. Ich habe ein ARIMA-Modell mit der Funktion auto.arima aus dem Prognosepaket an diesen Datensatz angepasst. Zu meiner Überraschung wendet die Funktion keine saisonalen Operationen an -...
Mein Design ist wie folgt. yyy ist Bernoulli Antwort x1x1x_1 ist eine stetige Variable x2x2x_2 ist eine kategoriale (Faktor-) Variable mit zwei Ebenen Das Experiment ist vollständig in Probanden. Das heißt, jedes Subjekt erhält jede Kombination von und .x 2x1x1x_1x2x2x_2 Dies ist eine...
Ich möchte die Standardfehler einer angepassten hyperbolischen Verteilung berechnen. In meiner Notation ist die Dichte gegeben durch Ich verwende das HyperbolicDistr-Paket in R. Ich schätze die Parameter über den folgenden Befehl:H.( l ; α , β, μ , δ)= α2- β2- -- -- -- -- -- -√2 α δK.1( δα2- β2- --...
Ich habe einige grundlegende Fragen zu PCA (Hauptkomponentenanalyse) und LDA (lineare Diskriminanzanalyse): In PCA gibt es eine Möglichkeit, den erklärten Varianzanteil zu berechnen. Ist es auch für LDA möglich? Wenn das so ist, wie? ldaEntspricht der von der Funktion (in der R MASS-Bibliothek)...
http://cran.r-project.org/web/packages/quadprog/quadprog.pdf Das R-Paket quadprogscheint in der Lage zu sein, das quadratische Programmierproblem nur zu lösen, wenn die Matrix positiv definit ist.DDD Es gibt jedoch einen Fall, in dem die Matrix nicht eindeutig positiv ist. sowieDDD...