Next Article in Journal
The Effect of Empirical Log Yield Observations on Carbon Storage Economics
Next Article in Special Issue
Low Population Differentiation but High Phenotypic Plasticity of European Beech in Germany
Previous Article in Journal
Assessment of Forest Ecosystem Development in Coal Mine Degraded Land by Using Integrated Mine Soil Quality Index (IMSQI): The Evidence from India
Previous Article in Special Issue
Improving Environmental and Energy Efficiency in Wood Transportation for a Carbon-Neutral Forest Industry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Algorithm Selection on Modeling Ozone Pollution: A Perspective on Box and Tiao (1975)

1
Faculty of Administration and Business, Research Institute of the University of Bucharest–ICUB, University of Bucharest, 050663 Bucharest, Romania
2
College of Science and Theoretical Studies, Saudi Electronic University, Jeddah 23442, Saudi Arabia
3
College of Forestry, Oregon State University, Corvallis, OR 97331, USA
4
Faculty of Silviculture and Forest Engineering, Transilvania University, 500123 Brasov, Romania
*
Author to whom correspondence should be addressed.
Forests 2020, 11(12), 1311; https://doi.org/10.3390/f11121311
Submission received: 29 October 2020 / Revised: 3 December 2020 / Accepted: 4 December 2020 / Published: 9 December 2020

Abstract

:
Estimation using a suboptimal method can lead to imprecise models, with cascading effects in complex models, such as climate change or pollution. The goal of this study is to compare the solutions supplied by different algorithms used to model ozone pollution. Using Box and Tiao (1975) study, we have predicted ozone concentration in Los Angeles with an ARIMA and an autoregressive process. We have solved the ARIMA process with three algorithms (i.e., maximum likelihood, like Box and Tiao, conditional least square and unconditional least square) and the autoregressive process with four algorithms (i.e., Yule–Walker, iterative Yule–Walker, maximum likelihood, and unconditional least square). Our study shows that Box and Tiao chose the appropriate algorithm according to the AIC but not according to the mean square error. Furthermore, Yule–Walker, which is the default algorithm in many software, has the least reliable results, suggesting that the method of solving complex models could alter the findings. Finally, the model selection depends on the technical details and on the applicability of the model, as the ARIMA model is suitable from the AIC perspective but an autoregressive model could be preferred from the mean square error viewpoint. Our study shows that time series analysis should consider not only the model shape but also the model estimation, to ensure valid results.

1. Introduction

Environmental investigations rely on evidence, which, since the empirical movement of Hume [1] and Locke [2] and the creation of the hypotheses testing by Fisher [3] and Pearson [4], requires formal analyses. Currently, the foundation of most of the analyses are procedures implemented using various software, such as R [5], SAS [6], or Matlab [7]. Essentially, numerical algorithms are pivotal in solving any problem involving computations. In most instances, one problem can be solved with more than one algorithm, but because the aim of most studies is usually to find the solution not the influence of the procedural details, selection of the best algorithm among the available options is customarily ignored. However, depending on the problem under investigation, an algorithm can be more accurate than other. Unfortunately, the entrenched approach of data analysis is to use the default options implemented in the software selected for analysis. Therefore, from a procedural perspective, the significance of studying the impact of different numerical algorithms on the results becomes of major importance. Many studies have been focused on the topic of algorithm impact on the findings [8], most of them deciding on which algorithm to choose based upon on a predefined criteria, such as bias, root mean square error, Akaike Information Criterion (AIC), or Schwartz Bayesian Information Criterion (BIC).
The wealth of data characterizing the last two decades promoted the usage of time series methods in studying the environmental processes. Currently there are many other processes than can be used in modeling time series data, the most popular one being the autoregressive (AR) and moving average (MA) models, as well as their combination (ARMA) [9,10]. Furthermore, several model enhanced ARMA [11], particularly ARIMA (autoregressive integrated moving average) [12] and ARFIMA (autoregressive fractionally integrated moving average) [13]. The stochastic approach to time series analysis based on ARMA was enhanced with artificial neural network (ANN) [14,15,16]. However, ANN can provide solution in an unfeasible amount of time, because of the iterative training phase. To alleviate the sequential requirements associated with ANNs, non-iterative supervised learning procedures were developed, such as the one based on the Ito decomposition and the neural-like structure of the successive geometric transformations model [17]. ANNs algorithms use more parameters than ARIMA, and therefore are more sensitive to the procedure used for estimation. Therefore, the most common approaches of studying time series are based on ARMA processes [18].
The plethora of software dedicated to the analysis of repeated measurements made the investigation seemingly challenge less to the untrained investigator. Nevertheless, time series analyses are not only difficult to model because of the inherent dependencies within the data [19,20], but also because they rely on technicalities [8]. Finding the appropriate model seems to be paramount when autocorrelated data are used [21,22], particularly for attributes with significant impact on land management, such as site productivity. In this paper we focused on time series modeling not only because it is challenging [23] but also because it plays a significant role in many scientific fields [24]. The most complicated implementation of time series models are argued to be in physical and environmental sciences [10]. Therefore, the goal of our study is to compare the impact of numerical algorithms on the model development rather than comparing the performances of the models. Among the possible models, we have focused on ozone pollution [25], which within the current climate change paradigm, adds undesired stresses to terrestrial ecosystems, particularly forests [26,27]. Therefore, our study will recommend an approach for modeling ozone pollution that is tailored to data availability, preset conditions (e.g., minimum AIC or mean square error), and parsimony [28].

2. Materials and Methods

To substantiate the idea of this study we have used a time series data, namely the famous intervention model for ozone of Box and Tiao [29]. We have focused on the Box and Tiao [29] study, as it is selected as a sample dataset in many teaching environments and software. The dataset represents the oxidant pollution in downtown Los Angeles from 1955 to 1972. Ozone is a chemical substance, referred as O3, which in high concentration leads to pollution [30]. Box and Tiao [29] used an autoregressive integrated moving average (ARIMA) process to model the oxidant pollution. The ARFIMA process are not suitable for modeling the downtown Los Angeles data used in the present study because it is appropriate for stationary processes whose autocorrelation function decay slowly; therefore a large number of observation is needed, which is not the case of the Box and Tiao [29] example.
The algorithm chosen by Box and Tiao to estimate the parameters of the ARIMA model was the maximum likelihood (ML). However, other algorithms are available for solving an ARIMA process. Therefore, we have chosen two other algorithms to solve ARIMA besides the ML, namely the conditional least square (CLS) and the unconditional least square (ULS) [31,32,33]. Besides ARIMA, the autoregressive process is commonly used to model time series data. Similarly to ARIMA, there are multiple algorithms that can be used in estimation, and in this study, we have used four [34]: Yule–Walker (YW), Iterative Yule–Walker (ITYW), maximum likelihood (ML), and unconditional least square (ULS).

2.1. Data

To ensure consistency of the assessment of the impact of the algorithms on the estimates, we have used the same ozone concentration data from downtown Los Angeles that was used by Box and Tiao [29]. The data contains the ozone concentration from 1955 to 1972. A downward trend is clearly present (Figure 1), with a peak ozone level of 8.7, observed in September of 1956 and the lowest is 1.3 in January and December of 1975, and December of 1977.
To model the ozone concentration, Box and Tiao considered two other data sources, one driven by the regulation and one by the climate. In early 1960s the Golden State Freeway opened, triggering new laws aiming at reduction of the fraction of hydrocarbons in the gasoline sold. The new laws were created to reduce the allowable proportion of hydrocarbons in the gasoline, which in turn affected the ozone level. A major change in the regulations governing the Los Angeles traffic, consequently with possible impact on ozone, occurred in 1960 when the rule 63 was enacted. To climate data was simply the separation of the calendar 12 month into warm and cold seasons, which conveniently were labeled summer and winter. To account for the major change in the design of the gasoline engine that occurred in the mid 1960s, which imposed a further reduction in the production of O3, Box and Tiao identified the seasonal variable only after 1966. Therefore, the winter variable selects the months from November to May between 1966 and 1972, and the summer variable the months from June to October between 1966 and 1972. The introduction of ancillary data is justified not only by the enforcement of new regulations but also by the summary statistics, which shows a change from 1955 to 1972 (Table 1).

2.2. Models and Algorithms Description

Box and Tiao introduced three additional binary variables to model the ozone concentration: (1) X1 (Equation (1)), which is called the intervention variable, as it represents the major changes that occurred in the auto regulations, (2) winter (Equations (2) and (3)) summer (Equation (3)).
X 1 = { 1 b e f o r e   J a n u a r y   1960 0 a f t e r   J a n u a r y   1960          
s u m m e r = { 1 t h e   m o n t h s   J u n e O c t o b e r   b e g i n n i n g   1966 0 o t h e r w i s e
w i n t e r = { 1 t h e   m o n t h s   N o v e m b e r M a y   b e g i n n i n g   1966 0 o t h e r w i s e
The summer and winter are also considered intervention variables because the temperature and sunlight intensity are affected by the presence of the oxidant pollution. Therefore, the change in regulations, and consequently in technology, influenced not only the pollution, but also the microclimate.
To remove the seasonal variation from the ozone concentration, the ozone series was differenced at lag 12, as was done in the study by Box and Tiao [29]. The value 12 was selected because the data was recorded monthly and the seasonality was annual. The time-series of the ozone level during the time of 1955 to 1973 exhibits a clear seasonal component (Figure 1), which after differencing, disappeared (Figure 2).
Our study starts with the work of Box and Tiao, which used the ARIMA procedure in combination with maximum likelihood (ML) to develop a model for the ozone concentration. However, we will expand their work by including besides ML, parameters of the same model estimated with two additional algorithms: CLS and ULS. We compared the performances of the algorithms with the Akaike Information Criteria (AIC).
We have tested the necessity of using ARIMA by developing the simpler autoregressive model without the moving average term. The autoregressive model explains the current value of the series based only of number p past values, AR (p). Similarly to ARIMA, we investigated the performances of four algorithms in modeling the ozone concentration. The results supplied by each algorithm were compared using the AIC and coefficient of determination (R2).

2.2.1. ARIMA Model and Estimation Algorithms

Box and Tiao [29] reached the conclusion that an ARIMA process (Equation (4)) is suitable for modeling the ozone concentration in downtown Los Angeles:
o z o n e = β 0 + β 1 x 1 + β 2 S u m m e r + β 3 W i n t e r + ε t
where ozone is the concentration of zone in pphm, X1, Summer, and Winter are the intervention variables, and ε is the noise or the stochastic variation.
Box and Tiao used the autocorrelation function to identify the lag needed to eliminate the seasonal behavior of ozone concentration. They found that lag 12 of the ozone leads to a significant correlation at lags 1 and 12 for noise, which prompted the model:
( 1 B 12 ) ε t = ( 1 θ 1 B ) ( 1 θ 2 B 12 ) e t
where B is the backshift operator of lag 1, such that Byt = yt−1, et is the white noise, namely an independently distributed normal variable with mean 0 and variance σa, regardless the time t, θ(B) = 1 − θ1B − θ2B2 − … −θqBq is the moving average polynomial of order q. The model fitted by Box and Tiao was
o z o n e t = β 1 x 1 t + β 2 s u m m e r t ( 1 B 12 ) + β 3 w i n t e r t ( 1 B 12 ) + ( 1 θ 1 B ) ( 1 θ 2 B 12 ) ( 1 B 12 ) e t .
The details of each estimation algorithm are presented in multiple sources [31,32,35,36,37], but in this paper we will just summarize the main aspect of each algorithm. The CLS estimation is conditional on the assumption that past unobserved errors are equal to 0 and it produces estimates by minimizing the Equation (7):
t = 1 n a t 2 = t = 1 n ( x t i = 1 π i x t i ) 2
where π is computed from the estimates of φ and θ at each iteration, which are the AR and MA parameters for the series x 1 , , x t 1 , x t , and t = 1 n a t 2 is e × e where e is the white noise, and e’ is the transpose of e.
The ULS estimates start from the estimates computed by using CLS algorithm, which are further minimized (Equation (8)):
t = 1 n a t 2 = t = 1 n ( x t C t V t 1 ( x 1 , , x t 1 ) ) 2
where Ct is the covariance matrix and V t is the variance matrix for the series x 1 , , x t 1 , x t .
The ML algorithm maximizes the likelihood function through nonlinear least square Marquardt’s method (Equation (9)). The initial estimates are computed using the CLS.
| H | 1 n e e | H | 1 n
where e is white noise vector and H is the lower triangular matrix with positive elements on the diagonal, such that H H = Ω where Ω is the determinant of the regression equation.
The three algorithms were assessed, with AIC (Equation (10)), which is a measure of the goodness of fit of the model as well as the complexity of that model.
AIC = − 2ln (L) + 2k
where L is the maximized value of the likelihood function for the estimated model and k is the number of estimated parameters.

2.2.2. Autoregressive Model and Estimation Algorithms

The autoregressive process used to model the time series is:
Y t = X t β + ε t
where Yt is the response variable, X’t is the matrix of predictor variables with β slope, and εt is the error term.
For consistency with the ARIMA model, we have considered an autoregression model that includes as predictor X1, summer, and winter (Equation (12)):
o z o n e = β 1 x 1 + β 2 s u m m e r + β 3 w i n t e r + ε t
where ozone is the predicted and the interventions X1, summer and winter are the predictors.
Cleveland [38] suggested that autocorrelation, inverse auto correlation and partial auto correlation functions should be used to determine the autoregressive order. Using the three autocorrelation functions we found 12 as the lag fulfilling the modeling requirements, confirming the ARIMA approach. To remove the possible correlation among residuals still present after differentiation, we have further differentiated the residuals:
ε t = ϕ 1 ε t 1 + ϕ 2 ε t 2 + + ϕ l ε t l + e t
where l is the lag and et is white noise.
The final model became:
12 o z o n e = β 1 12 x 1 + β 2 s u m m e r + β 3 w i n t e r + e t
where 12 o z o n e and 12 x 1 are seasonally differenced at lag 12 for Ozone and X1, respectively.
To estimate the parameters of the autoregressive model 12 we have considered four algorithms: Yule–Walker (YW), Iterative Yule–Walker (ITYW), maximum likelihood (ML), and unconditional least square (ULS).
The Yule–Walker algorithm starts from the ordinary least square estimates of β, which is used to estimates Φ, the parameters of an AR (p) process, from the sample autocorrelation function of the residuals. The equation supplying the Φ ^ is:
Φ ^ = R 1 r
where Φ ^ is the vector of autoregressive parameters, r is the vector (r1, …,rp), where ri is the lag i sample autocorrelation, R is the Toeplitz matrix with the elements ri,j=r|i-j|.
The Iterative Yule–Walker algorithm uses the resulted residuals stemmed from the YW algorithm to create new estimators of ϕ.
The ML algorithm is based on the likelihood function 16.
L ( ε , β , σ 2 ) = i = 1 n f ( ε i ) = 1 ( 2 π ) n 2 σ n e ( 1 2 σ 2 ε ε )
where εi are the residuals, β are the parameter, and σ2 is the variance of ε.
The likelihood function for an AR(1) model is
L ( ε , β , σ 2 ) = ( 2 π σ 2 ) n 2 1 ϕ 2 e 1 2 σ 2 S ( ϕ , μ )
S ( ϕ , μ ) = t = 2 n [ ( y t μ ) ϕ ( y t 1 μ ) ] 2 + ( 1 ϕ 2 ) ( y 1 μ ) = S c ( ϕ , μ ) + ( 1 ϕ 2 ) ( y 1 μ ) .
To maximize the likelihood function, it is customary to execute the computation on the logarithm of L rather than on L, after which partial derivatives are computed in respect with σ2. The derivation of likelihood function is in general considerably more involved, which requires estimations with numerical methods.
The ULS algorithm, is a compromise between conditional least square estimates and maximum likelihood, and is focused on minimization of the S(ϕ,μ) (Equation (18)):
Besides AIC, we have assessed the model using the mean square error (MSE) and the coefficient of determination (R2) [39].
Mean   Square   Error   M S E = i = 1 n ( Y i Y ^ i ) 2 n p
Coefficient   of   determination   R 2 = 1 i = 1 n ( Y i Y ^ i ) 2 / i = 1 n ( Y i Y ¯ ) 2
where Y i , Y ^ i , Y ¯ are the actual, the predicted, and the average predicted variable, n is the total number of observations and p is the number of parameters in each model.
We considered the analysis concluded when the residuals of the time series model are white noise, which was tested with the Durbin–Watson test. All computations were executed using SAS 9.4 [6].

3. Results

Our results confirm the findings of Box and Tiao, as all variables included in the ARIMA model were significant (p < 0.001), irrespective the algorithm except for the Winter variable, which was significant only for α = 0.1 (Table 2).
The results suggest that all variables have negative impact on ozone concentration. The highest AIC and SE values, therefore the worse performances, were supplied by the CLS algorithm, which is the default option in SAS. This finding suggests that the default algorithm is not always the best choice for estimation of the parameters defining a model. The best result, according to AIC, was achieved when ML algorithm was used, which explains the choice of Box and Tiao in solving the ARIMA process. However, ULS supplied the smallest mean square error among the tree, suggesting that different criteria could lead to different findings.
We forecasted the year 1973 with 95% confidence interval (Figure 3). For consistency, we predicted the monthly ozone concentration using the same three algorithms: ML, CLS, and ULS (Table 3). We noticed that the MSE for the forecasted months with the CLS algorithm were the highest among the three algorithms, which is an additional verification that the default option should not always be the preferred estimation choice.
The autoregressive model described by Equation (14) and solved with the four algorithms also revealed the importance of the algorithm in estimation of the parameters (Table 4). Consequently, the model performance depends on the algorithm chosen to estimates the parameters. The results show that YW algorithm, which is the default algorithm in SAS for the autoregressive processes, had the largest MSE and AIC, whereas ULS the smallest MSE and ML the smallest AIC (Table 4). The AR findings provides another supporting evidence that the ML algorithm could be the best choice, and that Box and Tiao selected the appropriate approach for solving the ozone concentration model. The Durbin–Watson test [40] revealed that the errors produced by the autoregressive model 14 are white noise (Table 4). Visual investigation and Durbin–Watson test on the residual produced by the autoregressive model reveals no pattern, and supports their distribution as white noise (Figure 4).
Autocorrelation Function (ACF), Partial Autocorrelation Function (PACF), and Inverse Autocorrelation Function (IACF) support the white noise distribution of the residuals produced by the autoregressive model (Figure 5) and revealed by the Durbin–Watson test (Table 4), as all correlation functions where less than two standard errors. Therefore, the ACF variation suggest that there is no need for additional autoregressive terms whereas the PACF and IACF indicate that there is no need to reduce the model parsimony by including moving average terms in modeling the ozone pollution.
Irrespective of the algorithm, the autoregressive model of Equation 14 produced similar results with ARIMA, in the sense that the winter intervention is not significant (p > 0.5). However, the autoregressive model also shows that the summer intervention is not significant, as the p-value > 0.25. In fact, the autoregressive model shows that the differenced intervention variable is the only significant variable.

4. Discussion

The two modeling approaches of ozone concentration (i.e., ARIMA and autoregressive) lead to models that cannot be challenged based on the achievement of white noise, as both produced the expected distribution of the residuals. However, if AIC was used as a criterion, then the ARIMA supplied superior results then the autoregressive (i.e., 501 vs 527), whereas if MSE is the guiding criterion then the opposite conclusion is reached (i.e., 0.763 vs 0.768). The conflicting findings should be considered in the context of the desired objective, as if the best overall model is needed then the AIC criterion should prevail, while if precise predictions are the focus, then the MSE criterion should be emphasized.
Irrespective of the approach on modeling the ozone concentration, ARIMA or autoregressive, the ML algorithm yielded the best results, according to AIC. However, if MSE is used as criterion, then the ML is not the best, but the least square method is (i.e., ULS for autoregressive models and CLS for ARIMA model), which suggests that for complex models, algorithms that are computational more efficient than ML could be used. Nevertheless, the default algorithms, Yule–Walker for the autoregressive approach and CLS for ARIMA, supplied the worse results, which points toward the necessity of trying various algorithms before reaching a conclusion.
A natural question arises whether ARIMA, a less parsimonious model than AR in term of structure, should be used instead of AR, given that not only similar results are obtained but also ARIMA is an enhancement of AR. The answer lies partially in the criteria used for model selection (i.e., AIC vs. MSE), partially in the overall parsimony, and partially in the predictive ability of the models. Smaller MSE suggests a more performant model, given the central limit theorem, but in the presence of outliers the least square estimator can produce biased results. Furthermore, selection of a model based on MSE offers no warranty that the model is not overfitted [41]. The AIC on the other hand, being based on the likelihood function, not only that is more robust to outliers, but is also an unbiased estimator of the Klluback discrepancy [41], a measure of model agreements. Therefore, in selection of a model, AIC criteria should be preferred over the MSE, as information based criteria consistently supply correct results [42]. In the ozone pollution case the AIC was more than 5% smaller for ARIMA model than the AR model, which suggest the appropriateness of the more parsimonious model (i.e., ARIMA). From the parsimony perspective, ARIMA model uses six parameters whereas the AR model uses five, which suggest the AR as the appropriate model. However, parsimony on itself is not sufficient, most of the time having more a discriminating power rather than being a selective criterion. From the predictive ability perspective, ARIMA is more flexible than AR, because it examines simultaneously the autoregressive part (i.e., variables) and the autocorrelation part (i.e., errors), whereas the AR considers only the autoregressive properties of the time series. Therefore, Box and Tiao selected ARIMA as the model because it meets several desired properties, namely robustness, predictability, and model agreement. However, a stationary process that does not present outliers could be modeled according to the MSE and parsimony, which is the case of AR. Finally, the choice of the model and the algorithm to estimate the parameters defining the model is not unique to ozone pollution or the dataset used in this study. Other studies pointed the importance of the algorithm selection [8,43,44], but not in forestry settings.
The main finding of our study suggests that careless usage of an estimation algorithm, even for simple problems, such as the ozone concentration in downtown Los Angeles that has 216 observations, could lead to questionable results. The conclusion holds irrespective of the approach, simpler, such as autoregressive, or more complex, such as ARIMA. Furthermore, depending on the objective, complex models even superior from some perspectives, could not be appropriate, particularly when considered from the parsimony perspective. Therefore, besides technical details that define the model, operational constraints should be considered, as the AIC should not be applied blindly in model selection.

5. Conclusions

Modeling is considered an area where subjective evaluations are not present, where objectivity rules, which is based on repeatable results. Objectivity and repeatability is of paramount importance for complex models, when not only the accuracy is important but also the precision. However, a model can be developed in multiple ways. Considering that usually models are part of a larger research question, the aim of such studies consists in finding the answer to the question of interest rather than on procedural details. Therefore, technical detail, such as selection of the most appropriate algorithm among the available options, is customarily ignored. Estimation with an inappropriate algorithm could lead to imprecise models, which have cascading effects in complex models such as climate change or pollution. The goal of this study is to compare the solutions supplied by different algorithms used to model a time series, with application to ozone pollution, of significant importance in the current climate change driven paradigm. We have focused on ozone as pollutant to be modeled, as it is important not only for forestry but also for climate change overall. We have used the ozone concentration data from the 1975 study of Box and Tiao, as it is well studied and provides the groundwork to modeling times series data. In addition to the maximum likelihood, which was the algorithm used by Box and Tiao, we have considered two other algorithms to solve the ARIMA process, namely conditional least square and unconditional least square. As an alternative to the ARIMA process, we have modeled the ozone concentration data with a more parsimonious approach, namely an autoregressive process. We have estimated the parameters of the autoregressive model with four numerical algorithms: Yule–Walker, Iterative Yule–Walker, maximum likelihood, and unconditional least square. Irrespective the approach, the ML and ULS produced the most suitable models, within 1% from each other, for MSE, and 5% for AIC. The larger difference in AIC compared with MSE, as well as the robustness to outliers, recommends ARIMA over AR in modeling the ozone pollution. Nevertheless, for a stationary process that is identified using data without outliers an AR model could be more appropriate, as MSE supplied by the ULS algorithm was always the smallest, regardless the model (i.e., 0.763 against 0.764 of ML for AR and 0.768 vs 0.797 of ML for ARIMA).
Our study proves that Box and Tiao’s choice of ML algorithm is the most suitable selection for working on the ozone data from downtown Los Angeles, according to AIC but not according to mean square error criterion. Furthermore, we found that Yule–Walker, which is the default algorithm in many software used for time series modeling, has the least reliable results, suggesting that the method of solving complex models could alter the findings. Finally, the model selection depends on the technical details as well as on the applicability of the model, as ARIMA model is suitable from AIC perspective but an autoregressive model could be preferred from mean square error viewpoint.
The present study proposes indirectly a procedure of selecting the algorithm that would produce the most appropriate results. However, the main limitation of the presented approach is the requirement to run multiple algorithms for multiple models, which, in eventuality of large dataset as inputs, could become time consuming. Therefore, future research should focus on the search of the algorithm that fits the model, the data, and the evaluation criteria.

Author Contributions

Conceptualization, B.M.S.; methodology, B.M.S. and N.G.; software, N.G. and M.P.; validation, N.G. and M.P.; formal analysis, N.G. and B.M.S.; investigation, M.P.; resources, M.P.; writing—original draft preparation, B.M.S. and N.G.; writing—review and editing, M.P. and N.G.; supervision, B.M.S.; project administration, M.P.; funding acquisition, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

The study was partially funded by the National Institute of Food and Agriculture, U.S. Department of Agriculture, McIntire Stennis Project OREZ-FERM-875, and by the Romanian ANCSI Project POC P-37–257.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hume, D. Enquiries Concerning the Human Understanding: And Concerning the Principles of Morals, 2nd ed.; Oxford University Press: Oxford, UK, 1902. [Google Scholar]
  2. Locke, J. An Essay Concerning Human Understanding: And a Treatise on the Conduct of the Understanding; Hayes &Zell: Philadelphia, PA, USA, 1860. [Google Scholar]
  3. Fisher, R.A. Statistical Methods for Research Workers; Oliver and Boyd: Edinburgh, UK, 1925. [Google Scholar]
  4. Pearson, K. On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1900, 50, 157–175. [Google Scholar] [CrossRef] [Green Version]
  5. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  6. SAS Institute. SAS; SAS Institute: Cary, NC, USA, 2017. [Google Scholar]
  7. The MathWorks Inc. Matlab; The MathWorks Inc: Natick, MA, USA, 2017. [Google Scholar]
  8. Seppelt, R.; Richter, O. “It was an artefact not the result”: A note on systems dynamic model development tools. Environ. Model. Softw. 2005, 20, 1543–1548. [Google Scholar] [CrossRef]
  9. Brockwell, P.J.; Davis, R.A. An Introduction to Time Series and Forecasting; Springer texts in statistics (Springer): New York, NY, USA, 1996. [Google Scholar]
  10. Schumway, R.H.; Stoffer, D.S. Time Series Analysis and Its Applications, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  11. Wold, H. A Study in the Analysis of Stationary Time Series, 2nd ed.; Almqvist and Wiksell Book Co: Uppsala, Sweden, 1954. [Google Scholar]
  12. Box, G.E.P.; Jenkins, G.M.; Reinsel, G.C. Time Series Analysis: Forecasting and Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1994. [Google Scholar]
  13. Granger, C.W.J.; Joyeux, R. An Introduction to Long-Memory Time Series Models and Fractional Differencing. J. Time Ser. Anal. 1980, 1, 15–29. [Google Scholar] [CrossRef]
  14. Izonin, I.; Ml, M.G.; Tkachenko, R.; Logoyda, M.; Mishchuk, O.; Kynash, Y. SGD-Based Wiener Polynomial Approximation for Missing Data Recovery in Air Pollution Monitoring Dataset. Adv. Comput. Intell. 2019, 11506, 781–793. [Google Scholar]
  15. Izonin, I.; Kryvinska, N.; Tkachenko, R.; Zub, K. An Approach towards Missing Data Recovery within IoT Smart System. Procedia Comput. Sci. 2019, 155, 11–18. [Google Scholar] [CrossRef]
  16. Izonin, I.; Kryvinska, N.; Tkachenko, R.; Zub, K.; Vitynskyi, P. An Extended-Input GRNN and its Application. Procedia Comput. Sci. 2019, 160, 578–583. [Google Scholar] [CrossRef]
  17. Tkachenko, R.; Izonin, I.; Vitynskyi, P.; Lotoshynska, N.; Pavlyuk, O. Development of the Non-Iterative Supervised Learning Predictor Based on the Ito Decomposition and SGTM Neural-Like Structure for Managing Medical Insurance Costs. Data 2018, 3, 46. [Google Scholar] [CrossRef] [Green Version]
  18. Kirchgässner, G.; Wolters, J.; Hassler, U. Introduction to Modern Time Series Analysis; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  19. Diggle, P.; Heagerty, P.; Liang, K.Y.; Zeger, S.L. Analysis of Longitudinal Data; Oxford University Press: New York, NY, USA, 2002. [Google Scholar]
  20. Fitzmaurice, G.M.; Laird, N.M.; Ware, J.H. Applied Longitudinal Analysis; Wiley: Hoboken, NJ, USA, 2004; ISBN 978-0-471-21487-8. [Google Scholar]
  21. Cieszewski, C.J. Developing a Well-Behaved Dynamic Site Equation Using a Modified Hossfeld IV Function Y3 = (axm)/(c + xm – 1), a Simplified Mixed-Model and Scant Subalpine Fir Data. For. Sci. 2003, 49, 539–554. [Google Scholar] [CrossRef]
  22. Cieszewski, C.J.; Strub, M. Comparing properties of self-referencing models based on nonlinear-fixed-effects versus nonlinear-mixed-effects modeling approaches. MCFNS 2018, 10, 46–57. [Google Scholar]
  23. Gujarati, D.N. Basic Econometrics; McGraw-Hill Companies: New York, NY, USA, 1995. [Google Scholar]
  24. Meek, C.; Chickering, D.; Heckerman, D. Autoregressive Tree Models for Time-Series Analysis. In Proceedings of the Second SIAM International Conference on Data Mining, Arlington, VA, USA, 11–13 April 2002; pp. 229–244, ISBN 978-0-89871-517-0. [Google Scholar]
  25. Krzyzanowski, J.; McKendry, I.G.; Innes, J.L. Evidence of Elevated Ozone Concentrations on Forested Slopes of the Lower Fraser Valley, British Columbia, Canada. Water Air Soil Pollut. 2006, 173, 273–287. [Google Scholar] [CrossRef]
  26. Laurence, J.A.; Retzlaff, W.A.; Kern, J.S.; Lee, E.H.; Hogsett, W.E.; Weinstein, D.A. Predicting the regional impact of ozone and precipitation on the growth of loblolly pine and yellow-poplar using linked TREGRO and ZELIG models. For. Ecol. Manag. 2001, 146, 247–263. [Google Scholar] [CrossRef]
  27. Subramanian, N.; Karlsson, P.-E.; Bergh, J.; Nilsson, U. Impact of Ozone on Sequestration of Carbon by Swedish Forests under a Changing Climate: A Modelling Study. For. Sci. 2015, 61, 445–457. [Google Scholar] [CrossRef] [Green Version]
  28. Strimbu, B.M.; Amarioarei, A.; Paun, M. A parsimonious approach for modeling uncertainty within complex nonlinear relationships. Ecosphere 2017, 8, e01945. [Google Scholar] [CrossRef] [Green Version]
  29. Box, G.E.P.; Tiao, G.C. Intervention Analysis with Applications to Economic and Environmental Problems. J. Am. Stat. Assoc. 1975, 70, 70–79. [Google Scholar] [CrossRef]
  30. Karlsson, P.E.; Pleijel, H.; Belhaj, M.; Danielsson, H.; Dahlin, B.; Andersson, M.; Hansson, M.; Munthe, J.; Grennfelt, P. Economic Assessment of the Negative Impacts of Ozone on Crop Yields and Forest Production. A Case Study of the Estate Östads Säteri in Southwestern Sweden. Ambi 2005, 34, 32–40. [Google Scholar] [CrossRef]
  31. Jacob, C. Conditional least squares estimation in nonstationary nonlinear stochastic regression models. Ann. Statist. 2010, 38, 566–597. [Google Scholar] [CrossRef] [Green Version]
  32. Ansley, C.F.; Newbold, P. Finite sample properties of estimators for autoregressive moving average models. J. Econom. 1980, 13, 159–183. [Google Scholar] [CrossRef]
  33. Davidson, J.E.H. Problems with the estimation of moving average processes. J. Econom. 1981, 16, 295–310. [Google Scholar] [CrossRef]
  34. Harvey, A.C. Time Series Models; MIT Press: Cambridge, MA, USA, 1981. [Google Scholar]
  35. Hartley, H.O.; Rao, J.N.K. Maximum Likelihood Estimation for Mixed Analysis Variance Model. Ann. Math. Stat. 1965, 36, 1610. [Google Scholar]
  36. Dolby, G.R.; Lipton, S. Maximum Likelihood Estimation of General Nonlinear Functional Relationship with Replicated Observations and Correlated Errors. Biometrika 1972, 59, 121–129. [Google Scholar] [CrossRef]
  37. Klimko, L.A.; Nelson, P.I. On Conditional Least Squares Estimation for Stochastic Processes. Ann. Stat. 1978, 6, 629–642. [Google Scholar] [CrossRef]
  38. Cleveland, W.S. Inverse Autocorrelations of a Time Series and Their Applications. Technometrics 1972, 14, 277–293. [Google Scholar] [CrossRef]
  39. Kobayashi, K.; Salam, M.U. Comparing Simulated and Measured Values Using Mean Squared Deviation and its Components. Agron. J. 2000, 92, 345–352. [Google Scholar] [CrossRef]
  40. Vinod, H.D. Generalization of the Durbin–Watson Statistic for Higher Order Autoregressive Process. Commun. Stat. 1973, 2, 115–144. [Google Scholar] [CrossRef]
  41. Cavanaugh, J.E.; Neath, A.A. The Akaike information criterion: Background, derivation, properties, application, interpretation, and refinements. WIREs Comput. Stat. 2019, 11, e1460. [Google Scholar] [CrossRef]
  42. Zhu, L.; Li, L.; Liang, Z. Comparison of six statistical approaches in the selection of appropriate fish growth models. Chin. J. Ocean. Limnol. 2009, 27, 457. [Google Scholar] [CrossRef]
  43. Strimbu, B.M.; Petritan, I.C.; Montes, C.; Biris, I.A. An assessment of the O-ring methodology using virgin stands of mixed European beech—Sessile oak. For. Ecol. Manag. 2017, 384, 378–388. [Google Scholar] [CrossRef]
  44. Tayyebi, A.; Pijanowski, B.C.; Linderman, M.; Gratton, C. Comparing three global parametric and local non-parametric models to simulate land use change in diverse areas of the world. Environ. Model. Softw. 2014, 59, 202–221. [Google Scholar] [CrossRef]
Figure 1. Ozone concentration in downtown Los Angeles from January 1955 to December 1972.
Figure 1. Ozone concentration in downtown Los Angeles from January 1955 to December 1972.
Forests 11 01311 g001
Figure 2. The 12-months differenced ozone series from January 1955 to December 1972.
Figure 2. The 12-months differenced ozone series from January 1955 to December 1972.
Forests 11 01311 g002
Figure 3. The forecast for year 1973 by three different algorithms: ML, ULS, and CLS.
Figure 3. The forecast for year 1973 by three different algorithms: ML, ULS, and CLS.
Forests 11 01311 g003
Figure 4. White noise residuals produced by the autoregressive model.
Figure 4. White noise residuals produced by the autoregressive model.
Forests 11 01311 g004
Figure 5. Autocorrelation, inverse autocorrelation, and partial autocorrelation function of the residuals produced by the autoregressive function.
Figure 5. Autocorrelation, inverse autocorrelation, and partial autocorrelation function of the residuals produced by the autoregressive function.
Forests 11 01311 g005
Table 1. Summary of statistics of the ozone concentration. The units are in parts per hundred million (pphm).
Table 1. Summary of statistics of the ozone concentration. The units are in parts per hundred million (pphm).
PeriodNumber of ValuesMeanStandard DeviationMinimumMaximum
1955–19722163.771.491.28.7
≥19601563.311.211.26.1
Winter ≥ 1966492.470.841.24.4
Summer ≥ 1966354.040.972.15.7
Table 2. Parameter estimates and fit statistics (mean square error (MSE), and Akaike Information Criterion (AIC)) of the ARIMA (autoregressive integrated moving average) model with three different algorithms of the ARIMA process. CLS—conditional least square; ULS—unconditional least square; ML—maximum likelihood.
Table 2. Parameter estimates and fit statistics (mean square error (MSE), and Akaike Information Criterion (AIC)) of the ARIMA (autoregressive integrated moving average) model with three different algorithms of the ARIMA process. CLS—conditional least square; ULS—unconditional least square; ML—maximum likelihood.
AlgorithmOzone lag (1)Ozone lag (12)X1 lag(12)SummerWinterMSEAIC
CLS−0.299830.59234−1.2624−0.2615−0.0820.856520.52
ULS−0.252641−1.3742−0.2027−0.0840.768510.78
ML−0.266840.76665−1.3306−0.2394−0.0800.797501.77
Table 3. Forecasting the year 1973 using three ARIMA estimation algorithms: CLS, ULS, and ML.
Table 3. Forecasting the year 1973 using three ARIMA estimation algorithms: CLS, ULS, and ML.
1973CLSULSML
MonthsForecastSEForecastSEForecastSE
January1.45380.85621.27380.76771.42050.7966
February1.87830.89381.58710.79191.84460.8244
March2.57050.89382.24820.79192.45670.8244
April2.87880.89382.76490.79192.85900.8244
May3.15530.89382.94820.79193.15010.8244
June2.80740.89382.79520.79192.72110.8244
July3.33550.89383.33410.79193.31470.8244
August3.44000.89383.35070.79193.47870.8244
September2.82530.89383.29520.79192.04050.8244
October2.07410.89382.75070.79192.35870.8244
November1.64650.89382.39820.79191.85880.8244
December1.21490.89381.56630.79191.28980.8244
Table 4. Analysis of the errors produced by the autoregressive model with the four algorithms.
Table 4. Analysis of the errors produced by the autoregressive model with the four algorithms.
AlgorithmMSEAICR2Durbin Watson at Lag 12
Pr < DWPr > DW
ML0.76425527.8620.51530.46470.5353
ULS0.76311528.2020.5160.32870.6713
ITYW0.77253529.0440.510.74130.2587
YW0.77311529.1530.50970.75260.2474
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paun, M.; Gunaime, N.; Strimbu, B.M. Impact of Algorithm Selection on Modeling Ozone Pollution: A Perspective on Box and Tiao (1975). Forests 2020, 11, 1311. https://doi.org/10.3390/f11121311

AMA Style

Paun M, Gunaime N, Strimbu BM. Impact of Algorithm Selection on Modeling Ozone Pollution: A Perspective on Box and Tiao (1975). Forests. 2020; 11(12):1311. https://doi.org/10.3390/f11121311

Chicago/Turabian Style

Paun, Mihaela, Nevine Gunaime, and Bogdan M. Strimbu. 2020. "Impact of Algorithm Selection on Modeling Ozone Pollution: A Perspective on Box and Tiao (1975)" Forests 11, no. 12: 1311. https://doi.org/10.3390/f11121311

APA Style

Paun, M., Gunaime, N., & Strimbu, B. M. (2020). Impact of Algorithm Selection on Modeling Ozone Pollution: A Perspective on Box and Tiao (1975). Forests, 11(12), 1311. https://doi.org/10.3390/f11121311

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop