#'S3 class methods for providing summary plots
#'@export
plot.stockID <- function(sID,.THEME=NULL){
with(sID,{
if(exists("idx")){
# -----------------------------------------------------------------#
# Vulnerable biomass.
bio.df <- data.frame(Year=data$year,
Biomass=t(ps.bt)[,idx])
bio.df <- gather(bio.df,key,bt,-Year)
p <- ggplot(bio.df,aes(Year,y=bt))
p <- p + stat_summary(geom="ribbon",fun.ymin="min",fun.ymax="max",
alpha=0.1)
p <- p + stat_summary(geom="ribbon", fun.data="median_hilow",
alpha=0.2,fill="blue")
p <- p + stat_summary(geom="ribbon", fun.data="mean_cl_boot",
alpha=0.2,fill="red")
p <- p + labs(x="Year",y="Vulnerable Biomass (lbs)")
print(p+.THEME)
# -----------------------------------------------------------------#
# -----------------------------------------------------------------#
# Spawning biomass depletion
ssb.df <- data.frame(Year=data$year,
Biomass=t(ps.dt)[,idx])
ssb.df <- gather(ssb.df,key,dt,-Year)
p <- ggplot(ssb.df,aes(Year,y=dt))
p <- p + stat_summary(geom="ribbon",fun.ymin="min",fun.ymax="max",
alpha=0.1)
p <- p + stat_summary(geom="ribbon", fun.data="median_hilow",
alpha=0.2,fill="blue")
p <- p + stat_summary(geom="ribbon", fun.data="mean_cl_boot",
alpha=0.2,fill="red")
p <- p + labs(x="Year",y="Spawning Biomass Depletion")
print(p+.THEME)
# -----------------------------------------------------------------#
# -----------------------------------------------------------------#
# Parameter distributions
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")
print(p+.THEME)
# -----------------------------------------------------------------#
}
})
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.