Re: [R] Question about plotting related to roll-up

From: Adaikalavan Ramasamy <ramasamy_at_cancer.org.uk>
Date: Thu 01 Jul 2004 - 02:45:17 EST

grp <- rep(1:5, each=3)
resp <- rnorm(15)

mu <- tapply(resp, grp, mean)
s <- tapply(resp, grp, sd)
stopifnot( identical( names(mu), names(s) ) )

LCL <- mu - 2*s # lower confidence limit UCL <- mu + 2*s
Here I choose 2 as we expect 95% of the data to fall under 4 sd.

# Type 1
plot(names(mu), mu, type="l", ylim=c( min(LCL), max(UCL) )) lines(names(mu), UCL, lty=3)
lines(names(mu), LCL, lty=3)

Your group must contain only numeric values. Otherwise, you will need to use a numerical coding followed by mtext() with proper characters.

# Type 2
plot(names(mu), mu, type="p", ylim=c( min(LCL), max(UCL) )) arrows( as.numeric(names(mu)), LCL, as.numeric(names(mu)), UCL, code=3, angle=90, length=0.1 )

On Wed, 2004-06-30 at 17:33, Coburn Watson wrote:
> Hello R'ers,
>
> I have a large set of data which has many y samples for each unit x. The data
> might look like:
>
> Seconds Response_time
> ---------- ----------------
> 0 0.150
> 0 0.202
> 0 0.065
> 1 0.110
> 1 0.280
> 2 0.230
> 2 0.156
> 3 0.070
> 3 0.185
> 3 0.255
> 3 0.311
> 3 0.120
> 4
> .... and so on
>
> When I do a basic plot with type=l or the default of points it obviously plots
> every point. What I would like to do is generate a line plot where the
> samples for each second are rolled up, averged and plotted with a bar which
> represents either std dev or some other aspect of variance. Can someone
> recommend a plotting mechanism to achieve this? I have adding lines using
> some of the smoothing functions but seem unable to remove the original plot
> line which is drawn (is there a way to just plot the data as feed through the
> smoothing function without the original data?).
>
> Please remove "_nospam" from the email address to reply directly to me.
>
> Thanks,
>
> Coburn Watson
> Software Performance Engineering
> DST Innovis
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>



R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Jul 01 02:49:58 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:54:38 EST