- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
The Scientific World Journal
Volume 2013 (2013), Article ID 362854, 7 pages
To Identify the Important Soil Properties Affecting Dinoseb Adsorption with Statistical Analysis
1State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China
2Institute of Meteorology and Climate Research, Karlsruhe Institute of Technology (KIT), 82467 Garmisch-Partenkirchen, Germany
3National Center for Computational Hydroscience and Engineering, University of Mississippi, Oxford, MS 38655, USA
4Nanjing Hydraulic Research Institute, Nanjing 210029, China
Received 8 March 2013; Accepted 3 April 2013
Academic Editors: H. Filik and A. Hursthouse
Copyright © 2013 Yiqing Guan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Investigating the influences of soil characteristic factors on dinoseb adsorption parameter with different statistical methods would be valuable to explicitly figure out the extent of these influences. The correlation coefficients and the direct, indirect effects of soil characteristic factors on dinoseb adsorption parameter were analyzed through bivariate correlation analysis, and path analysis. With stepwise regression analysis the factors which had little influence on the adsorption parameter were excluded. Results indicate that pH and CEC had moderate relationship and lower direct effect on dinoseb adsorption parameter due to the multicollinearity with other soil factors, and organic carbon and clay contents were found to be the most significant soil factors which affect the dinoseb adsorption process. A regression is thereby set up to explore the relationship between the dinoseb adsorption parameter and the two soil factors: the soil organic carbon and clay contents. A 92% of the variation of dinoseb sorption coefficient could be attributed to the variation of the soil organic carbon and clay contents.
Dinoseb (2-sec-butyl-4,6-dinitrophenol) is a member of the dinitrophenol family of pesticides, commonly used for controlling the growth of annual grassy and broadleaf weeds. It has long persistence, which leads to an accumulation in soil. It has been found in many areas of the world [1–3]. Many countries have prohibited the usage of the dinoseb. In USA, EPA banned dinoseb usage in 1986. Much research focuses on dinoseb’s toxic effects on human beings, animals, and microorganisms [4–6]. And the measuring technique has also well been studied [7–10].
After being applied to soil, the transport and fate of herbicides are controlled by many complicated mechanisms, including sorption to soil, uptake by plants, transport vial runoff and leaching, biodegradation, photodegradation, volatilization, and chemical degradation [11, 12]. Sorption is one of the most important mechanisms that influence the presence of herbicides in soil . To evaluate the sorption property of herbicides, the popular methods include batch equilibrium technique, column experiments, and field experiments. Comparing with other two techniques, batch experiments are easy and fast to perform, and the cost is low . The results of the batch sorption equilibrium experiments are usually fitted with linear sorption model or Freundlich sorption model to derive the sorption parameters.
Since the sorption of herbicides to a large extent was determined by soil properties, such as organic carbon (OC) content [14, 15], clay content [16, 17], and pH value , the multiple regression is usually applied to explore the relationships between the sorption parameters and soil environmental factors. However, in multilinear regression, if the predictor variables are not independent, multicollinearity will be a common statistical phenomenon. Multicollinearity may not affect the goodness of the multiple regression prediction, but it reflects the determination of the importance of each environmental factor. Stepwise regression can be used to correct for multicollinearity. It has been frequently applied in educational and psychological research, both to select useful subsets of variables and to evaluate the order of importance of variables . Path analysis was developed around 1918 by Sewall Wright. It was usually used to decompose correlations into different pieces for interpretation of effects. The two methods have been applied to many studies, including biology, sociology, and econometrics [20, 21]. The objectives of this study, therefore, are (1) to use the bivariable correlation analysis and path analysis to investigate the extent of these influences and explicitly explain direct and indirect effects of soil characteristic factors on dinoseb adsorption, (2) to use stepwise regression analysis for excluding the factors which have little influence on the adsorption parameter, and (3), based on these results, to set up regression equation of adsorption values with the most important soil factors.
2. Materials and Methods
2.1. Soil Sample Collection
The soil samples were collected in the upper Rhone river valley in Southwest Switzerland. Dinoseb had been found in the groundwater of the Rhone plain. This alluvial plain is well cultivated and of great economic and ecological importance, but it was alleged a high vulnerability of the groundwater to contamination . Along a transect, pits A, B, C, and D were excavated; the distances between the pits were 6.5 m, 8.1 m, and 6.5 m, respectively. At these four sites, altogether 55 small disturbed soil samples were collected at 15, 30, 55, 70, and 85 cm depths.
2.2. Characteristics of Dinoseb
The purity of the dinoseb production used in the study was 93% (Dr. Ehrenstorfer, Germany). The properties of dinoseb are summarized in Table 1.
2.3. Experimental Design
First of all, basic soil properties such as bulk density, porosity, particle size distribution, pH, cation exchange capacity (CEC), and organic carbon content were determined.
Soil samples collected from the field were air-dried at room temperature and sieved at 2 mm. Then 3 g of dry soil was mixed with 6 mL of dinoseb solutions at different concentrations into a 9 mL polypropylene centrifuge tube. The concentrations are 0, 1.5, 4.5, 9, and 15 mg/L, respectively. The tubes were shaken for 24 hrs on a rotary tumbler at 20°C. This duration was sufficient to achieve sorption equilibrium, but not long enough for chemical or biological transformations to significantly affect the results, as attested by sorption kinetics tests. The aqueous phase was separated from the solid phase by centrifugation for 15 min at 7000 rpm. The supernatant was filtered through a disposable 0.45 μm cellulous filter. The filtrates were analyzed by injection into a high performance liquid chromatography with diode array detection (Hewlett Packard series 1050) using a C18 column of 25 cm length (VYDAC). The light absorption wavelength used for the detection of dinoseb is 265 nm. The flow was set as 1 mL/min; the solvents used for HPLC were distilled water, acetonitrile (purity ≥ 99.8%), and 0.05% trifluoroacetic acid. The original composition of flow was 40% acetonitril, 40% water, and 20% trifluoroacetic acid. During each measuring, trifluoroacetic acid was kept as constant as 20%; acetonitrile increased from 40% to 65% over the first 15 min, then increased to 80% over the next 5 min after this measuring started; acetonitrile then decreased back to 40% over 5 min and sustained at this level until the end of this measuring, which was 30 min after measuring starts.
The adsorbed dinoseb mass was calculated with the difference between the initial dinoseb concentration and that measured in the supernatant. All batch sorption experiments were conducted in triplicate.
2.4. Statistical Analysis
Correlation analysis and path analysis in this study were used to demonstrate the degrees of the variables’ interactions or interferences with each other and the exact variable with the most exerting influence. Stepwise multiple-linear regression analysis was used for identifying the linear relationship between dinoseb absorption coefficients with soil properties. Significance of differences was either tested by using a parametric -test or -statistics in ANOVA (analysis of variance).
Stepwise multiple-linear regression  is one method in multiple linear regressions that used to analyze the linear relationship between single dependent variable with several independent variables. It was selected for this research because (1) multiple-linear regression makes use of the most of the directly observed and experimental information that has been available ; (2) the number of controlled variables (OC, CEC, pH, Clay) is fairly small so that it could be easily performed to analyze including all of them; (3) the bivariate correlations among soil properties with the dinoseb adsorption values are not explicitly fixed especially with the influence of multicollinearity; (4) the problem of overfitting could be avoided by adding or deleting variable with the specific criteria. Therefore, backward elimination  is applied to build up the final regression equation describing a predicted variable as a function of several independent variables. It follows these procedures: firstly adding all the independent variables into regression, secondly analyzing significance of difference about the partial coefficient of each independent variable and deleting the one with lowest significant contribution to the regression equation compared with the removing criteria (alpha-to-remove value), and finally repeating the regression modeling and testing with remaining variables and removing until all the remaining variables have significant contribution to the regression equation. But some issues of stepwise regression still exist such as that it cannot explicitly interpret the multicollinearity between controlled variables .
Due to the problem of multicollinearity in regression [28, 29], before setting up a stepwise multilinear regression, bivariate correlation analysis and path analysis  based on the causal relationship were adopted to make explicit the rational of conventional regression calculations. Path analysis have special usefulness in decomposing the soil property effects on the dinoseb adsorption into direct and indirect effects and quantifying the collinearity in the regression model. Note that the direct and indirect effects importantly depend on how the model is built . In this study, only the regression model which includes all variables was applied to path analysis to capture the overall direct and indirect effect from four soil properties on the dinoseb adsorption values.
3. Results and Discussion
3.1. Sorption Isotherms
The physical and chemical properties of the soil samples collected at four sits over 5 depth were summarized in Table 2.
Fifty-five sorption isotherms of dinoseb are determined. The isotherms have been fitted by Freundlich model: where is the adsorbed chemical concentration (g/g), is the Freundlich partition coefficient, (cm3/g), is an empirical coefficient and is the equilibrium concentration (mg/L).
Thirty out of the 55 dinoseb sorption isotherm fittings have a , more than 0.95. The values of 19 fittings are from 0.90 to 0.95. The left 6 fittings have a value from 0.87 to 0.90. Therefore, the model can well describe most of the dinoseb sorption at concentrations less than 15 mg/L. The derived Freundlich distribution coefficients are listed in Table 3.
3.2. Bivariate Correlation Analysis Results
Sorption behavior of dinoseb is believed to significantly depend on organic carbon content, and Clay content and pH have also been reported to affect the sorption of dinoseb . Correlation between the four soil properties and dinoseb adsorption capacity coefficients were assessed. Pearsons correlation coefficients stranded for the bivariate correlation among the dinoseb values and four soil properties (Table 4). On the basis of these data, the two-tailed parametric -test was performed to investigate the significance of differences for the relations between each two variable pair . Table 4 shows that correlations were all significant except that between CEC and the dinoseb values. The dinoseb values had highest positive correlation value with OC , followed by Clay and furthermore highest negative correlation value with pH , while there was no significant relationship between and CEC. These results indicated that the related soil properties with the dinoseb values were soil OC, Clay content, and pH. Moreover, it is noteworthy that correlation matrices among soil properties show several sets of relationships. The amount of OC was significantly and positively correlated with Clay content and CEC and negatively with pH at the significance level of 0.01. Similar to OC, the relationship between CEC, and Clay content was also significantly high. However, generally pH values were weakly and not significantly correlated with CEC and Clay. The results showed that not only the is positively correlated to soil OC, CEC and Clay content, and negatively with pH, but also relationships between two soil properties are still fairly high. With the limitation of bivariate correlation, the Pearson correlation coefficients cannot demonstrate the real relationships when multicollinearity exists.
3.3. Path Analysis Results
With path analysis, we can decompose the correlations into direct and indirect effects. The effects are quantified with the path coefficients (Table 5). According to the path coefficients, the sequence of direct effects to is OC > Clay > pH > CEC. Both zero-order correlation and path analysis show OC content has a significant positive effect on , and the direct effect on is much higher than the other three factors (path coefficient 1.056). In the zero-order correlation matrix, pH is significantly correlated with (correlation coefficient −0.659). The path analysis shows that this correlation is mainly due to the correlation of pH with OC (path coefficient −0.662). The direct effect of pH on is low (path coefficient −0.066). For CEC, with almost zero direct effect on , it can be considered that the moderate correlation (correlation coefficient 0.436) with is mainly due to the contribution of collinearity between OC content and CEC. Clay content has negative direct effect on (path coefficient −0.216), although the indirect effect due to correlation with OC is more obvious (path coefficient 0.746). Contrast to that, the correlation coefficient shows that Clay has a positive relationship with . Dinoseb is a weak acid with a pH of 4.4–4.62  and is mainly in anionic form at the pH of the studied soils . Therefore, it is more reasonable that its affinity to soil was negatively correlated with the content of the negatively charged clays.
3.4. Stepwise Multiple-Linear Regression Results
Based on the correlation matrix in Table 4 and path analysis coefficients in Table 5, it is obvious that it is not independent between pairs of the soil properties and that makes the interpretation of multiple linear regression equations between the dinoseb values and soil properties unreliable. The problem of multicollinearity among soil properties in linear model has been generally recognized in many studies [35, 36]. In order to overcome multicollinearity, stepwise regression, one of several standard procedures  for variable selection, was applied for multiple linear regression in this study. Due to the small number of correlated variables (OC, pH, CEC, Clay), the backward elimination was performed starting with all four soil properties as controlled variables and successively eliminates one at a time. And the criteria based on -statistics is to remove the lowest -to-remove statistic which is bigger than 0.05.
The regression coefficients and statistics summary of each prediction model of dinoseb values depending on soil properties as developed using stepwise multiple linear regression are presented in Tables 6 and 7.
In Table 6, the standardized coefficients (beta values) indicate the strength of the effect of the respective soil properties on dinoseb values; that is, the larger absolute value shows the stronger effect. Zero-order correlations have been discussed in correlation analysis. Partial correlations reveal the relationship between residualized dinoseb values and residualized soil properties, and part correlations express the correlations between residualized dinoseb values and unaltered soil properties.
The model 1, containing all four soil properties, explains 96.1% of the variation in dinoseb values. However, the significant levels of CEC, pH, and Clay content indicate that some of the soil properties can be removed from the model (significant levels are 0.999, 0.497, and 0.344, resp.). According to the removal principle, the soil property with highest significant level, which is CEC, should be removed and then Model 2 is built up with the remaining soil properties; in the same way, sequential stepwise regressions eliminated pH from model 2 since pH shows the highest significant level which is bigger than 0.05 (0.000 and 0.028, resp.). In model 3, both of the remaining variables show a significant level less than 0.05, thus elimination stops .
The statistics summary of each regression model is illustrated in Table 7. In addition to the three models in Table 6, Model 4 which uses only OC as a predictor is analysed. In all four models, the multiple correlations between the dinoseb values and predictors are strong ( varies from 0.961 to 0.945) and decrease slightly while one specific soil property is removed from the previous model. The changes from model 1 to model 2 and from model 2 to model 3 are not significant ( and 0.481, resp.). That means that removal of CEC and pH consecutively has minor effect on the goodness of the regression, whereas removal of Clay content from model 3 results in a significant change to . That also implies the clay factor is important for dinoseb sorption in soil.
3.5. Model Development
Combining the results from correlation analysis, path analysis, and stepwise regression, we can conclude that the soil OC and clay contents are the most important factors affecting the dinoseb sorption in soil. Therefore, the two factors are selected as the predictors of to build up the regression equation: in which OC is the soil organic content, and Clay is the clay content. The square is 0.92, that is, the variation of OC and clay in soil accounts for the 92% variation, in which 89% variation can be explained directly by OC variation, and the other 3% can be explained by the clay content variation. The -statics of the regression is 98.09, and the regression is found to be significant at .
A good multilinear regression was made, using all possible factors, including OC, Clay, CEC, and pH as the explanatory variables of the values. The sequence of correlation to was found to be OC > pH > Clay > CEC. However, The explanatory variables were not independent from each other, thus the multicollinearity may make the conclusion suspicious. With bivariate and correlation analysis and path analysis, it was found that the direct effects on should be OC > Clay > pH > CEC. Clay, pH, and CEC are mainly correlated with through the correlations with OC. The direct effects of pH and CEC on are very low. The zero-order correlation matrix shows that clay was positively correlated with , but the path analysis shows that the correlation is negative. The latter is more reasonable according to the dinoseb chemical properties. The backward stepwise regression showed that pH and CEC can be removed from the prediction model. Based on these results, a more efficient regression using OC and Clay as predictors is built up.
The authors gratefully acknowledge the financial support from the Special Fund of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering (2009585312), and the MOE Programme of Introducing Talents of Discipline to Universities (“111” Project B08048).
- M. Soutter and Y. Pannatier, “Groundwater vulnerability to pesticide contamination on a regional scale,” Journal of Environmental Quality, vol. 25, no. 3, pp. 439–444, 1996.
- H. J. O'Neill, T. L. Pollock, H. S. Bailey, P. Milburn, C. Gartley, and J. E. Richards, “Dinoseb presence in agricultural subsurface drainage from potato fields in Northwestern New Brunswick, Canada,” Bulletin of Environmental Contamination and Toxicology, vol. 43, no. 6, pp. 935–940, 1989.
- R. H. Kaake, D. J. Roberts, T. O. Stevens, R. L. Crawford, and D. L. Crawford, “Bioremediation of soils contaminated with the herbicide 2-sec-butyl-4,6- dinitrophenol (dinoseb),” Applied and Environmental Microbiology, vol. 58, no. 5, pp. 1683–1689, 1992.
- M. Matsumoto, T. Furuhashi, C. Poncipe, and M. Ema, “Combined repeated dose and reproductive/developmental toxicity screening test of the nitrophenolic herbicide dinoseb, 2-sec-butyl-4,6-dinitrophenol, in rats,” Environmental Toxicology, vol. 23, no. 2, pp. 169–183, 2008.
- N. Chèvre, A. R. Brazzale, K. Becker-van Slooten, R. Behra, J. Tarradellas, and H. Guettinger, “Modeling the concentration-response function of the herbicide dinoseb on Daphnia magna (survival time, reproduction) and Pseudokirchneriella subcapitata (growth rate),” Ecotoxicology and Environmental Safety, vol. 62, no. 1, pp. 17–25, 2005.
- K. L. Takahashi, H. Hojo, H. Aoyama, and S. Teramoto, “Comparative studies on the spermatotoxic effects of dinoseb and its structurally related chemicals,” Reproductive Toxicology, vol. 18, no. 4, pp. 581–588, 2004.
- M. Sreedhar, T. M. Reddy, K. R. irisha, and S. R. J. Reddy, “Differential pulse adsorptive stripping voltammetric determination of dinoseb and dinoterb at a modified electrode,” Analytical Sciences, vol. 19, no. 4, pp. 511–516, 2003.
- M. Pedrero, F. J. M. de Villena, J. M. Pingarrَn, and L. M. Polo, “Determination of dinoseb by adsorptive stripping voltammetry,” Electroanalysis, vol. 3, pp. 419–422, 1991.
- M. R. Viant, C. A. Pincetich, D. E. Hinton, and R. S. Tjeerdema, “Toxic actions of dinoseb in medaka (Oryzias latipes) embryos as determined by in vivo 31 NMR, HPLC-UV and 1H NMR metabolomics,” Aquatic Toxicology, vol. 76, no. 3-4, pp. 329–342, 2006.
- J. A. Arancibia, G. M. Delfa, C. E. Boschetti, G. M. Escandar, and A. C. Olivieri, “Application of partial least-squares spectrophotometric-multivariate calibration to the determination of 2-sec-butyl-4,6-dinitrophenol (dinoseb) and 2,6-dinitro-p-cresol in industrial and water samples containing hydrocarbons,” Analytica Chimica Acta, vol. 553, no. 1-2, pp. 141–147, 2005.
- A. W. Warrick and D. R. Nielsen, “Spatial variability of soil physical properties in the field,” in Application of Soil Physics, pp. 319–344, Academic Press, NewYork, NY, USA, 1980.
- L. Guo, W. A. Jury, R. J. Wagenet, and M. Flury, “Dependence of pesticide degradation on sorption: nonequilibrium model and application to soil reactors,” Journal of Contaminant Hydrology, vol. 43, no. 1, pp. 45–62, 2000.
- M. A. Maraqa, X. Zhao, R. B. Wallace, and T. C. Voice, “Retardation coefficients of nonionic organic compounds determined by batch and column techniques,” Soil Science Society of America Journal, vol. 62, no. 1, pp. 142–152, 1998.
- B. Thompson, “Stepwise regression and stepwise discriminant analysisi need not apply here: a guidelines editorial,” Education and Psychological Measurement, vol. 55, no. 4, pp. 525–534, 1995.
- Y. Dodge, The Oxford Dictionary of Statical Terms, 2003.
- M. Montanaro Gauci, T. F. Kruger, K. Coetzee, K. Smith, J. P. Van Der Merwe, and C. J. Lombard, “Stepwise regression analysis to study male and female factors impacting on pregnancy rate in an intrauterine insemination programme,” Andrologia, vol. 33, no. 3, pp. 135–141, 2001.
- M. Soutter and Y. Pannatier, “Groundwater vulnerability to pesticide contamination on a regional scale,” Journal of Environmental Quality, vol. 25, no. 3, pp. 439–444, 1996.
- A. G. Hornsby, R. Don Wauchope, and A. E. Herner, Pesticide Properties in the Environment, Springer, New York, NY, USA, 1996.
- I. M. M. Ghani and S. Ahmad, “Stepwise multiple regression method to forecast fish landing,” Procedia Social and Behavioral Sciences, vol. 8, pp. 549–554, 2010.
- K. L. Findell and E. A. B. Eltahir, “An analysis of the soil moisture-rainfall feedback, based on direct observations from Illinois,” Water Resources Research, vol. 33, no. 4, pp. 725–735, 1997.
- D. C. Montgomery, E. A. Peck, and G. G. Vining, Introduction to Linear Regression Analysis, John Wiley and Sons, Hoboken, NJ, USA, 4th edition, 2006.
- I. G. Chong and C. H. Jun, “Performance of some variable selection methods when multicollinearity is present,” Chemometrics and Intelligent Laboratory Systems, vol. 78, no. 1, pp. 103–112, 2005.
- M. S. Lachniet and W. P. Patterson, “Use of correlation and stepwise regression to evaluate physical controls on the stable isotope values of Panamanian rain and surface waters,” Journal of Hydrology, vol. 324, no. 1–4, pp. 115–140, 2006.
- P. S. Petraitis, A. E. Dunham, and P. H. Niewiarowski, “Inferring multiple causality: the limitations of path analysis,” Functional Ecology, vol. 10, no. 4, pp. 421–431, 1996.
- D. J. de Rodríguez, J. L. Angulo-Sánchez, and R. Rodríguez-García, “Correlation and path coefficient analyses of the agronomic trait of a native population of guayule plants,” Industrial Crops and Products, vol. 14, pp. 93–103, 2001.
- A. G. G. Vasconcelos, R. M. V. Almeida, and F. F. Nobre, “The path analysis approach for the multivariate analysis of infant mortality data,” Annals of Epidemiology, vol. 8, no. 4, pp. 262–271, 1998.
- S. B. Haderleln and R. P. Schwarzenbach, “Adsorption of substituted nitrobenzenes and nitrophenols to mineral surfaces,” Environmental Science and Technology, vol. 27, no. 2, pp. 316–326, 1993.
- H. C. Cho and S. Abe, “Is two-tailed testing for directional research hypotheses tests legitimate?” Journal of Business Research, 2012.
- L. Cox, M. C. Hermosin, and J. Cornejo, “Distribution coefficient of methomyl in soils from different depths,” Fresenius Environmental Bulletin, vol. 1, pp. 445–449, 1992.
- D. T. Chiou, “Theoretical considerations of partition uptake of non-ionic organic compounds by soil organic matter. In Reactions and movement of organic chemicals in soils,” SSSA Special Publication, no. 22, pp. 1–29, 1989.
- M. C. Hermosin, J. Cornejo, and J. L. Perez-rodriguez, “Adsorption and desorption of maleic hidrazide as a function of soil properties,” Soil Science, vol. 144, pp. 250–256, 1987.
- M. C. Hermosin and J. Cornejo, “Assessing factors related to pesticide adsorption by soils,” Environmental Toxicology and Chemistry, vol. 25, pp. 45–55, 1989.
- W. Lertpaitoonpan, S. K. Ong, and T. B. Moorman, “Effect of organic carbon and pH on soil sorption of sulfamethazine,” Chemosphere, vol. 76, no. 4, pp. 558–564, 2009.
- A. Mermoud, J. M. F. Martins, D. Zhang, and A. C. Favre, “Small-scale spatial variability of atrazine and dinoseb adsorption parameters in an alluvial soil,” Journal of Environmental Quality, vol. 37, no. 5, pp. 1929–1936, 2008.
- S. Jagadamma, R. Lal, R. G. Hoeft, E. D. Nafziger, and E. A. Adee, “Nitrogen fertilization and cropping system impacts on soil properties and their relationship to crop yield in the central Corn Belt, USA,” Soil and Tillage Research, vol. 98, no. 2, pp. 120–129, 2008.
- B. L. Bowerman and R. T. O' Connell, Linear Statistical Models: An Applied Approach, Duxbury Press, Belmont, Ca, USA, 2nd edition, 1990.