tmpar <- par(mfrow=c(1,2),mar=c(1.5,0,0,0),pty="s")

diam <- .2
tmp <- diam / sqrt(5)

plot(0:1,0:1,type="n",axes=F,xlab="",ylab="")
mtext("(a) ANOVA",1)
circle(1-diam,.5,diam,cen=0)

circle(diam,1-diam,diam,cen=0)
text(diam,1-diam,"score")
arrow(c(diam+2*tmp,1-diam-2*tmp),c(1-diam-tmp,.5+tmp),2)

circle(diam,diam,diam,cen=0)
text(diam,diam,"residual")
arrow(c(diam+2*tmp,1-diam-2*tmp),c(diam+tmp,.5-tmp),2)

tmp1 <- 180*atan(.5)/pi
text(1-diam,.5,"yield")
text(.5,.3+.5*(1-diam),"0.51 (3)",srt=-tmp1)
text(.5,.3+.5*diam,"0.86 (16)",srt=tmp1)

##################################

diam <- .15
path3(c("yield","nonlinear","linear","residual"),"(b) Linearity",diam=diam)
tmp1 <- 180*atan(.5)/pi

text(.5,.3+.5*(1-diam),"0.30 (2)",srt=-tmp1)
text(.5,.55,"0.41")
text(.5,.3+.5*diam,"0.86 (16)",srt=tmp1)

par(tmpar)