feather.plot<-function(r,theta,xpos,yref=0,use.arrows=TRUE,
col.refline="lightgray",fp.type="s",main="",xlab="",ylab="",
xlabels=NULL,...) {
if(missing(xpos)) xpos<-1:length(theta)
# reverse the angles and start at 12 o'clock
if(fp.type == "m") theta<-5*pi/2-theta
x<-r*cos(theta)
y<-r*sin(theta)
xmult<-diff(range(xpos))/(diff(range(y))*2)
x<-x*xmult
xlim<-range(c(xpos,x+xpos))
ylim<-range(c(y,yref))
oldpin<-par("pin")
xdiff<-xlim[2]-xlim[1]
ydiff<-ylim[2]-ylim[1]
plot(0,xlim=xlim,ylim=ylim,type="n",main=main,xlab=xlab,ylab=ylab,axes=FALSE)
box()
if(is.null(xlabels)) axis(1)
else axis(1,at=xpos,labels=xlabels)
abline(h=yref,col=col.refline)
if(use.arrows) arrows(xpos,yref,xpos+x,y,length=0.1,...)
else segments(xpos,yref,xpos+x,y,...)
par(pin=oldpin)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.