Comment puis-je activer la sortie numérique de la boîte à moustaches (avec plot=FALSE) en quelque chose d'utilisable?

Je suis avec succès à l'aide de la boxplot de fonction pour générer... boxplots. Maintenant, j'ai besoin de générer des tableaux contenant les statistiques boxplot calcule afin de créer des parcelles.

- Je le faire en utilisant le plot=FALSE option.

Le problème est que cela génère des données dans un assez bizarre format qui j'ai simplement ne peut pas faire n'importe quoi avec. Voici un exemple:

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"))

Ce que je veux, c'est un tableau pour chacune des stats -- min, max, médiane et les quartiles -- et leurs valeurs pour chaque groupe (celles de "noms").

Quelqu'un pourrait me donner un coup de main? Je suis très bien dans un R débutant.

Merci d'avance!

Gil, avez-vous besoin des numéros correspondant à la fin de la moustache ou le minimum et le maximum des données?

OriginalL'auteur Gil Williams | 2012-01-13

boxplot r

boxplot retourne une structure de R appelé list.

Une liste est plus ou moins un conteneur de données où vous pouvez vous référer à des éléments par nom.
Si vous ne A <- boxplot(...), vous pouvez accéder à la names avec A$names, le conf avec A$conf, etc.

Donc, en regardant la boxplot helpfile ?boxplot sous Value: (qui vous dit ce boxplot retours), nous voyons qu'il renvoie une liste avec les composants suivants:

   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.

De sorte que le tableau pour chacun des stats est en A$stats, chaque colonne appartient à un groupe et contient le minimum, le premier quartile, la médiane, le troisième quartile et max.

Que vous pouvez faire:

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

qui renvoie (pour 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

Alors vous pouvez vous référer à elle comme mytable['min','U'].

Stupéfiante, une réponse digne d'un bien écrit manuel de l'utilisateur en quelques minutes! Votre solution a fonctionné parfaitement. Un million de mercis!
Pour une raison quelconque, j'avais besoin d'utiliser row.names au lieu de rownames

OriginalL'auteur mathematical.coffee

4

Si vous voulez vraiment les quantiles de vos données au lieu d'une boîte à moustaches des chiffres, à l'aide de quantile directement serait mon choix (il est beaucoup plus facile à lire si vous regarder à travers ce que vous avez fait plus tard).
```
quantile (x, probs = c (0, .25, .5,.75, 1))
```
quantile lui-même ne pas travailler avec des groupes, mais vous pouvez également les combiner avec aggregate il est appelé pour chacun des groupes donnée en argument by (doit être une liste, de sorte que vous pouvez combiner plusieurs groupement des facteurs):
```
aggregate (chondro$x, by = list (chondro$clusters), 
           FUN = quantile, probs = c (0, .25, .5,.75, 1))
```
avec le résultat de:
```
   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
```
Si vous voulez vraiment avoir une boîte à moustaches des nombres (par exemple, dans quelle mesure les moustaches), d'avoir un regard sur ? fivenum et ? boxplot.stats.

Merci pour les suggestions. Depuis que j'ai besoin de la même stats de la boîte à moustaches est basée sur une table, en utilisant la même boîte à moustaches de la fonction pour calculer assure que je ne fais pas une erreur quelque part et de générer incorrect numéros.
Gil, puis aller avec boxplot.stats. Veuillez noter que la boîte à moustaches moustaches ne couvrent pas nécessairement l'ensemble des données, c'est à dire la première et la dernière ligne de la stats matrice ne peut pas être min et max de données.
Tout à fait raison! Merci, celeites.

OriginalL'auteur cbeleites
1

D'autres ont répondu à la question spécifique à propos de l'objet de retour pour la boxplot fonction, je voudrais simplement ajouter que si vous voulez vous renseigner sur le retour des objets en général, alors vous devriez vraiment apprendre plus sur les listes et comment utiliser le str fonction qui sera généralement vous donner un beaucoup plus de sens de la vue d'un objet alors que vous montrer ci-dessus. Il y a aussi le TkListView fonction dans le TeachingDemos package qui donne un plus à l'exploration interactive de la liste et d'autres objets. À l'aide de str et names et. (voir l'aide("[")) vous permettra d'avoir une idée de ce qui est dans un retour de l'objet (la page d'aide de la fonction qui a créé l'objet est également un bon endroit pour commencer), et comment accéder à la pièces que vous voulez.

Merci, Greg. Oui, obtenir une poignée sur la structure interne de ces objets sera en effet utile.

OriginalL'auteur Greg Snow

Vous devez vous connecter pour publier un commentaire.