Mike Livermore – CSI 773 – Week 2 Assignment – Lecture Examples – September 2010

1:

2:

3.1:

3.2:

windows()

cex = .9

xAxisTick = -.03

xAxisTickSep = .3 

yAxisTick = -.02      # These are for the y-axis

yAxisTickSep = .5  

# panel layout

panels = panelLayout(nrow=2,ncol=1,rowSep=c(0,.2,0),leftMar=.42)

# plot two curves and difference

rx = c(-.04,1.04)

ry = c(-.04,1.04)

xgrid = seq(0,1,by=.2)

ygrid = c(0,.25,.5,.75,1)

panelSelect(panels,1,1)

panelScale(rx,ry)

panelFill()

panelGrid(x=xgrid,y=ygrid,col="white",lwd=2)

axis(side=2,at=ygrid,tck=yAxisTick,mgp=c(2,yAxisTickSep,0),las=2,cex=cex)

     # mtext() could be used instead

     # mtext(side=2,line=.3,text=format(ygrid),at=ygrid,las=2,cex=cex)

lines(x,y1,lwd=3)

lines(x,y2,lwd=3)

lines(x,curveDif,lwd=3,col="blue")

panelOutline()

mtext(side=3,line=2.4,'Top Panel: Two Curves and Their Difference',cex=1.3)

mtext(side=3,line=.7,'Bottom Panel: Zoomed View of the Difference',cex=1.3)

# Rescale and plot difference

ry = c(.02,.14)

ygrid=c(.02,.05,.08,.11,.14)

panelSelect(panels,2,1)

panelScale(rx,ry)

panelFill()

panelGrid(x=xgrid,y=ygrid,col="white",lwd=2)

axis(side=1,at=xgrid,mgp=c(2,xAxisTickSep,0),tck=xAxisTick,las=1,cex=cex)

axis(side=2,at=ygrid,mgp=c(2,yAxisTickSep,0),tck=yAxisTick,las=2,cex=cex)

lines(x,curveDif,lwd=3,col="blue")

panelOutline()

4: