12.07.2015 Views

Fraser River sockeye salmon: data synthesis and cumulative impacts

Fraser River sockeye salmon: data synthesis and cumulative impacts

Fraser River sockeye salmon: data synthesis and cumulative impacts

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

as likely to be important by each of the other Cohen Commission Technical Reports (e.g., seasurface temperature). Regression analysis entails specifying a mathematical model that describesthe functional form of the relationship between the covariates <strong>and</strong> the response variable <strong>and</strong>using the observed <strong>data</strong> to estimate the parameters in the model.Many different models are possible. For example, models may include different covariates,linear <strong>and</strong> non-linear covariates, <strong>and</strong>/or interactions among different covariates. As long as thereare sufficient <strong>data</strong>, parameters for any model can be estimated but just because parameters can beestimated does not mean the model is sensible. Not surprisingly there is a vast amount ofliterature dedicated to the subject of model selection <strong>and</strong> comparison. We use the Burnham <strong>and</strong>Anderson (1998) hypothesis-driven approach to model selection <strong>and</strong> inference. In hypothesisdrivenanalyses, the only factors that would be allowed to enter the analyses would be those thatare connected to a logical, <strong>and</strong> in this case, biologically justified hypothesis. This reduces thepotential that some variables will emerge as significant simply by chance <strong>and</strong> not as a result ofany underlying mechanism, which is quite likely to happen in a project where there are largenumbers of covariates <strong>and</strong> hence potential models. St<strong>and</strong>ard practice is to select multiplefeasible c<strong>and</strong>idate models, fit each model (i.e., estimate the parameters), <strong>and</strong> then compare theperformance of each model. There are many approaches for comparing among models; we usedthe small sample size corrected Akaike’s information criterion (AIC c ) (Burnham <strong>and</strong> Anderson,1998).This project is unusual in its scope. While the response variable, ln(R/S) is available for 19stocks across B.C. <strong>and</strong> approximately 50 years of <strong>data</strong> are available for each stock, the number ofpotential covariates is very large. A total of 126 quantitative <strong>and</strong> 5 qualitative <strong>data</strong> sets wereprovided to us from the other technical reports (Table A3.4-9). We then calculated an additional32 <strong>data</strong> sets from the originals. It is possible for a single <strong>data</strong> set to be linked to (i.e.,hypothesized to impact) multiple life stages of <strong>Fraser</strong> <strong>River</strong> <strong>sockeye</strong>. In addition, there are up to4 different age types (i.e., 4sub2, 5sub2, 4sub1, <strong>and</strong> 3sub1). These links are described in TableA3.4-5 through Table A3.4-9 <strong>and</strong> result in a total of 1058 possible covariates to include in theanalysis. However, not all covariates are available for all years <strong>and</strong> stocks. Models can only becompared using AICs when the models are fit using the same <strong>data</strong>. The implication of this isthat we cannot compare all models of interest on the full <strong>data</strong> set but instead must identify timeperiods with complete <strong>data</strong> for different subsets of the covariates. For example, there is a smallsubset of the covariates (e.g., sea surface salinity for the Strait of Georgia) that have <strong>data</strong>extending back to 1950, but there are other covariates that only have <strong>data</strong> starting in 1996 (e.g.,chlorophyll a). If we wish to compare models with these two covariates (i.e. salinity <strong>and</strong>chlorophyll), we would have to either reduce the <strong>data</strong> set to those years with <strong>data</strong> for bothcovariates (i.e. limit the model to 1996-present <strong>and</strong> sacrifice the earlier <strong>data</strong> for salinity), or209

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!