Modifier les couleurs de la barre d'histogramme supérieure à une certaine valeur

Jusqu'à présent, j'ai réussi à changer la couleur d'une seule barre de l'histogramme suivant l'exemple ici

test <- rnorm(100);
h <- hist(test);
b <- cut(1, h$breaks);
clr <- rep("grey", length(h$counts));
clr[b] <- "red";
plot(h, col=clr);

Je veux être en mesure de modifier la couleur de l'histogramme des bacs qui sont au-dessus d'un certain axe des x de la valeur - par exemple, qui sont au-dessus de 1 dans la fonction de distribution dans l'exemple. Partie de la raison pour laquelle j'ai de la difficulté, c'est que je n'ai pas de comprendre exactement le facteur qui cut() retourne.

source d'informationauteur MattLBeck

  1. 8

    Fondamentalement, vous voulez une logique sélecteur sur test pas sur le cuts.

    Voici ce que votre couper l'objet ressemble:

    > bks <- cut(test,10)

    Les niveaux sont de type caractère:

    niveaux(bks)
    Un "(-2.53,-2.01]" "(-2.01,-1.5]" "(-1.5,-0.978]" "(-0.978,-0.459]"
    [5] "(-0.459,0.0596]" "(0.0596,0.578]" "(0.578,1.1]" "(1.1,1.62]"

    [9] "(1.62,2.13]" "(2.13,2.65]"

    Les données sont de type numérique:

    > head(as.numeric(bks))
[1] 5 6 6 6 3 5

    Voici une solution à l'aide de ggplot2 plutôt que de faire les coupes et ainsi de suite par la main:

    test <- rnorm(100)
dat <- data.frame( x=test, above=test>1 )
library(ggplot2)
qplot(x,data=dat,geom="histogram",fill=above)

    Modifier les couleurs de la barre d'histogramme supérieure à une certaine valeur

  2. 7

    Changer votre couleur de vecteur, clrde sorte qu'il affiche en rouge à chaque fois que la barre est plus grand que 1 et gris sinon.

    clr <- ifelse(h$breaks < 1, "grey", "red")[-length(h$breaks)]

    Puis parcelle avant.