G. Tremblay et al., THE USE OF POLYNOMIAL REGRESSION-ANALYSIS WITH INDICATOR VARIABLES FOR INTERPRETATION OF MERCURY IN FISH DATA, Biogeochemistry, 40(2-3), 1998, pp. 189-201
Mercury levels in fish in reservoirs and natural lakes have been monit
ored on a regular basis since 1978 at the La Grande hydroelectric comp
lex located in the James Bay region of Quebec, Canada. The main analyt
ical tools historically used were analysis of covariance (ANCOVA), lin
ear regression of the mercury-to-length relationship and Student-Newma
n-Keuls (SNK) multiple comparisons of mean mercury levels. Inadequacy
of linear regression (mercury-to-length relationships are often curvil
inear) and difficulties in comparing mean mercury levels when regressi
ons differ lead us to use polynomial regression with indicator variabl
es. For comparisons between years, polynomial regression models relate
mercury levels to length (L), length squared (L-2), binary (dummy) in
dicator variables (B-n), each representing a sampled year, and the pro
ducts of each of these explanatory variables (L x B-1, L-2 x B-1, L x
B-2, etc.). Optimal transformations of the mercury levels (for normali
ty and homogeneity) were found by the Box-Cox procedure. The models so
obtained formed a partially nested series corresponding to four situa
tions: (a) all years are well represented by a single polynomial model
; (b) the year-models are of the same shape, but the means may differ;
(c) the means are the same, but the year-models differ in shape; (d)
both the means and shapes may differ among years. Since year-specific
models came from the general one, rigorous statistical comparisons are
possible between models. Polynomial regression with indicator variabl
es allows rigorous statistical comparisons of mercury-to-length relati
onships among years, even when the shape of the relationships differ.
It is simple to obtain accurate estimates of mercury levels at standar
dized length, and multiple comparisons of these estimations are simple
to perform. The method can also be applied to spatial analysis (compa
rison of sampling stations), or to the comparison of different biologi
cal forms of the same species (dwarf and normal lake whitefish).