6
Gibt es eine Möglichkeit, einem qqplot Konfidenzintervalle hinzuzufügen?Hinzufügen von Konfidenzintervallen zu einem qq-Plot?
Ich habe einen Datensatz von Genexpressionswerten, die ich sichtbar gemacht habe mit PCA:
pca1 = prcomp(data, scale. = TRUE)
Ich suche jetzt Ausreißer durch die Verteilung der Daten gegen die Normalverteilung durch Überprüfung:
qqnorm(pca1$x,pch = 20, col = c(rep("red", 73), rep("blue", 33)))
qqline(pca1$x)
Das ist meine Daten:
data = [2,48 104 4,25 219 0,682 0,302 1,09 0,586 90,7 344 13,8 1,17 305 2,8 79,7 3,18 109 0,932 562 0,958 1,87 0,59 114 391 13,5 1,41 208 2,37 166 3,42]
würde ich Zeichnen Sie nun 95% -Konfidenzintervalle, um zu überprüfen, welche Datenpunkte außerhalb liegen. Irgendwelche Tipps, wie man das macht?
So mögen Sie Ihre Probenverteilung von der theoretischen Normalverteilung subtrahieren? Klingt wie das, was Sie tun möchten, ist 'nls', um Ihre Daten an eine normale dist-Funktion anzupassen und die Konfidenzdaten aus der Ausgabe von' nls' zu holen. –
Sie erhalten viel eher eine hilfreiche Antwort, wenn Sie einen [minimalen, reproduzierbaren Datensatz] bereitstellen (http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)/5963610 # 5963610) zusammen mit dem Code, den Sie ausprobiert haben. Vielen Dank! – Henrik
Ich habe meinen ersten Post mit einigen Daten bearbeitet. Ist es möglich, die Konfidenzdaten von der Ausgabe von qqnorm abzurufen? – user2846211