依照计量对比药物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)
