#' @export
CatchEffortPlot = function(CatchEffortData,yrs,graphic='R',col=rgb(0,1,0,0.8),ht=8,wd=6,path=file.path( project.datadirectory("bio.surfclam"), "figures")){
# CatchEffortPlots
if(graphic=='pdf')pdf(file.path(path,"CatchEffortPlot.pdf"),height=ht,width=wd)
if(graphic=="R")x11(ht,wd)
if(missing(yrs))yrs = CatchEffortData$Year
d = subset(CatchEffortData,Year%in%yrs)
par(mfrow = c(2,1), mar = c(0, 4, 0, 1), omi = c(0.5, 0.3, 0.3, 0.3),las=1)
#browser()
plot(yrs, d$Grand.Catch/1000, type = 'b', lwd = 2, ylim = c(0, max(d$Ban.Catch/1000,na.rm=T)), col='grey50', pch=17,cex.axis=1.2, ylab = "", xlab = "", xaxt='n')
lines(yrs, d$Ban.Catch/1000, type = 'b', lwd = 2,pch=16)
axis(1, lab = F)
mtext("Catch (kt)",2,3,las=0,cex=1.2)
legend('topleft',c("Banqureau","Grand Bank"),pch=16:17,col=c('black','grey50'),cex=1,bty='n')
plot(yrs, d$Grand.Effort, type = 'b', lwd = 2, ylim = c(0, max(d$Ban.Effort,na.rm=T)), col='grey50', pch=17, ylab = "",cex.axis=1.2, xlab = "", xaxt='n')
lines(yrs, d$Ban.Effort, type = 'b', lwd = 2,pch=16)
axis(1,cex.axis=1.2)
mtext(expression(Effort ( km^2 )),2,3,las=0,cex=1.2)
if(graphic!="R")dev.off()
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.