créer un boxplot dans R qui étiquette une boîte avec la taille de l'échantillon (N)

Est-il un moyen de créer une boîte à moustaches dans R qui affiche la boîte de (quelque part) "N=(taille de l'échantillon)"? Le varwidth logique ajuste la largeur de la zone sur la base de la taille de l'échantillon, mais qui ne permet pas de comparaisons entre les différentes parcelles.

FWIW, je suis l'aide de la boîte à moustaches de commande de la manière suivante, où f1 est un facteur:

boxplot(xvar ~ f1, data=frame, xlab="input values", horizontal=TRUE)

source d'informationauteur J Miller

  1. 24

    Voici quelques ggplot2 code. Il va afficher la taille de l'échantillon à la moyenne de l'échantillon, faisant de l'étiquette multifonctionnel!

    Tout d'abord, une fonction simple pour fun.data

    give.n <- function(x){
   return(c(y = mean(x), label = length(x)))
}

    Maintenant, pour démontrer, avec les diamants de données

    ggplot(diamonds, aes(cut, price)) + 
   geom_boxplot() + 
   stat_summary(fun.data = give.n, geom = "text")

    Vous pourriez avoir à jouer avec la taille du texte pour le rendre bien, mais maintenant vous avez une étiquette de la taille de l'échantillon qui donne également un sens de l'inclinaison.

  2. 9

    Vous pouvez utiliser le names paramètre pour écrire les n en regard de chaque facteur nom.

    Si vous ne voulez pas de calculer la n vous-même, vous pouvez utiliser cette petite astuce:

    # Do the boxplot but do not show it
b <- boxplot(xvar ~ f1, data=frame, plot=0)
# Now b$n holds the counts for each factor, we're going to write them in names
boxplot(xvar ~ f1, data=frame, xlab="input values", names=paste(b$names, "(n=", b$n, ")"))
  3. 4

    Pour obtenir le n sur le dessus de la barre, vous pouvez utiliser text avec le stat les détails fournis par la boîte à moustaches comme suit

    b <- boxplot(xvar ~ f1, data=frame, plot=0)
text(1:length(b$n), b$stats[5,]+1, paste("n=", b$n))

    Les stats de champ de b est
    une matrice, chaque colonne contient l'extrême de la partie inférieure de la moustache, la partie inférieure de la charnière, la médiane, la charnière supérieure et l'extrémité de la moustache supérieure pour un groupe/de la parcelle.

  4. 0

    La gplots paquet fournit boxplot.nqui, selon la documentation produit une boîte à moustaches annoté avec le nombre d'observations.

  5. 0

    J'ai trouvé une solution de contournement à l'aide de la Envstats paquet. Ce paquet doit être téléchargé, chargé et activé à l'aide de: 

    library(Envstats)

    La stripChart (différent de stripchart) ne s'ajouter à la graphique de certaines valeurs telles que la n valeurs. D'abord j'ai tracé ma boîte à moustaches. Ensuite, j'ai utilisé l'add=T dans le stripChart. De toute évidence, beaucoup de choses ont été cachés dans les stripChart code de sorte qu'ils n'apparaissent pas sur la boîte à moustaches. Voici le code que j'ai utilisé pour la stripChart de cacher la plupart des articles.

    Une boîte à moustaches intégré stripChart pour afficher les n valeurs:

    stripChart(data.frame(T0_G1,T24h_G1,T96h_G1,T7d_G1,T11d_G1,T15d_G1,T30d_G1), show.ci=F,axes=F,points.cex=0,n.text.line=1.6,n.text.cex=0.7,add=T,location.scale.text="none")

    Donc une boîte à moustaches

    boxplot(data.frame(T0_G1,T24h_G1,T96h_G1,T7d_G1,T11d_G1,T15d_G1,T30d_G1),main="All Rheometry Tests on Egg Plasma at All Time Points at 0.1Hz,0.1% and 37 Set 1,2,3", names=c("0h","24h","96h","7d ", "11d", "15d", "30d"),boxwex=0.6,par(mar=c(8,4,4,2)))

    Puis stripChart

    stripChart(data.frame(T0_G1,T24h_G1,T96h_G1,T7d_G1,T11d_G1,T15d_G1,T30d_G1), show.ci=F,axes=F,points.cex=0,n.text.line=1.6,n.text.cex=0.7,add=T,location.scale.text="none")

    Vous pouvez toujours ajuster la hauteur des numéros (n valeurs) pour qu'elles correspondent à l'endroit où vous voulez.