Wie schalte ich die numerische Ausgabe von boxplots (mit plot=FALSE) in etwas brauchbar?

Ich bin erfolgreich mit dem boxplot Funktion zu generieren... boxplots. Nun muss ich das generieren von Tabellen mit Statistiken, die boxplot berechnet, um plots.

Dafür benutze ich die plot=FALSE option.

Das problem ist, dass dieses produziert Daten in eine eher seltsame format, dass ich einfach nichts tun kann mit. Hier ist ein Beispiel:

structure(list(stats = structure(c(178.998262143545, 182.227431564442, 
202.108456373209, 220.375358994654, 221.990406228232, 216.59986775699, 
217.054997032148, 228.509462713206, 267.070720949859, 284.832378859975, 
189.864120937198, 201.876421960518, 219.525439081472, 234.260088973545, 
279.343359793024, 209.472617639903, 209.526516071858, 214.785213079737, 
230.027361556731, 240.0647114578, 202.057148813419, 207.375619207685, 
220.093663781351, 226.246698737471, 240.343646265795), .Dim = c(5L, 
5L)), n = c(4, 6, 8, 4, 8), conf = structure(c(171.971593703341, 
232.245319043076, 196.247705331772, 260.771220094641, 201.435457751239, 
237.615420411705, 198.589545146688, 230.980881012787, 209.552007821332, 
230.635319741371), .Dim = c(2L, 5L)), out = numeric(0), group = numeric(0), 
names = c("U", "UM", "M", "LM", "L")), .Names = c("stats", "n", "conf", "out", "group", 
"names"))

Was ich will ist eine Tabelle, die für jeden von der Statistik -- min, max, median und die Quartile -- und deren Werte für jede Gruppe (die, die im "Namen").

Könnte jemand mir eine hand mit diesem? Ich bin sehr viel ein R Anfänger.

Vielen Dank im Voraus!

  • Gil, benötigen Sie die entsprechenden zahlen an das Ende der Bärte oder die minimum und maximum der Daten?
InformationsquelleAutor Gil Williams | 2012-01-13
  1. 23

    boxplot gibt eine Struktur in R genannt list.

    Einer Liste ist mehr oder weniger ein Daten-container, wo Sie können, finden Sie Elemente nach name.
    Wenn Sie A <- boxplot(...) können Sie den Zugriff auf die names mit A$names, die conf mit A$conf usw.

    So, Blick auf die boxplot helpfile ?boxplot unter Value: (die Ihnen sagt, was boxplot gibt), sehen wir, dass es gibt eine Liste mit den folgenden Komponenten:

       stats: a matrix, each column contains the extreme of the lower
          whisker, the lower hinge, the median, the upper hinge and the
          extreme of the upper whisker for one group/plot.  If all the
          inputs have the same class attribute, so will this component.
       n: a vector with the number of observations in each group.    
    conf: a matrix where each column contains the lower and upper
          extremes of the notch.    
     out: the values of any data points which lie beyond the extremes
          of the whiskers.    
   group: a vector of the same length as ‘out’ whose elements indicate
          to which group the outlier belongs.    
   names: a vector of names for the groups.

    Damit die Tabelle für jeden der stats ist in A$stats jede Spalte gehört zu einer Gruppe und enthält Minimum, unteres Quartil, median, Oberes Quartil und max.

    Könnten Sie tun:

    A <- boxplot(...)
mytable <- A$stats
colnames(mytable)<-A$names
rownames(mytable)<-c('min','lower quartile','median','upper quartile','max')
mytable

    gibt (für mytable):

                          U       UM        M       LM        L
min            178.9983 216.5999 189.8641 209.4726 202.0571
lower quartile 182.2274 217.0550 201.8764 209.5265 207.3756
median         202.1085 228.5095 219.5254 214.7852 220.0937
upper quartile 220.3754 267.0707 234.2601 230.0274 226.2467
max            221.9904 284.8324 279.3434 240.0647 240.3436

    Dann finden Sie auf es wie mytable['min','U'].

    • Erstaunlich -- eine Antwort verdient, eine gut geschriebene Bedienungsanleitung in wenigen Minuten! Deine Lösung funktioniert perfekt. Dank einer million!
    • Aus irgendeinem Grund, ich brauchte, um row.names statt rownames
    InformationsquelleAutor mathematical.coffee
  2. 4

    Wenn Sie wirklich wollen, Quantile Ihrer Daten statt boxplot zahlen, mit quantile direkt wäre meine Wahl (ist es viel einfacher zu Lesen, wenn man sich durch das, was Sie Taten später).

    quantile (x, probs = c (0, .25, .5,.75, 1))

    quantile sich nicht für die Arbeit mit Gruppen, aber Sie können kombinieren Sie es mit aggregate so heißt es für jede der Gruppen gegeben in argument by (muss eine Liste, so können Sie kombinieren hier verschiedene Gruppierung der Faktoren):

    aggregate (chondro$x, by = list (chondro$clusters), 
           FUN = quantile, probs = c (0, .25, .5,.75, 1))

    mit dem Ergebnis:

       Group.1   x.0%  x.25%  x.50%  x.75% x.100%
1  matrix -11.55  -6.55   5.45  14.45  22.45
2  lacuna -11.55  -2.55   4.45  10.45  22.45
3    cell  -8.55  -1.55  11.45  15.45  20.45

    Wenn Sie wirklich wollen, um ein boxplot-Nummern (z.B., wie weit der Backenbart gehen), haben Sie einen Blick auf ? fivenum und ? boxplot.stats.

    • Vielen Dank für die Anregungen. Da ich die selben stats der boxplot basiert auf einer Tabelle, mit dem gleichen boxplot-Funktion zu berechnen, die Ihnen versichert, dass ich das nicht machen irgendwo einen Fehler und erzeugen falsche zahlen.
    • Gil, dann gehen Sie mit boxplot.stats. Bitte beachten Sie, dass der boxplot-Whisker nicht unbedingt decken das gesamte Spektrum der Daten, d.h. die erste und die Letzte Zeile in der stats matrix kann nicht sein min und max der Daten.
    • Ganz Recht! Danke, celeites.
    InformationsquelleAutor cbeleites
  3. 1

    Andere haben beantwortet die konkrete Frage über die Rückgabe-Objekt für die boxplot Funktion, ich möchte nur hinzufügen, dass, wenn Sie möchten, erfahren Sie mehr über die Objekte zurückgeben im Allgemeinen, dann sollten Sie wirklich lernen, über die Listen und die Verwendung der str Funktion, die in der Regel geben Sie einem viel mehr sinnvolle Sicht auf ein Objekt, dann das, was du oben zeigen. Es gibt auch die TkListView Funktion in der TeachingDemos - Paket, das gibt eine mehr interaktive exploration der Liste und andere Objekte. Mit str und names - und Untergruppen (siehe auch help("[")) lassen Sie ein Gefühl dafür bekommen, was in einem Objekt zurück (die Hilfeseite für die Funktion, die das Objekt erstellt hat, ist auch ein guter Ort, um zu starten) und wie Sie auf die Stücke, die Sie wünschen.

    • Danke, Greg. Ja, in den Griff zu bekommen, die interne Struktur dieser Objekte wird in der Tat nützlich sein.
    InformationsquelleAutor Greg Snow
Schreibe einen Kommentar