Message-ID: <434A5FC0.7030400@sms.ed.ac.uk>
Date: 2005-10-10T12:34:08Z
From: Iain Gallagher
Subject: plot - no main title and missing abscissa value
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