Comparing Different Estimators for Parameters of Kumaraswamy Distribution

Similar documents
Available online at J. Math. Comput. Sci. 2 (2012), No. 4, ISSN:

Supplementary Information

On a Z-Transformation Approach to a Continuous-Time Markov Process with Nonfixed Transition Rates

The Alpha-Logarithmic Series Distribution of Order k and Some of Its Applications

هقارنت طرائق تقذير هعلواث توزيع كاها ري املعلوتني

S, we call the base curve and the director curve. The straight lines

BAYESIAN ESTIMATION METHOD FOR PARAMETER OF EPIDEMIC SIR REED-FROST MODEL. Puji Kurniawan M

The Central Limit Theorems for Sums of Powers of Function of Independent Random Variables

Degree of Approximation of Fourier Series

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Statistical Estimation

ECSE Partial fraction expansion (m<n) 3 types of poles Simple Real poles Real Equal poles

The Central Limit Theorem

The Non-Truncated Bulk Arrival Queue M x /M/1 with Reneging, Balking, State-Dependent and an Additional Server for Longer Queues

Comparison between Fourier and Corrected Fourier Series Methods

Cameras and World Geometry

Suppose we have observed values t 1, t 2, t n of a random variable T.

Inference of the Second Order Autoregressive. Model with Unit Roots

A Weighted Moving Average Process for Forcasting

Relations on the Apostol Type (p, q)-frobenius-euler Polynomials and Generalizations of the Srivastava-Pintér Addition Theorems

On The Estimation of Two Missing Values in Randomized Complete Block Designs

Capítulo. of Particles: Energy and Momentum Methods

F D D D D F. smoothed value of the data including Y t the most recent data.

STUDY OF THE STRESS-STRENGTH RELIABILITY AMONG THE PARAMETERS OF GENERALIZED INVERSE WEIBULL DISTRIBUTION

Pure Math 30: Explained!

PRICING AMERICAN PUT OPTION WITH DIVIDENDS ON VARIATIONAL INEQUALITY

On imploding cylindrical and spherical shock waves in a perfect gas

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

Professor Wei Zhu. 1. Sampling from the Normal Population

Outline. Review Homework Problem. Review Homework Problem II. Review Dimensionless Problem. Review Convection Problem

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP)

The analysis of the method on the one variable function s limit Ke Wu

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Big O Notation for Time Complexity of Algorithms

Dynamic h-index: the Hirsch index in function of time

ON GENERALIZED FRACTIONAL INTEGRAL OPERATORS. ( ρ( x y ) T ρ f(x) := f(y) R x y n dy, R x y n ρ( y )(1 χ )

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

July 24-25, Overview. Why the Reliability Issue is Important? Some Well-known Reliability Measures. Weibull and lognormal Probability Plots

GENERALIZED FRACTIONAL INTEGRAL OPERATORS AND THEIR MODIFIED VERSIONS

A Weighted Moving Average Process for Forecasting

Extended Laguerre Polynomials

Moment Generating Function

Page 1. Before-After Control-Impact (BACI) Power Analysis For Several Related Populations. Richard A. Hinrichsen. March 3, 2010

[ ] 0. = (2) = a q dimensional vector of observable instrumental variables that are in the information set m constituents of u

Spectrum of The Direct Sum of Operators. 1. Introduction

Mathematical Statistics. 1 Introduction to the materials to be covered in this course

Generalized Fibonacci-Type Sequence and its Properties

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 4 9/16/2013. Applications of the large deviation technique

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

The Nehari Manifold for a Class of Elliptic Equations of P-laplacian Type. S. Khademloo and H. Mohammadnia. afrouzi

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

F.Y. Diploma : Sem. II [AE/CH/FG/ME/PT/PG] Applied Mathematics

ECE 350 Matlab-Based Project #3

Some Newton s Type Inequalities for Geometrically Relative Convex Functions ABSTRACT. 1. Introduction

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

S n. = n. Sum of first n terms of an A. P is

INVESTMENT PROJECT EFFICIENCY EVALUATION

arxiv: v4 [math.pr] 20 Jul 2016

( ) ( ) Weibull Distribution: k ti. u u. Suppose t 1, t 2, t n are times to failure of a group of n mechanisms. The likelihood function is

A Weighted Moving Average Process for Forecasting. Shou Hsing Shih Chris P. Tsokos

Generalized Fibonacci-Lucas Sequence

THE ANALYSIS OF SOME MODELS FOR CLAIM PROCESSING IN INSURANCE COMPANIES

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b),

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

MATH Midterm Solutions

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

Transistor configurations: There are three main ways to place a FET/BJT in an architecture:

This web appendix outlines sketch of proofs in Sections 3 5 of the paper. In this appendix we will use the following notations: c i. j=1.

Sampling Example. ( ) δ ( f 1) (1/2)cos(12πt), T 0 = 1

BINOMIAL THEOREM OBJECTIVE PROBLEMS in the expansion of ( 3 +kx ) are equal. Then k =

Statistical Optics and Free Electron Lasers

Research Design - - Topic 2 Inferential Statistics: The t-test 2010 R.C. Gardner, Ph.D. Independent t-test

Research Article Estimating the Reliability Function for a Family of Exponentiated Distributions

Exercise 3 Stochastic Models of Manufacturing Systems 4T400, 6 May

Lecture 8 April 18, 2018

Greatest term (numerically) in the expansion of (1 + x) Method 1 Let T

On Some Fractional Integral Operators Involving Generalized Gauss Hypergeometric Functions

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

On the Validity of the Pairs Bootstrap for Lasso Estimators

LESSON 15: COMPOUND INTEREST

Using Difference Equations to Generalize Results for Periodic Nested Radicals

O & M Cost O & M Cost

1 Notes on Little s Law (l = λw)

ABSOLUTE INDEXED SUMMABILITY FACTOR OF AN INFINITE SERIES USING QUASI-F-POWER INCREASING SEQUENCES

On a Problem of Littlewood

B. Maddah INDE 504 Simulation 09/02/17

INF 5460 Electronic noise Estimates and countermeasures. Lecture 13 (Mot 10) Amplifier Architectures

Numerical Solution of Sine-Gordon Equation by Reduced Differential Transform Method

ON POINTWISE APPROXIMATION OF FUNCTIONS BY SOME MATRIX MEANS OF FOURIER SERIES

Reinforcement learning

Progression. CATsyllabus.com. CATsyllabus.com. Sequence & Series. Arithmetic Progression (A.P.) n th term of an A.P.

A New Statistical competitive Data Clustering Approaches: SRPCL and SRPCCL

Institute of Actuaries of India

OLS bias for econometric models with errors-in-variables. The Lucas-critique Supplementary note to Lecture 17

By the end of this section you will be able to prove the Chinese Remainder Theorem apply this theorem to solve simultaneous linear congruences

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier

The shortest path between two truths in the real domain passes through the complex domain. J. Hadamard

A Note on Prediction with Misspecified Models

Using Counting Techniques to Determine Probabilities

Transcription:

Compaig Diffee Esimaos fo Paamees of Kumaaswamy Disibuio ا.م.د نذير عباس ابراهيم الشمري جامعة النهرين/بغداد-العراق أ.م.د نشات جاسم محمد الجامعة التقنية الوسطى/بغداد- العراق Absac: This pape deals wih compaig diffee mehods o esimae scale paamee ( ) ad shape paamee ( ) of Kumaaswamy Disibuio, hese esimaos ae mome esimaos,maximum likelihood,pobabiliy weighed mome, mehod of maximum likelihood of odeed obsevaio ad L-mome esimaos. he compaiso is doe hough simulaio usig diffee sample size ad diffee se of values of (, ) ad he compaig he esuls usig saisical measues mea squae eo (MSE) keywods:kumaaswamy Disibuio, mome mehod,maximum likelihood mehod,pobabiliy weighed mome mehod,l-mome mehod. Ioducio: This disibuio is oe of he coiuous pobabiliy disibuios which have may applicaios,such as heigh of idividuals scoes obaied i a es which lies bewee (-) ad ca be asfomed o ieval (,).May eseaches,nadaajah discussed his disibuio as a esul of bea family,also Fle Che sudied he disibuio ad show is applicaio i hydology daa daily seam flow. also Gupa ad Nadaajah (24),wie abou he esimaio of is paamees.i (23) Tabassum Sidhu ad Navid Feoze ad Muhammed Aslam focus o o-bayesia esimaio fo he shape paamees of he Kumaaswamy disibuio. Hee we coiue he wok abou his disibuio o compae diffee esimaos fo he wo paamees by usig fou mehods such as :mome mehod,maximum likelihood mehod, pobabiliy weighed momes,maximum likelihood of odeed obsevaio, he compaiso is doe hough simulaio The pobabiliy desiy fucio of he disibuio is: f T ( ) ( ),,, () 399

Ad he cumulaive disibuio fucio is: F ( ) ( ) T.. (2) Ad he eliabiliy fucio is: R( ) ( ) (3) The hazad fucio is: H( ) (4) Which is ime depede. - Mehod of mome This mehod based o piciple of equaig populaio momes by hei coespodig sample momes. he populaio momes is: E( ) ( ) d (5) Le z / ( z) d dz dz ( z) / Afe some algebic opeaios,the.. (6) / ) ) E( ) (7) / ) h The fomula of populaio mome abou oigi. whe, / ) ) E ( ) / ) Fis Sample mome is i i The he fis equaio is: / ) ). (8) / ) Ad whe 2 2 / ) ) E ( 2 ) 2 / ) 4

Secod Sample mome 2 i i he 2 / ) ) E ( 2 ) = 2 / ) 2 i i (9) Fom hese wo equaio we obai a momes esimaos ( mom, mom ) by solvig i umeically hough Newo-Raphso Mehod. 2- Mehod of Maximum Likelihood The mehod of maximum likelihood selecs he se of values of he model paamees ha maximizes he likelihood fucio. I ideed maximizes he pobabiliy of he obseved daa ude he esulig disibuio. Maximum-likelihood esimaio gives a uified appoach o esimaio, which is well-defied i he case of he omal disibuio ad may ohe poblems. So he likelihood of he disibuio is : L i ( i )... () i i The logaihm of likelihood fucio is: LogL log( ) log( ) ( ) log( i ) ( ) log( ) () i To obai he maximum likelihood esimaos we deive equaio (),wih espec o ad,he log L log L i log( i ) i i log( i ) ( ) i i ( i ) i.. (2) log( ) i (3) Which ca be solved umeically o obai a maximum likelihood esimaos ( mle, mle ) by solvig i umeically hough Newo-Raphso Mehod. 3- Mehod of Pobabiliy Weighed Momes This mehod depeds o equaig weighed momes k wih populaio momes obaied, whee 4

k i k ( pi ). (4) i Whee p i ca be esimaed by ay o-paameic esimao, he commo oe is ( i.35) p i. (5) By equaig he heoeical mome by coespodig weighed momes ' M ),we ge: ( k / ) ) / ) 2 / ) ) 2 / ) = i = i k i ( pi ).. (6) i ( p 2 i )... (7) These equaio's Which ca be solved umeically o obai a weighed mome esimaos ( pwm, pwm ) by solvig i umeically hough Newo- Raphso Mehod. 4- Mehod of L-Mome The esimaio of paamees by his mehod deped o equaig populaio L-momes by coespodig sample momes.i.e Le populaio L-mome defied by: F ( ) f ( ) d. (8) Ad coespodig sample l-mome is: b C i C i ( i)..(9) Now fo he p.d.f give i equaio () ad c.d.f give i equaio (2) we mus solve he iegal ( ) ( ) d ( ) ( ) d Le: z ( ) / / z Afe usig above asfomaio we fid i i ( ) Ci Bea, (2) i Fom his fomula Bea (/, ) Bea (/,/ ) (2) Ad 42

2 Bea (/,2/ ) 2Bea (/,/ ) Bea (/, ). (22) Equaig wih sample momes obaied fom i b C ( i).. (23) C i These equaio's Which ca be solved umeically o obai a L-mome esimaos ( Lmom, Lmom ) by solvig i umeically hough Newo- Raphso Mehod. Resuls I his simulaio sudy, we have chose,5,25,5, o epese small, modeae ad lage sample size, seveal values of paamee.5,,.5 ad.5,,. 5.The umbe of eplicaio used was (L=).The simulaio pogam was wie by usig malab-r2b pogam. Afe esimaig he paamees he mea squae eo (MSE) was calculaed o compae he mehods of esimaio, Whee: ^ 2 ( ) ^ Mse( ) l.. (24) ^ 2 ( ) ^ Mse( ) l (25) The esuls of he simulaio sudy ae summaized ad abulaed i able() ad able(2) of he fou esimaos fo all sample size ad (, ) values especively. I each ow of able() ad able(2) we have fou values of esimaos ha is he mome mehod, maximum likelihood mehod,pobabiliy weighed mome mehod,l-mome mehod. he bes mehod is he L- mome ha gives smalles value of (MSE). So we ecommeded o use i i esimaio of he paamees of he disibuio. Table : MSE fo diffee esimaos of fo Kummaasumay disibuio sample size () Mom mle Pwm L-mom Bes.5.34.453.673.375 L-mom.987.987 2..9545 L-mom.5.23.342.792.95 L-mom.5.44.32.752.5 L-mom 5.983.97.42.955 L-mom.5..98.42.775 L-mom 25.5.22..532.9895 L-mom 43

5.932.92.9632.8995 L-mom.5.22..532.9895 L-mom.5.976.964.72.9435 L-mom.876.864.972.8435 L-mom.5.987.975.82.9545 L-mom.5.954.942.9852.925 L-mom.654.6.6433.63 Mle.5.546.27.342.275 Mle Table 2: MSE fo diffee esimaos of fo Kummaasumay disibuio sample size () mom mle Pwm L-mom Bes 5 25 5.5.483 2.24 2.6.25 L-mom.9523.94.933.94 L-mom.5.373.463.95.9 L-mom.5.9973.53.2383.8923 L-mom.9363.93.773.833 L-mom.5.633.273.343.9583 L-mom.5.9753.293.263.873 L-mom.8853.393.263.783 L-mom.5.9753.293.263.873 L-mom.5.9293.833.73.8243 L-mom.8293.9833.73.7243 L-mom.5.943.943.83.8353 L-mom.5.973.792.483.823 L-mom.5654.794.864.464 L-mom.5.2363.393.4773.33 L-mom Refeces - Al-Ahai F.M.(2),"Paamee Esimaio fo he Double Paeo Disibuio"Joual of Mahemaics ad Saisics,7(4):289-294. 2- Sidhu T.N., Feoze N. ad Aslam M.,(23)" Bayesia Aalysis of he Kumaaswamy Disibuio ude Failue Cesoig Samplig Scheme ".Ieaioal Joual of Advaced Sciece ad Techology Vol. 5. 3- P. Kumaaswamy, (98) A Geealized Pobabiliy Desiy Fucio fo Double-Bouded Radom Pocesses, Joual of Hydology, vol. 46, pp. 79-88. 4- S. Nadaajah,(28) O he Disibuio of Kumaaswamy, J Hydol, vol. 348, pp. 568-569. 5- M. C. Joes,(29) Kumaaswamy s Disibuio: A Bea Type Disibuio wih Some Tacabiliy Advaages, Sa Mehodol, vol. 6, pp. 7-8. 44

6-M.Gag,(29) O Geealized Ode Saisics fom Kumaaswamy Disibuio, Tamsui Oxfod J Mah Sci., vol. 25, o. 2, pp.53-66. 45