Using SAS/IML to Calculate L-Moments for the Univariate

Total Page:16

File Type:pdf, Size:1020Kb

Using SAS/IML to Calculate L-Moments for the Univariate USE OF PROC IML TO CALCULATE L-MOMENTS FOR THE UNIVARIATE DISTRIBUTIONAL SHAPE PARAMETERS SKEWNESS AND KURTOSIS Michael A. Walega Covance, Princeton, New Jersey Introduction and how it pertains to the assumption of normality. Exploratory data analysis statistics, such as those As discussed by Glass et al. (1972), incorrect generated by the SAS® procedure PROC conclusions may be reached when the normality UNIVARIATE (1990), are useful tools to assumption is not valid, especially when one-tail characterize the underlying distribution of data prior tests are employed or the sample size or to more rigorous statistical analyses. Assessment of significance level are very small. Hopkins and the distributional shape of data is usually Weeks (1990) also discuss the effects of highly non- accomplished by careful examination of the values normal data on hypothesis testing of variances. of the third and fourth central moments, skewness Thus, the exam-ination of the skewness (departure and kurtosis. However, when the sample size is from symmetry) and kurtosis (deviation from a small or the underlying distribution is non-normal, normal curve) is an important component of the information obtained from the sample skewness exploratory data analyses. and kurtosis can be misleading. Various methods to estimate skewness and kurtosis One alternative to the central moment shape have been proposed (MacGillivray and Balanda, statistics is the use of a linear combination of order 1988). For many years, the conventional statistics (L-moments) to examine the distributional coefficients of skewness and kurtosis, ϒ and κ shape characteristics of data. L-moments have (Hosking, 1990), have been used to describe the several theoretical advantages over the central shape characteristics of distributions. However, as moment shape statistics: Characterization of a pointed out by Hosking (1990) and Royston (1992), wider range of distributions, robustness to outliers these coefficients are not without limitations. Both and more accurate estimates in small sample sizes. are sensitive to minute changes in the tails of a distribution, susceptible to moderate outliers and This paper focuses on the development of a macro biased in small to moderately-sized samples from skew distributions. Also, the information conveyed program that uses SAS/IML® (1989) to generate the by the third and fourth central moments with regards central moment and L-moment distributional shape to the shape of a distribution can be difficult to parameters. In addition, the results of simulations, assess. Thus, it would be appropriate to determine conducted with various sample sizes and if other, more robust measures of skewness and distributions, will be presented. kurtosis can be used to assess the shape of a distribution. Background L-moments Largely through the influence of John Tukey’s work (1977), statisticians have increasingly emphasized One more robust measure are linear combinations the exploratory analysis of data prior to more formal of order statistics, or L-moments. In theory, L- statistical analyses (t-tests, ANOVA, etc.). Tukey moments are less prone to the effects of sampling has suggested that to fully understand the nature of variability as compared to conventional moments. a variable and its measurement, characteristics Hosking (1990) provides an excellent overview of other than the central tendency (mean) and the theory behind the derivation and application of variability (standard deviation) need to be L-moments as summary statistics for univariate examined. Many classical statistical tests rely on probability distributions. Royston (1992) compares the assumption that the underlying distribution of the the prop-erties of the conventional shape data (or residuals) is Gaussian. Bickel (1988) and Van Der Laan and Verdooren (1987) discuss the parameters to their L-moment counterparts for two concept of robustness lognormal dis-tributions. Rather than discuss the detailed theory behind L-moments, the reader is referred to the two aforementioned papers. Instead, a brief overview of the development of the outer extremes and the central portion samples of equations necessary to apply L-moments is size 4 (Royston, 1992). described below. As with the paper by Royston, the notation of Hosking (1990) will be employed. Scale-free versions of the L-moments for skewness, τ τ 3, and kurtosis, 4, can be written as For the random variables X1, …, Xn of sample size τ λ λ n drawn from the distribution of a random variable X 3 = 3 / 2 and µ σ2 ≤ ≤ with the mean and variance of , let X1:1 … τ λ λ X1:n be the order statistics such that the L-moments 4 = 4 / 2 . of X are defined by , An alternative measure of skewness, τ 3, is defined τ τ r-1 as (1 + 3) / (1 - 3). This measure is the ratio of the expected length of the upper tail to that of the lower -1 r - 1 tail in samples of size 3, and as such may be easier r k λr ≡ ∑ (-1) ( ) EX r = 1,2, … , r-k:r, to interpret than τ . λ , τ , τ and λ are subject to k 3 2 3 3 4 k=0 the constraints λ τ τ, ∞ where r is the rth L-moment of a distribution and 2 > 0 , -1 < 3 < 1 , 0 < 3 < and th EXi:r is the expected value of the i smallest 1/4(5τ2 - 1) ≤ τ < 1. observation in a sample of size r. 3 4 The first four central moments of a random variable If a random sample of size n is drawn from a X can be written as distribution of the random variable X and ≤ ≤ x1:n … xn:n are the ordered sample values then µ λ λ λ λ = E(X), estimates of the L-moments 1, 2, 3 and 4, namely I1, I2, I3, I4, can be calculated as follows. σ2 = E(X - µ)2 , First, define w2, w3, and w4 as ϒ = E(X - µ)3 /σ3 and 1 n w = ------- ∑ (i- 1)x , κ = E(X - µ)4 / σ4. 2 i:n n(n-1) i=2 In a similar fashion, the first four L-moments of a random variable X can be written as 1 n ∑ w3 = -------------- (i - 1)(i -2)xi:n and λ 1 = E(X), n(n-1)(n-2) i=3 λ 2 = 1/2E(X2:2 - X1:2), 1 n w = ---------------------------- ∑(i - 1)(i - 2)(i - 3)x λ 4 i:n 3 = 1/3E (X3:3 - 2X2:3 + X1:3) and n(n-1)(n-2)(n-3) i=4 λ 4 = 1/4E(X4:4 - 3X3:4 + 3X2:4 - X1:4) Then the L-moments and the corresponding shape λ It can be seen that 1 is equivalent to the usual statistics can be estimated as µ λ σ2 measure of central tendency, . 2 is similar to I = ∑ / n, in that both measure the difference between two 1 xi randomly selected values of X; however, by its I = 2w - I , nature σ2 assigns more weight to extreme sample 2 2 1 values than does λ . λ is a scale-dependent 2 3 I = 6w - 6w + I , λ 3 3 2 1 measure of skewness for a sample of size 3, and 4 is proportional to a weighted difference between I4 = 20w4 - 30w3 + 12w2 - I1, t3 = I3 / I2 and If BY variable processing is requested, the data are sorted before submission to PROC UNIVARIATE for analysis. PROC PRINTTO is used to capture the t4 = I4 / I2. usual output and send it to the file ‘UNI.DAT’. An output data set from PROC UNIVARIATE is used to t and t are the sample L-skewness and 3 4 store the number of non-missing observations for L-kurtosis, respectively, The sample estimate of the , each analysis. alternative measure of skewness, t 3, is defined as If the user chooses to generate a hardcopy of the (1 + t3) / (1 - t3). results, a DATA step is used to process the UNI.DAT file. The functions PUT and SUBSTR are The Program used in conjunction with the $HEX16. format to search for pagebreaks and set a flag that will be The macro program L_MOMENTS was originally ® used to fire a PUT _PAGE_ in a DATA_NULL_ step written using SAS v6.08 under the VMS operating at the end of the program. Next, a flag is set if BY environment. It has been modified to run under variable box plots are created. For each page of ® V6.09 and V6.12 on HP-UNIX . With slight output that does not contain BY variable box plots, a modification (detailed below), the program should counter is incremented. The counter is used to run on any operating system. The user is required facilitate direct read access of the shape statistics to provide the name of the SAS data set (macro data set created by PROC IML for use in the variable INDS) to be used in the analyses and the DATA_NULL_ that generates the hardcopy output. name(s) of the variables (macro variable VARS), Next, that part of the output line that displays the separated by spaces, to be analyzed. There is no values for skewness and kurtosis is removed. limit as to the number of variables that can be Finally, a flag is set that indicates the last line of the analyzed. Options available to the user include: tabular portion of the PROC UNIVARIATE output. • Specify the location of the SAS data library For each analysis variable, the raw data are sorted, (macro variable LIB). Default is current user then merged and transposed. The output data set location. from PROC UNIVARIATE that contains the number of non-missing observations is also transposed. • BY group processing (macro variable BYVAR). PROC IML is then used to calculate the central No limit of number of BY variables, delimited by moments and L-moments for skewness and a space.
Recommended publications
  • Use of Proc Iml to Calculate L-Moments for the Univariate Distributional Shape Parameters Skewness and Kurtosis
    Statistics 573 USE OF PROC IML TO CALCULATE L-MOMENTS FOR THE UNIVARIATE DISTRIBUTIONAL SHAPE PARAMETERS SKEWNESS AND KURTOSIS Michael A. Walega Berlex Laboratories, Wayne, New Jersey Introduction Exploratory data analysis statistics, such as those Gaussian. Bickel (1988) and Van Oer Laan and generated by the sp,ge procedure PROC Verdooren (1987) discuss the concept of robustness UNIVARIATE (1990), are useful tools to characterize and how it pertains to the assumption of normality. the underlying distribution of data prior to more rigorous statistical analyses. Assessment of the As discussed by Glass et al. (1972), incorrect distributional shape of data is usually accomplished conclusions may be reached when the normality by careful examination of the values of the third and assumption is not valid, especially when one-tail tests fourth central moments, skewness and kurtosis. are employed or the sample size or significance level However, when the sample size is small or the are very small. Hopkins and Weeks (1990) also underlying distribution is non-normal, the information discuss the effects of highly non-normal data on obtained from the sample skewness and kurtosis can hypothesis testing of variances. Thus, it is apparent be misleading. that examination of the skewness (departure from symmetry) and kurtosis (deviation from a normal One alternative to the central moment shape statistics curve) is an important component of exploratory data is the use of linear combinations of order statistics (L­ analyses. moments) to examine the distributional shape characteristics of data. L-moments have several Various methods to estimate skewness and kurtosis theoretical advantages over the central moment have been proposed (MacGillivray and Salanela, shape statistics: Characterization of a wider range of 1988).
    [Show full text]
  • On the Meaning and Use of Kurtosis
    Psychological Methods Copyright 1997 by the American Psychological Association, Inc. 1997, Vol. 2, No. 3,292-307 1082-989X/97/$3.00 On the Meaning and Use of Kurtosis Lawrence T. DeCarlo Fordham University For symmetric unimodal distributions, positive kurtosis indicates heavy tails and peakedness relative to the normal distribution, whereas negative kurtosis indicates light tails and flatness. Many textbooks, however, describe or illustrate kurtosis incompletely or incorrectly. In this article, kurtosis is illustrated with well-known distributions, and aspects of its interpretation and misinterpretation are discussed. The role of kurtosis in testing univariate and multivariate normality; as a measure of departures from normality; in issues of robustness, outliers, and bimodality; in generalized tests and estimators, as well as limitations of and alternatives to the kurtosis measure [32, are discussed. It is typically noted in introductory statistics standard deviation. The normal distribution has a kur- courses that distributions can be characterized in tosis of 3, and 132 - 3 is often used so that the refer- terms of central tendency, variability, and shape. With ence normal distribution has a kurtosis of zero (132 - respect to shape, virtually every textbook defines and 3 is sometimes denoted as Y2)- A sample counterpart illustrates skewness. On the other hand, another as- to 132 can be obtained by replacing the population pect of shape, which is kurtosis, is either not discussed moments with the sample moments, which gives or, worse yet, is often described or illustrated incor- rectly. Kurtosis is also frequently not reported in re- ~(X i -- S)4/n search articles, in spite of the fact that virtually every b2 (•(X i - ~')2/n)2' statistical package provides a measure of kurtosis.
    [Show full text]
  • The Normal Probability Distribution
    BIOSTATISTICS BIOL 4343 Assessing normality Not all continuous random variables are normally distributed. It is important to evaluate how well the data set seems to be adequately approximated by a normal distribution. In this section some statistical tools will be presented to check whether a given set of data is normally distributed. 1. Previous knowledge of the nature of the distribution Problem: A researcher working with sea stars needs to know if sea star size (length of radii) is normally distributed. What do we know about the size distributions of sea star populations? 1. Has previous work with this species of sea star shown them to be normally distributed? 2. Has previous work with a closely related species of seas star shown them to be normally distributed? 3. Has previous work with seas stars in general shown them to be normally distributed? If you can answer yes to any of the above questions and you do not have a reason to think your population should be different, you could reasonably assume that your population is also normally distributed and stop here. However, if any previous work has shown non-normal distribution of sea stars you had probably better use other techniques. 2. Construct charts For small- or moderate-sized data sets, the stem-and-leaf display and box-and- whisker plot will look symmetric. For large data sets, construct a histogram or polygon and see if the distribution bell-shaped or deviates grossly from a bell-shaped normal distribution. Look for skewness and asymmetry. Look for gaps in the distribution – intervals with no observations.
    [Show full text]
  • 2.8 Probability-Weighted Moments 33 2.8.1 Exact Variances and Covariances of Sample Probability Weighted Moments 35 2.9 Conclusions 36
    Durham E-Theses Probability distribution theory, generalisations and applications of l-moments Elamir, Elsayed Ali Habib How to cite: Elamir, Elsayed Ali Habib (2001) Probability distribution theory, generalisations and applications of l-moments, Durham theses, Durham University. Available at Durham E-Theses Online: http://etheses.dur.ac.uk/3987/ Use policy The full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that: • a full bibliographic reference is made to the original source • a link is made to the metadata record in Durham E-Theses • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders. Please consult the full Durham E-Theses policy for further details. Academic Support Oce, Durham University, University Oce, Old Elvet, Durham DH1 3HP e-mail: [email protected] Tel: +44 0191 334 6107 http://etheses.dur.ac.uk 2 Probability distribution theory, generalisations and applications of L-moments A thesis presented for the degree of Doctor of Philosophy at the University of Durham Elsayed Ali Habib Elamir Department of Mathematical Sciences, University of Durham, Durham, DHl 3LE. The copyright of this thesis rests with the author. No quotation from it should be published in any form, including Electronic and the Internet, without the author's prior written consent. All information derived from this thesis must be acknowledged appropriately.
    [Show full text]
  • Tests of Normality and Other Goodness-Of-Fit Tests
    Statistics 2 7 3 TESTS OF NORMALITY AND OTHER GOODNESS-OF-FIT TESTS Ralph B. D'Agostino, Albert J. Belanger, and Ralph B. D' Agostino Jr. Boston University Mathematics Department, Statistics and Consulting Unit Probability plots and goodness-of-fit tests are standardized deviates, :z;. on the horizontal axis. In table useful tools in detennining the underlying disnibution of I, we list the formulas for the disnibutions we plot in our a population (D'Agostino and Stephens, 1986. chapter 2). macro. If the underlying distribution is F(x), the resulting Probability plotting is an informal procedure for describing plot will be approximately a straight line. data and for identifying deviations from the hypothesized disnibution. Goodness-of-fit tests are formal procedures TABLE 1 which can be used to test for specific hypothesized disnibutions. We will present macros using SAS for Plotting Fonnulas for the six distributions plotted in our macro. creating probability plots for six disnibutions: the uniform. (PF(i-0.5)/n) normal, lognonnai, logistic, Weibull. and exponential. In Distribution cdf F(x) Vertical Horizontal addition these macros will compute the skewness <.pi;> . Axis Axis 7, kurtosis (b2), and D'Agostino-Pearson Kz statistics for testing if the underlying disttibution is normal (or Unifonn x-11 for 1'<%<1' +o ~ P... ~ ·lognormal) and the Anderson-Darling (EDF) A2 statistic G II for testing for the normal, log-normal, and exponential disnibutions. The latter can be modified for general ··I( ~;j i-3/8) distributions. 11+1(4 PROBABILITY PLOTS In(~~ ·-1( 1-3/8) 11+1{4 Say we desire to investigate if the underlying cumulative disnibution of a population is F(x) where this Weibull 1-exp( -(.!.t) ln(lw) In( -In( 1-p,)) disnibution depends upon a location parameter )l and scale 6 plll31lleter a.
    [Show full text]
  • Univariate Analysis and Normality Test Using SAS, STATA, and SPSS
    © 2002-2006 The Trustees of Indiana University Univariate Analysis and Normality Test: 1 Univariate Analysis and Normality Test Using SAS, STATA, and SPSS Hun Myoung Park This document summarizes graphical and numerical methods for univariate analysis and normality test, and illustrates how to test normality using SAS 9.1, STATA 9.2 SE, and SPSS 14.0. 1. Introduction 2. Graphical Methods 3. Numerical Methods 4. Testing Normality Using SAS 5. Testing Normality Using STATA 6. Testing Normality Using SPSS 7. Conclusion 1. Introduction Descriptive statistics provide important information about variables. Mean, median, and mode measure the central tendency of a variable. Measures of dispersion include variance, standard deviation, range, and interquantile range (IQR). Researchers may draw a histogram, a stem-and-leaf plot, or a box plot to see how a variable is distributed. Statistical methods are based on various underlying assumptions. One common assumption is that a random variable is normally distributed. In many statistical analyses, normality is often conveniently assumed without any empirical evidence or test. But normality is critical in many statistical methods. When this assumption is violated, interpretation and inference may not be reliable or valid. Figure 1. Comparing the Standard Normal and a Bimodal Probability Distributions Standard Normal Distribution Bimodal Distribution .4 .4 .3 .3 .2 .2 .1 .1 0 0 -5 -3 -1 1 3 5 -5 -3 -1 1 3 5 T-test and ANOVA (Analysis of Variance) compare group means, assuming variables follow normal probability distributions. Otherwise, these methods do not make much http://www.indiana.edu/~statmath © 2002-2006 The Trustees of Indiana University Univariate Analysis and Normality Test: 2 sense.
    [Show full text]
  • Typical Distribution Functions in Geophysics, Hydrology and Water Resources
    Chapter 6 Typical distribution functions in geophysics, hydrology and water resources Demetris Koutsoyiannis Department of Water Resources and Environmental Engineering Faculty of Civil Engineering, National Technical University of Athens, Greece Summary In this chapter we describe four families of distribution functions that are used in geophysical and engineering applications, including engineering hydrology and water resources technology. The first includes the normal distribution and the distributions derived from this by the logarithmic transformation. The second is the gamma family and related distributions that includes the exponential distribution, the two- and three-parameter gamma distributions, the Log-Pearson III distribution derived from the last one by the logarithmic transformation and the beta distribution that is closely related to the gamma distribution. The third is the Pareto distribution, which in the last years tends to become popular due to its long tail that seems to be in accordance with natural behaviours. The fourth family includes the extreme value distributions represented by the generalized extreme value distributions of maxima and minima, special cases of which are the Gumbel and the Weibull distributions. 5.1 Normal Distribution and related transformations 5.1.1 Normal (Gaussian) Distribution In the preceding chapters we have discussed extensively and in detail the normal distribution and its use in statistics and in engineering applications. Specifically, the normal distribution has been introduced in section 2.8, as a consequence of the central limit theorem, along with two closely related distributions, the χ2 and the Student (or t), which are of great importance in statistical estimates, even though they are not used for the description of geophysical variables.
    [Show full text]
  • Inference for Unreplicated Factorial and Fractional Factorial Designs. William J
    Louisiana State University LSU Digital Commons LSU Historical Dissertations and Theses Graduate School 1994 Inference for Unreplicated Factorial and Fractional Factorial Designs. William J. Kasperski Louisiana State University and Agricultural & Mechanical College Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_disstheses Recommended Citation Kasperski, William J., "Inference for Unreplicated Factorial and Fractional Factorial Designs." (1994). LSU Historical Dissertations and Theses. 5731. https://digitalcommons.lsu.edu/gradschool_disstheses/5731 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Historical Dissertations and Theses by an authorized administrator of LSU Digital Commons. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand corner and continuing from left to right in equal sections with small overlaps.
    [Show full text]
  • 5 Checking Assumptions
    5 CHECKING ASSUMPTIONS 5 Checking Assumptions Almost all statistical methods make assumptions about the data collection process and the shape of the population distribution. If you reject the null hypothesis in a test, then a reasonable conclusion is that the null hypothesis is false, provided all the distributional assumptions made by the test are satisfied. If the assumptions are not satisfied then that alone might be the cause of rejecting H0. Additionally, if you fail to reject H0, that could be caused solely by failure to satisfy assumptions also. Hence, you should always check assumptions to the best of your abilities. Two assumptions that underly the tests and CI procedures that I have discussed are that the data are a random sample, and that the population frequency curve is normal. For the pooled variance two sample test the population variances are also required to be equal. The random sample assumption can often be assessed from an understanding of the data col- lection process. Unfortunately, there are few general tests for checking this assumption. I have described exploratory (mostly visual) methods to assess the normality and equal variance assump- tion. I will now discuss formal methods to assess these assumptions. Testing Normality A formal test of normality can be based on a normal scores plot, sometimes called a rankit plot or a normal probability plot or a normal Q-Q plot. You plot the data against the normal scores, or expected normal order statistics (in a standard normal) for a sample with the given number of observations. The normality assumption is plausible if the plot is fairly linear.
    [Show full text]
  • Probability Plots for Exploratory Data Analysis
    SESUG Paper 172-2019 Probability Plots for Exploratory Data Analysis Dennis J. Beal, Leidos, Oak Ridge, Tennessee ABSTRACT Probability plots are used in statistical analysis to check distributional assumptions, visually check for potential outliers, and see the range, median, and variability of a data set. Probability plots are an important statistical tool to use for exploratory data analysis. This paper shows SAS® code that generates normal and lognormal probability plots using the Output Delivery System (ODS) on a real environmental data set using PROC UNIVARIATE and interprets the results. This paper is for beginning or intermediate SAS users of Base SAS® and SAS/GRAPH®. Key words: normal probability plots, lognormal probability plots, PROC UNIVARIATE INTRODUCTION When the data analyst receives a data set, the first step to understanding the data is to conduct an exploratory data analysis. This first step might include a PROC CONTENTS to list the variables that are in the data set and their data types (numeric, character, date, time, etc.). The next step might be to run PROC FREQ to see the values of all the variables in the data set. Another common practice is to run PROC UNIVARIATE to see summary statistics for all numeric variables. These steps help guide the data analyst to see whether the data needs cleaning or standardizing before conducting a statistical analysis. In addition to summary statistics, PROC UNIVARIATE also generates histograms and normal probability plots for numeric variables. In SAS v. 9.3 and 9.4 PROC UNIVARIATE generates generic histograms and normal probability plots for the numeric variables with the ODS.
    [Show full text]
  • Normality Tests
    NCSS Statistical Software NCSS.com Chapter 194 Normality Tests Introduction This procedure provides seven tests of data normality. If the variable is normally distributed, you can use parametric statistics that are based on this assumption. If a variable fails a normality test, it is critical to look at the histogram and the normal probability plot to see if an outlier or a small subset of outliers has caused the non-normality. If there are no outliers, you might try a transformation (such as, the log or square root) to make the data normal. If a transformation is not a viable alternative, nonparametric methods that do not require normality may be used. Always remember that a reasonably large sample size is required to detect departures from normality. Only extreme types of non-normality can be detected with samples less than fifty observations. There is a common misconception that a histogram is always a valid graphical tool for assessing normality. Since there are many subjective choices that must be made in constructing a histogram, and since histograms generally need large sample sizes to display an accurate picture of normality, preference should be given to other graphical displays such as the box plot, the density trace, and the normal probability plot. Normality tests generally have small statistical power (probability of detecting non-normal data) unless the sample sizes are at least over 100. Technical Details This section provides details of the seven normality tests that are available. Shapiro-Wilk W Test This test for normality has been found to be the most powerful test in most situations.
    [Show full text]
  • Analysis of the Behaviour of Stocks of Dar Es Salaam Stock Exchange (DSE)
    Mathematical Theory and Modeling www.iiste.org ISSN 2224-5804 (Paper) ISSN 2225-0522 (Online) Vol.4, No.10, 2014 Analysis of the Behaviour of Stocks of Dar es Salaam Stock Exchange (DSE) Phares Kaboneka1*, Wilson C. Mahera1,2, Silas Mirau1 1 School of Computational and Communication Science and Engineering, Nelson Mandela African institution of science and technology PO box 447, Arusha, Tanzania 2 Department of Mathematics, Dar es Salaam University Tel: +255752456180 *E-mail: [email protected] The research is financed by the Nelson Mandela African Institution of Science and Technology. Abstract A stock market is a place where investors trade certificates that indicate partial ownership in businesses for a set price. Different countries in the world have stock markets where other countries started their stock markets long time ago like the USA and they have investigated the trend of their market if it is normally distributed or not. Also they have strong models that assist them in making predictions and also help the investors on the choice of the stocks to invest so as to gain the profit in the future. On the other hand other countries just started few years ago. Tanzania is among the countries where stock markets has just started recently and hence there is a need to study the nature of the stocks distribution and see whether the Dar-Es-Salaam Stock of Exchange (DSE) market do follow the theoretical conclusions or not. Thus in this study we adapt the Markowitz modern portfolio theory (MPT) and using the mean variance analysis theory together with the DSE data to investigate if the DSE stock market follows a normal distribution or not.
    [Show full text]