3 图形初阶

dose<-c(20,30,40,45,60)
drugA<-c(16,20,27,40,60)
drugB<-c(15,18,25,31,40)
plot(dose,drugA,type="b")
opar<-par(no.readonly=TRUE) # 复制一份图形参数设置
par(pin=c(2,3))#修改图形参数
par(lwd=2,cex=1.5)
par(cex.axis=.75,font.axis=3)
plot(dose,drugA,type="b",pch=19,lty=2,col="red")#直接在绘图函数中制定对应的参数
plot(dose,drugB,type="b",pch=23,lty=6,col="blue",bg="green")
plot(opar)#应用新图形参数
#pch 绘制点使用的符号，cex 指定符号的大小  #lty 指定线条类型，lwd指定线条宽度
例子：3-3 # 成功绘制目标图形


dose<-c(20,30,40,45,60)
drugA<-c(16,20,27,40,60)
drugB<-c(15,18,25,31,40)


opar<-par(no.readonly=TRUE) 
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(h=c(30),lwd=1.5,lty=2,col="gray")
library(Hmisc)#加载函数
minor.tick(nx=3,ny=3,tick.ratio=0.5)# 添加次要刻度线
legend("topleft",inset=0.05,title="Drug Type",c("A","B"),lty=c(1,2),pch=c(15,17),col=c("red","blue")) # 添加图例
par(opar)


3-2 #自定义坐标轴
x<-c(1:10) #生成数据
y<-x
z<-10/x
opar<-par(no.readonly=TRUE)
par(mar=c(5,4,4,8)+0.1)#增加边界大小
plot(x,y,type="b",   # 绘制x对y的图形
        pch=21,col="red",yaxt="n",lty=3,ann=FALSE)
lines(x,z,type="b",pch=22,col="blue",lty=2)# 绘制x对10/x的图形
axis(2,at=x,labels=x,col.axis="red",las=2)
axis(4,at=z,labels=round(z,digits=2),col.axis="blue",las=2,cex.axis=0.7,tck=-0.01)
mtext("y=1/x",side=4,line=3,cex.lab=1,las=2,col="blue")
title("An Example of Creative Axes",xlab="X values",ylab="Y=X")
par(opar)