Handling Missing Data on Asymmetric Distribution
|
|
- Kevin Henry
- 5 years ago
- Views:
Transcription
1 International Matematical Forum, Vol. 8, 03, no. 4, Handling Missing Data on Asymmetric Distribution Amad M. H. Al-Kazale Department of Matematics, Faculty of Science Al-albayt University, Al-Mafraq-Jordan Abstract Te problem of imputation of missing observations emerges in many areas. Data usually contained missing observations due to many factors, suc as macine failures and uman error. Incomplete dataset usually causes bias due to differences between observed and unobserved data. Tis paper proposed Neyman allocation metod to estimate asymmetric winsorizing mean for andling missing observations wen te data follow te exponential distribution. Different values of te exponential distribution parameters were used to illustrate. A set of data from exponential distribution were generated to compare te performance of te proposed metods suc as regression trend, average of te wole data, naive forecast and average bound of te oles and te proposed Neyman allocation metod. Te goodness-of-fit criterions used were te mean absolute error (MAE) and te mean squared error (MSE). It was found tat te proposed metod gave te best fit in te sense of aving smaller error, in particular for a large percentage of missing observations. Keywords: Neyman Allocation, Missing Data, Single imputation metods, Winsorized mean. INTRODUCTION Te attendance of missing observations in statistical survey data is an important issue to deal wit (Little and Rubin, 987). Besides missing observations, te time series data are possibly contaminated by outliers or are eterogeneous. Incomplete datasets may lead to results tat are different from tose tat would ave been obtained from a complete dataset (Hawtorne and Elliot, 005).
2 54 A. M. H. Al-Kazale Tis paper proposed a metod of substituting tese missing observations by using asymmetric winsorized mean. In tis metod, we prepared te wole data by dividing it into groups by using te stratified sampling to determine te boundaries among te strata. Stratification is te process of grouping members of te population into relatively omogeneous groups before sampling. Te strata sould be mutually exclusive were every element in te population must be assigned to only one stratum (Cyert, and Davidson, 96). One of te main objectives of stratified sampling is to reduce te variance of te estimator and to get more statistical precision tan wit te simple random sampling (Cocran, 977; Hanif, 000; Hess et al., 966; Seik and Amad, 00). In tis paper, Neyman allocation metod is used to determine te boundaries of te groups. Te missing observations problem is an old one for analysis tasks. Te waste of te data wic can result from casewise deletion of missing values, oblige to propose alternative approaces. A problem frequently encounters in data collection is missing observations or observations may be virtually impossible to obtain, eiter because of time or cost constrains. In order to replace tose observations, tere are several different options available to te researcers (Pankratz, 983; Patricia, 994). Firstly, replace te missing observations wit te mean of te series. Secondly; replace te missing observations wit te naive forecast. Also; replace te missing observations wit a simple trend forecast. Finally replace te missing observations wit an average of te last two known observations tat bound te missing observations.. NEYMAN ALLOCATION Neyman allocation is a sample allocation metod tat may be used wit stratified samples. We used te Neyman allocation to determine te stratum boundaries, wen te population is igly skewed te Neyman allocation sould be used (Cocran,977; Amad Mair et al., 007; Samiuddin et al., 998). Te suffix denotes te stratum and i te unit witin te stratum. Wen a population of N units are being stratified into L strata and te samples from eac stratum are selected wit simple random sampling an unbiased estimate of te population mean for te estimation variable x, is x st (st for stratified), were x L N x = = = st N = L W x (.) For te stratified random sampling, te variance of te estimate x st is
3 Handling missing data on asymmetric distribution 55 L L N S W S N n n = = V ( x ) = = st (.) Te variance depend on Neyman allocation is reduced to L VNey ( xst ) = ( WS ) n (.3) = Let x 0, x L be te smallest and largest values of x in te population. Te problem is to find intermediate stratum boundaries y, y,..., y L and by differentiating wit respect to te stratum boundaries y. We get te minimum variance VNey ( x st ) wit respect to y, since y appears in te sum only in te terms of W, S and W+, S+ (Cocran,977). Hence we ave te formula for finding te optimum stratum boundaries. ( y μ) + S ( y μ + ) + S + = were =,,3,..., L (.6) S S + Suppose tat te distribution of x be continuous wit te density function, f ( x ), a < x < bin order to make L (groups) strata, te range of x is to be cut up at points y < y < y3 <... < y L. Te relative frequency W, te population mean μ and te population variance S at te -t stratum are given by. y W = f ( x) dx (.7) y y μ = x f ( x) dx W (.8) y y = ( ) μ W y S x f x dx (.9) Tese equations are difficult to be solved, since μ and S depend on y. We must use te iterative metod to solve tem by using computer program C++ (Amad Mair et al. 007). We divide te wole data into groups by putting ( y0, y ) into group one were y0 = x0, ( y, y ) in group two and so on. 3. WINSORIZED MEAN Winsorized mean is a winsorized statistical measurement of central tendency, like te mean and te median and even more similar to te truncated mean (Mingxin and
4 56 A. M. H. Al-Kazale Yijun, 009; Barnett and Lewis 994), te winsorized mean eliminates te outliers at bot ends of an ordered set of observations. Unlike te trimmed mean, te winsorized mean replaces te outliers wit observed values, rater tan discarding tem. It involves te calculation of te mean after replacing te given parts of a probability distribution at te ig and low ends wit te most extreme remaining values. Te winsorized mean is defined by n s xw = [( r + ) x r+ + xi ( s ) x n s ] n + + (3.) i = r+ wic means tat te winsorized mean is te average of te observations were te first r smallest values are replaced by te ( r + ) -t smallest value x r +, and te first s largest values are replaced by te ( s + ) -t largest value x n s (Barnett and Lewis 994). Wen, we determine te suitable number of groups (must be more tan two groups) depending on te wole data by using Equation (.6). We obtain te boundaries of tese groups at te first and te last groups. Consequently we count te number of observations inside tese two groups. Te number of observations in te first group is r and te number of observations in te last group is s. Subsequently, te asymmetrical winsorized mean can be calculated using Equation (3.). Te data were generated from exponential distribution. Ten, te wole data ordered from te smallest to te largest, te smallest value is zero and te largest value is 7.. Te data were assumed aving 5%, 0%, and 0% randomly missing observations in tis study. 4 THE BOUNDARIES OF GROUPS IN EXPONENTIAL DISTRIBUTION Suppose tat te time series x, x,..., x n of n observations come from an exponential distribution wit probability density function f ( x) = λe λx, λ > 0, x > 0, (4.) were λ is te parameter of exponential function. In order to make L (groups) strata, te domain of f ( x ) was truncated from, (0, ) to (0, b ] were b is te largest value of x. From equations (.7), (.8) and (.9), te relative frequency W, mean μ and variance S of te -t stratum can be computed for te exponential distribution as te following λy λy W = e e (4.)
5 Handling missing data on asymmetric distribution 57 S λx λx ( y + ) e ( y + ) e μ = λ λ λ y λ y e e ( y + y + ) e ( y + y + ) e λx λx = λ λ λ λ ( μ ) λy λy e e (4.3) (4.4) Table Stratum boundaries for groups size of tree to eigt in te case of exponential distribution wen λ = 0., Boundaries 3G 4G 5G 6G 7G 8G y y y y y y y 7.83 Table Stratum boundaries for groups size of tree to eigt in te case of exponential distribution wen λ = 0.4, Boundaries 3G 4G 5G 6G 7G 8G y y y y y y y Table 3 Stratum boundaries for groups size of tree to eigt in te case of exponential distribution wen λ = 0.6 Boundaries 3G 4G 5G 6G 7G 8G y y y y y y y
6 58 A. M. H. Al-Kazale Te number of groups or strata depends on te size of te population. If we want to divide te wole data into tree groups, we need to calculate two boundaries y and y since x 0 = y 0 = 0 and y 3 = 7. wic is te largest value. Te boundary y depend on te μ and μ and S and S, also y depends on μ and μ 3 and S and S 3. We calculated tese values by applying equation (.6) depending on equations (.8) and (.9) and in te calculation of te boundaries we assumed te parameter of exponential distribution are λ = 0., 0.4, and 0.6. If we divide te wole data into four groups, we need to calculate tree boundaries y, y and y 3. Similarly we find te boundaries of te strata (group) for te five, six, seven and eigt groups in te same way by using te C++ program. Tables,, and 4 ave sown te boundaries for all groups for te different values of te parameter of exponential distribution for λ = 0., 0.4, and 0.6 respectively. 5 SUMMARY AND DISCUSSION Te metod to determine te two parameters of asymmetric winsorized mean ( rs, ) wen te wole data was divided into tree to eigt groups as follow. We assumed tat te data wit different percentage of missing observations suc as 5%, 0%, and 0% wic as been selected randomly. In te first case wen te data as 5% missing observations, te asymmetric winsorized mean was calculated as te following: Suppose te stratum boundaries for te wole data was divided into tree groups wit 5% missing data wit λ = 0.. Te first group as te interval (0, 3.95), te second group as te interval (3.95, ) and te tird group as te interval (7.9365, 7.) (see Table ). Since we are interested in te first and te last group, ten we count te number of observations in bot groups. Te number of observations in te first one was 9, and te number of observations in tird group was 55 as sown in Table 5a. Finally te asymmetric winsorized mean was calculated by taking r= 9 and s= 55 and by applying equation (3.), wic is equal to Te same procedure to compute te asymmetric winsorized mean was followed for different groups for 5% missing observations as sown in tables 5b, and 5c.
7 Handling missing data on asymmetric distribution 59 Table 5a Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0., and 5%missing observations for tree groups to eigt groups r s W ( r, s ) Moreover, for te data wit 0% missing observations, wen te wole data was divided into tree groups wit λ = 0.,, te number of observations in te first group was 9 wic was te value of r and te number of observations in te last one 50 wic was s, ten te asymmetric winsorized mean was equal to In te same way, te asymmetric winsorized mean was calculated for different groups. Tables 6a, 6b, and 6d ave sown te numbers of observations in te first and last groups for different values of exponential distributions parameter and te values of asymmetric winsorized mean. Table 5b Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0.4 and 5%missing observations for tree groups to eigt groups r s W ( r, s ) Table 5c Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0.6 and 5%missing observations for tree groups to eigt groups r s W ( r, s ) Table 6a Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0., and 0% missing observations for tree groups to eigt groups r s W ( r, s )
8 60 A. M. H. Al-Kazale Table 6b Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0.4 and 0%missing observations for tree groups to eigt groups r s W ( r, s ) Table 6c Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0.6 and 0%missing observations for tree groups to eigt groups r s W ( r, s ) Finally, for 0% missing observations, te data was divided into tree groups. Te number of te observations of te first group was 03 wic was te value of te parameter r and te number of observations in te tird group was 5 wic was te value of parameter s, ten te asymmetric winsorized mean computed was found to be , as sown in Table 7a. Similarly te asymmetric winsorized mean computed for te oter groups wit 0% percentage of missing observations as obtained in Tables 7b, and 7c. Table 7a Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0., and 0% missing observations for tree groups up to eigt groups r s W ( r, s ) Table 7b Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0.4 and 0% missing observations for tree groups up to eigt groups r s W ( r, s )
9 Handling missing data on asymmetric distribution 6 Table 7c Te values of te two parameters and te asymmetric winsorizing mean wen λ = 0.6 and 0% missing observations for tree groups up to eigt groups r s W ( r, s ) We compared te result of using asymmetric winsorized mean metod for missing observations wit oter estimation metods for missing observations. Te first metod is to replace te missing observations wit te mean of te series. Tis mean can be calculated over te entire range of te sample. Second metod is to replace te missing observations wit te naive forecast. Naive model is te simplest form of a univariate forecast model. Tis model uses te current time period's value for te next time period, tat is X ˆt + = X t. Also, we can replace te missing observations wit a simple trend forecast. Tis is accomplised by estimating te regression equation of te form X t = a+ bt (were t is te time) for te periods prior to te missing value. Ten use te equation to fit te time periods missing. Finally, we can replace te missing observations wit an average of te last two known observations tat cover te missing observations (Patricia,994). Te accuracy of estimating te missing observations wit asymmetric winsorized mean depends on ow close te estimating values to te actual values. In practice, we define te difference between te actual and te estimating values as an error. If te estimation is doing a good job, te error will be relatively small. Tis means tat te error for eac time period is purely random fluctuation around original value. So we sould get a value equal 0 or near 0. We tested asymmetric winsorized mean approac to estimating missing data points on time series data sets wit respect to te oter metods. In tis paper te following statistical measures were used as te estimators of accuracy (Patricia, 994). - Te mean absolute errors - Te mean square errors MAE = MSE = n t = n n e t = t n e t (5.) (5.)
10 6 A. M. H. Al-Kazale were n is te number of imputations, e i is difference between te actual and te estimating values. In order to evaluate te accuracy of te estimation procedure for missing observations, we used mean absolute error (MAE) as in Equation (5.) and mean square error (MSE) as in Equation (5.). Tables 8a to 9c sown te amount of te MAE and MSE for te estimation of te missing observations by using te proposed metod, and Tables and for te oter metods. Te MAE values as sown in Tables 8a to 8c for te proposed metods and Table 0 for te oter metods. At te 5% missing observations case, te simple average metod does better tan te asymmetric winsorized metods wen λ = 0., in tree and four groups, but for te oter values, it is clear tat te proposed metod is better tan te oter metods. For te 0% missing observations, te average coverage metod does better tan te oter metods, wile te simple average metod does better tan te asymmetric metod only wen λ = 0., at tree and four groups, λ = 0.4 at tree groups and wen λ = 0.6 at tree, four and five groups. Finally, for te 0% missing observations, tere is only one case; te simple average metod does better tan te asymmetric winsorized metod wen λ = 0., in tree groups. For te second accuracy measures, te MSE as sown in Tables 9a to 9c for te proposed metods and Table for te oter metods. For te 5% missing observations case in general, te asymmetric winsorized mean metods do better tan te oters, but te simple average metod do better tan te proposed metod wen λ = 0., in tree and four groups. For te 0% missing observations, te asymmetric winsorized metods do better tan te naive and trend metods, wile te average bound metod do better tan te proposed metod, te mixed results between te proposed metods and te simple average metod. Finally, for te 0% missing observations, te simple average metod does better tan te proposed metod in two cases wen te λ = 0., in tree groups and λ = 0.6 in tree and four groups. In tis paper, 8 asymmetric winsorized mean metods were proposed to andle missing observations wit 5%, 0%, and 0% missing observations. General results sown tat te proposed metods do better tan te oter metods, in particular wen te data aving 0% missing observations. We recommend using te asymmetric winsorized metod wen te data ave more missing observations and dividing te wole data into more tan four groups. Table 8a Te values of mean absolute error of asymmetric winsorized mean for 5% missing data, λ = 0., 0.4, and 0.6 for different groups MAE 3G 4G 5G 6G 7G 8G λ = 0., λ = λ =
11 Handling missing data on asymmetric distribution 63 Table 8b Te values of mean absolute error of asymmetric winsorized mean for 0% missing data, λ = 0., 0.4, and 0.6 for different groups MAE 3G 4G 5G 6G 7G 8G λ = 0., λ = λ = Table 8c Te values of mean absolute error of asymmetric winsorized mean for 0% missing data, λ = 0., 0.4, and 0.6 for different groups MAE 3G 4G 5G 6G 7G 8G λ = 0., λ = λ = Table 9a Te values of mean square error for asymmetric winsorized mean for 5% missing data, λ = 0., 0.4, and 0.6 for different groups MSE 3G 4G 5G 6G 7G 8G λ = 0., λ = λ = Table 9b Te values of mean square error for asymmetric winsorized mean for 0% missing data, λ = 0., 0.4, and 0.6 for different groups MSE 3G 4G 5G 6G 7G 8G λ = 0., λ = λ = Table 9c Te values of mean square error for asymmetric winsorized mean for 0% missing data, λ = 0., 0.4, and 0.6 for different groups MSE 3G 4G 5G 6G 7G 8G λ = 0., λ = λ =
12 64 A. M. H. Al-Kazale Table 0 Te values of mean absolute error for te oter metods for 5, 0 and 0%. missing MAE Simple Bound Trend Naive Average Average Regression 5% % % Table Te values of mean square error for te oter metods for 5, 0 and 0%. missing MSE Simple Bound Trend Naive Average Average Regression 5% % % REFERENCES [] A.K. Seik, and M. Amad,. Statistical Models of Accelerated Life Testing. Pak. J. Statist. 7 (): (00) [] A. Pankratz,. Forecasting wit Univariate Box - Jenkins Models. Jon Wiley, USA, 983. [3] E. G. Patricia, Introduction to Time Series Modeling and Forecasting in Business and Economics. McGraw-Hill, Inc. New York, 994. [4] G. Hawtorne, and P. Elliot,. Imputing Cross-Sectional Missing Data: Comparison of Common Tecniques. Australian and New Zealand Journal of Psyciatry. 39: (005) [5] I. Hess, V.K. Seti, and and T.R. Balakrisnan, Stratification: a practical investigation. J. Amer. Statist. Assoc. 6, (966) [6] M. Hanif,. Design and Model-based sampling inference. Pak. J. Statist. 6(3): (000) [7] M. Samiuddin, M. Hanif, and A.K.A Kattan, An optimum form for te design based ratio estimator. Pak. J. Statist. 4 () (998) [8] R. Amad Mair, A.M.H. Al-Kazale and M.M.T. Al-Kassab,. Approximation Metod in Finding Optimum Stratum Depending on Neyman Allocation Applied on Beta Distribution. Proceedings of te t WSEAS International Conference on Applied Matematics. Cairo, Egypt. WSEAS Press. (007)
13 Handling missing data on asymmetric distribution 65 [9] R.M. Cyert, and H.J. Davidson, Statistical Sampling for Accounting Information. Prentice-Hall, Englewood Cliffs, NJ, pp. (96) 6 7,. [0] R. J. A. Little, and B.B. Rubin,. Statistical Analysis wit Missing Data. Jon Wiley: New York, 987. [] V. Barnett, and T. Lewis Outliers in Statistical Data. 3 rd Edition, Jon Wiley, New York, USA,994. [] W.G. Cocran. Sampling Tecniques. 3 rd Edition. Jon Wiley, New York, USA, 977. [3] W. Mingxin, and Z. Yijun,. Trimmed and winsorised means based on a scaled deviation. J. of Stat. Planning and Inference. 39: (009) Received: September, 0
New Method to Estimate Missing Data by Using the Asymmetrical Winsorized Mean in a Time Series
Applied Mathematical Sciences, Vol. 3, 2009, no. 35, 1715-1726 New Method to Estimate Missing Data by Using the Asymmetrical Winsorized Mean in a Time Series Ahmad Mahir R. and A. M. H. Al-Khazaleh 1 School
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationVARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR
Sankyā : Te Indian Journal of Statistics 1995, Volume 57, Series B, Pt. 1, pp. 85-92 VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR By SANJAY KUMAR SAXENA Central Soil and Water Conservation Researc
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
More informationEFFICIENT REPLICATION VARIANCE ESTIMATION FOR TWO-PHASE SAMPLING
Statistica Sinica 13(2003), 641-653 EFFICIENT REPLICATION VARIANCE ESTIMATION FOR TWO-PHASE SAMPLING J. K. Kim and R. R. Sitter Hankuk University of Foreign Studies and Simon Fraser University Abstract:
More informationRatio estimation using stratified ranked set sample
METRON - International Journal of Statistics 003, vol. LXI, n. 1, pp. 75-90 HANI M. SAMAWI MAHMOUD I. SIAM Ratio estimation using stratified ranked set sample Summary - Ratio estimation metod is used to
More informationEstimation Approach to Ratio of Two Inventory Population Means in Stratified Random Sampling
American Journal of Operational Researc 05, 5(4): 96-0 DOI: 0.593/j.ajor.050504.03 Estimation Approac to Ratio of Two Inventor Population Means in Stratified Random Sampling Subas Kumar Yadav, S. S. Misra,*,
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationESTIMATION OF A POPULATION MEAN OF A SENSITIVE VARIABLE IN STRATIFIED TWO-PHASE SAMPLING
Pak. J. Stati. 06 Vol. 3(5), 393-404 ESTIMATION OF A POPUATION MEAN OF A SENSITIVE VARIABE IN STRATIFIED TWO-PHASE SAMPING Nadia Muaq, Muammad Noor-ul-Amin and Muammad Hanif National College of Business
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More informationThese errors are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationEFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS
Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationTHE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA EXAMINATION MODULE 5
THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA EXAMINATION NEW MODULAR SCHEME introduced from te examinations in 009 MODULE 5 SOLUTIONS FOR SPECIMEN PAPER B THE QUESTIONS ARE CONTAINED IN A SEPARATE FILE
More informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
More informationIEOR 165 Lecture 10 Distribution Estimation
IEOR 165 Lecture 10 Distribution Estimation 1 Motivating Problem Consider a situation were we ave iid data x i from some unknown distribution. One problem of interest is estimating te distribution tat
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationInvestigation of Tangent Polynomials with a Computer Algebra System The AMATYC Review, Vol. 14, No. 1, Fall 1992, pp
Investigation of Tangent Polynomials wit a Computer Algebra System Te AMATYC Review, Vol. 14, No. 1, Fall 199, pp. -7. Jon H. Matews California State University Fullerton Fullerton, CA 9834 By Russell
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationRegularized Regression
Regularized Regression David M. Blei Columbia University December 5, 205 Modern regression problems are ig dimensional, wic means tat te number of covariates p is large. In practice statisticians regularize
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationChapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1
Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 969) Quick Review..... f ( ) ( ) ( ) 0 ( ) f ( ) f ( ) sin π sin π 0 f ( ). < < < 6. < c c < < c 7. < < < < < 8. 9. 0. c < d d < c
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationLogistic Kernel Estimator and Bandwidth Selection. for Density Function
International Journal of Contemporary Matematical Sciences Vol. 13, 2018, no. 6, 279-286 HIKARI Ltd, www.m-ikari.com ttps://doi.org/10.12988/ijcms.2018.81133 Logistic Kernel Estimator and Bandwidt Selection
More informationThe Derivative The rate of change
Calculus Lia Vas Te Derivative Te rate of cange Knowing and understanding te concept of derivative will enable you to answer te following questions. Let us consider a quantity wose size is described by
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationTeaching Differentiation: A Rare Case for the Problem of the Slope of the Tangent Line
Teacing Differentiation: A Rare Case for te Problem of te Slope of te Tangent Line arxiv:1805.00343v1 [mat.ho] 29 Apr 2018 Roman Kvasov Department of Matematics University of Puerto Rico at Aguadilla Aguadilla,
More informationDe-Coupler Design for an Interacting Tanks System
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationWeek #15 - Word Problems & Differential Equations Section 8.2
Week #1 - Word Problems & Differential Equations Section 8. From Calculus, Single Variable by Huges-Hallett, Gleason, McCallum et. al. Copyrigt 00 by Jon Wiley & Sons, Inc. Tis material is used by permission
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationMath 1241 Calculus Test 1
February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationSchool of Geomatics and Urban Information, Beijing University of Civil Engineering and Architecture, Beijing, China 2
Examination Metod and Implementation for Field Survey Data of Crop Types Based on Multi-resolution Satellite Images Yang Liu, Mingyi Du, Wenquan Zu, Scool of Geomatics and Urban Information, Beijing University
More informationA Generalization of the Lavallée and Hidiroglou Algorithm for Stratification in Business Surveys
Survey Metodology, December 00 191 Vol. 8, No., pp. 191-198 Statistics Canada A Generalization of te Lavallée and Hidiroglou Algoritm for Stratification in Business Surveys LOUIS-PAUL RIVEST 1 ABSTRACT
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationEstimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution
International Journal of Clinical Medicine Researc 2016; 3(5): 76-80 ttp://www.aascit.org/journal/ijcmr ISSN: 2375-3838 Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationWYSE Academic Challenge 2004 State Finals Mathematics Solution Set
WYSE Academic Callenge 00 State Finals Matematics Solution Set. Answer: c. We ave a sstem of tree equations and tree unknowns. We ave te equations: x + + z 0, x + 6 + 7z 9600, and 7x + + z 90. Wen we solve,
More informationSolution for the Homework 4
Solution for te Homework 4 Problem 6.5: In tis section we computed te single-particle translational partition function, tr, by summing over all definite-energy wavefunctions. An alternative approac, owever,
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationPolynomials 3: Powers of x 0 + h
near small binomial Capter 17 Polynomials 3: Powers of + Wile it is easy to compute wit powers of a counting-numerator, it is a lot more difficult to compute wit powers of a decimal-numerator. EXAMPLE
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationFast optimal bandwidth selection for kernel density estimation
Fast optimal bandwidt selection for kernel density estimation Vikas Candrakant Raykar and Ramani Duraiswami Dept of computer science and UMIACS, University of Maryland, CollegePark {vikas,ramani}@csumdedu
More informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More informationOutline. MS121: IT Mathematics. Limits & Continuity Rates of Change & Tangents. Is there a limit to how fast a man can run?
Outline MS11: IT Matematics Limits & Continuity & 1 Limits: Atletics Perspective Jon Carroll Scool of Matematical Sciences Dublin City University 3 Atletics Atletics Outline Is tere a limit to ow fast
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationFinancial Econometrics Prof. Massimo Guidolin
CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationFundamentals of Concept Learning
Aims 09s: COMP947 Macine Learning and Data Mining Fundamentals of Concept Learning Marc, 009 Acknowledgement: Material derived from slides for te book Macine Learning, Tom Mitcell, McGraw-Hill, 997 ttp://www-.cs.cmu.edu/~tom/mlbook.tml
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationRightStart Mathematics
Most recent update: January 7, 2019 RigtStart Matematics Corrections and Updates for Level F/Grade 5 Lessons and Workseets, second edition LESSON / WORKSHEET CHANGE DATE CORRECTION OR UPDATE Lesson 7 04/18/2018
More informationBootstrap confidence intervals in nonparametric regression without an additive model
Bootstrap confidence intervals in nonparametric regression witout an additive model Dimitris N. Politis Abstract Te problem of confidence interval construction in nonparametric regression via te bootstrap
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More informationNew Distribution Theory for the Estimation of Structural Break Point in Mean
New Distribution Teory for te Estimation of Structural Break Point in Mean Liang Jiang Singapore Management University Xiaou Wang Te Cinese University of Hong Kong Jun Yu Singapore Management University
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationThe Krewe of Caesar Problem. David Gurney. Southeastern Louisiana University. SLU 10541, 500 Western Avenue. Hammond, LA
Te Krewe of Caesar Problem David Gurney Souteastern Louisiana University SLU 10541, 500 Western Avenue Hammond, LA 7040 June 19, 00 Krewe of Caesar 1 ABSTRACT Tis paper provides an alternative to te usual
More informationStrati cation by Size Revisited
Journal of Of cial Statistics, Vol. 16, No. 2, 2000, pp. 139±154 Strati cation by Size Revisited Alan H. Dorfman 1 and Ricard Valliant 2 Strati cation by size is used in nite population sampling as a means
More informationFast Exact Univariate Kernel Density Estimation
Fast Exact Univariate Kernel Density Estimation David P. Hofmeyr Department of Statistics and Actuarial Science, Stellenbosc University arxiv:1806.00690v2 [stat.co] 12 Jul 2018 July 13, 2018 Abstract Tis
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationarxiv: v1 [stat.me] 5 Nov 2008
arxiv:0811.0659v1 [stat.me] 5 Nov 2008 Estimation of missing data by using the filtering process in a time series modeling Ahmad Mahir R. and Al-khazaleh A. M. H. School of Mathematical Sciences Faculty
More informationVariance Estimation in Stratified Random Sampling in the Presence of Two Auxiliary Random Variables
International Journal of Science and Researc (IJSR) ISSN (Online): 39-7064 Impact Factor (0): 3.358 Variance Estimation in Stratified Random Sampling in te Presence of Two Auxiliary Random Variables Esubalew
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationAdditional Lecture Notes
Additional Lecture Notes Lecture 1: Introduction Overview Te purposes of tis lecture are (i) to take care of introductions; (ii) to transact course business; (iii) to consider te value of models; (iv)
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More information