R绘图学习笔记(二)－

依照计量对比药物A和药物B的响应情况

＃计量向量

dose <- c(20,30,40,45,60)

＃药物A的反应向量数据

drugA <- c(16,20,27,40,60)

＃药物B的反应向量数据

drugB <-c(15,18,25,31,40)

＃par函数用于设定或询问绘图参数。参数设定可通过par(参数名 = 取值)或par(赋值参数列表)的形式进行。

opar <-par(no.readonly=TRUE)

par(lwd=2,cex=1.5,font.lab=2)

plot(dose,drugA,type=“b”,pch=15,lty=1,col=“red”,ylim=c(0,60),main=“drug a vs drug b”,xlab=“drug dosage”,ylab=“drug response”)

lines(dose,drugB,type=“b”,pch=17,lty=2,col=“blue”)

＃基准线abline设置

abline(h=c(30),lwd=1.5,lty=2,col=“blue”)

abline(h=c(30),lwd=1.5,lty=2,col=“gray”)

＃添加小刻度

library(Hmisc)

minor.tick(nx=3,ny=3,tick.ratio=0.5)

＃添加标签信息

legend(“top left”,inset=.001,title=“drug type”,c(“A”,”B”),lty=c(1,2),pch=c(15,17),col=c(“red”,”blue”))

par(opar)