[R] help with histogram

Martin Maechler maechler at stat.math.ethz.ch
Fri May 7 17:20:28 CEST 2004


>>>>> "Guazzetti" == Guazzetti Stefano <Stefano.Guazzetti at ausl.re.it>
>>>>>     on Fri, 7 May 2004 14:45:22 +0200 writes:

    Guazzetti> Maybe you want something like:

    Guazzetti> x<-rnorm(1000)
    Guazzetti> hist(x, breaks=100,
    Guazzetti>    col=ifelse(abs((hist(x, breaks=100, main=""))$breaks) < 1.669,
    Guazzetti>    4,2))

    Guazzetti> see also the density argument in ?hist

very good!

Even slightly better {not calling  hist() twice} is

 x <- rnorm(1000)
 hx <- hist(x, breaks=100,plot=FALSE)
 plot(hx, col=ifelse(abs(hx$breaks) < 1.669, 4, 2))

Regards, Martin




More information about the R-help mailing list