Next Article in Journal
Novel Linguistic Steganography Based on Character-Level Text Generation
Next Article in Special Issue
Generalized Quasi-Einstein Manifolds in Contact Geometry
Previous Article in Journal
Measuring Arithmetic Word Problem Complexity through Reading Comprehension and Learning Analytics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling Tap Water Consumer Ratio

Department of Statistics, Faculty of Science, University of Gazi, Ankara 06560, Turkey
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(9), 1557; https://doi.org/10.3390/math8091557
Submission received: 4 August 2020 / Revised: 5 September 2020 / Accepted: 6 September 2020 / Published: 10 September 2020

Abstract

:
Increasing population and the rising air temperatures are known as factors that cause water depletion in the watersheds. Therefore, it is important to accurately predict the future ratios of tap water consumers using the same watershed to the population living in the specified area, to produce better water policies and to take the necessary measures. Predictions can be made by a growth curve model (GCM). Parameter estimations of the GCM are usually based on the ordinary least square (OLS) estimator. However, the outlier presence affects the estimations and the predictions, which are obtained by using the estimated model. The present article attempts to construct first- and third-order GCMs with robust least median square (LMS) and M estimators to make short-term predictions of ratios of tap water consumers. According to the findings, parameter estimations of the models, the outliers, and the predictions vary with respect to the estimators. The M estimator for short-term predictions is suggested for use, due to its robustness against outlier points.

1. Introduction

Water has vital importance for the survival of living beings. It is the basic element in terms of maintaining life in nature and human activities. The fact that 75% of the Earth’s surface is covered with water creates a perception that water scarcity would never be an issue to be discussed; more than 97% of this water is seawater, 2% is a mass of ice, and the large part of the remaining 1% is groundwater which is difficult to reach [1]. Only a tiny fraction of the water that forms the large part of our planet we live consists of healthy drinkable water. Thankfully, this water is renewed by nature’s solar-powered water cycle. With the evaporation fueled by the sun’s energy, water vapor is carried to the atmosphere. Of this evaporation, 86% occurs from the sea and 14% from land [1]. Even though an equal amount falls back to the Earth as rain, sleet, or snow, the distribution of water is more on continents than in oceans. With the transfer of water from the sea to land in this repeatable process, the renewable local drinking water resources occur. However, increasing population and air temperatures are effective factors in reducing drinking water availability per capita. The demand for drinking water will increase due to increasing population. A reduction in drinking water availability with increasing water demand, will equal an increase in stress on water supplies, as water supplies cannot be replenished quickly enough to meet demand. Moreover, agriculture already accounts for approximately 70% of the drinking water withdrawals in the world and is commonly seen as one of the main factors behind the growing global scarcity of freshwater [2].
Although a lot of international agreements and declarations such as the International Convention on the Rights of the Child, [3], the Human Rights Council Resolution, [4], and the Water Framework Directive, [5], stated the right to safe drinking water and sanitation is an internationally recognized human right, a great part of the world’s population cannot benefit from this right [6].
Turkey, a country with land both in Asia and Europe and a population of approximately 81 million people, is exposed to dry seasons. European Environment Agency reported that Turkey will encounter moderate and high-level water scarcity in many areas [7]. Thus, it is obvious that Turkey is a candidate country to experience problems with water scarceness. Considering that the population is predicted to be close to 100 million in 2030, [8], it is crucial to take precautions to avoid water shortages and to produce better water policies. Hence, it is important for a country to both follow and also predict the tap water consumer ratios (TWCRs) particularly to take measures to decrease the negative effects of the reduction in tap water that could occur very soon.
Registered water use in a country is important to ensure planned and economical water use. Therefore, monitoring the TWCRs of watersheds and predictions of their near future values would make it easier to establish a basis for water-related precautions to form. In the literature, there are studies on water and its future predictions [9,10,11,12,13]. Information on the assessment of sustainable water consumption perception, the evaluation of direct and indirect water consumption through the water footprint indicator, and the link between urban services and water uses are examined elsewhere [14]. Water footprint is described as an indicator of water use in relation to human consumption [15]. To incorporate the advances of life cycle assessment and water footprint analysis, an associated indicators set has been developed (see Reference [16]). However, ignoring outliers which are points that differ from the bulk of the data or fit a different distribution could cause biases in the findings [17].
This work aims to construct growth curve models (GCMs) to predict Turkey’s TWCRs. It is the ratio of the number of households using tap water in a particular region to the total number of households in that region. In this study, there are 26 grouped cities corresponding to 26 particular regions. Grouping of the cities is according to the watershed use of households. In the construction of the models, the non-robust ordinary least square estimator (OLS), which is employed in general, and robust least median square (LMS) and M estimators are used. This study demonstrates that detected outlying points differ according to which estimator is employed. Hence, every estimated model and hence the short-term predictions that are produced by the estimated model vary with respect to estimators. Thus, making better predictions for TWCRs to take convenient water policies depends on employing robust estimators that handle outliers very well during the parameter estimations of the GCMs. It is highlighted that the presence of outlying points has an undue impact on the model’s parameter estimations and future predictions.
The paper is organized as follows. In Section 2 we introduce GCMs based on non-robust OLS and, robust LMS and M estimators. To the best of our knowledge, this is the first time that robust LMS and M estimators are studied in this context. We show the differences in the outlying points, estimated GCMs, and hence the predictions by using the data originated in Section 3. Here, the TWCRs of Turkey are estimated by considering that the straight-line growth model, in other words, the first-order GCM, of time could be fitted to data. In the estimation procedure, we use OLS, LMS, and M estimators as mentioned above. As well as detection of outlying points, predictions to TWCRs of chosen years are obtained. Furthermore, to follow the estimations and predictions obtained from the estimated first-order GCMs based on OLS, LMS, and M estimators, separately, the curves are plotted on a single figure. Considering that the data could match more a third-order polynomial the progress is repeated for third-order GCMs.

2. GCMs Based on OLS, LMS, and M Estimators

The GCM usually expressed as
Y = X B Z + ε ,
is the change in a growth that corresponds to the response variable Y . This model indicates analytically how the parameters ( B 0 , B 1 ,   etc . ) and their standard errors ( ε 0 , ε 1 ,   etc . ) behave in a deterministic procedure for varying points of time [18]. X and Z are the design matrices. Here, Z used for grouped repeated measures is not taken into account since only the growth of Turkey’s TWCRs in watersheds on different time points is the subject to be researched. At this point, the vector of unknown parameters, the error, and the design matrix are denoted as B m × 1 , ε p × n , and X p × m , respectively. Each column of ε is assumed to be distributed as p -variate normal with 0 _ the mean vector and Σ > 0 the unknown covariance matrix. Additionally, Y is distributed as N p × n ( X B Z ,   Σ ,   I n ) where X B Z is the expected value, Σ and   I n are the covariance matrices of Y i j ( i fixed and j = 1 , , p ) and Y i ( i = 1 , , n ), respectively [18]. The number of time points examined on each of n observations is denoted by p and m 1 is the degree of the polynomial in time.
The OLS estimator of B , which is defined as B ^ O L S , is obtained from B ^ O L S = ( X X ) 1 X Y Z ( Z Z ) 1 . The B ^ 0 is the expected value of Y at time point 0 and called as the estimation of coefficient B 0 . The B ^ 1 is the expected value of Y when a one-unit change in time has occurred for observation i and called the estimation of coefficient B 1 . In addition, the OLS estimation of Σ , described as Σ ^ O L S , is based on B ^ O L S and is calculated from [18].
Σ ^ O L S = 1 n ( ( Y X B ^ O L S Z ) ( Y X B ^ O L S Z ) ) ,
Regarding the detection of outliers, the sum of squares of residuals of the ith observation is calculated from e i 2 ( B ^ O L S , Σ ^ O L S ) = ( y i ( X B ^ O L S Z ) i ) Σ ^ O L S 1 ( y i ( X B ^ O L S Z ) i ) . Since e i 2 ( B ^ O L S , Σ ^ O L S ) is chi-square distributed with p degrees of freedom, the calculated value of it is compared with the critical value determined from χ p , 1 α 2 , where α denotes the significance level. If the sum of squares of a suspicious observation is larger than the critical value, it would be appraised as an outlier [19,20]. The definition and explanations for e i 2 ( B ^ O L S , Σ ^ O L S ) mentioned above are also valid for LMS and M estimators.
The estimation procedure of B and Σ with robust LMS and M estimators depends on the weighted least square (WLS) estimator. Thus, the estimation procedure for WLS is based on minimizing i = 1 n w i e i 2 ( B ^ W L S , Σ ^ W L S ) , [21]. The WLS estimator of B , which is denoted as B ^ W L S , is computed from
B ^ W L S = ( X Σ ^ W L S 1 X ) 1 X Σ ^ W L S 1 Y W Z ( Z W Z ) 1 ,
and the estimation of the weighted covariance matrix is computed from
Σ ^ W L S = Y H Y t r ( H ) ,
where H = W W Z ( Z W Z ) 1 Z W . The notation “ t r ” denotes the trace. The elements of the diagonal weight matrix W that is used in Equation (3) and the calculation of H vary according to which estimator will be obtained. For instance, the ith element, w i , of the diagonal weight matrix, W t , is defined as
w i = { 1 , if   i th   observation   exist   in   i th   combination 0 , otherwise ,
when employing the LMS estimator and t is a value that ranges from 1 to C ( n , h ) . Here, C ( n , h ) denotes the number of h-combinations from a given data set of n elements and h is calculated from
h = n / 2 + ( p + 1 ) / 2 .
The notation . means rounding to a lower integer. Then, the LMS estimators of B and Σ , defined as B ^ L M S and Σ ^ L M S , can be easily obtained by regarding the minimization problem of the objective function
min B ^ W L S t   ( median ( e i 2 ( B ^ W L S t , Σ ^ W L S t ) ) ) ,
where i = 1 , , n [19].
The M estimator, B ^ M , is obtained by solving the objective function
min B ^ W L S k   i = 1 n ρ ( e i 2 ( B ^ W L S k , Σ ^ W L S k ) ) ,
where k = 0 , 1 , 2 , . Here, the value of B ^ L M S is used for the initial point B ^ W L S 0 . In Equation (8), ρ indicates a function which has a minimum at 0 for all e i 2 ( B ^ W L S k , Σ ^ W L S k ) and k shows the iteration number. In this instance, Tukey’s ρ function defined as
w i ( e i 2 ( B ^ W L S k , Σ ^ W L S k ) ) = ρ ( e i 2 ( B ^ W L S k , Σ ^ W L S k ) ) e i 2 ( B ^ W L S k , Σ ^ W L S k ) ,
is used to compute B ^ M [20]. In Equation (9), ρ ( e i 2 ( B ^ W L S k , Σ ^ W L S k ) ) means the derivative of ρ ( e i 2 ( B ^ W L S k , Σ ^ W L S k ) ) . The ith diagonal element of the diagonal weight matrix W k is obtained from:
w i ( e i 2 ( B ^ W L S k + 1 , Σ ^ W L S k + 1 ) ) = { ( 1 ( e i 2 ( B ^ W L S k , Σ ^ W L S k ) / c ) 2 ) 2 , 0 e i 2 ( B ^ W L S k , Σ ^ W L S k ) c 0       , e i 2 ( B ^ W L S k , Σ ^ W L S k )   > c   .
Here, c = 18.0909 and is calculated as the value that provides E χ p 2 [ ρ ( e i 2 ( B ^ M , Σ ^ M ) ) ] / ρ ( c ) = h / n , where E χ p 2 is the expected value obtained from the chi-squared distribution with p degrees of freedom.

3. Application to the GCMs on Predictions of TWCRs in Turkey for Deciding Which Estimator Is the Best

3.1. The Dataset

Turkey consists of eighty-one cities. These cities are categorized into twenty-six groups according to which local watershed they benefit from [8]. Table 1 summarizes TWCRs for each of them from 2001 to 2004 and at two-year intervals from 2006 to 2016.
In this research work, the response variable Y 26 × 10 represents a matrix and the observed value of Y i j denoted by y i j is the TWCR of group i in the year j , with i = 1 , , 26 and j = 1 , , 10 . The design matrix X is employed in two ways since the data could specify a functional form of linear or a cubic growth. Firstly, it is a 10 × 2 dimensional matrix where the first column consists of 1′s and the second column is the numbers 1 to 4 and 6, 8, 10, 12, 14, 16 used for the chosen years 2001 to 2004 and 2006, 2008, 2010, 2012, 2014, 2016, respectively. Hence, it is possible to preserve years as the unit of time. As it has been explained, Z does not affect the estimations. Thus, it is taken as an 1 × 26 dimensional vector consisting of 1′s. With the benefit of this design, the parameters of GCMs denoted as B 0 and B 1 would be estimated. Three different methods including OLS, LMS, and M are used to reconstruct GCMs, separately. This makes it possible to show that the differences in the results of identified outlier observations vary regarding the methods.
In the second part of the study, the design matrix
X 10 × 4 = [ 1 1 1 1 1 2 2 2 2 3 1 3 3 2 3 3 . . . . . . . . . . . . 1 16 16 2 16 3 ] ,
where the first column consists of ones and the other columns consists of the numbers that correspond to the chosen years [18], is employed. The reason for using the first, second, and third power of these numbers in the design matrix, respectively, is to build third-order GCMs. [17]. Here, Z is employed as defined previously in the construction of the first-order GCMs.

3.2. Detection of Outliers

The results of detecting outliers in the data, which consists of TWCRs in Turkey’s watersheds at different time points, and the parameter estimations of the GCMs according to methods mentioned above are summarized in Table 2, when α = 0.05 . This table summarizes the findings that are observed for first- and third-order GCMs. Watershed number 9 has been identified as outlier with both non-robust and robust estimators. This is strong evidence that there is an outlier in the data. However, when applying robust LMS and M estimators, watersheds numbered 4, 19, 25, and 1, 5, 22, respectively, are detected as outlying points, besides the watershed numbered 9. Thus, it is safe to infer that the predictions obtained from the estimated GCM based on OLS can be adversely affected by the undetected outliers. By definition, LMS and M estimators are more resistant to outliers compared to the OLS estimator [18,19,20]. Therefore, it is suggested to consider the predictions obtained from the estimated GCM based on these estimators.

3.3. Results

To show the differences in the estimated GCMs depended on OLS, LMS, and M estimators, Figure 1a,b are plotted. The horizontal line denotes the numbers corresponding to years and the vertical line denotes the predictions of ratios. Regarding the estimated first-order GCMs in Figure 1a, it is observed that the GCMs based on the OLS, LMS, M estimators are different. The observed predictions from the OLS appear to be larger when compared with the predictions from the LMS and M estimators. Knowing that the OLS estimator is being influenced by the outlier points, it is better to evaluate the predictions obtained from LMS and M estimators due to their robustness to outlier points [18,19,20]. Even, in general, the predictions of TWCRs obtained from the OLS and LMS estimators are higher than the predictions of TWCRs obtained from the M estimator. The M estimator is more resistant to outliers than the OLS and LMS estimators [18,19,20]. Therefore, it is recommended to evaluate the results obtained from the M estimator. For instance, the predictions for 2021 seem to be approximately 86%, 89%, and 90% in the case of using M, LMS, and OLS estimators, respectively. Thus, it is highlighted that the outlying observations affect the results. The predictions of the ratios of watersheds at years 2017 to 2021 (corresponds to 17 to 21, respectively) could be seen from these graphs as well. The values of the predictions based on the three methods tend to increase over the years.
In addition, estimating and predicting procedures are repeated as the data could be more appropriate for a third-order GCM. Figure 1b illustrates the estimated third-order GCMs after using OLS, LMS, and M estimators and the predictions of TWCRs for watersheds in Turkey from 2017 to 2021 (corresponds to 17 to 21, respectively).
The differences between the estimated third-order polynomials can be seen clearly in this figure. Predictions on TWCRs of Turkey’s watersheds have risen steadily, particularly after 2015. Moreover, even the results for M estimator are much lower than the results for OLS and LMS estimators. In addition, the predictions based on the M estimator are observed below 100. Thus, predictions obtained from this estimator are said to be acceptable since the vertical line in Figure 1b denotes the ratios. Consequently, it is proposed to consider the observed predictions from the M estimator because of its robustness against outlying points.

4. Conclusions

GCMs, as statistical growth models for short-term predictions, are used for various studies. Based on tap water consumer data in Turkey recorded between the years 2001 and 2016, this study investigated the TWCRs of Turkey’s watersheds with first- and third-order GCMs for short-time predictions. To estimate the parameters of the models, OLS, LMS, and M estimators are used. Usage of both robust and non-robust estimators allowed us to remark on the differences in parameter estimates (Table 2) and short-term predictions for TWCRs (Figure 1). A legitimate clarification for these findings seems to be the existence of outliers in the data. The predictions obtained through the M estimator are assumed to be the best, due to its robustness against outlier points [18,19,20]. According to these predictions, the TWCRs for Turkey’s watersheds will constantly increase. Furthermore, a prediction based on the estimated third-order GCM for the year 2021 is expected to be approximately 5% more than in 2020. Hence, making short-term predictions with the robust M estimator means a better view of the truth, which will lead us to produce better improvements on water policies.

Author Contributions

Conceptualization, M.E. and O.U.E.; methodology, M.E.; software, O.U.E.; validation, M.E. and O.U.E.; formal analysis, M.E. and O.U.E.; investigation, M.E.; resources, M.E. and O.U.E.; data curation, O.U.E.; writing—original draft preparation, M.E.; writing—review and editing, M.E. and O.U.E.; visualization, M.E.; supervision, M.E. and O.U.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors wish to thank the anonymous reviewers for their helpful suggestions and comments which led to the improvement of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Postel, S. The Last Oasis: Facing Water Scarcity; Earthscan: London, UK, 2013. [Google Scholar]
  2. Alexandratos, N.; Bruinsma, J. World Agriculture Towards 2030/2050: The 2012 Revision; ESA Working Papers No. 12-03; FAO: Rome, Italy, 2012. [Google Scholar]
  3. UNGA. Convention on the Rights of the Child; United Nations General Assembly: New York, NY, USA, 1989. [Google Scholar]
  4. UNGA. Human Rights and Access to Safe Drinking Water and Sanitation: Resolution Adopted by the Human Rights Council; United Nations General Assembly: New York, NY, USA, 2010. [Google Scholar]
  5. WFD (2000/60/EC). Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 Establishing a Framework for Community Action in the Field of Water Policy; L 327/3; Official Journal of the European Communities: Strasbourg, France, 2000. [Google Scholar]
  6. WWAP. The United Nations World Water Development Report: Water and Jobs; United Nations World Water Assessment Programme, UNESCO: Paris, France, 2016. [Google Scholar]
  7. EEA. The European Environment—State and Outlook 2005; EEA: Copenhagen, Denmark, 2005. [Google Scholar]
  8. Turkish Statistical Institute. Available online: https://biruni.tuik.gov.tr/medas/?kn=121&locale=tr (accessed on 18 February 2018).
  9. Nieswiadomy, M.L.; Molina, D.J. Urban water demand estimates under increasing block rates. Growth Chang. 1988, 19, 1–12. [Google Scholar] [CrossRef] [Green Version]
  10. Shukla, A.K.; Ojha, C.S.P.; Mijic, A.; Buytaert, W.; Pathak, S.; Garg, R.D.; Shukla, S. Population growth, land use and land cover transformations, and water quality nexus in the upper Ganga river basin. Hydrol. Earth Syst. Sci. 2018, 22, 4745–4770. [Google Scholar] [CrossRef] [Green Version]
  11. Wong, J.S.; Zhang, Q.; Chen, Y.D. Statistical modeling of daily urban water consumption in Hong Kong: Trend, changing patterns, and forecast. Water Resour. Res. 2010, 46, 1–10. [Google Scholar] [CrossRef]
  12. Zhou, S.L.; McMahon, T.A.; Walton, A.; Lewis, J. Forecasting daily urban water demand: A case study of Melbourne. J. Hydrol. 2000, 236, 153–164. [Google Scholar] [CrossRef]
  13. Chiu, C.C.; Shiang, W.J.; Lin, C.J.; Wang, C.H.; Chang, D.M. Water footprint analysis of second-generation bioethanol in Taiwan. J. Clean. Prod. 2015, 101, 271–277. [Google Scholar] [CrossRef]
  14. Gómez-Llanos, E.; Durán-Barroso, P.; Robina-Ramírez, R. Analysis of consumer awareness of sustainable water consumption by the water footprint concept. Sci. Total Environ. 2020, 721, 137743. [Google Scholar] [CrossRef] [PubMed]
  15. Chapagain, A.K.; Hoekstra, A.Y. Water Footprints of Nations; IHE Institute for Water Education, UNESCO: Delft, The Netherlands, 2004. [Google Scholar]
  16. Zhu, Y.; Liang, J.; Yang, Q.; Zhou, H.; Peng, K. Water use of a biomass direct-combustion power generation system in China: A combination of life cycle assessment and water footprint analysis. Renew. Sustain. Energy Rev. 2019, 115, 109396. [Google Scholar] [CrossRef]
  17. Ekiz, M.; Ekiz, O.U. Modelling tap water consumer ratio. In Proceedings of the International Conference on Mathematics and its Applications on Science and Engineering, Ankara, Turkey, 9–10 July 2020. [Google Scholar]
  18. Pan, J.X.; Fang, K.T. Growth Curve Models and Statistical Diagnostics, 1st ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  19. Rousseeuw, P.J.; Leroy, A.M. Robust Regression and Outlier Detection; Wiley: New York, NY, USA, 1987. [Google Scholar]
  20. Maronna, R.A.; Martin, D.R.; Yohai, V.J. Robust Statistics: Theory and Methods; Wiley: New York, NY, USA, 2006. [Google Scholar]
  21. Pendergast, J.F.; Broffitt, J.D. Robust estimation in growth curve models. Commun. Stat. Theory Methods 1985, 14, 1919–1939. [Google Scholar] [CrossRef]
Figure 1. Estimated growth curve models (GCMs) to TWCRs of Turkey’s watersheds and short-term predictions; (a) first-order GCMs; (b) third-order GCMs.
Figure 1. Estimated growth curve models (GCMs) to TWCRs of Turkey’s watersheds and short-term predictions; (a) first-order GCMs; (b) third-order GCMs.
Mathematics 08 01557 g001
Table 1. Grouped cities of Turkey according to which watershed they benefit from and tap water consumer ratios (TWCRs) over some years (all of the numerical values presented in the table are percentages).
Table 1. Grouped cities of Turkey according to which watershed they benefit from and tap water consumer ratios (TWCRs) over some years (all of the numerical values presented in the table are percentages).
No. of the Watershed CitiesYears
2001200220032004200620082010201220142016
1. Adana, Mersin82828285878887869898
2. Ankara919192969798969899100
3. Antalya, Isparta, Burdur79798081838384859294
4. Aydın, Denizli, Muğla67686869737374759599
5. Ağrı, Kars, Iğdır, Ardahan45474848515148505052
6. Balıkesir, Çanakkale59596062676768709090
7. Bursa, Eskişehir, Bilecik76848485909090919999
8. Erzurum, Erzincan, Bayburt66676869686869719293
9. Gaziantep, Adıyaman, Kilis76777779828183738192
10. Hatay, Kahramanmaraş, Osmaniye72707474767776779394
11. Kastamonu, Çankırı, Sinop52525253555558606063
12. Kayseri, Sivas, Yozgat73737576818282838889
13. Kocaeli, Sakarya, Düzce, Bolu, Yalova73737576828283839394
14. Konya, Karaman85858586888888899697
15. Kırıkkale, Aksaray, Niğde, Nevşehir, Kırşehir79797979838382847879
16. Malatya, Elazığ, Bingöl, Tunceli70717373777876788587
17.Manisa, Afyonkarahisar, Kütahya, Uşak72727273767778778384
18. Mardin, Batman, Şırnak, Siirt65676769707068707373
19. Samsun, Tokat, Çorum, Amasya61616265697071737684
20. Tekirdağ, Edirne, Kırklareli71737474818181848990
21. Trabzon, Ordu, Giresun, Rize, Artvin, Gümüşhane53565759646464667474
22. Van, Muş, Bitlis, Hakkari51525458565753557978
23. Zonguldak, Karabük, Bartın56565658626364656768
24. İstanbul9597979899999999100100
25. İzmir8588889190939292100100
26. Şanlıurfa, Diyarbakır65646568666966669596
Table 2. Parameter estimations and outlier detections with ordinary least squares (OLS), least median square (LMS), and M estimators.
Table 2. Parameter estimations and outlier detections with ordinary least squares (OLS), least median square (LMS), and M estimators.
Type of GCMMethodOutlier’s Number B ^ 0 B ^ 1 B ^ 2 B ^ 3
First-order GCMOLS968.51841.0376--
LMS4, 9, 1967.81851.0066--
M5, 9, 2268.67490.8559--
Third-order GCMOLS9, 1967.52542.29540.25420.0120
LMS9, 2565.00822.34700.26950.0128
M1, 9, 2261.36791.24840.09800.0054

Share and Cite

MDPI and ACS Style

Ekiz, M.; Ekiz, O.U. Modelling Tap Water Consumer Ratio. Mathematics 2020, 8, 1557. https://doi.org/10.3390/math8091557

AMA Style

Ekiz M, Ekiz OU. Modelling Tap Water Consumer Ratio. Mathematics. 2020; 8(9):1557. https://doi.org/10.3390/math8091557

Chicago/Turabian Style

Ekiz, Meltem, and Osman Ufuk Ekiz. 2020. "Modelling Tap Water Consumer Ratio" Mathematics 8, no. 9: 1557. https://doi.org/10.3390/math8091557

APA Style

Ekiz, M., & Ekiz, O. U. (2020). Modelling Tap Water Consumer Ratio. Mathematics, 8(9), 1557. https://doi.org/10.3390/math8091557

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