Application of Univariate Probability Distributions Fitting using Monte Carlo Simulation

  • Muhammad Ilyas Laboratory for Applied Mathematics and Data Analysis (LAMDA) Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.
  • Shaheen Abbas Laboratory for Applied Mathematics and Data Analysis (LAMDA) Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.
  • Afzal Ali Laboratory for Applied Mathematics and Data Analysis (LAMDA) Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.
Keywords: Anderson-Darling Test (ADT), Adequate Fitting Tests (AFT), Kolmogorov-Smirnov D-test (KST), Maximum Likelihood Estimation (MLE), probability distribution, Probability of Correct Selection (PCS)

Abstract

Abstract Views: 280

In this study, we showcased a univariate probability distribution through the application of the three sub- and super-exponential heavier-longer and lighter-shorter tails fitting. The univariate family includes the Lognormal, Gamma, and Weibull distributions. The adequacy of distribution tails was determined using Adequate Fitting Tests (ADTs) and a descriptive criterion. For this purpose, this study investigated the logarithm population data sets of the Karachi region from 1729 to 1946 and again from 1951 to 2018. The data contained both irregular and regular lengths and peaks. We considered three Lognormal, Gamma, and Weibull distributions for tails fitting, statistical methods for validation, and normality tests to verify stochastic behavior in both data sets. Weibull and Lognormal distribution tails were found to possess heavier distribution using two validation tests (maximum likelihood estimation and probability of correct selection). Furthermore, the univariate probability distributions were applied to Monte Carlo simulation in order to generate the actual population data. The results indicated that heavy-tailed Lognormal and Weibull distributions were adequate as compared to the more commonly used lighter tailed Gamma distribution. Hence, it was determined that the Monte Carlo Simulation performs appropriate Lognormal and Weibull distributions for irregular and regular data. Our results indicated that Lognormal-Weibull distributions are suitable for the prediction of both long-term and short-term forecasting.

Downloads

Download data is not yet available.

Author Biographies

Muhammad Ilyas, Laboratory for Applied Mathematics and Data Analysis (LAMDA) Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.

Assistant Professor (IPFP) at Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.

Afzal Ali, Laboratory for Applied Mathematics and Data Analysis (LAMDA) Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.

M.Phil. student at Mathematical Sciences Research Centre, Federal Urdu University of Arts, Sciences and Technology, Karachi, Pakistan.

References

Maio C, Schexnayder C. Simulation model probability distribution functions: Relationships between input data and goodness-of-fit tests. InIAARC/IFAC/IEEE. International symposium 1999 (pp. 103-108).

Silveira FV, Gomes-Silva F, Brito CC, et al. Modelling wind speed with a univariate probability distribution depending on two baseline functions. arXiv preprint arXiv:2101.03622. 2021. https://arxiv.org/abs/2101.03622

Leemis LM, McQueston JT. Univariate Distribution Relationships. Am Stat. 2008;62(1):45-53. https://doi.org/10.1198/000313008X270448

Johnson NL, Kotz S, Kemp AW. Univariate discrete distributions. John Wiley & Sons; 2005.

Johnson NL, Kotz SI, Balakrishnan N. Beta distributions. Continuous univariate distributions. 2nd ed. New York, NY: John Wiley and Sons. 1994:221-35.

Kalbfleisch JD, Prentice RL. The statistical analysis of failure time data. John Wiley & Sons; 2011.

Kale MD, Butar FB. Fractal analysis of time series and distribution properties of Hurst exponent (Doctoral dissertation, Sam Houston State University).

Hassan D, Abbas S, Ansari MR, Jan B. The study of sunspots and K-index data in the perspective of probability distributions. Int J Phys Social Sci. 2014;4(1):23-41.

Vose, David. Risk analysis: a quantitative guide. John Wiley & Sons;2008.

Leemis and M. Lawrence. Relationships among common univariate distributions. Am Stat. 1986;40:143-46.

Protassov, Rostislav S. EM-based maximum likelihood parameter estimation for multivariate generalized hyperbolic distributions with fixed λ. Stat Comput. 2004;14:67-77. https://doi.org/10.1023/B:STCO.0000009419.12588.da

Johnson, Norman L., Samuel Kotz, and Narayanaswamy Balakrishnan. Continuous univariate distributions, John wiley & sons; 1995.

Kim, Jin Seon, and Bong-Jin Yum. Selection between Weibull and lognormal distributions: A comparative simulation study. Comput Stat Data Analy. 2008;53:477-485. https://doi.org/10.1016/j.csda.2008.08.012

Metzler, Ralf, and Joseph Klafter. The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics. J Phys A: Math Gen. 2004;37:R161.

Papalexiou SM, Koutsoyiannis D, Makropoulos C. How extreme is extreme? An assessment of daily rainfall distribution tails. Hydrol Earth System Sci. 2013;17(2):851-62.

Gibbons JD, Chakraborti S. Nonparametric statistical inference. CRC press; 2014.

Michael JR. The stabilized probability plot. Biometrika. 1983;70(1):11-7.

Werner T, Upper C. Time variation in the tail behavior of Bund future returns. J Futures Markets: Futures, Options, and Other Derivative Products. 2004;24(4):387-98.

Kundu D, Manglick A. Discriminating between the Weibull and log‐normal distributions. Naval Res Logist (NRL). 2004;51(6):893-905. https://doi.org/10.1002/nav.20029

Prabhu SS, Montgomery DC, Runger GC. A combined adaptive sample size and sampling interval X control scheme. J Quality Technol. 1994; 26(3):164-76. https://doi.org/10.1080/00224065.1994.11979524

Corrado CJ, Su T. Implied volatility skews and stock index skewness and kurtosis implied by S&P 500 index option prices. J Derivat. 1997; 4(4):8-19.

Shpilberg, David C. The probability distribution of fire loss amount. J Risk Insur. 1977;103-115. https://doi.org/10.2307/251860

Améndola C, Kohn K, Reichenbach P, Seigal A. Invariant theory and scaling algorithms for maximum likelihood estimation. SIAM J Appl Algebr Geom. 2021;5(2):304-37.

Klüppelberg, Claudia. Subexponential distributions and integrated tails. J Appl Prob. 1988;25:132-141. https://doi.org/10.2307/3214240

Law, Averill MW. David Kelton, and W. David Kelton. Simulation modeling and analysis. New York: McGraw-Hill; 2000.

AbouRizk, Simaan M, Daniel W, Halpin, James R, Wilson. Fitting beta distributions based on sample data. J Construct Eng Manag. 1994;120:288-305. https://doi.org/10.1061/(ASCE)0733-9364(1994)120:2(288)

AbouRizk, Simaan M., and Daniel W. Halpin. Probabilistic simulation studies for repetitive construction processes. J Construct Eng Manag. 1990;116:575-594.https://doi.org/10.1061/(ASCE)0733-9364(1990)116:4(575)

Bakhtiari, Hamed, Jin Zhong, and Manuel Alvarez. Predicting the stochastic behavior of uncertainty sources in planning a stand-alone renewable energy-based microgrid using Metropolis–coupled Markov chain Monte Carlo simulation. Appl Energy. 2021;290:116719. https://doi.org/10.1016/j.apenergy.2021.116719

Published
2022-03-30
How to Cite
1.
Ilyas M, Abbas S, Ali A. Application of Univariate Probability Distributions Fitting using Monte Carlo Simulation. Sci Inquiry Rev. [Internet]. 2022Mar.30 [cited 2024Dec.22];6(1). Available from: https://journals.umt.edu.pk/index.php/SIR/article/view/1948
Section
Orignal Article