From: Adaikalavan Ramasamy <ramasamy_at_cancer.org.uk>

Date: Thu 01 Jul 2004 - 02:45:17 EST

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

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
*