Model for the substantiation of road safety improvement measures on the roads of Lithuania/Eismo saugumo gerinimo priemoniu pagrindimo modelis Lietuvos automobiliu keliams/Celu satiksmes drosibas uzlabosanas pasakumu ekonomiska pamatojuma modelis Lietuva/Liiklusohutuse parendamise meetmete pohjendusmudel Leedu teedele.
Ratkeviciute, Kornelija
State of the problem
After Lithuania joined the European Union (EU), due to the
expanding economic relations with EU member-states, the flows of transit
traffic on the roads of Lithuania and the vehicle ownership have been
increasing every year. Based on data of international accident
statistics Lithuania is one of those few countries where the number of
people killed per one million inhabitants is the highest compared to
other EU countries (Cygas et al. 2004). In 1994-2008 on the roads of
Lithuania 10 681 people were killed and 102 850 were injured. Most often
the victims are young people, road accidents cause large moral and
material losses for the public (Ratkeviciute et al. 2008).
In September 2001 the White Book, approved by the European
Commission, confirmed an ambitious objective to reduce the number of
accident victims by 50% by 2010 (from 50 000 to 25 000 per year). This
can be implemented by seeking to unify the amount of penalties, to
implement additional measures to ensure road traffic safety and to
introduce new modern technologies. Having joined the EU Lithuania also
undertook this objective: until 2010 to reduce the number of people
killed during road accidents by 50%. Unfortunately, statistical data
shows that safety situation on the roads of Lithuania is still one of
the worst in the European Union. The number of black spots is high,
therefore, it is necessary to develop an effective model for the
substantiation of road safety improvement measures (further--road safety
measures) which would allow to select road safety measures first of all
for those road sections where accident indices exceed the limit values.
The results of accident data research showed that the currently
used methodology for the substantiation of road safety measures in
Lithuania does not allow to accurately select road safety measures for
high-accident road sections, to forecast accidents and their preventive
measures (Ratkeviciute et al. 2008). Thus, it is necessary to correct it
taking into consideration the experience of Lithuania and other foreign
countries, to use mathematical statistics and mathematical optimisation models for preventing road accidents. Improvement of this methodology
describes the topicality of this dissertation.
Research object covers the injury accident concentration places,
black spots, methodology for the substantiation of road safety
improvement measures, mathematical statistics and optimisation models
for preventing road accidents. Experimental object of the work is
1995-2008 accident data on Lithuanian roads.
Aim and tasks of the work are to improve methodology for the
substantiation of road safety measures in order to eliminate black
spots, to construct a draft mathematical model for predicting fatal and
injury accidents and to form the optimization problem which enables to
model the optimal selection of road safety measures for the main roads
of Lithuania.
The following tasks must be solved to achieve the aim of the work:
--to make the analysis of the currently used Lithuanian methodology
for the substantiation of road safety measures;
--to make the analysis and evaluation of road safety indices after
implementation of special measures on the roads of national significance
of Lithuania and to determine the effect of methodology used;
--to study and evaluate similar methodologies for the
substantiation of road safety measures used in other foreign countries;
--having evaluated the experience of other countries and the
shortcomings of methodology currently used in Lithuania, to construct
the effective model for the substantiation of road safety measures;
--to carry out testing of a new computer program and experimental
calculations;
--to study the use of mathematical methods for predicting fatal and
injury accidents;
--to construct a draft mathematical model for predicting fatal and
injury accidents and to form the optimization problem which enables to
model the selection of road safety measures for the main roads of
Lithuania.
Scientific novelty of the work consists of the development of the
model for evaluating the effect of road safety measures and its
adaptation to the roads of national significance of Lithuania, based on
the existing traffic volume and traffic conditions.
For the first time the constructed evaluation model was based on
the results of analysis of engineering road safety measures implemented
on the high-accident road sections of Lithuania in 1997-2008, and on the
experience of foreign countries (Ratkeviciute et al. 2008). A new
computer program Evaluation of the Road Safety Measures was created, its
testing was carried out as well as experimental calculations.
For the first time in Lithuania a draft mathematical model is
presented for predicting road accidents and modelling the optimal
selection of road safety measures on the main roads of Lithuania.
The constructed model for the evaluation of road safety measures
and the developed computer program which is used for the economic
evaluation of road safety measures before their implementation on
high-accident road sections and on the black spots of Lithuania is a
practical value. The constructed draft mathematical stochastic model and
the optimization problem is used for accident forecasting, selection of
the optimal road safety measures to be implemented on the black spots of
the main roads of Lithuania.
[FIGURE 1 OMITTED]
Road safety situation in Lithuania
Though the number of vehicles and the average annual daily traffic
(AADT, vpd) has been increasing every year, the number of people killed
on the roads is decreasing (Fig. 1). Due to a continuous improvement of
vehicles (air bags, modern safety belts, anti-lock braking systems,
reinforced vehicle body, etc.) accident severity becomes more slight.
From 2001 to 2009 the number of people killed on the roads of
Lithuania was reduced by 26%. In 2008 the number of people killed during
road accidents in Lithuania was one of the least in the recent 40 years.
In the year 2008 on the roads of 27 EU member-states 39 thousand people
were killed, however, this number has decreased by 15.4 thousand since
2001. Such a rapid decrease in the number of accidents was achieved
merely by the united efforts of all the institutions concerned and their
joint actions in saving human lives. A perfect result was achieved by a
tightened road user control, gradual implementation of the engineering
road safety measures for the number of years, active educational
activities, therefore, the road users' culture and responsibility
on the road has notably changed, and the executed scientific works
allowed to more effectively work and adopt the best foreign practice
(Ratkeviciute et al. 2008).
In order to select the most suitable road safety measures the
causes of road accidents and their influencing factors were determined.
This is the most important task for road safety specialists which is
attempted to be strategically solved all over the EU in a way of
creating common methodologies. The basis of these
methodologies--identification of the accident-influencing factors and
implementation of appropriate road safety measures on the high-accident
road sections.
Analysis of methodologies for the substantiation of road safety
measures used in Lithuania and foreign countries
When comparing methodologies for the substantiation of road safety
measures used in Lithuania (TARVAL abbreviation from Finnish
Turvallisuusvaikutusten Arviointi Vaikutuskertoimilla TARVA Lithuanian
version), Finland (TARVA), Belarus and Poland the following main
differences could be distinguished:
--Belarusian and Polish methodologies are more complicated than
TARVAL and TARVA (Ratkeviciute et al. 2008);
--Polish methodology enables to predict the reduction of collisions
(%) with vehicles, pedestrians and bicyclists (Ratkeviciute et al.
2008);
--Belarusian methodology enables to predict the probability for the
reduction of accident number in parts of a unit from the total number of
accidents as well as the probability for the reduction of fatal or
injury accidents (Kapski et al. 2007; 2008; Ratkeviciute et al. 2008);
--TARVAL and TARVA versions enables to predict the impact
coefficients for the accidents with motor vehicles, pedestrians/cyclists
and animals;
--the impact coefficients for animal-involved accidents are used
only in TARVAL and TARVA methodologies (but not in Belarusian or Polish
methodologies).
The Finnish methodology for the substantiation of road safety
measures TARVA is very close to the Lithuanian version TARVAL, however,
at present it is much more new and differs from the initial TARVA
version (Cygas et al. 2004). The specialists of the Technical Research
Centre of Finland (VTT) have been continuously improving it and adapting
to the current traffic situation on the national roads (Peltola 2000;
2007). The latst version of TARVA program also allows to forecast road
accidents on the newly constructed or reconstructed road sections, to
select road safety measures and to carry out accident prevention
(Peltola 2007).
Sweden has no methodology for the substantiation of road safety
measures, similar to that of TARVAL used in Lithuania. However, the
Swedish methodology Traffic Conflict Technique (TCT) enables to identify
conflict situations between different road users on the road section or
junction before the occurrence of road accident. Having this type of
data it is possible to select those road safety measures which would
prevent from possible accidents (Jonsson 2005).
Research and evaluation of accident indices on the roads of
Lithuania
For the evaluation of the effect of road safety measures
implemented on the roads of Lithuania in 1999-2002 fatal and injury
accidents were analyzed of 1995-2006. Accident data is compared before
implementing road safety measures and 4 years after their implementation
on the selected 48 high-accident road sections.
Economic evaluation was carried out using the cost-benefit analysis
(Cygas et al. 2004; Hauer 2005). Based on this method benefits,
consisting of savings in accident losses, were compared to the costs
consisting of the implementation costs of road safety measures. Economic
evaluation of safety measures was conducted for each black spot.
Sensitivity analysis was carried out based on three sensitivity tests
(growth of traffic volume, costs of safety measures and change in
evaluation period). All 48 cases showed that investments into road
safety measures will pay back (Ratkeviciute et al. 2007; 2008).
Analysis of the effect of road safety measures implemented on 48
high-accident road sections showed that the expected results of the
improvement of safety situation were not achieved on 24 road sections,
i. e. 50%; the expected reduction in the accident rate was not achieved
on 7 road sections, i. e. 15%; the implemented road safety measures have
justified themselves on 24 road sections, i. e. 50% (Ratkeviciute et al.
2007).
In order to make an analysis of only partly justified road safety
measures the visual investigations of high-accident road sections were
carried out and a detail analysis of safety situation.
The use of research data in the models for the substantiation of
road safety measures
Lithuania has a comprehensive accident databases and makes an
accurate recording of accident locations. The amount and accuracy of
other data required for the accident analysis has highly increased:
traffic volume and its change, road safety measures, precise locations
for their implementation, etc. All this allow collecting reliable data
and based on this data to carry out a detail accident analysis, to
correct the list of road safety measures and their impact coefficients,
to more accurately forecast accident indices and to execute economic
evaluation.
Analysis of accident research data showed that the currently used
TARVAL software is not sufficiently accurate to evaluate safety measures
to be implemented and to forecast accident indices. Therefore, it was
necessary to correct it taking into consideration the experience in
Lithuania and other countries and the current safety situation on the
roads of Lithuania (Ratkeviciute et al. 2007).
Fig. 2 gives stages for the improvement of methodology used for the
substantiation of road safety measures in Lithuania. For the improvement
purposes the results of the analysis of effects of road safety measures
implemented in Lithuania in 1997-2008 were taken into consideration,
also the foreign experience and methodologies (Elvik 2009; Hakkert,
Gitelman 2004; Kapski et al. 2007; 2008; Peltola 2000; Ratkeviciute et
al. 2008; Yannisa et al. 2008), the road safety measures and their
effect on road safety. The list of road safety measures has been
supplemented with 59 new measures not earlier used.
Improvement of the method for the substantiation of road safety
measures under Lithuanian conditions
The computer program Saugaus eismo priemoniu vertinimas (Evaluation
of the Road Safety Measures) for the evaluation of road safety measures
consists of two parts:
--normative tables;
--data and calculations.
Before starting working with the program it is necessary to input
the newest initial data into the already made-up normative tables. The
cost of road accident is determined according to the unit prices
approved by the Director General of the Lithuanian Road Administration
under the Ministry of Transport and Communications of the Republic of
Lithuania. The latest unit prices were approved by the order No. V-410
of 19 November 2008. In order to make the analysis of measures
implemented in the previous years and the effect of those measures the
unit prices of that period must be used. Road safety measures data
stored in a table of the impact of road safety measures is the main
input data of the road safety management system. Impact coefficients
divided into three groups: vehicle-involved accidents,
pedestrian-involved accidents and animal-involved accidents. Road safety
measures are standard or individual. Individual road safety measures not
yet included into the system, but there is a possibility to enter new
safety measures.
[FIGURE 2 OMITTED]
Having entered data into normative tables the calculations are
carried out. Calculations consist of three stages:
--data input;
--forecast calculations of road accidents and accident losses;
--economic calculations.
Data input is in detail described in the User's Manual
presented in the dissertation. A very important feature of data input of
a new program is that the data input area is joined to the Lietuvos
valstybines reiksmes keliu informacine sistema (LAKIS) (The Lithuanian
State Road Information System of the Lithuanian Road Administration.
Economic calculations show if the project is expedient to be
implemented from the economic point of view. If several alternatives are
studied it is determined which of them is the most acceptable.
Summary table of results gives data on the road section where a
certain road safety measure is planned to be implemented, data on the
traffic volume and accident rate of this road section, the planned road
safety measures, summary data on the forecasted effect of measures after
their implementation. Based on economic indices and the forecasted
accident rate the software itself makes the evaluation of the
project's attractiveness which can be:
I--unsatisfactory;
II--satisfactory;
III--good;
IV--very good.
A test of the improved model for the substantiation of road safety
measures
Experimental calculations were carried out on those road sections
where after implementation of road safety measures or their complex the
expected decrease in the accident rate was not achieved. Based on the
results of experimental calculations the impact coefficients of road
safety measures were once again corrected.
Accident modelling
The developed computer program Evaluation of the Road Safety
Measures is used for the evaluation of safety measures to be implemented
on black spots. In order to make forecasts of the safety situation on
the newly built or reconstructed road sections and to prevent the
occurrence of black spots it is necessary to carry out a very
comprehensive analysis of historical data on the change of road
accidents and traffic volume on the roads of Lithuania, to divide the
roads of Lithuania into plenty of homogenous sections and to adapt the
methods of mathematical statistics and game theory.
Forecasts of accident number
In order to forecast the number of road accidents on the main roads
of Lithuania a statistical mathematical model of time series was used.
Since in the analysis of statistical data of the number of accidents in
a time series the clear seasonal variations were observed (Figs 3, 4) a
search for regression curve was made which would well reflect the
seasonal variations. A random quantity--the number of accidents in 3
months (in a quarter of the year)--is marked by Y:
[??] = [a.sub.0] + at + [a.sub.1][t.sub.1] + [a.sub.2][t.sub.2] +
[a.sub.3][t.sub.3], (1)
where [??]--forecasted average Y value; t--trend variable;
[t.sub.i]--variable taking the value 1 in the quarter i of the year and
the value 0 in other quarters (the fourth quarter of the year is
corresponded by the values of variables [t.sup.1] = [t.sub.2] =
[t.sub.3] = +).
[FIGURE 3 OMITTED]
Seasonal variations are described by the regression equations
(obtained using the Microsoft Excel tool Data
Analysis) (Table 1):
[??] = 190.19 +0.69t - 76.64[t.sub.1] - 81.63[t.sub.2] -
22.11[t.sub.3], (1997-2006) (2)
[??] = 187.32 + 0.60t - 72.92[t.sub.1] - 74.70[t.sub.2] -
16.49[t.sub.3], (1997-2007) (3)
[??] = 193.52+ 0.02t - 67.78[t.sub.1] - 69.54[t.sub.2,] -
14.48[t.sub.3]. (1997-2008) (4)
Fig. 3 gives the dependency of the observed (1997-2008) and
forecasted (based on 1997-2006 data) Y values in time. Accident data
studied in 1997-2007 showed that in the quarters III and IV the number
of recorded accidents was significantly larger compared to the quarters
I and II. Accident statistics of 2008 breaks even this regularity. The
observed dispersion of the random quantity Y in the quarters III and IV
is also larger compared to the quarters I and II.
[FIGURE 4 OMITTED]
The decrease in the number of accidents only in one year, i. e.
2008, is not statistically important to obtain a downward linear trend
or a certain non-linear trend. This requires further observations of the
accidents on the roads of Lithuania.
A detail analysis of the dependency of the number of accidents on
the length of road section and its traffic volume is given. A detail
statistical analysis of data about the black spots is presented too.
With the help of a hig-volume sample (n = 7095) the main numerical characteristics of people killed and injured in one road accident were
found. Forecasts were based on the fatal and injury accidents on the
main roads of Lithuania in 1997-2008.
When analyzing the number of accidents on the black spots of the
main roads in 2004-2007 data was grouped according to the type t of a
black spot:
--[t.sub.1]--group I--black spots having no junctions;
--[t.sub.2]--group II--black spots having at-grade junctions;
--[t.sub.3]--group III--black spots having two-level junctions.
The number of black spots of the group I was 26, of the group
II--70, of the group III--14. Such grouping enabled to rather reliably
forecast the average number of accidents of each group (marked with
[X.sub.1], [X.sub.2], [X.sub.3]), separately of vehicle-involved
accidents and accidents with pedestrians and cyclists (marked with
[X.sub.1m], [X.sub.2m], [X.sub.3m] and [X.sub.1dp], [X.sub.2dp],
[X.sub.3dp]) and the average number of people killed or injured in a
4-year period (marked with [Z.sub.1], [Z.sub.2], [Z.sub.3] and
[S.sub.1], [S.sub.2], [S.sub.3]). Regression equations of the
statistical models of these parameters are:
[X.sub.i] = a[N.sup.b][L.sup.c], (5)
[X.sub.im] = a[N.sup.b][L.sup.c], (6)
[X.sub.idp] = a[N.sup.b][L.sup.c], (7)
[Z.sub.i] = a[N.sup.b][L.sup.c], (8)
[S.sub.i] = a[N.sup.b][L.sup.c], (9)
where a = 1; N--AADT on the studied road section, vpd; L--length of
the studied road section, km; b--linear regression coefficient depending
on traffic volume; c--linear regression coefficient depending on the
length of road section; i = 1, 2, 3.
Regression equations, obtained for the black spots of all types (t
= 1, 2, 3), forecast the averages of the studied random variables with
the errors in the interval (1; 2.3), determination coefficients
[R.sup.2] get the values in the interval (0.69; 0.99), thus, it could be
stated that the dependency of random variables on the traffic volume and
the length of the black spot gives a fairly good explanation of the
dispersion of random variables [X.sub.1], [X.sub.2], [X.sub.3],
[X.sub.1m], [X.sub.2m], [X.sub.3m], [X.sub.1dp], [X.sub.2dp],
[X.sub.3dp], [Z.sub.1], [Z.sub.2], [Z.sub.3] and [S.sub.1], [S.sub.2],
[S.sub.3]
The optimization of the selection of measures which was conducted
to reduce the number of black spots on the roads also is described. In
the result different mathematical models were constructed. One of them
allows the specialists to optimally select road safety measures to
maximally reduce the number of people killed under unrestricted amount
of funds (10), the other--under restricted financial possibilities (11).
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)
and
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)
where [bar.[[alpha].sub.m]], [bar.[[alpha].sub.dp]]--average number
of people killed in one vehicle-involved road accident
([bar.[[alpha].sub.m]]) in one road accident with pedestrians and
cyclists ([bar.[[alpha].sub.dp]]); [kappa]--number of preventive
measures; [d.sub.imj], [d.sub.idoj]--average decrease in the number of
vehicle-involved accidents ([d.sub.imj],) and in the number of accidents
with pedestrians and cyclists ([d.sub.idoj]) after implementation of the
measure j in the black spot i;--cost of implementation of the measure
' in the black spot i. Parameter C is the amount of money allocated
for the implementation of measures.
Conclusions
Having made the analysis of the effect of road safety measures
implemented it was determined that from the studied 48 high-accident
road sections the expected improvement results in safety situation were
not achieved on 50% of road sections; the expected reduction in accident
rate was not achieved on 15% of road sections; the implemented road
safety measures were justified on 50% of road sections.
Examination of 1995-2006 results of the substantiation of road
safety measures implemented on the roads of national significance of
Lithuania, the analysis of accident indices after implementation of
special measures and the detail analysis of only partly justified road
safety measures showed that methodology for the substantiation of road
safety measures, currently used in Lithuania, does not allow to
accurately enough select road safety measures for the high-accident road
sections.
Having made the analysis and evaluation of foreign methodologies
for the substantiation of road safety measures and based on the
Microsoft Office Access program the database was created operating as a
separate program or as a sub-program Evaluation of the Road Safety
Measures.
The created and tested calculation program Evaluation of the Road
Safety Measures gives a possibility to determine the effect of the
engineering road safety measures planned to be implemented on a certain
high-accident location or black spot, using the suggested forecast
coefficients of the impact and elasticity of road safety measures, of
traffic volume and the growth in prices of preventive measures and
taking into consideration the country's economic situation the
program enables to forecast road accidents and accident losses.
In the result of research activities the list of road safety
measures has been expanded by 59 new measures. A new program Evaluation
of the Road Safety Measures uses already 131 safety measure. It gives a
possibility to enter individual road safety measures.
With the use of high-volume sample n = 7095, data on the fatal and
injury accidents on the main roads in 1997-2008, the main numerical
characteristics of the number of people killed and injured in one road
accident were found. The performed accident modelling allows
constructing mathematical models for the optimum selection of preventive
measures for the main roads.
In order to make forecasts of the safety situation on the newly
built or reconstructed road sections and to prevent the occurrence of
black spots it is necessary to carry out a very comprehensive analysis
of historical data on the change of road accidents and traffic volume on
the roads of Lithuania, to divide the roads into plenty of homogenous
sections and to adapt the methods of mathematical statistics and
optimization methods.
DOI: 10.3846/bjrbe.2010.17
Received 21 December 2009; accepted 12 April 2010
References
Cygas, D.; Laurinavicius, A.; Ratkeviciute, K. 2004. Economic
Estimation of Traffic Safety Improvement Measures in Lithuania, in Proc
of the 3rd European Transport Congress (ETC) "Current Changes on
the European Transport Market": collection of papers, 22-23 April
2004, Opatia, Croatia, 321-326.
Elvik, R. 2009. An Exploratory Analysis of Models for Estimating
the Combined Effects of Road Safety Measures, Accident Analysis &
Prevention 41(4): 876-880. doi:10.1016/j.aap.2009.05.003
Hakkert, A. S.; Gitelman, V. 2004. Road Safety Measures on Rural
Roads in a Cost-Benefit Context, in Proc of the 17th International
Co-operation on Theories and Concepts in Traffic Safety (ICTCT) Workshop
"Effective Solution for Improving Road Safety in Rural Areas
Integrating the 4 Es: Education, Enforcement, Engineering and
Electronics". October 28-29, 2004, Tartu, Estonia, 9-28.
Hauer, E. 2005. Observational Before-After Studies in Road Safety.
Estimating the Effect of Highway and Traffic Engineering Measures on
Road Safety. Dept of Civil Engineering, University of Toronto. Canada.
289 p.
Jonsson, T. 2005. Predictive Models for Accidents on Urban links--a
Focus on Vulnerable Road Users. PhD thesis. Lund University, Lund
Institute of Technology. Bulletin 226. 109 p.
Kapski, D.; Leonovich, I.; Ratkeviciute, K. 2007. Theoretical
Principles of Forecasting Accident Rate on the Conflict Sections of the
Cities by the Method of Potential Danger, The Baltic Journal of Road and
Bridge Engineering 2(3): 133-140.
Kapski, D.; Leonovich, I.; Ratkeviciute, K.; Miskinis, D. 2008.
Implementation of Experimental Research in Road Traffic: Theory and
Practice, The Baltic Journal of Road and Bridge Engineering 3(2):
101-108. doi:10.3846/1822-427X.2008.3.101-108
Peltola, H. 2000. Background and Principles of the Finnish Safety
Evaluation Tool, TARVA, in Proc of the 13th International Co-operation
on Theories and Concepts in Traffic Safety (ICTCT) Workshop
"Evaluation of Traffic Safety Measures". October5-7, 200,
Corfu, Greece, 67-79.
Peltola, H. 2007. Evaluating Measures in Order to Achieve Safety
Targets, 14th International Conference "Road Safety on Four
Continents" (RS4C). November 14-16, 2007, Bangkok, Thailand. 8 p.
Ratkeviciute, K.; Cygas, D.; Laurinavicius, A.; Maciulis, A. 2007.
Analysis and Evaluation of the Efficiency of Road Safety Measures
Implemented on Lithuanian Roads, The Baltic Journal of Road and Bridge
Engineering 2(2): 81-87.
Ratkeviciute, K.; Cygas, D.; Bernotaite, I. 2008. Analysis of
Methodologies for the Evaluation of Effects of Road Safety Measures, in
Proc of the 7th International Conference "Environmental
Engineering": selected papers, vol. 3. Ed. by Cygas, D.; Froehner,
K. D. May 22-23, 2008, Vilnius, Lithuania. Vilnius: Technika, 1214-1222.
Yannisa, G; Gitelmanb, V; Papadimitrioua, E.; Hakkert, A. S.;
Winkelbauerc, M. 2008. Testing a Framework for the Efficiency Assessment
of Road Safety Measures, Transport Reviews 28(3): 281-301.
doi:10.1080/01441640701615120
Kornelija Ratkeviciute
Technological Science, Civil Engineering 02T Dept of Roads, Vilnius
Gediminas Technical University, Sauletekio al. 11, 10223 Vilnius,
Lithuania
E-mail: kornelija.ratkeviciute@vgtu.lt
Table 1. Reliability and adequacy of forecasting results
Number of
Statistical observations Standard
data (quarters) R [R.sup.2] deviation p
1997-2006 40 0.92 0.85 16.59 8.41E-14
1997-2007 44 0.90 0.80 18.19 2.41E-13
1997-2008 48 0.79 0.63 25.31 8.23E-09