Enhancing the Performance of the Mont Carlo Simulation for Generating Correlated Data: An Application to Climate Change

Document Type : Research Paper

Authors

Abstract

The Mont Carlo simulation is widely used for data generation based on specified distribution. However, generation of correlated data is only possible for specific distributions (bivariate normal), while it is required to deal with non-normal distributions, e.g. in regional frequency analysis. To tackle the problem, an innovative approach was proposed by using the genetic algorithm. Since it is possible only to change the order of numbers in an array, our problem is not completely in agreement with the classic genetic algorithm. Therefore, we reordered the algorithm; the internal of each chromosome was changed instead of crossing two different ones. The proposed method was applied for one random variable, and it was shown that the genetic algorithm produces multi-solutions, while unique solution is obtained using the Mont Carlo simulation procedure. We presented an appropriate objective function to include any bivariate distribution and showed its hither versatility as compared with the Mont Carlo simulation (only for the bivariate normal distribution). The multi-solutions in our proposed methodology is due to the fact that any array can be reordered in different possible combinations, each of which having a unique correlation coefficient with another fixed array. It was shown that these different combinations increase as the correlation coefficient decreases.

Keywords


  1. Ang, A.H-S., and W.H. Tang. 1984. Probability concepts in engineering and design. Volume II- Decision, Risk, and Reliability. John Wiley and Sons, Inc. 562 p.
  2. Ang, A.H-S., and W.H. Tang. 2007. Probability concepts in engineering. Empasis on applications to civil and environmental engineering. John Wiley and Sons, Inc. 2nd edition. 406 p.
  3. Casterllarin, A., D.H., Burn, and A. Brath, 2008. Homogeniety testing: how homogeneous do heterogeneous cross-correlated regions seem? J. Hydrol.360,67-76.
  4. Hosking, J.R.M. 1994. The four- parameter kappa distribution. IBM J. Res. Develop.38, 251-258.
  5. Hosking, J.R.M., and J.R. Wallis. 1997. Regional frequency alalysis. An approach based on L-moments. Cambridge University Press. 224 p.
  6. Kleidon, A. 2010. Non-equilibrium thermodynamics, maximum entropy production and earth-system evolution. Philos. Tran. Royal Soc. A. 368(1910), 181-196. doi: 10.1098/ rsta. 2009. 01188.
  7. Rodriguez-Iturbe, I., A. Rinaldo, R. Rigon, R.L. Bras, A. Marani, and E. Ijjasz-Vasquez. 1992. Fractal structure as least energy patterns: The case of river networks. Geophys. Res. Lett. 19:889-892.