The object fr.data.lpue stores generated Observed Landings Per Unit Effort which could be derived from commercial bottom trawl fleets. As real data are confidential, data presented in the example are generated data.
The object fr.data.lpue is a data.frame made of 6 columns named year (calendar year), fishing.season (the fishing season during which the cuttlefish is exploited; a fishing season starts in July of a calendar year Y and ends in June of a calendar year Y+1), month (calendar month from January (1) to December (12)), rectangle (ICES rectangle code), power.class (power class of the trawler concerned).
Each of the 9110 rows of the data frame refer to the fishing effort carried out during one cruise in one ICES rectangle.
The development of the two-stage biomass model to assess the English Channel cuttlefish stock was carried out in the framework of the EU funded project CRESH (under the Interreg IV A France-Manche-England programme). The development of the R package to perform the routine assessment of the cuttlefish stock was co-funded by France Filiere Peche and by the Departement des Peches Maritimes et de l'Aquaculture.
Michael Gras and Jean-Paul Robin
Gras, M., Roel, B. A., Coppin, F., Foucher, E. and Robin, J.-P. (2014). A two-stage biomass model to assess the English Channel cuttlefish (Sepia officinalis L.) stock. Submitted to ICES Journal of Marine Science.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
data(fr.data.lpue) fr.delta.glm<-delta.glm(input.data=fr.data.lpue) par(mfrow = c(2,2)) #Histogram of the binomial error GLM residuals hist(fr.delta.glm$binomial.residuals) #Plot with the fitted data on the x axis and and the re plot(fr.delta.glm$binomial.fit, fr.delta.glm$binomial.residuals) #QQplot of the residuals from the binomial error GLM qqnorm(fr.delta.glm$binomial.residuals) qqline(fr.delta.glm$binomial.residuals) par(mfrow = c(2,2)) #Histogram of the residuals from the Gaussian error GLM hist(fr.delta.glm$gaussian.residuals) #Plot of fitted values vs residuals from the Gaussian error GLM plot(fr.delta.glm$gaussian.fit,fr.delta.glm$gaussian.residuals) qqnorm(fr.delta.glm$gaussian.residuals) qqline(fr.delta.glm$gaussian.residuals) #Aggregation of the standardised LPUE per year. Aggregation #can be done on the 3 other factors in the same way. fr.yearly.lpue<-aggregate(fr.delta.glm$predicted.lpue$st.lpue, list(fr.delta.glm$predicted.lpue$fishing.season), FUN="mean") fr.yearly.lpue<-data.frame(c(1900:1905), fr.yearly.lpue) colnames(fr.yearly.lpue)<-c("year","fishing.season","fr.st.lpue")