** Atmospheric and Climate Sciences** Vol.4 No.2(2014), Article ID:45310,11 pages DOI:10.4236/acs.2014.42033

UV Index Modeling by Autoregressive Distributed Lag (ADL Model)

Alexandre Boleira Lopo^{1,2}, Maria Helena Constantino Spyrides^{1,3}, Paulo Sérgio Lucio^{1,3}, Javier Sigró^{4}

^{1}Doctoral Program in Climate Science (PPGCC/UFRN), Natal, Brazil

^{2}State University of Bahia (UNEB), Salvador, Brazil

^{3}Department of Statistics, Federal University of Rio Grande do Norte (UFRN), Natal, Brazil

^{4}Department of Geography, Centre for Climate Change (C3), University Rovira i Virgili, Tarragona, Spain

Email: alexandrelopo@hotmail.com, spyrides@ccet.ufrn.br

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 14 January 2014; revised 10 February 2014; accepted 17 February 2014

ABSTRACT

The objective of this work is to model statistically the ultraviolet radiation index (UV Index) to make forecast (extrapolate) and analyze trends. The task is relevant, due to increased UV flux and high rate of cases non-melanoma skin cancer in northeast of Brazil. The methodology utilized an Autoregressive Distributed Lag model (ADL) or Dynamic Linear Regression model. The monthly data of UV index were measured in east coast of the Brazilian Northeast (City of Natal-Rio Grande do Norte). The Total Ozone is single explanatory variable to model and was obtained from the TOMS and OMI/AURA instruments. The Predictive Mean Matching (PMM) method was used to complete the missing data of UV Index. The results mean squared error (MSE) between the observed UV index and interpolated data by model was of 0.36 and for extrapolation was of 0.30 with correlations of 0.90 and 0.91 respectively. The forecast/extrapolation performed by model for a climatological period (2012-2042) indicated a trend of increased UV (Seasonal Man-Kendall test scored τ = 0.955 and p-value < 0.001) if the Total Ozone remain on this tendency to reduce. In those circumstances, the model indicated an increase of almost one unit of UV index to year 2042.

**Keywords:**UV Flux; Dynamic Linear Regression Model; Seasonal Man-Kendall Test; Mean Squared Error; Residuals

1. Introduction

The modeling of Ultraviolet Radiation (UV) and its index (UV Index or UVI) is important, among other reasons; due to lack of information on these variables and small number of stations [1] -[3] . This fact leads researchers to invest in development of computational, physical, statistical and stochastic models to estimate or forecast/projections of UV/UVI [4] [5] .

In Brazil, some statistical models for UV-UV Index were validated and performed by researchers [6] built statistical models and artificial neural networks to estimate UV and obtained a mean square error of less than 5%. Researchers of Center for Weather Forecasting and Climate Studies, National Institute for Space Research (INPE/CPTEC) [7] conducted a study on the estimates of the radiative transfer model UVSIM (Ultra-Violet SImplified Model) to analyze the effect of cloud in UV Index. The results showed a high correlation model (0.8) with observed values. In this model, the input parameters are the coordinates of UVSIM time (hour, day), geographic (latitude, longitude), cloud cover and Ozone concentration estimated by instrument OMI (Ozone Monitoring Instrument) of satellite AURA/NASA. Finally, the model UVGAME (UV Global Atmospheric Model) was used in research on the UV Index variations and regional and seasonal distribution of the number of cases of skin cancer in the Brazilian’s skin color [8] .

Internationally, the historical anthropogenic changes in the surface all-sky UV-B radiation through 1850-2005 are evaluated by analyzing the CMIP5 transient historical simulations performed with MI-ROC-ESM-CHEM were the researchers studied [5] . The results of the study have indicated that changes in ozone transport in the lower stratosphere, which is induced by increasing greenhouse gas concentrations, increase ozone concentration in the extra tropical upper troposphere and lower stratosphere. These transient changes work to decrease the amount of UV-B reaching the Earth’s surface, counteracting the well-known effect increasing UV-B due to stratospheric ozone depletion, which developed rapidly after 1980.

The research [9] conducted an analysis by 13 models used in prediction schemes for UV Index, including simple regression. The models had parameters like location data, Total Ozone column and solar zenith angle. The authors considered that the differences between the models were derived from treatment of unknown input parameters, especially aerosols. Another research [4] used data from the period 1976-2006 and a regression model to establish a statistical relationship between UV and content ozone, global solar radiation and cloud cover at different scales of time and explanatory variables in order to reduce the standard error. Others Researchers [10] established a statistical relationship between Total Ozone, solar zenith angle and UV irradiance. Modeling of UV and problems on the radiative transfer in clear and turbid environments were described by [11] . The researchers of European Communities [12] emphasizes that prediction methods vary in simple statistical methods used to certain times and local until robust methods to forecast some hours to several days, either for all sky or clear conditions sky. The accuracy of forecasts of UVI is mainly limited by the quantity and quality of input data. Finally, [12] explains that in future the data assimilation of large-scale ground-based observations of Total Ozone, aerosol and cloud cover through satellite should considerably improve accuracy for models.

This study presents the proposal to apply an autoregressive distributed lag model (ADL) or Dynamic Linear Regression model. This is a linear regression model involving time series that includes current and past values of the variable under study and explanatory variables with or without time lags [13] [14] . It is used when there is a dependency structure between these variables. These models have been used in the past in the environmental field [15] and dynamic regression models are widely used because they express and model the behaviour of a system over time [16] . As the variables of the model are indexed by time and since there are lags in both exogenous and dependent variables, then the ADL model should be used [15] .

The ADL model can be well utilized to make forecasts/projections (extrapolation). In this study, UV Index is response variable and Total Ozone is explanatory variable. It is worth noting that the more current information is of greater importance to improve forecasts [17] .

The modeling of UV/UV Index and its use for forecasts/projections is important due to increase in Flow of Ultraviolet Radiation (FUV) on the surface of the Earth. Researches [18] -[20] have warned of an increase in biologically active ultraviolet radiation (erythemal dose) due to decreasing levels of stratospheric ozone.

Ozone is the component responsible in absorbing the UV in earth’s atmosphere, preventing it to come fully to the surface of the Earth. The reduction in the quantity (Total Ozone column or Total Ozone) is the main cause for the increase of UV and brings impacts on nature and human health [8] [20] [21] . The Global Total Ozone average from 2005 to 2009 is around 3.5% below that recorded in the period 1964-1980. The changes in this variable occurred since 1970 (base year) to 2010, in which there is a reduction of 3% [20] [21] .

The researchers [22] analyzed the effects of climate change on the ozone layer using climate models. The models indicated an accelerated circulation in the stratosphere with changes in the spatial distribution of tropospheric and stratospheric ozone. The authors showed that, in future scenarios of the IPCC (until 2095), the UV Index would change with clear skies, with a 9% reduction in high northern latitudes and increases of 4% in the tropics and up to 20% at high latitudes of the south, during late spring and early summer. The results suggest that climate change will alter the balance of tropospheric ozone and UV Index, which would have consequences for radiative forcing in the troposphere, air quality and human health and ecosystems, however the amount of ozone in the atmosphere is recuperating as the result of the Montreal Protocol, in force since 1989 [21] -[23] .

In Northeast Brazil (NEB) there is an increase of Flow Ultraviolet Radiation (FUV) [16] . This fact has motivated this study, whose goal is to build a forecasting model of UV Index in city of Natal, capital of state Rio Grande do Norte, located in the NEB.

This article aims to conduct modeling, using ADL model, for variability of UV Index in city of Natal, east coast of the Brazilian Northeast, a function of Total Ozone and perform predictions (interpolation), forecasts/ projections (extrapolation) and trend analysis, collaborating to understand the issues presented.

2. Material and Methods

2.1. Study Area

Natal is a tourist town with beautiful beaches and 853,928 inhabitants [24] located in east coast of NEB between the sea (South Atlantic) and the right bank of the Potengi River near the equator (5˚45'54''S and 35˚12'05''W). The city is called by its inhabitants of “Sun City” [25] due to the great sunlight related to the intensity of solar radiation and total annual insolation of 2968.4 hours [26] . The “Sun City” has problems of non-melanoma skin cancer rates above the regional average for women (54%) and men (87%) [27] .

In this city the UV Index is classified as “extreme” from October to April and “very high” from May to August. The annual variability has a characteristic in September and October that consists of stabilization and reduction in UV Index because of highest concentration of Total Ozone and increased presence of marine aerosols [28] . The annual UV Index average observed for the period was 11 (±1.0). In the daily cycle, the maximum UV Index occurs around 11:20 am, classified as “high” starting from 9:00 am local time (time zone GMT-03) to 10:00 am achieves intensity considered “too high” [29] .

2.2. UV Index-Formulation

The UV Index describes the intensity of UV in relation to its photo-biological effect [1] being defined by equation (1):

(1)

in which E_{λ} is the spectral irradiance expressed in W×m^{−2}×nm^{−1} to the wavelength λ and dλ it is the wavelength range used in the integral calculus. S_{er} is the reference action spectrum erythema and K_{er} is a constant equal to 40 m^{2}×W^{−1}.

The UV Index was proposed by the World Health Organization (WHO) [1] based on the reference spectrum of action erythema of [30] .

2.3. Data for UV Index and Total Ozone Column

The daily data (2001-2012) of UV Index were measured in the surface by: 1) Radiometer GUV (Ground-based Ultraviolet Radiometer) [31] installed in Laboratory of Tropical Environmental Variables of National Institute for Space Research/Center Regional of Northeast (INPE/CRN/LAVAT) and 2) Sensor Model UV-6490 of Meteorological Station Davis installed at the Laboratory of Hydraulic Machines and Solar Energy in Technology Center of University Federal of Rio Grande of Norte (UFRN/LMHES). The daily maximum of UV Index was collected on interval of 11 h - 13 h, independent of sky conditions A set of daily data (2001-2012) of Total Ozone (DU units) was collected at TOMS (Total Ozone Mapping Spectrometer) e OMI (Ozone Monitoring Instrument), available in http://avdc.gsfc.nasa.gov/. The spectrometer TOMS is atmospheric sensor has been flying on different missions within NASA’s Earth Probes Program. The objective is to extend the global ozone data set that began in 1978 with the flight of TOMS on NIMBUS-7. The end of the operation occurred in 2005 when he worked on the platform Terra Probe [32] . The instrument OMI was launched in July 2004 on board EOS-Aura. OMI monitors the recovery of the ozone layer in response to the phase out of chemicals, such as CFCs. Together with its companion instruments MLS and HIRDLS it will measure criteria pollutants such as O_{3}, NO_{2}, SO_{2} and aerosols [33] .

2.4. Completing the Missing Data

The presence of missing data on UV Index series and the need to complete time series for stochastic models was applied multiple imputation technique for each group of same month, using Predictive Mean Matching method (PMM) [34] [35] placed in MICE package (Multivariate imputation by Chained Equations) in the R software free, available at http://www.r-project.org/ [36] . The software MICE used allowed program their own imputation function, while at the same time it supports a variety of imputation methods [37] .

The PMM method is an imputation method that combines parametric and nonparametric techniques. It imputes missing values by means of the Nearest Neighbour Donor where the distance is computed on the expected values of the missing variables conditional on the observed covariates, instead of directly on the values of the covariates [38] . The PMM is a variant of linear regression that determines an imputed value calculated by the regression model closest to the observed value [39] [40] . The PMM consider the following formulation (equation (2)) for each i missing in Y

(2)

Consider X a variable without missing data, the group of observed values; and , one obtains Ŷ^{obs} found as the nearest observation .

There are a number of imputation methods: simple imputation methods, regression imputation, hot deck imputation methods and distance function matching. PPM method is considered to be the most accurate, since it combines elements of regression, nearest-neighbor and hot deck imputation [41] and is characterized as a general purpose method [42] .

PPM can overcome the difficulties of both parametric and nonparametric imputation techniques, given the fact that parametric techniques may fail when the model is not suitable for the available data and nonparametric techniques require high amount of observations [38] .

The PMM application results were promising for monthly scale when compared with the original data [43] [44] , although performance of the predictive mean method varies considerably with the predictive power of the imputation regression model and the percentage of cases with missing data on income [45] .

The procedures for filling the missing values using the method of PMM must follow the criteria: data are monthly and the regression model was applied to each group of same months (same name) of full data series [39] [40] , because the use of locally adjusted PMM method provides reduced bias [46] .

2.5. Autoregressive Distributed Lag Model (ADL)

The ADL is a parametric model that combines the dynamics of time series and the effect of explanatory variables. It consists of stochastic regression involving time series that includes current and past values of the variable under study and explanatory variables, including lags [13] [14] [47] . This model uses the notation ADL (p_{j}, p), wherein p_{j} and p indicate the lag order of the variables or variable, respectively, dependent and independent [14] . General representation of the model is defined by equation (3):

(3)

in which: y_{t}: dependent variable in time t; β_{0}: a constant; y_{t‒i} the dependent variable in t ‒ i; x_{jt‒i} is the j^{th} independent variable in t ‒ i, by i = {1, ×××, p_{j}} e j = {1, ×××, p}; β_{ji}: coefficient of the j^{th} independent variable in t ‒ i; f_{i}: coefficient of the dependent variable in t ‒ i; ε_{t} random residual [13] [14] [48] .

The method of Ordinary Least Squares was used to estimate the parameters in the model: β_{0} representing the intercept and β_{1},_{ }β_{2}, …, β_{n} are the angular coefficients [49] [50] . In the model proposed in this study, y_{t} refers to UV Index at time t and x_{t} was considered the Total Ozone at time t.

In applications of ADL model, the residual error must meet the following assumptions: the errors ε_{t} are random and independent variables following a normal distribution: ε_{t} ~ N (0, σ^{2}) with zero mean, and constant variance (σ^{2}) (homoscedasticity) [51] . The Kolmogorov-Smirnov and Shapiro-Wilk tests were used to verify the normality assumption [52] . The homoscedasticity assumption was measured using the Breusch-Pagan test [53] [54] .

The F test was used to test the significance of the regression equation and t-test to measure contribution of explanatory variable [55] [56] . In this study was applied 5% significance for the tests.

There are a number of distributed lag models and the selection of the appropriate one consists in specifying the lag length correctly [16] . The ADL model includes lagged values of both independent and dependent variables and was chosen due to the fact that: 1) it is not so strict, as geometric lag models or finite distributed lag models and 2) it is a general form that can capture the current and lagged effects of an independent variable over the dependent [16] .

In this study the Cross-Correlation Function (CCF) was used to determine the number of lags between UV Index and Total Ozone for considered on the ADL model [57] [58] . The autocorrelation function (ACF) was applied to the original series of variability of UV Index in order to verify the seasonality and non-stationarity this series [58] .

Box-Pierce test [59] and Durbin-Watson test [60] [61] were applied to assess the independence of errors.

The calculations and results were obtained through the procedures performed by the R software free [62] [63] and package dynlm (Dynamic Linear Models and Time Series Regression) authored by [64] .

2.6. Trend Analysis

The Mann-Kendall test nonparametric [65] was applied with 95% statistical significance to analyze the trend of the number and intensity of UV Index. This test compares each value of the temporal series with the remaining values in sequential order, counting the number of times that the remaining terms are greater than the analyzed value [66] . This test is the most appropriate method to analyze weather tendencies in climatological series and has been used to calculate climatic tendencies [67] -[69] . In applying the Mann-Kendall test, we used the package SeasonalMannKendall (Mann-Kendall trend test for monthly environmental time series) in library of software R [70] [71] .

3. Results and Discussions

3.1. Descriptive Study

The variability and monthly and annual average of UV Index and Total Ozone are presented in Figure 1 and Table 1. The average annual of UV Index is 11 and for Total Ozone are 264.8 UD. The colors in the graph of UV Index are associated with the categorization of risk of the WHO. The UV index reaches the value classified as “extreme” (color violet) in seven months of the year, between spring and summer.

The annual variability of the UV Index in city of Natal has a stabilization/reduction in September e October as associated with a higher ozone concentration [28] . This feature was observed in this data series, as shown in Figure 1.

3.2. Completing the Data Daily of UV Index.

There missing data in 13 months of the total of 144, corresponding to 9.03% of the total number of observations. The process of multiple imputations through PMM was used to completing the missing data. The Figure 2 shows the time series of UV Index with data filled by method PMM and Total Ozone (simultaneously).

In the series very high UV Index values were observed in three periods, February 2005 and 2007 with 13.7 (252.6 DU) and 13.5 (250.5 DU) respectively and March 2010 with 13.6 (248.4 DU). The high values were coincident with low levels of Total Ozone. This study did not assess the causes of the extreme values of UV Index, however in Figure 2 identified that the low values of Total Ozone was associated with these measures.

3.3. Modeling in UV Index

The ADL model final was adjusted with historical values (period 2001-2011) of UV Index considering a lag 1, 4, 7 and 12 to capture seasonality. The signal of Total Ozone (explanatory variable) was used without and with lag 2. This ADL (4, 1) is represented by coefficients in Table 2.

The significance of regression equation was confirmed through of F test (F = 80.61 and p-value < 0.001).

A regression model is accepted if the residuals are standardized and homoscedastic. The Figure 3 shows graph

Figure 1. Boxplot by variability of the monthly UV Index and Total Ozone to Natal-RN for period 2001-2012. Source: INPE/CRN-LMHES/UFRN and TOMS-OMI/AURA.

Table 1. UV Index and Total Ozone in the city of Natal (monthly and annual average for 2001-2012).

Source: INPE/CRN/LAVAT-LMHES/UFRN and TOMS-OMI/AURA.

Figure 2. Time series of UV Index X Ozone total (vertical lines in green) for period 2001-2012 in city of Natal-RN.

Table 2. Coefficient regression (Coeff), standard error (SE) and p-value for ADL model.

of residuals (on top) and graphs for analysis of variance and normality. The normality of residuals was confirmed by Kolmogorov-Smirnov test (KS) (p-value = 0.76) and by Shapiro-Wilk test (SW) (p-value = 0.40).

The Breusch-Pagan test (p-value = 0.86) confirmed the hypothesis of homoscedasticity of residuals. The Box-Pierce test indicated independence of residuals (p-value = 0.07) together with Durbin-Watson test (p-value < 0.001).

Figure 3. Graphs for diagnostic of residuals of model ADL (4, 1).

3.4. Validation of Model ADL

The Figure 4(a) shows the UV Index observations and predicted values (interpolated) for the years 2002‒2011 by the model ADL.

The ADL model presented good results for the interpolation with a mean squared error (MSE) of 0.36 and correlation of 0.90 between the interpolated data and observed of UV Index and MSE of 0.30 and correlation of 0.91 for extrapolation (forecast to 2012) in Figure 4(b). These results were considered appropriate and validated the model to make forecasts.

3.5. Trend Analysis of UV Index and Total Ozone.

The trend study for UV Index and Total Ozone in observed period occurred through the Seasonal Man-Kendall test. Was identified stability in UV levels (τ = −0.009 and p-value = 0.896) and a downward trend in the levels of Total Ozone (τ = −0.302 and p-value < 0.001). The Figure 5 shows the time series with lowess smooth.

The projection performed by model for a period climatological (2012-2042) indicated a trend of increased of

(a) (b)

Figure 4. (a) UV Index data observed (black) and prediction or interpolation (blue) for the ADL model for the city of Natal for period 2001-2011; (b) Observations of UV Index (red), interpolation data by model (black), forecast values for 2012 (blue) and confidence interval (orange).

Figure 5. The time series plot with lowess smooth for UV Index and Total Ozone for the period observed.

UV (τ = 0.955 and p-value < 0.001) whereas Total Ozone remain on this trend. The current average annual of UV Index is 11.0, however the model predicts a rise in 2042 for average of 11.8, an increase of almost one unit of the UV Index. At present the amount of ozone in the atmosphere is recuperating [2] [5] and must change this trend.

4. Final Considerations

In constructing the model ADL (4, 1) have been found MSE and correlation index between UV index data (interpolated and extrapolated X observed) considered appropriate to validate the model. Furthermore, in model we confirmed that residues were random variables with zero mean, independents, normalized and with homoscedasticity, fact that can provide errors in estimates of regression coefficients and the general applicability of the model if not confirmed.

In the study, the UV Index model was built with a single explanatory variable, the Total Ozone along with the autoregressive signal of UV index. This is characteristic of ADL model which can combine the autoregressive signal and explanatory variables with the possibility of introducing lags.

The forecast/extrapolation performed by model for a climatological period (2012-2042) indicated a trend of increased UV (Seasonal Man-Kendall test scored τ = 0.955 and p-value < 0.001) if the Total Ozone remain on this tendency to reduce. In those circumstances, the model indicated an increase of almost one unit of UV index to year 2042.

Finally, it is noteworthy that there is a scenario of increased the UV flux at Earth’s surface, being relevant the development of forecast models to UV index which can collaborate on preventive actions against high rate of cases skin cancer in northeast of Brazil.

Acknowledgements

This work was supported by the CNPq, National Council for Scientific and Technological Development, Brazil: via the Doctor Sandwich—SWE/CsF (246611/2012-0) and projects “ATMANTAR”-MCTI/PROANTAR/CNPq (520182/2006-5) and “A Atmosfera Antártica, Conexões e Impactos Ambientais na América do Sul”, INCT-APA, MCTI/PROANTAR/CNPq (574018/2008-5).

The author thanks the State University of Bahia by financial support during the PhD of PPGCC/UFRN and to Federal Institute of Education, Science and Technology of Bahia by study leave.

References

- World Health Organization and International Commission on Non-Ionizing Radiation Protection (2002) Global Solar UV Index: A Practical Guide.
- WMO, World Meteorological Organization (1998) Report of the WMO-WHO Meeting of Experts on Standardization of UV Indices and Their Dissemination to the Public.
- Lee-Taylor, J. and Madronich, S. (2007) Climatology of UV-A, UV-B, and Erythemal Radiation at the Earth’s Surface, 1979-2000. National Center for Atmospheric Research, Boulder.
- Borkowski, J.L. (2008) Modelling of UV Radiation Variations at Different Time Scales. Annales Geophysicae, 26, 441-446. http://dx.doi.org/10.5194/angeo-26-441-2008
- Watanabe, S., Takemura, T., Sudo, K., Yokohata, T. and Kawase, H. (2012) Anthropogenic Changes in the Surface All-Sky UV-B Radiation through 1850-2005 Simulated by an Earth System Model. Atmospheric Chemistry and Physics, 12, 5249-5257. http://dx.doi.org/10.5194/acp-12-5249-2012
- Leal, S.S., Tíba, C. and Piacentini, R. (2011) Daily UV Radiation Modeling with the Usage of Statistical Correlations and Artificial Neural Networks. Renewable Energy, 36, 3337-3344. http://dx.doi.org/10.1016/j.renene.2011.05.007
- Costa, S.S., Rodrigues, M., Ceballos, J. and Corrêa, M.P. (2010) Índices de Radiação Ultravioleta: Estudo Comparativo entre Modelo de Transferência Radiativa e Observações à Superfície. Proceedings of the 16th Brazilian Congress of Meteorology, Belém, 13-17 September 2010.
- Corrêa, M.P. and Plana-Fattori, A. (2006) Uma análise das variações do índice ultravioleta em relação às observações de conteúdo de ozônio e da espessura óptica dos aerossóis sobre a cidade de São Paulo. Revista Brasileira de Meteorologia, 21, 24-32.
- De Backer, H., Koepke, P., Bais, A., de Cabo, X., Frei, T., Gillotay, D. and Vanicek, K. (2001) Comparison of Measured and Modelled UV Indices for the Assessment of Health Risks. Meteorological Applications, 8, 267-277. http://dx.doi.org/10.1017/S1350482701003024
- Fioletov, V.E., Kerr, J. and Wardle, D.I. (1997) The Relationship between Total Ozone Column and Spectral UV Irradiance from Brewer Observations and Its Use for Derivation of Total Ozone Column from UV Measurements. Geophysical Research Letters, 24, 2997-3000. http://dx.doi.org/10.1029/97GL53153
- Zerefos, C. (1997) Factors Influencing the Transmission of Solar Ultraviolet Irradiance through the Earth’s Atmosphere. In: Solar Ultraviolet Radiation, Springer, Berlin, 133-141. http://dx.doi.org/10.1007/978-3-662-03375-3_9
- Vanicek, K., Frei, T., Litynska, Z. and Schmalwieser, A. (2000) UV-Index for the Public. Publication of the European Communities, Brussels.
- Pankratz, A. (2012) Forecasting with Dynamic Regression Models. Vol. 935, John Wiley & Sons, Hoboken.
- Kiviet, J.F. and Dufour, J.M. (1997) Exact Tests in Single Equation Autoregressive Distributed Lag Models. Journal of Econometrics, 80, 325-353. http://dx.doi.org/10.1016/S0304-4076(97)00048-1
- Ravines, R.R., Schmidt, A.M. and Migon, H.S. (2006) Revisiting Distributed Lag Models through a Bayesian Perspective. Applied Stochastic Models in Business and Industry, 22, 193-210. http://dx.doi.org/10.1002/asmb.628
- Diniz, C.A.R., Rodrigues, C.P., Leite, J.G. and Pires, R.M. (2014) A Bayesian Estimation of Lag Lengths in Distributed Lag Models. Journal of Statistical Computation and Simulation, 84, 415-427. http://dx.doi.org/10.1080/00949655.2012.712972
- Buss, G. (2009) Economic Forecasts with Bayesian Autoregressive Distributed Lag Model: Choosing Optimal Prior in Economic Downturn. http://epp.eurostat.ec.europa.eu/cache/ITY_PUBLIC/EWP-2011-006/EN/EWP-2011-006-EN.PDF
- Madronich, S., McKenzie, R.L., Björn, L.O. and Caldwell, M.M. (1998) Changes in Biologically Active Ultraviolet Radiation Reaching the Earth’s Surface. Journal of Photochemistry and Photobiology B: Biology, 46, 5-19. http://dx.doi.org/10.1016/S1011-1344(98)00182-1
- Herman, J.R. (2010) Global Increase in UV Irradiance during the Past 30 Years (1979-2008) Estimated from Satellite Data. Journal of Geophysical Research: Atmospheres (1984-2012), 115, Article ID: D04203. http://dx.doi.org/10.1029/2009JD012219
- Herman, J.R., Labow, G., Hsu, N.C. and Larko, D. (2009) Changes in Cloud and Aerosol Cover (1980-2006) from Reflectivity Time Series Using SeaWiFS, N7-TOMS, EP-TOMS, SBUV-2, and OMI Radiance Data. Journal of Geophysical Research: Atmospheres (1984-2012), 114. http://dx.doi.org/10.1029/2007JD009508
- Fahey, D.W. and Hegglin, M.I. (2011) Twenty Questions and Answers about the Ozone Layer: 2010 Update. Scientific Assessment of Ozone Depletion, World Meteorological Organization, Geneva, 72 p.
- Hegglin, M.I. and Shepherd, T.G. (2009) Large Climate-Induced Changes in Ultraviolet Index and Stratosphere-toTroposphere Ozone Flux. Nature Geoscience, 2, 687-691. http://dx.doi.org/10.1038/ngeo604
- World Meteorological Organization (2011) Scientific Assessment of Ozone Depletion: 2010. Global Ozone Research and Monitoring Project-Report No. 52, 516 p, World Meteorological Organization, Geneva.
- Instituto Brasileiro de Geografia e Estatística (2013) BRASIL, IBGE, Estimativa Populacional 2013, Instituto Brasileiro de Geografia e Estatística, 2013. http://www.ibge.gov.br/
- Furtado, E.M. (2005) A onda do turismo na cidade do sol: A reconfiguração urbana de Natal. Dissertação, UFRN, Natal, Brasil.
- Instituto Nacional de Meteorologia-INMET (2009) Normais Climatológicas do Brasil 1961-1990. Brasília.
- Instituto Nacional de Câncer-INCA (2011) Estimativa 2012-Incidência de Câncer no Brasil, Ministério da Saúde do Brasil.
- Lopo, A.B., Spyrides, M.H.C., Lucio, P.S. and Sigró, J. (2013) Ozone and Aerosol Influence on Ultraviolet Radiation on the East Coast of the Brazilian Northeast. Atmospheric and Climate Sciences, 4, 92-99. http://dx.doi.org/10.4236/acs.2014.41012
- Da Silva, F.R. (2008) Estudo da radiacão ultravioleta na cidade de Natal-RN. Doctoral Dissertation, Universidade Federal do Rio Grande do Norte.
- McKinlay, A.F. and Diffey, B.L. (1987) A Reference Action Spectrum for Ultraviolet Induced Erythema in Human Skin. CIE Journal, 6, 17-22.
- Biospherical Instruments (2013) Application Note: GUV Diffey Dose Calculations. http://www.biospherical.com/
- Bramstedt, K., Gleason, J., Loyola, D., Thomas, W., Bracher, A., Weber, M. and Burrows, J.P. (2003) Comparison of Total Ozone from the Satellite Instruments GOME and TOMS with Measurements from the Dobson network 1996- 2000. Atmospheric Chemistry and Physics, 3, 1409-1419. http://dx.doi.org/10.5194/acp-3-1409-2003
- Team O.M.I. (2013) Ozone Monitoring Instrument (OMI) Data User’s Guide. http://avdc.gsfc.nasa.gov/
- Rubin, D.B. (2009) Multiple Imputation for Nonresponse in Surveys. Vol. 307, John Wiley & Sons, Hoboken.
- Rubin, D.B. and Little, R.J. (2002) Statistical Analysis with Missing Data. John Wiley & Sons, Hoboken.
- van Buuren, S. and Groothuis-Oudshoorn, K. (2011) Mice: Multivariate Imputation by Chained Equations in R. Journal of Statistical Software, 45.
- Horton, N.J. and Kleinman, K.P. (2007) Much Ado about Nothing. The American Statistician, 61, 79-90. http://dx.doi.org/10.1198/000313007X172556
- Di Zio, M. and Guarnera, U. (2009) Semiparametric Predictive Mean Matching. AStA Advances in Statistical Analysis, 93, 175-186. http://dx.doi.org/10.1007/s10182-008-0081-2
- Li, K.H., Raghunathan, T.E. and Rubin, D.B. (1991) Large-Sample Significance Levels from Multiply Imputed Data Using Moment-Based Statistics and an F Reference Distribution. Journal of the American Statistical Association, 86, 1065-1073. http://dx.doi.org/10.1080/01621459.1991.10475152
- Schafer, J.L. (1997) Analysis of Incomplete Multivariate Data. CRC Press, Boca Raton. http://dx.doi.org/10.1201/9781439821862
- Durrant, G.B. (2005) Imputation Methods for Handling Item-Nonresponse in the Social Sciences: A Methodological Review. ESRC National Centre for Research Methods and Southampton Statistical Sciences Research Institute. NCRM Methods Review Papers NCRM/002.
- Vink, G., Frank, L.E., Pannekoek, J. and Buuren, S. (2014) Predictive Mean Matching Imputation of Semicontinuous Variables. Statistica Neerlandica, 68, 61-90. http://dx.doi.org/10.1111/stan.12023
- Horton, N.J. and Lipsitz, S.R. (2001) Multiple Imputation in Practice: Comparison of Software Packages for Regression Models with Missing Variables. The American Statistician, 55, 244-254. http://dx.doi.org/10.1198/000313001317098266
- Kenward, M.G. and Carpenter, J. (2007) Multiple Imputation: Current Perspectives. Statistical Methods in Medical Research, 16, 199-218. http://dx.doi.org/10.1177/0962280206075304
- Landerman, L.R., Land, K.C. and Pieper, C.F. (1997) An Empirical Evaluation of the Predictive Mean Matching Method for Imputing Missing Values. Sociological Methods & Research, 26, 3-33. http://dx.doi.org/10.1177/0049124197026001001
- Okamoto, M. (2006) Bias-Reduced Multivariate Computation: Use of the Locally-Adjusted Predictive Mean Matching Method. http://www.amstat.org/sections/srms/proceedings/y2006/Files/JSM2006-000044.pdf
- Makridakis, S., Wheelwright, S.C. and Hyndman, R.J. (1998) Forecasting: Methods and Applications. 3rd Edition, Wiley, New York.
- Song, H., Witt, S.F. and Jensen, T.C. (2003) Tourism Forecasting: Accuracy of Alternative Econometric Models. International Journal of Forecasting, 19, 123-141. http://dx.doi.org/10.1016/S0169-2070(01)00134-0
- Faraway, J.J. (2004) Linear Models with R. CRC Press, Boca Raton.
- Seber, G.A. and Lee, A.J. (2012) Linear Regression Analysis. Vol. 936, John Wiley & Sons, Hoboken.
- Cook, R.D. and Weisberg, S. (1982) Residuals and Influence in Regression. Chapman and Hall.
- Verzani, J. (2004) Using R for Introductory Statistics. CRC Press, Boca Raton.
- Breusch, T.S. and Pagan, A.R. (1979) A Simple Test for Heteroscedasticity and Random Coefficient Variation. Econometrica: Journal of the Econometric Society, 47, 1287-1294. http://dx.doi.org/10.2307/1911963
- Conover, W.J., Johnson, M.E. and Johnson, M.M. (1981) A Comparative Study of Tests for Homogeneity of Variances, with Applications to the Outer Continental Shelf Bidding Data. Technometrics, 23, 351-361. http://dx.doi.org/10.1080/00401706.1981.10487680
- Neter, J., Kutner, M., Nachtssheim, C.J. and Wasserman, W. (1996) Applied Linear Statistical Models. 4th Edition, WCB/McGraw-Hill, Boston, 791 p.
- Kutner, M.H. (1996) Applied Linear Statistical Models. Vol. 4, Irwin, Chicago.
- Helleseth, T. (1976) Some Results about the Cross-Correlation Function between Two Maximal Linear Sequences. Discrete Mathematics, 16, 209-232. http://dx.doi.org/10.1016/0012-365X(76)90100-X
- Shumway, R.H., Stoffer, D.S. and Stoffer, D.S. (2000) Time Series Analysis and Its Applications. Vol. 3, Springer, New York. http://dx.doi.org/10.1007/978-1-4757-3261-0
- Box, G.E. and Pierce, D.A. (1970) Distribution of Residual Autocorrelations in Autoregressive-Integrated Moving Average Time Series Models. Journal of the American Statistical Association, 65, 1509-1526. http://dx.doi.org/10.1080/01621459.1970.10481180
- Durbin, J. and Watson, G.S. (1950) Testing for Serial Correlation in Least Squares Regression. I. Biometrika, 37, 409- 428. http://dx.doi.org/10.1093/biomet/37.3-4.409
- Durbin, J. and Watson, G.S. (1951) Testing for Serial Correlation in Least Squares Regression. II. Biometrika, 38, 159-177.
- R Development Core Team. R (2005) A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org
- Crawley, M.J. (2005) Statistics: An Introduction Using R. Wiley, New York, 342 p.
- Zeileis, A. (2009) Dynlm: Dynamic Linear Regression. R Package Version 0.2-3. http://CRAN.R-project.org/package= dynlm
- Mann, H.B. (1945) Non-Parametric Test against Trend. Econometrica, 13, 245-259. http://dx.doi.org/10.2307/1907187
- Kendall, M.G. (1975) Rank Correlation Methods. Charles Griffin, London, 120 p.
- Oliveira, P.T., Lima, K.C. and Silva, S. (2013) Synoptic Environment Associated with Heavy Rainfall Events on the Coastland of Northeast Brazil. Advances in Geosciences, 35, 73-78. http://dx.doi.org/10.5194/adgeo-35-73-2013
- dos Santos, D.N., de PR da Silva, V., de AS Sousa, F. and Silva, R.A. (2010) Estudo de alguns cenários climáticos para o Nordeste do Brasil. Revista Brasileira de Engenharia Agricola e Ambiental-Agriambi, 14. http://dx.doi.org/10.1590/S1415-43662010000500006
- Hirsch, R.M., Slack, J.R. and Smith, R.A. (1982) Techniques of Trend Analysis for Monthly Water Quality Data. Water Resources Research, 18, 107-121. http://dx.doi.org/10.1029/WR018i001p00107
- McLeod, A.I. (2005) Kendall Rank Correlation and Mann-Kendall Trend Test. R Package “Kendall”. http://www.stats.uwo.ca/faculty/aim
- Belle, G. and Hughes, J.P. (1984) Nonparametric Tests for Trend in Water Quality. Water Resources Research, 20, 127-136. http://dx.doi.org/10.1029/WR020i001p00127