Re: [R] plot - no main title and missing abscissa value

From: Iain Gallagher <i.j.gallagher_at_sms.ed.ac.uk>
Date: Mon 10 Oct 2005 - 22:34:08 EST


Hi. Sorry (esp to Uwe for the repeated messages!)

Here is the data and my code in full. Thanks for the help.

Data.

Day	Ym1Imp	Ym1sham	Semimp	Semsham
0	5.78	5.78	1.22	1.36
1	44.36	42.1	16.26	18.83
3	38.39	14.66	18.02	2.86
5	57.76	1.03	15.28	0.29
7	72.93	2.71	18.6	1.06
10	48.57	4.61	11.26	5.21
14	74.08	1.53	9.66	0.11
21	73.86	0.14	7.2	0.02


Code

ym<- read.table("ym1expression.csv", header=T, sep="\t", quote="\"") #read in ym data
attach(ym)# make data visible to R
par(mar=c(5,5,4,5),las=1, xpd=NA)
x<- c(0,1,3,5,7,10,14,21)
plot(Day, Ym1Imp, ylim=c(0,100), type="b", bty="l", main="Ym1 Expression", cex=1.3, xaxt="n", yaxt="n") #plot implant data
axis(side=1, at=c(0,1,3,5,7,10,14,21),
labels=c(0,1,3,5,7,10,14,21)) # label x axis mtext("Day", side =1, at=10, line=3, cex=1.2) # title x axis
mtext("AU", side=2, at=50, line=3, cex=1.2)# y axis title
axis(side=2, at=c(0, 25, 50, 75, 100),
labels=expression("0", "25", "50", "75", "100")) # label y axis
arrows(x, Ym1Imp-Semimp, x, Ym1Imp+Semimp, code=3, angle=90, length=0.1)# place error bars
points(Day, Ym1sham, type="b", pch=16, cex=1.3)# plot sham data
arrows(x, Ym1sham-Semsham, x, Ym1sham+Semsham, code=3, angle=90, length=0.1)# plot sham error bars legend(20, 60, legend="Implant", pch=1, lty=1, bty="n")# implant legend
legend(20, 50, legend="Sham", pch=16, lty=1, bty="n")# sham legend

Iain

-- 
Iain Gallagher
Institute for Infection & Immunology Research
Ashworth Laboratories
Kings Buildings
University of Edinburgh
Edinburgh
EH9 3JT
UK

(+44) 0131 651 3630

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Mon Oct 10 22:40:48 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:40:41 EST