Stochastic modelling of skewed data exhibiting long range dependence
|
|
- David Boone
- 5 years ago
- Views:
Transcription
1 IUGG XXIV General Assembly 27 Perugia, Italy, 2 3 July 27 International Association of Hydrological Sciences, Session HW23 Analysis of Variability in Hydrological Data Series Stochastic modelling of skewed data exhibiting long range deendence S.M. Paalexiou Deartment of Water Resources, National Technical University of Athens
2 . Abstract Time series with long-range deendence aear in many fields including hydrology and there are several studies that have rovided evidence of long autocorrelation tails. Provided that the intensity of the long-range deendence in time series of a certain rocess, quantified by the self-similarity arameter, also known as the Hurst exonent H, could not be falsified, it is then essential that the variable of interest is modelled by a model reroducing long-range deendence. Common models of this category that have been widely used are the fractional Gaussian noise (FGN) and the fractional ARIMA (FARIMA). In case of a variable exhibiting skewness, the revious models can not be imlemented in a direct manner. In order to reserve skewness in the simulated series, a normalizing transformation is tyically alied in the real-life data at first. The models are then fitted to the normalized data and the roduced synthetic series are finally denormalized. In this aer, a different method is roosed, consisting of two arts. The first one regards the aroximation of the long-range deendence by an autoregressive model of high order AR(), while the second one regards the direct calculation of the main statistical roerties of the random comonent, that is mean, variance and skewness coefficient. The skewness coefficient calculation of the random comonent is done using joint samle moments. The advantage of the method is its efficiency and simlicity and the analytical solution.
3 2. Motivation Since Hurst (95) observed the long-term ersistence henomenon in the annual average streamflows of Nile, the same behaviour has been identified in numerous natural rocesses while, its imortance has been underlined by scientists in many controversial discilines. It seems that the Hurst henomenon is ubiquitous in nature and this makes it necessary to find adequate ways to model it. Many models have been roosed in the literature that reserve the Hurst behaviour, such as Fractional Gaussian Noise (FGN) (i.e. Mandelbrot, 969; Mandelbrot and Wallis 969), fast FGN (Mandelbrot, 97), broken line models (i.e. Ditlevsen, 97), fractional ARIMA (Hosking, 98), and recently symmetric moving average models (SMA) (Koutsoyiannis, 2; 22). If the Hurst behaviour aears in a rocess, it needs to modeled as it affects dramatically the time series structure. Another distinguished characteristic of hydrological rocesses, that needs to be modeled, is asymmetry. In this direction have been made many attemts to adat standard models to reserve the skewness (i.e. Matalas and Wallis, 976). Some of the revious models are not easy to aly as the arameters are not easy to estimate. while other can reserve the skewness but not the Hurst behaviour and vice versa. Other roblems are the narrow tye of autocorrelation functions that those model can simulate (excetion is the SMA model). In this study is roosed a general methodology to reserve both the Hurst behaviour and skewness. The framework of the methodology is simle: the Hurst henomenon is modeled from an autoregressive model of high order, AR(), while the skewness is reserved by evaluating the skewness coefficient of the random comonent of the model. The model should be easy to aly and suitable for any ractical uroses such as hydrologic design or water resources management.
4 3. Modelling Aroach In order to reserve the long-range deendence or the Hurst henomenon in the simulated time series, a high order autoregressive model is imlemented. The long-range deendence behaviour, is essentially the slow decay of the autocorrelation function with time. On the contrary, the AR() models are considered to be short-range deendence models. Nevertheless, as this study reveals, AR() models of high order can reroduce the Hurst henomenon sufficiently enough for any ractical modelling uroses. In the general case of order, the AR() model takes the following form: X t t + i= where ε t is the innovation or the random comonent and α i are coefficients. In order to fit the model to a dataset, the α i coefficients and the basic statistics (mean, standard deviation) of the ε t have to be estimated. The auto-covariance function γ k of the AR() model for lag k and for k> is given by g k a i g i-k i= X t-i a i The relacement of γ k with the samles estimates and the imlementation of the last equation times gives a linear system of equations that can be solved straightforwardly, evaluating therefore the α i coefficients. Finally, the mean and the variance of the ε t can be estimated using the following two equations. m t m Xt i j - y a i g i z k i= { s t 2 g - i= a i g i
5 4. Preserving the Skewness in an AR() Model To reserve asymmetry in the simulated time series, it is necessary to evaluate the skewness coefficient of the innovation, Csk ε. It can be shown that the third central moment of the innovation of the AR() model is Defining as multi-auto-covariance of order (m, m 2,,m n ) and lag (l, l 2,, l 3 ), E i j k i= it can be roven that the following equation is valid, y a i X t-i z { 3 i= (, ) - HL a3 i + 3 i= m H3L m 3 t m 3 X t - E i j k i= m H2,L j=i+ - H, j-il a 2 i a j + 3 i= m H,2L j=i+ -2 H, j-il a i a2 j + 6 i= Relacing the multi-auto-covariance terms in the revious equation with the samle estimates, given by y a i X t-i z { ( m ) ( l ) 2, l2 l m m n = E ( X m m m t l X) ( Xt l X) ( Xt l X) μ μ μ μ 2 n n - j=i+ k= j+ H, m j-i,k-il H,,L a i a j a k k max ( ) ( l, l2 l ), n + l l2 l m m2 m n n μ ( m, m2 m ) = ( ) ( x n i l μ X ) ( xi l μ 2 X ) ( xi l μ X ) k max l, l2 ln + i= it is then straightforward to estimate the μ 3ε in () and thus the μ 3ε C sk = ε 3 σ ε ( )
6 5. The Generalized AutoCorrelation Function (GACF) The major criticism of a high order AR() model would focus on the lack of arsimony, as estimation of the autocorrelation function u to lag is required to fit the model. Moreover, it is well known that the estimator of the ACF is highly variable and that it increases its variability with increasing lag (Bras and Rodriquez, 985). Consequently, the uncertainty in the estimation of the ACF would lead to uncertain validation of the model arameters. To overcome this disadvantage, it is roosed to fit a generalized ACF, ρ (G), to the first few emirical ACF values (where α, β, δ are ositive arameters and j is the lag). Subsequently, the fitted GACF can be used to extraolate ACF values for Emirical ACF high j. The figure deicts the emirical autocorrelation function of the Nilometer dataset (analysis follows) and the fitted GACF. The GACF has been fitted to the first emirical values of the ACF by minimizing the square error. ρ ( ) ( ) G = α β + j j δ β ACF Fitted GACF Lag
7 6. The Generalized Lambda Distribution (GLD) In order to reserve the skewness in the simulated series, the innovation ε t must be samled from a distribution with variable skewness. Such a flexible distribution is the GLD family. The GLD(λ, λ 2, λ 3, λ 4 ) family distributions originated from the one-arameter lambda distribution roosed by John Tukey (96) and was generalized for Monte Carlo simulation uroses by John Ramberg and Bruce Schmeiser (974). Although the GLD has been alied in many fields since the early 97s (Karian and Dudewicz, 2), it has never been used in hydrology. The GLD family with arameters λ, λ 2, λ 3 and λ 4, is defined in terms of its ercentile function, QH yl yl 3 - H - yl l 4 +l l 2 where < y <. The arameters λ and λ 2 are, resectively, location and scale arameters, while λ 3 and λ 4 determine the skewness and kurtosis of the distribution. The GLD robability density function is f HxL l 2 l 3 y l H - yl l 4 - l 4 The restrictions on λ, λ 2, λ 3 and λ 4, that yield a valid GLD distribution, the arameter sace and the skewness kurstosis sace are discussed in detail by Karian and Dudewicz (2). In the next figure GLD dfs are lotted with mean =, variance = and skewness coefficient Csk ranging from 4.5 at x QH yl PDF Csk Random Variable
8 7. Fitting the GLD and Samling If X is GLD(λ, λ 2, λ 3, λ 4 ) with λ 3 >-/4 and λ 4 >-/4, then its first four moments (Ramberg et al., 979), μ, μ 2, μ 3, μ 4 (mean, variance, skewness coefficient, and kurtosis coefficient), are given by m A l 2 where +l m 2 B-A2 l2 C sk m 3 s3 2 A3-3 BA+C s3 l2 3 C k m 4 s4-3 A4 +6 BA 2-4 CA+D s4 l2 4 A l l 4 + B -2 ΒHl 3 +, l 4 + L + 2 l l 4 + C 3 ΒHl 3 +, 2 l 4 + L - 3 ΒH2 l 3 +, l 4 + L + 3 l l 4 + D -4 ΒHl 3 +, 3 l 4 + L + 6 ΒH2 l 3 +, 2 l 4 + L - 4 ΒH3 l 3 +, l 4 + L + 4 l l 4 + and B is the Beta function defined as ΒHa, bl = t a- H - tl b- dt If we consider the innovation ε t as a random variable with known estimation of the mean, the variance, the skewness and the kurtosis coefficient, a GLD distribution can be fitted by solving numerically the revious nonlinear system. The mean, the variance and the skewness coefficient of ε t can be analytically estimated as described in slide four. At the moment there is no analytical way to estimate the kurtosis coefficient of ε t, but heuristically for this study was taken the minimum so as λ 3 < and λ 4 <, which imlies that the fitted GLD ranges form - to (Karian and Dudewicz, 2). Once the arameters λ, λ 2, λ 3, and λ 4 of the GLD are estimated, the samling is very easy as the ercentile function has a simle and analytical formulae.
9 8. Simulation Organogram Recorded data Standardized data Emirical Autocorrelation Function Generalized Autocorrelation Function Innovation arameters μ σ α i arameters 2 ε, ε, sk C ε Generalized Lambda Distribution Simulation of series for samle size equal to the recorded data. Autoregressive model X = a X + ε t i t i t i= Statistical analysis of the synthetic series
10 9. Original Data I: Nilometer Index Standardized Nilometer series indicating the annual minimum water level of the Nile river for the years 622 to 284 A.D. (663 years; Beran, 994) A data with small ositive skewness but with a large Hurst exonent value that verifies the multiscale fluctuations. Standardized Nilometer Index Annual scale 3 year scale year scale Year A.D. Frequency Nilometer Index Samle size Average St. deviation Skewness Hurst Exonent log Σ k H.5 R H Standardized Nilometer Index log k
11 . Original Data II: Annual Temeratures Northern Hemishere temerature anomalies in C with reference to mean (Standardized in the figure on the right). (2 years, Moberg et al., 25) A highly variable Northern Hemishere temerature reconstruction that reveals the large natural variability of the climate in multile scales. The multiscale variation is verified by the large Hurst exonent value. Standardized Average Temerature Annual scale 3 year scale year scale Year A.D. Frequency Annual Temeratures Samle size 979 Average St. deviation Skewness -.3 Hurst Exonent.9 log Σ k H.5 R H Standardized Annual Average Temerature log k
12 . Original Data III: Daily Average Temeratures Standardized average daily temeratures in July recorded at Den Helder station in Netherlands, from 9 to25. (source: Royal Netherlands Meteorological Institute) The histogram below deicts the ositive asymmetry of the dataset, while the long range deendence is manifested from the high Hurst exonent value. Standardized Average Temerature Daily scale day scale Monthly scale Time day Frequency Daily Temeratures Samle size 369 Average St. deviation Skewness.93 Hurst Exonent.72 log Σ k R H H Standardized Daily Average Temerature log k
13 2. Original Data IV: Daily Average Dew Points Standardized average daily dew oints in January recorded at Den Helder station in Netherlands, from 9 to 25. (source: Royal Netherlands Meteorological Institute) This dataset is suitable for the uroses of this study as it shows high negative asymmetry and a high Hurst exonent value. Standardized Average Dew Points Daily scale day scale Monthly scale Time day Frequency Daily Dew Points Samle size Average St. deviation Skewness Hurst Exonent log Σ k.5..5 R H H Standardized Daily Average Dew Points log k
14 3. Simulated Data I: Nilometer Index An auto regressive model of order 3, AR(3), was fitted to the original dataset in order to reserve the Hurst behaviour. The skewness coefficient of the innovation was evaluated, according to the methodology analysed in this study, to Csk ε =.46. A simulated series with statistics in close roximity to the observed ones is resented here. Standardized Nilometer Index Annual scale 3 year scale year scale Time Year Frequency Nilometer Index Samle size Average St. deviation Skewness Hurst Exonent log Σ k R H H Standardized Nilometer Index log k
15 4. Simulated Data II: Global Temeratures An auto regressive model of order 4, AR(4), was fitted to the original dataset in order to reserve the Hurst behaviour. The skewness coefficient of the innovation was evaluated, according to the methodology analysed in this study, to Csk ε =.48. A simulated series with statistics in close roximity to the observed ones is resented here. Standardized Average Temerature Annual scale 3 year scale year scale Time Year Frequency Standardized Annual Average Temerature Annual Temeratures Samle size 979 Average.4 St. deviation.7 Skewness -.2 Hurst Exonent.9 log Σ k.5..5 H.5 R 2.99 H log k
16 5. Simulated Data III: Daily Average Temeratures An auto regressive model of order 2, AR(2), was fitted to the original dataset in order to reserve the Hurst behaviour. The skewness coefficient of the innovation was evaluated, according to the methodology analysed in this study, to Csk ε =.94. A simulated series with statistics in close roximity to the observed ones is resented here. Standardized Average Temerature Daily scale day scale Monthly scale Time Day Frequency Daily Temeratures Samle size 369 Average.2 St. deviation Skewness.94 Hurst Exonent.7 log Σ k H.5 R H Standardized Daily Average Temerature log k
17 6. Simulated Data IV: Daily Average Dew Points An auto regressive model of order 2, AR(2), was fitted to the original dataset in order to reserve the Hurst behaviour. The skewness coefficient of the innovation was evaluated, according to the methodology analysed in this study, to Csk ε = A simulated series with statistics in close roximity to the observed ones is resented here. Standardized Average Dew Points Daily scale day scale Monthly scale Time Day Frequency Daily Dew Points Samle size Average St. deviation Skewness Hurst Exonent log Σ k.5..5 R H H Standardized Daily Average Dew Points log k
18 7. Simulated Mean and Standard Deviation The box lots below deict the estimated mean values of the simulated series. Blue dots reresent the observed means of each dataset ( as the series were standardised). The effect of the Hurst behaviour is clearly manifested by larger variability of the estimator in the series with the larger Hurst exonent values. The box lots below deict the estimated standard deviation values of the simulated series. Blue dots reresent the standardized values. The variability of the classic standard deviation estimator, is larger in the series with larger values of Hurst exonent, and it is also negative biased Simulated Mean Values Simulated Standard Deviation Values Nilometer Index Annual Temeratures Daily Temeratures Daily Dew Points Nilometer Index Annual Temeratures Daily Temeratures Daily Dew Points
19 8. Simulated Skewness and Hurst Exonent The box lots below deict the estimated skewness coefficient values of the simulated series. Blue dots reresent the observed values of each dataset. Given that the estimation of the Csk is highly uncertain, as its value is sensitive to outliers, it is encouraging that the model reroduces Csk values in roximity with the observed ones and with low variability..5 The box lots below deict the estimated Hurst exonent values of the simulated series. Blue dots reresent the observed values. Again the model, as the lot reveals, manages to reroduce sufficiently the Hurst behaviour, with the mean Hurst exonent value of the simulated series in agreement with the observed ones. Simulated Skewness Coefficient Values Simulated Hurst Exonent Values Nilometer Index Annual Temeratures Daily Temeratures Daily Dew Points Nilometer Index Annual Temeratures Daily Temeratures Daily Dew Points
20 9. Conclusions While the autoregressive models are considered to be short range ersistence models, it is concluded in this study that a higher order AR model reserves adequately the Hurst behaviour, for Hurst exonent values as high as.9. It seems that the model can reserve even more intense long-term ersistence but this needs to be further examined. To reserve the asymmetry, an analytical exression for the estimation of the skewness coefficient of the innovation is given. Subsequently, the innovation sequence is samled from a flexible skewed distribution, the so-called Generalized Lambda Distribution. The model manages to reserve sufficiently the skewness as the mean skewness coefficient of the simulated series is in roximity with the observed ones. As the simulated series are in accordance with the observed ones, the model can be used for any ractical modeling uroses. Overall, the roosed methodology is simle and robust.
21 2. References Beran, J., Statistics for Long Memory rocesses, Volume 6 of Monograhs on Statistics and Alied Probability, New York: Chaman and Hall, 994. Bras, R. L., and I. Rodriguez-Iturbe, Random Functions and Hydrology, 559., Addison-Wesley-Longman, Reading, Mass., 985. Ditlevsen, O. D., Extremes and first assage times, Doctoral dissertation, Tech. Univ. of Denmark, Lyngby, Denmark, 97. Hosking, J. R. M.,Fractional Differencing, Biometrika, vol. 68, , 98. Hurst, H., Long-term storage caacity of reservoirs, Transactions of the American Society of Civil Engineers 6, 77 88, 95. Karian, Z., Dudewicz, E., Fitting Statistical Distributions: The Generalized Lambda Distribution and Generalized Bootstra Methods, Boca Raton: CRC Press, 2. Koutsoyiannis, D., A generalized mathematical framework for stochastic simulation and forecast of hydrologic time series, Water Resources Research, 36(6), , 2. Koutsoyiannis, D., The Hurst henomenon and fractional Gaussian noise made easy, Hydrological Sciences Journal, 47(4), , 22. Mandelbrot, B. B., Une class de rocessus stochastiques homothetiques a soi: Alication a la loi climatologique de H. E. Hurst, C. R. Hebd. Seances Acad. Sci., 26, , 965. Mandelbrot, B. B., A fast fractional Gaussian noise generator, WaterResour. Res., 7(3), , 97. Mandelbrot, B. B., and J. R. Wallis, Comuter exeriments with fractional Gaussian noises,, Averages and variances, Water Resour. Res., 5(), , 969a. Mandelbrot, B. B., and J. R. Wallis, Comuter exeriments with fractional Gaussian noises, 2, Rescaled ranges and sectra, Water Resour. Res., 5(), , 969b. Mandelbrot, B. B., and J. R. Wallis, Comuter exeriments with fractional Gaussian noises, 3, Mathematical aendix, Water Resour. Res., 5(), , 969c. Matalas, N. C., and J. R. Wallis, Generation of synthetic flow sequences, in Systems Aroach to Water Management, edited by A. K. Biswas, McGraw-Hill, New York, 976. Moberg, A., D. M. Sonechkin, K. Holmgren, N. M. Datsenko, and W. Karlen, Highly variable Northern Hemishere temeratures reconstructed from low- and high-resolution roxy data, Nature, 433(726), 63 67, 25. Ramberg, J., Schmeiser, B., An aroximate method for generating asymmetric random variables, Communications of the ACM 7(2):78 82, 974. Ramberg, J. S., Dudewich, E. J., Tadikamalla, P. R., Mykytka, E. F., A robability distribution and its uses in fitting data. Communications in Statistics Theory and Methods 2(2):2 24, 979. Tukey, J. W., The Practical Relationshi Between the Common Transformations of Percentages of Counts and of Amounts, Technical Reort 36, Statistical Techniques Research Grou, Princeton University,96.
Statistics and Probability Letters. Variance stabilizing transformations of Poisson, binomial and negative binomial distributions
Statistics and Probability Letters 79 (9) 6 69 Contents lists available at ScienceDirect Statistics and Probability Letters journal homeage: www.elsevier.com/locate/staro Variance stabilizing transformations
More informationSINGLE SAMPLING PLAN FOR VARIABLES UNDER MEASUREMENT ERROR FOR NON-NORMAL DISTRIBUTION
ISSN -58 (Paer) ISSN 5-5 (Online) Vol., No.9, SINGLE SAMPLING PLAN FOR VARIABLES UNDER MEASUREMENT ERROR FOR NON-NORMAL DISTRIBUTION Dr. ketki kulkarni Jayee University of Engineering and Technology Guna
More informationON JARQUE-BERA TESTS FOR ASSESSING MULTIVARIATE NORMALITY
Journal of Statistics: Advances in Theory and Alications Volume, umber, 009, Pages 07-0 O JARQUE-BERA TESTS FOR ASSESSIG MULTIVARIATE ORMALITY KAZUYUKI KOIZUMI, AOYA OKAMOTO and TAKASHI SEO Deartment of
More informationConfidence Intervals for a Proportion Using Inverse Sampling when the Data is Subject to False-positive Misclassification
Journal of Data Science 13(015), 63-636 Confidence Intervals for a Proortion Using Inverse Samling when the Data is Subject to False-ositive Misclassification Kent Riggs 1 1 Deartment of Mathematics and
More informationA NOTE ON SKEW-NORMAL DISTRIBUTION APPROXIMATION TO THE NEGATIVE BINOMAL DISTRIBUTION
A NOTE ON SKEW-NORMAL DISTRIBUTION APPROXIMATION TO THE NEGATIVE BINOMAL DISTRIBUTION JYH-JIUAN LIN 1, CHING-HUI CHANG * AND ROSEMARY JOU 1 Deartment of Statistics Tamkang University 151 Ying-Chuan Road,
More informationQuantitative Aggregate Effects of Asymmetric Information
Quantitative Aggregate Effects of Asymmetric Information Pablo Kurlat February 2012 In this note I roose a calibration of the model in Kurlat (forthcoming) to try to assess the otential magnitude of the
More informationOnline Robustness Appendix to Are Household Surveys Like Tax Forms: Evidence from the Self Employed
Online Robustness Aendix to Are Household Surveys Like Tax Forms: Evidence from the Self Emloyed October 01 Erik Hurst University of Chicago Geng Li Board of Governors of the Federal Reserve System Benjamin
More informationDr. Maddah ENMG 625 Financial Eng g II 10/16/06. Chapter 11 Models of Asset Dynamics (1)
Dr Maddah ENMG 65 Financial Eng g II 0/6/06 Chater Models of Asset Dynamics () Overview Stock rice evolution over time is commonly modeled with one of two rocesses: The binomial lattice and geometric Brownian
More informationRevisiting the risk-return relation in the South African stock market
Revisiting the risk-return relation in the South African stock market Author F. Darrat, Ali, Li, Bin, Wu, Leqin Published 0 Journal Title African Journal of Business Management Coyright Statement 0 Academic
More informationSampling Procedure for Performance-Based Road Maintenance Evaluations
Samling Procedure for Performance-Based Road Maintenance Evaluations Jesus M. de la Garza, Juan C. Piñero, and Mehmet E. Ozbek Maintaining the road infrastructure at a high level of condition with generally
More informationModeling and Estimating a Higher Systematic Co-Moment Asset Pricing Model in the Brazilian Stock Market. Autoria: Andre Luiz Carvalhal da Silva
Modeling and Estimating a Higher Systematic Co-Moment Asset Pricing Model in the Brazilian Stock Market Autoria: Andre Luiz Carvalhal da Silva Abstract Many asset ricing models assume that only the second-order
More informationLinda Allen, Jacob Boudoukh and Anthony Saunders, Understanding Market, Credit and Operational Risk: The Value at Risk Approach
P1.T4. Valuation & Risk Models Linda Allen, Jacob Boudoukh and Anthony Saunders, Understanding Market, Credit and Operational Risk: The Value at Risk Approach Bionic Turtle FRM Study Notes Reading 26 By
More informationA GENERALISED PRICE-SCORING MODEL FOR TENDER EVALUATION
019-026 rice scoring 9/20/05 12:12 PM Page 19 A GENERALISED PRICE-SCORING MODEL FOR TENDER EVALUATION Thum Peng Chew BE (Hons), M Eng Sc, FIEM, P. Eng, MIEEE ABSTRACT This aer rooses a generalised rice-scoring
More informationA Multi-Objective Approach to Portfolio Optimization
RoseHulman Undergraduate Mathematics Journal Volume 8 Issue Article 2 A MultiObjective Aroach to Portfolio Otimization Yaoyao Clare Duan Boston College, sweetclare@gmail.com Follow this and additional
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationA Comparative Study of Various Loss Functions in the Economic Tolerance Design
A Comarative Study of Various Loss Functions in the Economic Tolerance Design Jeh-Nan Pan Deartment of Statistics National Chen-Kung University, Tainan, Taiwan 700, ROC Jianbiao Pan Deartment of Industrial
More information2/20/2013. of Manchester. The University COMP Building a yes / no classifier
COMP4 Lecture 6 Building a yes / no classifier Buildinga feature-basedclassifier Whatis a classifier? What is an information feature? Building a classifier from one feature Probability densities and the
More informationObjectives. 3.3 Toward statistical inference
Objectives 3.3 Toward statistical inference Poulation versus samle (CIS, Chater 6) Toward statistical inference Samling variability Further reading: htt://onlinestatbook.com/2/estimation/characteristics.html
More informationCapital Budgeting: The Valuation of Unusual, Irregular, or Extraordinary Cash Flows
Caital Budgeting: The Valuation of Unusual, Irregular, or Extraordinary Cash Flows ichael C. Ehrhardt Philli R. Daves Finance Deartment, SC 424 University of Tennessee Knoxville, TN 37996-0540 423-974-1717
More informationForecasting Exchange Rate between Thai Baht and the US Dollar Using Time Series Analysis
Forecasting Exchange Rate between Thai Baht and the US Dollar Using Time Series Analysis Kunya Bowornchockchai International Science Index, Mathematical and Computational Sciences waset.org/publication/10003789
More informationEffects of Size and Allocation Method on Stock Portfolio Performance: A Simulation Study
2011 3rd International Conference on Information and Financial Engineering IPEDR vol.12 (2011) (2011) IACSIT Press, Singaore Effects of Size and Allocation Method on Stock Portfolio Performance: A Simulation
More informationValue at Risk and Self Similarity
Value at Risk and Self Similarity by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) St. Andrews, March 17 th, 2009 Value at Risk and Self Similarity 1 1 Introduction The concept
More informationSupplemental Material: Buyer-Optimal Learning and Monopoly Pricing
Sulemental Material: Buyer-Otimal Learning and Monooly Pricing Anne-Katrin Roesler and Balázs Szentes February 3, 207 The goal of this note is to characterize buyer-otimal outcomes with minimal learning
More informationStochastic model of flow duration curves for selected rivers in Bangladesh
Climate Variability and Change Hydrological Impacts (Proceedings of the Fifth FRIEND World Conference held at Havana, Cuba, November 2006), IAHS Publ. 308, 2006. 99 Stochastic model of flow duration curves
More information***SECTION 7.1*** Discrete and Continuous Random Variables
***SECTION 7.*** Discrete and Continuous Random Variables Samle saces need not consist of numbers; tossing coins yields H s and T s. However, in statistics we are most often interested in numerical outcomes
More informationApplication of Monte-Carlo Tree Search to Traveling-Salesman Problem
R4-14 SASIMI 2016 Proceedings Alication of Monte-Carlo Tree Search to Traveling-Salesman Problem Masato Shimomura Yasuhiro Takashima Faculty of Environmental Engineering University of Kitakyushu Kitakyushu,
More informationAsian Economic and Financial Review A MODEL FOR ESTIMATING THE DISTRIBUTION OF FUTURE POPULATION. Ben David Nissim.
Asian Economic and Financial Review journal homeage: htt://www.aessweb.com/journals/5 A MODEL FOR ESTIMATING THE DISTRIBUTION OF FUTURE POPULATION Ben David Nissim Deartment of Economics and Management,
More informationThe Economic and Social BOOTSTRAPPING Review, Vol. 31, No. THE 4, R/S October, STATISTIC 2000, pp
The Economic and Social BOOTSTRAPPING Review, Vol. 31, No. THE 4, R/S October, STATISTIC 2000, pp. 351-359 351 Bootstrapping the Small Sample Critical Values of the Rescaled Range Statistic* MARWAN IZZELDIN
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of
More informationInformation and uncertainty in a queueing system
Information and uncertainty in a queueing system Refael Hassin December 7, 7 Abstract This aer deals with the effect of information and uncertainty on rofits in an unobservable single server queueing system.
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationThe Great Moderation Flattens Fat Tails: Disappearing Leptokurtosis
The Great Moderation Flattens Fat Tails: Disappearing Leptokurtosis WenShwo Fang Department of Economics Feng Chia University 100 WenHwa Road, Taichung, TAIWAN Stephen M. Miller* College of Business University
More informationJaime Frade Dr. Niu Interest rate modeling
Interest rate modeling Abstract In this paper, three models were used to forecast short term interest rates for the 3 month LIBOR. Each of the models, regression time series, GARCH, and Cox, Ingersoll,
More informationStatistical properties and Hurst- Kolmogorov dynamics in proxy data and temperature reconstructions
European Geosciences Union General Assembly Vienna, Austria 7 April May Session HS7 Change in climate, hydrology and society Statistical properties and Hurst- Kolmogorov dynamics in proxy data and temperature
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -26 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -26 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Hydrologic data series for frequency
More informationRobust Critical Values for the Jarque-bera Test for Normality
Robust Critical Values for the Jarque-bera Test for Normality PANAGIOTIS MANTALOS Jönköping International Business School Jönköping University JIBS Working Papers No. 00-8 ROBUST CRITICAL VALUES FOR THE
More informationMODELLING OF INCOME AND WAGE DISTRIBUTION USING THE METHOD OF L-MOMENTS OF PARAMETER ESTIMATION
International Days of Statistics and Economics, Prague, September -3, MODELLING OF INCOME AND WAGE DISTRIBUTION USING THE METHOD OF L-MOMENTS OF PARAMETER ESTIMATION Diana Bílková Abstract Using L-moments
More informationThe Skewed Generalized T Distribution Tree Package Vignette
The Skewed Generalized T Distribution Tree Package Vignette Carter Davis Setember, 05 Summary The Skewed Generalized T Distribution is a univariate 5-arameter distribution introduced by Theodossiou 998)
More informationTowards an advanced estimation of Measurement Uncertainty using Monte-Carlo Methods- case study kinematic TLS Observation Process
Towards an advanced estimation of Measurement Uncertainty using Monte-Carlo Methods- case study Hamza ALKHATIB and Hansjörg KUTTERER, Germany Key words: terrestrial laser scanning, GUM, uncertainty, Monte
More informationEmpirical Analysis of the US Swap Curve Gough, O., Juneja, J.A., Nowman, K.B. and Van Dellen, S.
WestminsterResearch http://www.westminster.ac.uk/westminsterresearch Empirical Analysis of the US Swap Curve Gough, O., Juneja, J.A., Nowman, K.B. and Van Dellen, S. This is a copy of the final version
More informationStochastic simulation of periodic processes with arbitrary marginal distributions
15 th International Conference on Environmental Science and Technology Rhodes, Greece, 31 August to 2 September 2017 Stochastic simulation of periodic processes with arbitrary marginal distributions Tsoukalas
More informationA Semi-parametric Test for Drift Speci cation in the Di usion Model
A Semi-arametric est for Drift Seci cation in the Di usion Model Lin hu Indiana University Aril 3, 29 Abstract In this aer, we roose a misseci cation test for the drift coe cient in a semi-arametric di
More informationIIntroduction the framework
Author: Frédéric Planchet / Marc Juillard/ Pierre-E. Thérond Extreme disturbances on the drift of anticipated mortality Application to annuity plans 2 IIntroduction the framework We consider now the global
More informationTrends in currency s return
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Trends in currency s return To cite this article: A Tan et al 2018 IOP Conf. Ser.: Mater. Sci. Eng. 332 012001 View the article
More informationSome Characteristics of Data
Some Characteristics of Data Not all data is the same, and depending on some characteristics of a particular dataset, there are some limitations as to what can and cannot be done with that data. Some key
More informationGENERATION OF STANDARD NORMAL RANDOM NUMBERS. Naveen Kumar Boiroju and M. Krishna Reddy
GENERATION OF STANDARD NORMAL RANDOM NUMBERS Naveen Kumar Boiroju and M. Krishna Reddy Department of Statistics, Osmania University, Hyderabad- 500 007, INDIA Email: nanibyrozu@gmail.com, reddymk54@gmail.com
More informationApproximate Variance-Stabilizing Transformations for Gene-Expression Microarray Data
Approximate Variance-Stabilizing Transformations for Gene-Expression Microarray Data David M. Rocke Department of Applied Science University of California, Davis Davis, CA 95616 dmrocke@ucdavis.edu Blythe
More informationCS522 - Exotic and Path-Dependent Options
CS522 - Exotic and Path-Deendent Otions Tibor Jánosi May 5, 2005 0. Other Otion Tyes We have studied extensively Euroean and American uts and calls. The class of otions is much larger, however. A digital
More informationNon-Inferiority Tests for the Ratio of Two Correlated Proportions
Chater 161 Non-Inferiority Tests for the Ratio of Two Correlated Proortions Introduction This module comutes ower and samle size for non-inferiority tests of the ratio in which two dichotomous resonses
More informationPARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS
PARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS Melfi Alrasheedi School of Business, King Faisal University, Saudi
More informationMonetary policy is a controversial
Inflation Persistence: How Much Can We Exlain? PAU RABANAL AND JUAN F. RUBIO-RAMÍREZ Rabanal is an economist in the monetary and financial systems deartment at the International Monetary Fund in Washington,
More informationCalculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the
VaR Pro and Contra Pro: Easy to calculate and to understand. It is a common language of communication within the organizations as well as outside (e.g. regulators, auditors, shareholders). It is not really
More informationOn accuracy of upper quantiles estimation
Hydrol. Earth Syst. Sci., 14, 2167 2175, 2010 doi:10.5194/hess-14-2167-2010 Author(s 2010. CC Attribution 3.0 License. Hydrology and Earth System Sciences On accuracy of upper quantiles estimation I. Markiewicz,
More information1 < = α σ +σ < 0. Using the parameters and h = 1/365 this is N ( ) = If we use h = 1/252, the value would be N ( ) =
Chater 6 Value at Risk Question 6.1 Since the rice of stock A in h years (S h ) is lognormal, 1 < = α σ +σ < 0 ( ) P Sh S0 P h hz σ α σ α = P Z < h = N h. σ σ (1) () Using the arameters and h = 1/365 this
More informationA data-adaptive maximum penalized likelihood estimation for the generalized extreme value distribution
Communications for Statistical Alications and Methods 07, Vol. 4, No. 5, 493 505 htts://doi.org/0.535/csam.07.4.5.493 Print ISSN 87-7843 / Online ISSN 383-4757 A data-adative maximum enalized likelihood
More informationLONG MEMORY IN VOLATILITY
LONG MEMORY IN VOLATILITY How persistent is volatility? In other words, how quickly do financial markets forget large volatility shocks? Figure 1.1, Shephard (attached) shows that daily squared returns
More informationInt. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS019) p.4301
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 0, Dublin (Session CPS09.430 RELIABILITY STUDIES OF BIVARIATE LOG-NORMAL DISTRIBUTION Pusha L.Guta Deartment of Mathematics and Statistics
More informationFeasibilitystudyofconstruction investmentprojectsassessment withregardtoriskandprobability
Feasibilitystudyofconstruction investmentrojectsassessment withregardtoriskandrobability ofnpvreaching Andrzej Minasowicz Warsaw University of Technology, Civil Engineering Faculty, Warsaw, PL a.minasowicz@il.w.edu.l
More informationEMH vs. Phenomenological models. Enrico Scalas (DISTA East-Piedmont University)
EMH vs. Phenomenological models Enrico Scalas (DISTA East-Piedmont University) www.econophysics.org Summary Efficient market hypothesis (EMH) - Rational bubbles - Limits and alternatives Phenomenological
More informationOvernight Index Rate: Model, calibration and simulation
Research Article Overnight Index Rate: Model, calibration and simulation Olga Yashkir and Yuri Yashkir Cogent Economics & Finance (2014), 2: 936955 Page 1 of 11 Research Article Overnight Index Rate: Model,
More informationBy choosing to view this document, you agree to all provisions of the copyright laws protecting it.
Coyright 2015 IEEE. Rerinted, with ermission, from Huairui Guo, Ferenc Szidarovszky, Athanasios Gerokostooulos and Pengying Niu, On Determining Otimal Insection Interval for Minimizing Maintenance Cost,
More informationIdiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective
Idiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective Alisdair McKay Boston University June 2013 Microeconomic evidence on insurance - Consumption responds to idiosyncratic
More informationAssessing the performance of Bartlett-Lewis model on the simulation of Athens rainfall
European Geosciences Union General Assembly 2015 Vienna, Austria, 12-17 April 2015 Session HS7.7/NP3.8: Hydroclimatic and hydrometeorologic stochastics Assessing the performance of Bartlett-Lewis model
More informationVolumetric Hedging in Electricity Procurement
Volumetric Hedging in Electricity Procurement Yumi Oum Deartment of Industrial Engineering and Oerations Research, University of California, Berkeley, CA, 9472-777 Email: yumioum@berkeley.edu Shmuel Oren
More informationA random variable X is a function that assigns (real) numbers to the elements of the sample space S of a random experiment.
RANDOM VARIABLES and PROBABILITY DISTRIBUTIONS A random variable X is a function that assigns (real) numbers to the elements of the samle sace S of a random exeriment. The value sace V of a random variable
More informationThe Supply and Demand for Exports of Pakistan: The Polynomial Distributed Lag Model (PDL) Approach
The Pakistan Develoment Review 42 : 4 Part II (Winter 23). 96 972 The Suly and Demand for Exorts of Pakistan: The Polynomial Distributed Lag Model (PDL) Aroach ZESHAN ATIQUE and MOHSIN HASNAIN AHMAD. INTRODUCTION
More information2.1 Properties of PDFs
2.1 Properties of PDFs mode median epectation values moments mean variance skewness kurtosis 2.1: 1/13 Mode The mode is the most probable outcome. It is often given the symbol, µ ma. For a continuous random
More informationWorst-case evaluation complexity for unconstrained nonlinear optimization using high-order regularized models
Worst-case evaluation comlexity for unconstrained nonlinear otimization using high-order regularized models E. G. Birgin, J. L. Gardenghi, J. M. Martínez, S. A. Santos and Ph. L. Toint 2 Aril 26 Abstract
More informationBasic Procedure for Histograms
Basic Procedure for Histograms 1. Compute the range of observations (min. & max. value) 2. Choose an initial # of classes (most likely based on the range of values, try and find a number of classes that
More informationUtility and the Skewness of Return in Gambling
The Geneva Paers on Risk and Insurance Theory, 9: 145 163, 004 c 004 The Geneva Association Utility and the Skewness of Return in Gambling MICHAEL CAIN School of Business, University of Wales, Hen Goleg,
More informationVolatility Clustering of Fine Wine Prices assuming Different Distributions
Volatility Clustering of Fine Wine Prices assuming Different Distributions Cynthia Royal Tori, PhD Valdosta State University Langdale College of Business 1500 N. Patterson Street, Valdosta, GA USA 31698
More informationLECTURE NOTES ON MICROECONOMICS
LECTURE NOTES ON MCROECONOMCS ANALYZNG MARKETS WTH BASC CALCULUS William M. Boal Part : Consumers and demand Chater 5: Demand Section 5.: ndividual demand functions Determinants of choice. As noted in
More informationRe-testing liquidity constraints in 10 Asian developing countries
Re-testing liquidity constraints in 0 Asian develoing countries Kuan-Min Wang * Deartment of Finance, Overseas Chinese University Yuan-Ming Lee Deartment of Finance, Southern Taiwan University Abstract
More informationAssicurazioni Generali: An Option Pricing Case with NAGARCH
Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance
More informationCausal Links between Foreign Direct Investment and Economic Growth in Egypt
J I B F Research Science Press Causal Links between Foreign Direct Investment and Economic Growth in Egyt TAREK GHALWASH* Abstract: The main objective of this aer is to study the causal relationshi between
More informationFinancial Econometrics
Financial Econometrics Volatility Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) Volatility 01/13 1 / 37 Squared log returns for CRSP daily GPD (TCD) Volatility 01/13 2 / 37 Absolute value
More informationThis homework assignment uses the material on pages ( A moving average ).
Module 2: Time series concepts HW Homework assignment: equally weighted moving average This homework assignment uses the material on pages 14-15 ( A moving average ). 2 Let Y t = 1/5 ( t + t-1 + t-2 +
More informationFinancial Economics (I) Instructor: Shu-Heng Chen Department of Economics National Chengchi University
Financial Economics (I) Instructor: Shu-Heng Chen Department of Economics National Chengchi University Lecture 7: Rescale Range Analysis and the Hurst Exponent Hurst exponent is one of the most frequently
More informationA Comparative Study of Various Forecasting Techniques in Predicting. BSE S&P Sensex
NavaJyoti, International Journal of Multi-Disciplinary Research Volume 1, Issue 1, August 2016 A Comparative Study of Various Forecasting Techniques in Predicting BSE S&P Sensex Dr. Jahnavi M 1 Assistant
More informationEstimating the Parameters of Closed Skew-Normal Distribution Under LINEX Loss Function
Australian Journal of Basic Applied Sciences, 5(7): 92-98, 2011 ISSN 1991-8178 Estimating the Parameters of Closed Skew-Normal Distribution Under LINEX Loss Function 1 N. Abbasi, 1 N. Saffari, 2 M. Salehi
More informationStock Market Risk Premiums, Business Confidence and Consumer Confidence: Dynamic Effects and Variance Decomposition
International Journal of Economics and Finance; Vol. 5, No. 9; 2013 ISSN 1916-971X E-ISSN 1916-9728 Published by Canadian Center of Science and Education Stock Market Risk Premiums, Business Confidence
More informationExample 1 of econometric analysis: the Market Model
Example 1 of econometric analysis: the Market Model IGIDR, Bombay 14 November, 2008 The Market Model Investors want an equation predicting the return from investing in alternative securities. Return is
More informationInformation Processing and Limited Liability
Information Processing and Limited Liability Bartosz Maćkowiak European Central Bank and CEPR Mirko Wiederholt Northwestern University January 2012 Abstract Decision-makers often face limited liability
More informationApplication of Sarima Models in Modelling and Forecasting Nigeria s Inflation Rates
American Journal of Alied Mathematics and Statistics, 4, Vol., No., 6-8 Available online at htt://ubs.scieub.com/ajams///4 Science and Education Publishing DOI:9/ajams---4 Alication of Sarima Models in
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (42 pts) Answer briefly the following questions. 1. Questions
More informationSELFIS: A Short Tutorial
SELFIS: A Short Tutorial Thomas Karagiannis (tkarag@csucredu) November 8, 2002 This document is a short tutorial of the SELF-similarity analysis software tool Section 1 presents briefly useful definitions
More information1. You are given the following information about a stationary AR(2) model:
Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4
More informationPower of t-test for Simple Linear Regression Model with Non-normal Error Distribution: A Quantile Function Distribution Approach
Available Online Publications J. Sci. Res. 4 (3), 609-622 (2012) JOURNAL OF SCIENTIFIC RESEARCH www.banglajol.info/index.php/jsr of t-test for Simple Linear Regression Model with Non-normal Error Distribution:
More informationGARCH Models. Instructor: G. William Schwert
APS 425 Fall 2015 GARCH Models Instructor: G. William Schwert 585-275-2470 schwert@schwert.ssb.rochester.edu Autocorrelated Heteroskedasticity Suppose you have regression residuals Mean = 0, not autocorrelated
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider
More informationTESTING THE CAPITAL ASSET PRICING MODEL AFTER CURRENCY REFORM: THE CASE OF ZIMBABWE STOCK EXCHANGE
TESTING THE CAPITAL ASSET PRICING MODEL AFTER CURRENCY REFORM: THE CASE OF ZIMBABWE STOCK EXCHANGE Batsirai Winmore Mazviona 1 ABSTRACT The Caital Asset Pricing Model (CAPM) endeavors to exlain the relationshi
More informationAbsolute Return Volatility. JOHN COTTER* University College Dublin
Absolute Return Volatility JOHN COTTER* University College Dublin Address for Correspondence: Dr. John Cotter, Director of the Centre for Financial Markets, Department of Banking and Finance, University
More information2.4 STATISTICAL FOUNDATIONS
2.4 STATISTICAL FOUNDATIONS Characteristics of Return Distributions Moments of Return Distribution Correlation Standard Deviation & Variance Test for Normality of Distributions Time Series Return Volatility
More informationBrownian Motion, the Gaussian Lévy Process
Brownian Motion, the Gaussian Lévy Process Deconstructing Brownian Motion: My construction of Brownian motion is based on an idea of Lévy s; and in order to exlain Lévy s idea, I will begin with the following
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationManagement of Pricing Policies and Financial Risk as a Key Element for Short Term Scheduling Optimization
Ind. Eng. Chem. Res. 2005, 44, 557-575 557 Management of Pricing Policies and Financial Risk as a Key Element for Short Term Scheduling Otimization Gonzalo Guillén, Miguel Bagajewicz, Sebastián Eloy Sequeira,
More informationINDEX NUMBERS. Introduction
INDEX NUMBERS Introduction Index numbers are the indicators which reflect changes over a secified eriod of time in rices of different commodities industrial roduction (iii) sales (iv) imorts and exorts
More informationHigh-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5]
1 High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5] High-frequency data have some unique characteristics that do not appear in lower frequencies. At this class we have: Nonsynchronous
More informationResearch Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and Its Extended Forms
Discrete Dynamics in Nature and Society Volume 2009, Article ID 743685, 9 pages doi:10.1155/2009/743685 Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and
More informationHyetosR: An R package for temporal stochastic simulation of rainfall at fine time scales
European Geosciences Union General Assembly 2012 Vienna, Austria, 22-27 April 2012 Session HS7.5/NP8.3: Hydroclimatic stochastics HyetosR: An R package for temporal stochastic simulation of rainfall at
More information