Boxplots pour les groupes?

J'ai un jeu de données de test () comme indiqué ci-dessous:

Type    Met1    Met2    Met3    Met4
TypeA   65  43  97  77
TypeA   46  25  76  77
TypeA   44  23  55  46
TypeA   46  44  55  77
TypeA   33  22  55  54
TypeB   66  8   66  47
TypeB   55  76  66  65
TypeB   55  77  88  46
TypeB   36  67  55  44
TypeB   67  55  76  65

J'ai vérifié un grand nombre de liens sur la boîte parcelles, mais je n'ai toujours pas réussi pour le type de diagramme en boîte que je veux. Je tiens à avoir une boîte à moustaches avec mon X-axis type A (jaune, orange) pour tous les Mets (Met1, Met2, Met3, Met4). En substance, je veux quelque chose comme ce qui suit (tiré de ici):

Boxplots pour les groupes?

Je suis en train d'essayer des choses comme,

boxplot(formula = len ~ Type , data = test, subset == "TypeA")
boxplot(formula = len ~ Type , data = test, subset == "TypeA", add=TRUE)
Legend(legend = c( "typeA", "typeB" ), fill = c( "yellow", "orange" ) )

Mais je ne suis pas en mesure de s'en sortir avec tout ça. Quelqu'un peut-il m'aider à savoir comment dois-je faire une telle boîte de parcelles sur mes données de test dans le corrigé de la sorte?

Vous avez besoin d'apprendre à lire le code avec soin, par exemple, subset == "TypeA" est clairement pas ce qu'ils montrent dans l'exemple que vous avez lié.

OriginalL'auteur Letin | 2013-01-30

  1. 19

    Une solution avec ggplot2.

    D'abord, transformer votre bloc de données test à long format à l'aide de melt:

    library(reshape2)
test.m <- melt(test)

    Tracer les données:

    library(ggplot2)
ggplot(test.m, aes(x = variable, y = value, fill = Type)) +
  geom_boxplot() +
  scale_fill_manual(values = c("yellow", "orange"))

    Boxplots pour les groupes?

    OriginalL'auteur Sven Hohenstein

  2. 5

    Comme d'autres l'ont dit, vous devez d'abord faire fondre vos données.

    df <- read.table(text="Type    Met1    Met2    Met3    Met4
TypeA   65  43  97  77
TypeA   46  25  76  77
TypeA   44  23  55  46
TypeA   46  44  55  77
TypeA   33  22  55  54
TypeB   66  8   66  47
TypeB   55  76  66  65
TypeB   55  77  88  46
TypeB   36  67  55  44
TypeB   67  55  76  65",header=TRUE)

library(reshape2)
df2 <- melt(df)

boxplot(
  formula = value ~ variable,
  data    = df2,
  boxwex  = 0.25,
  at      = 1:4 - 0.2,
  subset  = Type == "TypeA",
  col     = "yellow",
  main    = "blah",
  xlab    = "x",
  ylab    = "y",
  ylim    = c( 0, ceiling( max( df2$value ) ) + 1 ),
  yaxs    = "i" )


boxplot(
  formula = value ~ variable,
  data    = df2,
  boxwex  = 0.25,
  at      = 1:4 + 0.2,
  subset  = Type == "TypeB",
  col     = "orange",
  add     = TRUE )

    OriginalL'auteur Roland

  3. 3

    comme ça,

    test <- structure(list(Type = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 
    2L, 2L, 2L), .Label = c("TypeA", "TypeB"), class = "factor"), 
    Met1 = c(65L, 46L, 44L, 46L, 33L, 66L, 55L, 55L, 36L, 67L
    ), Met2 = c(43L, 25L, 23L, 44L, 22L, 8L, 76L, 77L, 67L, 55L
    ), Met3 = c(97L, 76L, 55L, 55L, 55L, 66L, 66L, 88L, 55L, 
    76L), Met4 = c(77L, 77L, 46L, 77L, 54L, 47L, 65L, 46L, 44L, 
    65L)), .Names = c("Type", "Met1", "Met2", "Met3", "Met4"), 
    class = "data.frame", row.names = c(NA, -10L))


# install.packages("ggplot2", dependencies = TRUE)
require(ggplot2)
require(reshape2)
df <- melt(test)
p <- ggplot(df, aes(factor(variable), value)) + geom_boxplot(aes(fill = Type))
p

    Boxplots pour les groupes?

    Vous jetez un oeil à le geom_boxplot page de manuel.

    OriginalL'auteur Eric Fail

  4. 2

    Vous pouvez utiliser le reshape fonction pour faire une colonne avec toutes les données de Met1, .., Met4 colonnes. Il crée également une colonne (temp$time) à identifier, à partir de la colonne les données proviennent de ce que vous pouvez l'utiliser pour stratifier votre diagramme en boîte (temp$Type*temp$time).

    df <- read.table(text=
"Type Met1 Met2 Met3 Met4
TypeA   65  43  97  77
TypeA   46  25  76  77
TypeA   44  23  55  46
TypeA   46  44  55  77
TypeA   33  22  55  54
TypeB   66  8   66  47
TypeB   55  76  66  65
TypeB   55  77  88  46
TypeB   36  67  55  44
TypeB   67  55  76  65",header=TRUE)

temp <- reshape(df, direction='long', varying = 2:5, sep='')

boxplot(temp$Met ~ temp$Type*temp$time, col=c("yellow", "orange"))

    OriginalL'auteur hector