Estimation of Quantiles and the Interquartile Range in Complex Surveys Carol Ann Francisco Iowa State University

Total Page:16

File Type:pdf, Size:1020Kb

Estimation of Quantiles and the Interquartile Range in Complex Surveys Carol Ann Francisco Iowa State University Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1987 Estimation of quantiles and the interquartile range in complex surveys Carol Ann Francisco Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Statistics and Probability Commons Recommended Citation Francisco, Carol Ann, "Estimation of quantiles and the interquartile range in complex surveys " (1987). Retrospective Theses and Dissertations. 11680. https://lib.dr.iastate.edu/rtd/11680 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. INFORMATION TO USERS While the most advanced technology has been used to photograph and reproduce this manuscript, the quality of the reproduction is heavily dependent upon the quality of the material submitted. For example: • Manuscript pages may have indistinct print. In such cases, the best available copy has been filmed. • Manuscripts may not always be complete. In such cases, a note will indicate that it is not possible to obtain missing pages. • Copyrighted material may have been removed from the manuscript. In such cases, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, and charts) are photographed by sectioning the original, beginning at the upper left-hand corner and continuing from left to right in equal sections with small overlaps. Each oversize page is also filmed as one exposure and is available, for an additional charge, as a standard 35mm slide or as a 17"x 23" black and white photographic print. Most photographs reproduce acceptably on positive microfilm or microfiche but lack the clarity on xerographic copies made from the microfilm. For an additional charge, 35mm slides of 6"x 9" black and white photographic prints are available for any photographs or illustrations that cannot be reproduced satisfactorily by xerography. Order Number 8721882 Estimation of quantiles and the interquartile range in complex surveys Francisco, Carol Ann, Ph.D. Iowa State University, 1987 UMI SOON.ZeebRd. Ann Arbor, MI 48106 Estimation of quantiles and the interquartile range in complex surveys by Carol Ann Francisco A Dissertation Submitted to the Graduate Faculty in Partial Fulfillment of the Requirements for the Degree of DOCTOR OF PHILOSOPHY Major: Statistics Approved: Signature was redacted for privacy. In Chaise of^Majot Work Signature was redacted for privacy. For the Major Department Signature was redacted for privacy. For the Grad e College Iowa State University Ames, Iowa 1987 il TABLE OF CONTENTS Page I. INTRODUCTION 1 II. PREVIOUS WORK IN QUANTILE ESTIMATION 8 A. Simple Random Sampling 8 1. Infinite Population Sampling 8 2. Finite Population Sampling 27 B. Stratified Random Sampling 46 C. Stratified Cluster Sampling 60 III. ESTIMATION OF QUANTILES AND THE INTERQUARTILE RANGE IN LARGE-SCALE SURVEYS 71 A. Properties of the Sample Mean 75 B. Properties of the Empirical Distribution Function 102 C. Confidence Intervals for Quantlles and the Interquartile Range 106 1. Confidence Intervals Based on Sample Quantlles 107 2. Quantile Confidence Intervals Based on Test Inversion 138 D. Extensions to More Complex Designs 156 IV. IMPLEMENTATION OF THE PROPOSED ESTIMATION PROCEDURES 158 A. Conputer Algorithm 159 B. Monte Carlo Simulations 167 1. Percentage of Urban Land Population 167 2. Soil Loss Population 171 3. Generated Lognormal Population 175 V. SUMMARY AND CONCLUSIONS 183 ill Page VI. REFERENCES 187 VII. ACKNOWLEDOIENTS 193 I. IHTBODUCTION Researchers are often interested in estimating cumulative distribution functions from survey data. For example, Sedransk and Sedransk (1979) examined the feasibility of using estimated cumulative distribution functions as one means of making comparisons among radi­ ation therapy facilities in a large-scale national survey of cancer patient medical records. Also of interest in many analytical surveys are estimators of functions of the cumulative distribution function such as the median and other quantiles, along with estimators of their standard errors and confidence intervals. Quantiles are of particular Interest when a variable is known to have a markedly skewed distribu­ tion. For example, in a survey of urban households, a point or interval estimate of median household Income may be desired. Studies which assume simple random sampling without replacement from an absolutely continuous distribution form the basis of most of the existing literature on quantile estimation. Estimators which have been proposed for this situation include the sample median and various linear combinations of order statistics (Karrell and Davis 1982, Kaigh and Lachenbruch 1982). Estimators of the standard error of the sample median which have been proposed Include a small sample variance estimator by Marltz and Jarrett (1978), a method which is essentially equivalent to the bootstrap "Method 1" estimator described by Efron (1979), and a generalized least squares estimator of the large sample standard error of the sample median (Sheather and Marltz 1983). 2 Procedures for constructing confidence intervals for the median under simple random sampling from a continuous distribution date back to the late 1930s (Nair 1940, Savur 1937, Thonçson 1936) and since then have been thoroughly investigated (David 1981). Sampling designs which are more complex than simple random sampling are generally used in survey sampling. Some information about the population of interest is usually available and often is Incorporated into the survey design and into the estimation process to Increase the precision of estimates. Extension of results derived under an assumption of simple random sampling from an absolutely continuous distribution to the complex sampling designs used in finite population sampling has met with limited success. Woodruff (1952) and Jonrup (1975) have proposed using a weighted sample median to estimate the population median, where the weight assigned to each observation is proportional to the inverse of its selection probability. Using the approach taken by Maritz and Jarrett (1978), Gross (1980) derived a small-sample estimator of the variance of the weighted sample median estimator for stratified sampling without replacement from a finite population. Current methods for constructing confidence intervals for quantiles vary depending upon the approach taken relative to making inferences from samples to the finite population. Broadly speaking, the finite population can be considered to be a fixed set of elements, or it can be considered to come from an infinite population, called a superpopula­ tion. 3 Under the classical or fixed population approach (Cassel, Sarndal, and Wretman 1977), each population value is regarded as a fixed but unknown quantity. The sample is viewed as random, and the units within the sample are considered fixed quantities. Inference from the sample to the finite population is based primarily upon the stochastic element introduced by the sampling plan. The superpopulation approach considers the observed values in a survey to be outcomes of an underlying process or superpopulation model. IMlike the classical approach, the superpopulation approach takes the view that associated with each population unit is a random variable which has a given stochastic structure. The superpopulation model plays a vital role in inference with this approach. From the observations in a sample, inferences are made about the superpopulation model. The superpopulation model then is used to make predictions concerning the population units not included in the sample and, ultimately, the population parameter of interest. A number of authors have investigated model-free procedures for constructing exact 100(1 - a) percent confidence intervals for quantiles in finite populations (0 < a < 1) . Inferences from the sang)le to the finite population are based upon confidence Intervals which take into account the sampling scheme. Thompson (1936), Wilks (1962), and Konijn (1973) have given design-based confidence intervals for the sample median when simple random sampling from a finite population is assumed. Meyer (1972) and Sedransk and Meyer (1978) investigated three exact confidence interval procedures for quantiles 4 when sampling is from a stratified population. Stratified random sang)ling when the population has been divided into two strata either randomly or in an ordered fashion (based on the known distribution of an increasing concomitant variable) was assumed. Confidence intervals for two of the proposed procedures are formed from pairs of order statistcs; a third procedure is derived from the empirical distribution function. For more than two strata, the confidence Interval procedures proposed by Meyer (1972) and Sedransk and Meyer (1978) become very complex and require substantial amounts of computation to implement. Theoretical results are tedious but can be extended to more than two strata. Meyer (1972) has given results for the case of three strata. Blesseos (1976) has extended Meyer's procedures for determining the confidence coefficient of confidence Intervals based on order statistics from the combined sample to the case where the finite population has up to five ordered or
Recommended publications
  • Lesson 7: Measuring Variability for Skewed Distributions (Interquartile Range)
    NYS COMMON CORE MATHEMATICS CURRICULUM Lesson 7 M2 ALGEBRA I Lesson 7: Measuring Variability for Skewed Distributions (Interquartile Range) Student Outcomes . Students explain why a median is a better description of a typical value for a skewed distribution. Students calculate the 5-number summary of a data set. Students construct a box plot based on the 5-number summary and calculate the interquartile range (IQR). Students interpret the IQR as a description of variability in the data. Students identify outliers in a data distribution. Lesson Notes Distributions that are not symmetrical pose some challenges in students’ thinking about center and variability. The observation that the distribution is not symmetrical is straightforward. The difficult part is to select a measure of center and a measure of variability around that center. In Lesson 3, students learned that, because the mean can be affected by unusual values in the data set, the median is a better description of a typical data value for a skewed distribution. This lesson addresses what measure of variability is appropriate for a skewed data distribution. Students construct a box plot of the data using the 5-number summary and describe variability using the interquartile range. Classwork Exploratory Challenge 1/Exercises 1–3 (10 minutes): Skewed Data and Their Measure of Center Verbally introduce the data set as described in the introductory paragraph and dot plot shown below. Exploratory Challenge 1/Exercises 1–3: Skewed Data and Their Measure of Center Consider the following scenario. A television game show, Fact or Fiction, was cancelled after nine shows. Many people watched the nine shows and were rather upset when it was taken off the air.
    [Show full text]
  • Normal Probability Distribution
    MATH2114 Week 8, Thursday The Normal Distribution Slides Adopted & modified from “Business Statistics: A First Course” 7th Ed, by Levine, Szabat and Stephan, 2016 Pearson Education Inc. Objective/Goals To compute probabilities from the normal distribution How to use the normal distribution to solve practical problems To use the normal probability plot to determine whether a set of data is approximately normally distributed Continuous Probability Distributions A continuous variable is a variable that can assume any value on a continuum (can assume an uncountable number of values) thickness of an item time required to complete a task temperature of a solution height, in inches These can potentially take on any value depending only on the ability to precisely and accurately measure The Normal Distribution Bell Shaped f(X) Symmetrical Mean=Median=Mode σ Location is determined by X the mean, μ μ Spread is determined by the standard deviation, σ Mean The random variable has an = Median infinite theoretical range: = Mode + to The Normal Distribution Density Function The formula for the normal probability density function is 2 1 (X μ) 1 f(X) e 2 2π Where e = the mathematical constant approximated by 2.71828 π = the mathematical constant approximated by 3.14159 μ = the population mean σ = the population standard deviation X = any value of the continuous variable The Normal Distribution Density Function By varying 휇 and 휎, we obtain different normal distributions Image taken from Wikipedia, released under public domain The
    [Show full text]
  • 2018 Statistics Final Exam Academic Success Programme Columbia University Instructor: Mark England
    2018 Statistics Final Exam Academic Success Programme Columbia University Instructor: Mark England Name: UNI: You have two hours to finish this exam. § Calculators are necessary for this exam. § Always show your working and thought process throughout so that even if you don’t get the answer totally correct I can give you as many marks as possible. § It has been a pleasure teaching you all. Try your best and good luck! Page 1 Question 1 True or False for the following statements? Warning, in this question (and only this question) you will get minus points for wrong answers, so it is not worth guessing randomly. For this question you do not need to explain your answer. a) Correlation values can only be between 0 and 1. False. Correlation values are between -1 and 1. b) The Student’s t distribution gives a lower percentage of extreme events occurring compared to a normal distribution. False. The student’s t distribution gives a higher percentage of extreme events. c) The units of the standard deviation of � are the same as the units of �. True. d) For a normal distribution, roughly 95% of the data is contained within one standard deviation of the mean. False. It is within two standard deviations of the mean e) For a left skewed distribution, the mean is larger than the median. False. This is a right skewed distribution. f) If an event � is independent of �, then � (�) = �(�|�) is always true. True. g) For a given confidence level, increasing the sample size of a survey should increase the margin of error.
    [Show full text]
  • L-Moments Based Assessment of a Mixture Model for Frequency Analysis of Rainfall Extremes
    Advances in Geosciences, 2, 331–334, 2006 SRef-ID: 1680-7359/adgeo/2006-2-331 Advances in European Geosciences Union Geosciences © 2006 Author(s). This work is licensed under a Creative Commons License. L-moments based assessment of a mixture model for frequency analysis of rainfall extremes V. Tartaglia, E. Caporali, E. Cavigli, and A. Moro Dipartimento di Ingegneria Civile, Universita` degli Studi di Firenze, Italy Received: 2 December 2004 – Revised: 2 October 2005 – Accepted: 3 October 2005 – Published: 23 January 2006 Abstract. In the framework of the regional analysis of hy- events in Tuscany (central Italy), a statistical methodology drological extreme events, a probabilistic mixture model, for parameter estimation of a probabilistic mixture model is given by a convex combination of two Gumbel distributions, here developed. The use of a probabilistic mixture model in is proposed for rainfall extremes modelling. The papers con- the regional analysis of rainfall extreme events, comes from cerns with statistical methodology for parameter estimation. the hypothesis that rainfall phenomenon arises from two in- With this aim, theoretical expressions of the L-moments of dependent populations (Becchi et al., 2000; Caporali et al., the mixture model are here defined. The expressions have 2002). The first population describes the basic components been tested on the time series of the annual maximum height corresponding to the more frequent events of low magnitude, of daily rainfall recorded in Tuscany (central Italy). controlled by the local morpho-climatic conditions. The sec- ond population characterizes the outlier component of more intense and rare events, which is here considered due to the aggregation at the large scale of meteorological phenomena.
    [Show full text]
  • Hand-Book on STATISTICAL DISTRIBUTIONS for Experimentalists
    Internal Report SUF–PFY/96–01 Stockholm, 11 December 1996 1st revision, 31 October 1998 last modification 10 September 2007 Hand-book on STATISTICAL DISTRIBUTIONS for experimentalists by Christian Walck Particle Physics Group Fysikum University of Stockholm (e-mail: [email protected]) Contents 1 Introduction 1 1.1 Random Number Generation .............................. 1 2 Probability Density Functions 3 2.1 Introduction ........................................ 3 2.2 Moments ......................................... 3 2.2.1 Errors of Moments ................................ 4 2.3 Characteristic Function ................................. 4 2.4 Probability Generating Function ............................ 5 2.5 Cumulants ......................................... 6 2.6 Random Number Generation .............................. 7 2.6.1 Cumulative Technique .............................. 7 2.6.2 Accept-Reject technique ............................. 7 2.6.3 Composition Techniques ............................. 8 2.7 Multivariate Distributions ................................ 9 2.7.1 Multivariate Moments .............................. 9 2.7.2 Errors of Bivariate Moments .......................... 9 2.7.3 Joint Characteristic Function .......................... 10 2.7.4 Random Number Generation .......................... 11 3 Bernoulli Distribution 12 3.1 Introduction ........................................ 12 3.2 Relation to Other Distributions ............................. 12 4 Beta distribution 13 4.1 Introduction .......................................
    [Show full text]
  • Descriptive Statistics
    NCSS Statistical Software NCSS.com Chapter 200 Descriptive Statistics Introduction This procedure summarizes variables both statistically and graphically. Information about the location (center), spread (variability), and distribution is provided. The procedure provides a large variety of statistical information about a single variable. Kinds of Research Questions The use of this module for a single variable is generally appropriate for one of four purposes: numerical summary, data screening, outlier identification (which sometimes is incorporated into data screening), and distributional shape. We will briefly discuss each of these now. Numerical Descriptors The numerical descriptors of a sample are called statistics. These statistics may be categorized as location, spread, shape indicators, percentiles, and interval estimates. Location or Central Tendency One of the first impressions that we like to get from a variable is its general location. You might think of this as the center of the variable on the number line. The average (mean) is a common measure of location. When investigating the center of a variable, the main descriptors are the mean, median, mode, and the trimmed mean. Other averages, such as the geometric and harmonic mean, have specialized uses. We will now briefly compare these measures. If the data come from the normal distribution, the mean, median, mode, and the trimmed mean are all equal. If the mean and median are very different, most likely there are outliers in the data or the distribution is skewed. If this is the case, the median is probably a better measure of location. The mean is very sensitive to extreme values and can be seriously contaminated by just one observation.
    [Show full text]
  • Measures of Variability
    CHAPTER 4 Measures of Variability Chapter 3 introduced descriptive statistics, the purpose of which is to numerically describe or summarize data for a Chapter Outline variable. As we learned in Chapter 3, a set of data is often 4.1 An Example From the described in terms of the most typical, common, or frequently Research: How Many “Sometimes” occurring score by using a measure of central tendency such in an “Always”? as the mode, median, or mean. Measures of central ten- 4.2 Understanding Variability dency numerically describe two aspects of a distribution of scores, modality and symmetry, based ondistribute what the scores 4.3 The Range have in common with each other. However, researchers • Strengths and weaknesses of also describe the amount of differences among the scores of a the range variable in analyzing and reportingor their data. This chapter 4.4 The Interquartile Range introduces measures of variability, which are designed to rep- • Strengths and weaknesses of resent the amount of differences in a distribution of data. In the interquartile range this chapter, we will present and evaluate different measures 4.5 The Variance (s2) of variability in terms of their relative strengths and weak- • Definitional formula for the nesses. As in the previous chapters, our presentation will variance depend heavily on the findings and analysis from a published • Computational formula for the research study.post, variance • Why not use the absolute value of the deviation in 4.1 An Example From the calculating the variance? Research: How Many • Why divide by N – 1 rather than “Sometimes” in an “Always”? N in calculating the variance? 4.6 The Standard Deviation (s) copy,Throughout our daily lives, we are constantly asked to fill out • Definitional formula for the questionnaires and surveys.
    [Show full text]
  • Measures of Spread
    community project encouraging academics to share statistics support resources All stcp resources are released under a Creative Commons licence stcp-gilchristsamuels-6 The following resources are associated: Descriptive Statistics – Measures of Middle Value Descriptive Statistics – Measures of Spread Research question type: Most What kind of variables: Nominal, ordinal and interval/scale Common Applications: Most quantitative studies Descriptive statistics are a way of summarising your data into a few values, often for comparative purposes. Commonly they are presented in tables, although appropriate charts are also used. They include measures of middle values (or central tendency), and measures of spread (or dispersion) of your data. The following is a rough guide to appropriate measures for level of data: Data level Middle value (central tendency) Dispersion (spread) Nominal Mode Ordinal Mode, median Range, interquartile range (IQR) Interval/scale Median, mean Variance, standard deviation Range The range is defined as the difference between the maximum and the minimum values of a set of data. Example Find the range for: 2, 6, 3, 9, 5, 6, 2, 6. The maximum value of the data is 9, and the minimum value is 2. Hence the range is 9 - 2 = 7. Interquartile range (IQR) Quartiles split a dataset into four quarters when the values are written in ascending order. !!! The lower quartile is 25% of the way through a data set. This is the th value, where n is the ! number of data values. (Note: Quarter values are found by starting with the lower value and adding either ¼ or ¾ of the difference. For example, the 2¼th value is the 2nd value plus ¼ of the 3rd value minus the 2nd).
    [Show full text]
  • Chapter 4 – Analyzing Skewed Quantitative Data Introduction: in Chapter 3, We Focused on Analyzing Bell Shaped (Normal) Data, but Many Data Sets Are Not Bell Shaped
    Chapter 4 – Analyzing Skewed Quantitative Data Introduction: In chapter 3, we focused on analyzing bell shaped (normal) data, but many data sets are not bell shaped. How do we analyze quantitative data when it is not bell shaped? The main issue is that the mean and standard deviations are not accurate and should not be used in the analysis. Then what statistics should we use? We will be introducing a new kind of graph that is specially designed for analyzing skewed data. It is called the “box and whisker plot” or “box plot” for short. When data sets are not bell shaped, we will focus on the median, quartiles, interquartile range and boxplots to measure center and spread. Quartiles are more accurate because they are based on the order of the numbers instead of distances and so are not as effected by the skewed shape and extremely unusual values. 143 Section 4A – Review of Shapes and Centers, Creating Histograms and Dot Plots with Technology Let us start by reviewing how to create dot plots and histograms with technology and determine the shape of a quantitative data set. Here are the directions for making dot plots and histograms in Statcato. Making a dot plot in Statcato: Graph => Dot plot => Pick a column => OK Making a histogram in Statcato: Graph => Histogram => Pick a column => Chose number of bins => OK Example 1 Use the women’s weight data (in pounds) from the Health data and Statcato to create a dot plot, histogram and boxplot. What is the shape of the data set? Plugging in the women’s weight data into Statcato, we can use the directions above to create the following graphs.
    [Show full text]
  • Outing the Outliers – Tails of the Unexpected
    Presented at the 2016 International Training Symposium: www.iceaaonline.com/bristol2016 Outing the Outliers – Tails of the Unexpected Outing the Outliers – Tails of the Unexpected Cynics would say that we can prove anything we want with statistics as it is all down to … a word (or two) from the Wise? interpretation and misinterpretation. To some extent this is true, but in truth it is all down to “The great tragedy of Science: the slaying making a judgement call based on the balance of of a beautiful hypothesis by an ugly fact” probabilities. Thomas Henry Huxley (1825-1895) British Biologist The problem with using random samples in estimating is that for a small sample size, the values could be on the “extreme” side, relatively speaking – like throwing double six or double one with a pair of dice, and nothing else for three of four turns. The more random samples we have the less likely (statistically speaking) we are to have all extreme values. So, more is better if we can get it, but sometimes it is a question of “We would if we could, but we can’t so we don’t!” Now we could argue that Estimators don’t use random values (because it sounds like we’re just guessing); we base our estimates on the “actuals” we have collected for similar tasks or activities. However, in the context of estimating, any “actual” data is in effect random because the circumstances that created those “actuals” were all influenced by a myriad of random factors. Anyone who has ever looked at the “actuals” for a repetitive task will know that there are variations in those values.
    [Show full text]
  • Problem Max. Points Your Points 1-10 10 11 10 12 3 13 4 14 18 15 8 16 7 17 14 Total 75
    M 140 Test 1 B Name__________________ (1 point) SHOW YOUR WORK FOR FULL CREDIT! Problem Max. Points Your Points 1-10 10 11 10 12 3 13 4 14 18 15 8 16 7 17 14 Total 75 Multiple choice questions (1 point each) For questions 1-3 consider the following: The distribution of scores on a certain statistics test is strongly skewed to the right. 1. How would the mean and the median compare for this distribution? a) mean < median b) mean = median c) mean > median d) Can’t tell without more information When a distribution skewed to the right, the “tail” is on the right. Since the tail “pulls” the mean, the mean is greater than the median. 2. Which set of measures of center and spread are more appropriate for the distribution of scores? a) Mean and standard deviation b) Median and interquartile range c) Mean and interquartile range d) Median and standard deviation For skewed distributions the median and the interquartile range are the appropriate measures of center and spread because they are resistant to outliers. 3. What does this suggest about the difficulty of the test? a) It was an easy test b) It was a hard test c) It wasn’t too hard or too easy d) It is impossible to tell If the distribution is skewed to the right, then most of the values are to the left, toward the lower values. Thus, more students got lower scores, which suggests that the test was hard. 4. Which one of the following variables is NOT categorical? a) whether or not an individual has glasses b) weight of an elephant c) religion d) social security number All the other variables are categorical—the average doesn’t make sense for those.
    [Show full text]
  • For Questions 1-4, Consider the Following Data Set: 4 3 6 2 8 8
    Math 225 Quiz 2 Name:_________________ For questions 1-4, consider the following data set: 4 3 6 2 8 8 1. What is the median of the data set? a. 4 b. 5 c. 2 d. 4.5 2. If I change one of the 8’s in the data set to 18, how would that change the median? a. The median would be higher b. The median would be lower c. The median would remain the same 3. If I change one of the 8’s in the data set to 18, how would that change the mean? a. The mean would be higher b. The mean would be lower c. The mean would remain the same 4. What is the mode of the data set? a. there is no mode b. 8 c. 2 5. For the distribution in the picture a. mean ≈ median b. mean < median c. mean > median d. Can’t tell the relationship of the mean and the median without looking at the data. 6. Which one of these is NOT true about the standard deviation? a. The standard deviation is zero only if the observations are all the same. b. The standard deviation has the same unit as the data. c. The standard deviation is used as a measure of spread when the mean is use as the measure of center. d. The standard deviation is resistant to outliers. 1 7. What percent of the data lies between the first quartile, Q 1, and the Median? a. About 50% b. About 25% c.
    [Show full text]