plot(sID, .THEME = NULL)

[johnrsibert]

This method seems to produce 3 plots, but the first two are obliterated by the third one. Does it need a call to x11() between plots?

[johnrsibert]

These changes produce nicely offset graphics windows with titles:

$ git diff catchMSY_graphics.R
diff --git a/R/catchMSY_graphics.R b/R/catchMSY_graphics.R
index 8bc9f1b..d3ac57d 100644
--- a/R/catchMSY_graphics.R
+++ b/R/catchMSY_graphics.R
@@ -4,6 +4,10 @@ plot.stockID <- function(sID,.THEME=NULL){

        with(sID,{
                if(exists("idx")){
+                       xinc = 50
+                       yinc = 50
+                       xpos = 2*xinc
+                       ypos = yinc
                        # ------------------------------------------------------
                        # Vulnerable biomass.
                        bio.df <- data.frame(Year=data$year,
@@ -18,11 +22,14 @@ plot.stockID <- function(sID,.THEME=NULL){
                        p <- p + stat_summary(geom="ribbon", fun.data="mean_cl_b
                                              alpha=0.2,fill="red")
                        p <- p + labs(x="Year",y="Vulnerable Biomass (lbs)")
+                       x11(title="Biomass",xpos=xpos,ypos=ypos)
                        print(p+.THEME)
                        # ------------------------------------------------------

                        # ------------------------------------------------------
                        #  Spawning biomass depletion
+                       xpos = xpos + xinc
+                       ypos = ypos + yinc
                        ssb.df <- data.frame(Year=data$year,
                                             Biomass=t(ps.dt)[,idx])
                        ssb.df <- gather(ssb.df,key,dt,-Year)
@@ -35,16 +42,20 @@ plot.stockID <- function(sID,.THEME=NULL){
                        p <- p + stat_summary(geom="ribbon", fun.data="mean_cl_b
                                              alpha=0.2,fill="red")
                        p <- p + labs(x="Year",y="Spawning Biomass Depletion")
+                       x11(title="Depletion",xpos=xpos,ypos=ypos)
                        print(p+.THEME)
                        # ------------------------------------------------------

                        # ------------------------------------------------------
                        # Parameter distributions
+                       xpos = xpos + xinc
+                       ypos = ypos + yinc
                        gS <- gather(as.data.frame(S[idx,]))
                        p <- ggplot(gS,aes(x=value)) 
                        p <- p + geom_histogram(aes(y=..density..),alpha=0.3) 
                        p <- p + geom_density(col=NA,fill="red",alpha=0.2) 
                        p <- p + facet_wrap(~key,scales="free")
+                       x11(title="Parameters",xpos=xpos,ypos=ypos)
                        print(p+.THEME)
                        # ------------------------------------------------------

@@ -54,4 +65,4 @@ plot.stockID <- function(sID,.THEME=NULL){

        })

-}
\ No newline at end of file
+}