## ----install-------------------------------------------------------------
# Install package
install.packages("devtools", repos = "http://cran.us.r-project.org")
library("devtools")
install_github("pfmc-assessments/nwfscDeltaGLM")
# Load package
library(nwfscDeltaGLM)
## ------------------------------------------------------------------------
data(SA3)
## ------------------------------------------------------------------------
data(Example_Species)
names(Example_Species)
## ----load----------------------------------------------------------------
# Load data and strata
data(Example_Species)
masterDat = Example_Species
strata.limits = data.frame("STRATA" = c("First", "Second", "Mid", "D", "South"),
NLat = c(49, 49, 49, 36, 34.5),
SLat = c(36, 36, 34.5, 34.5, 34.5),
MinDepth = c(55, 183, 549, 55, 55),
MaxDepth = c(183, 549, 700, 549, 549), stringsAsFactors = FALSE)
## ------------------------------------------------------------------------
names(masterDat)
species = "Example_Species"
names(masterDat)[which(names(masterDat)=="EXPANDED_WT_KG")] = "Example_Species"
## ------------------------------------------------------------------------
DataList = processData(masterDat = masterDat, strata.limits = strata.limits)
## ------------------------------------------------------------------------
mcmc.control = list(chains=1, thin=1, burnin=10, iterToSave=20)
modelStructure1 = list("StrataYear.positiveTows"="fixed", "VesselYear.positiveTows"="zero", "Vessel.positiveTows"="zero", "StrataYear.zeroTows"="fixed", "VesselYear.zeroTows"="zero", "Vessel.zeroTows"="zero", "Catchability.positiveTows"="one", "Catchability.zeroTows"="zero", "year.deviations"="fixed", "strata.deviations"="fixed")
## ------------------------------------------------------------------------
fitted_models = list()
fitted_models[[1]] = fitDeltaGLM(datalist = DataList, modelStructure=modelStructure1, mcmc.control=mcmc.control, Species=species)
## ----mcmcdiag------------------------------------------------------------
# Make sure that Data is attached prior to running
data(SA3)
doMCMCDiags(datalist=DataList, mods=fitted_models, strata.limits=strata.limits)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.