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

Similar documents
Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

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

Estimation of Lag Time Between Onset of and Death from an Occult. Tumor

Computer Propagation Analysis Tools

Lecture 17: Kinetics of Phase Growth in a Two-component System:

New sufficient conditions of robust recovery for low-rank matrices

On Control Problem Described by Infinite System of First-Order Differential Equations

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

Low-complexity Algorithms for MIMO Multiplexing Systems

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

Variance and Covariance Processes

Exponential and Logarithmic Equations and Properties of Logarithms. Properties. Properties. log. Exponential. Logarithmic.

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

336 ERIDANI kfk Lp = sup jf(y) ; f () jj j p p whee he supemum is aken ove all open balls = (a ) inr n, jj is he Lebesgue measue of in R n, () =(), f

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

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

Psychometric Methods: Theory into Practice Larry R. Price

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

6 PROBABILITY GENERATING FUNCTIONS

Design Guideline for Buried Hume Pipe Subject to Coupling Forces

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Representing Knowledge. CS 188: Artificial Intelligence Fall Properties of BNs. Independence? Reachability (the Bayes Ball) Example

( ) exp i ω b ( ) [ III-1 ] exp( i ω ab. exp( i ω ba

Circular Motion. Radians. One revolution is equivalent to which is also equivalent to 2π radians. Therefore we can.

envionmen ha implemens all of he common algoihmic deails of all nodal mehods u pemis he specic mehod o e used in any concee insance o e specied y he u

The Substring Search Problem

Determining solar characteristics using planetary data

Consider a Binary antipodal system which produces data of δ (t)

An Automatic Door Sensor Using Image Processing

Auchmuty High School Mathematics Department Advanced Higher Notes Teacher Version

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

7 Wave Equation in Higher Dimensions

Physics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!

The sudden release of a large amount of energy E into a background fluid of density

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

Comparing Different Estimators for Parameters of Kumaraswamy Distribution

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

Special Subject SC434L Digital Video Coding and Compression. ASSIGNMENT 1-Solutions Due Date: Friday 30 August 2002

Molecular Evolution and Phylogeny. Based on: Durbin et al Chapter 8

Chapter 6 Balanced Incomplete Block Design (BIBD)

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

Orthotropic Materials

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Degree of Approximation of a Class of Function by (C, 1) (E, q) Means of Fourier Series

Deflection of light due to rotating mass a comparison among the results of different approaches

LET a random variable x follows the two - parameter

Multiple Experts with Binary Features

Quantum Algorithms for Matrix Products over Semirings

Lecture 22 Electromagnetic Waves

IMPROVED DESIGN EQUATIONS FOR ASYMMETRIC COPLANAR STRIP FOLDED DIPOLES ON A DIELECTRIC SLAB

arxiv: v1 [math.co] 4 Apr 2019

Online Completion of Ill-conditioned Low-Rank Matrices

How to Obtain Desirable Transfer Functions in MIMO Systems Under Internal Stability Using Open and Closed Loop Control

Two-dimensional Effects on the CSR Interaction Forces for an Energy-Chirped Bunch. Rui Li, J. Bisognano, R. Legg, and R. Bosch

Reichenbach and f-generated implications in fuzzy database relations

Unsupervised Segmentation of Moving MPEG Blocks Based on Classification of Temporal Information

The Production of Polarization

KINEMATICS OF RIGID BODIES

Surveillance Points in High Dimensional Spaces

LINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE

Nuclear Medicine Physics 02 Oct. 2007

16 Modeling a Language by a Markov Process

Comparison of some one sample confidence intervals for estimating the mean of the weibull distribution INTRODUCTION.

AN EFFICIENT INTEGRAL METHOD FOR THE COMPUTATION OF THE BODIES MOTION IN ELECTROMAGNETIC FIELD

b Ψ Ψ Principles of Organic Chemistry lecture 22, page 1

Research Design - - Topic 17 Multiple Regression & Multiple Correlation: Two Predictors 2009 R.C. Gardner, Ph.D.

EFFECT OF PERMISSIBLE DELAY ON TWO-WAREHOUSE INVENTORY MODEL FOR DETERIORATING ITEMS WITH SHORTAGES

Numerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method

3.6 Applied Optimization

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2 " P 1 = " #P L L,

Nuclear and Particle Physics - Lecture 20 The shell model

On a quantity that is analogous to potential and a theorem that relates to it

n 1 Cov(X,Y)= ( X i- X )( Y i-y ). N-1 i=1 * If variable X and variable Y tend to increase together, then c(x,y) > 0

Encapsulation theory: the transformation equations of absolute information hiding.

Chapter 7. Interference

1 Spherical multipole moments

arxiv: v2 [stat.me] 13 Jul 2015

ADVANCED SUBSIDIARY (AS) General Certificate of Education Mathematics Assessment Unit F1. assessing. Module FP1: Further Pure Mathematics 1

Chapter Finite Difference Method for Ordinary Differential Equations

Pressure Vessels Thin and Thick-Walled Stress Analysis

Alternative Tests for the Poisson Distribution

Light Time Delay and Apparent Position

6 Matrix Concentration Bounds

Introduction to Arrays

-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.


Risk tolerance and optimal portfolio choice

156 There are 9 books stacked on a shelf. The thickness of each book is either 1 inch or 2

Recent Advances in Chemical Engineering, Biochemistry and Computational Chemistry

MEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION

On the local convexity of the implied volatility curve in uncorrelated stochastic volatility models

Relative and Circular Motion

C/CS/Phys C191 Shor s order (period) finding algorithm and factoring 11/12/14 Fall 2014 Lecture 22

Chapter 2: Introduction to Implicit Equations

b g [2] M N QP L N M

In statistical computations it is desirable to have a simplified system of notation to avoid complicated formulas describing mathematical operations.

Measures the linear dependence or the correlation between r t and r t-p. (summarizes serial dependence)

Appendix B The Relativistic Transformation of Forces

APPLICATION OF MAC IN THE FREQUENCY DOMAIN

Topic 4a Introduction to Root Finding & Bracketing Methods

Transcription:

Mahemaical Theoy and Modeling ISSN 45804 (Pape ISSN 505 (Online Vol.6, No.7, 06 www.iise.og On The Esimaion of Two Missing Values in Randomized Complee Bloc Designs EFFANGA, EFFANGA OKON AND BASSE, E. E. DEPARTMENT OF STATISTICS, UNIVERSIT OF CALABAR, CALABAR. NIGERIA. DEPARTMENT OF STATISTICS, UNIVERSIT OF CALABAR, CALABAR. NIGERIA. ABSTRACT: This pape eviews he wo of Bhad and Ahmed (0, idenify flaws in hei model fomulaion and mae coecions accodingly. Fuhemoe, he fomula fo he esimaion of wo missing values in andomized complee loc designs is deived and is applied o he example given in Bhad and Ahmed (0. The esuls oained fom ou fomula poduces a ee esimae of missing values han ha of Bhad and Ahmed (0. Key wods: Randomized Complee Bloc Design, Sum of Squaes of Eos, Missing Values, Esimaion, and Souces of Vaiailiy.. Inoducion Occasionally one o moe values ae missing in a andomized complee loc design (RCBD due o some easons. Such missing values would have o e esimaed efoe pefoming ANOVA ecause of loss of ohogonaliy. Reseaches on he mehodology fo he esimaion of missing values aound in lieaue. Fomula fo he esimaion of a single missing value in RCBD oiginaed fom he wo of ae (933. Mongomey (976 used he fomula fo single missing value ieaively o esimae wo o moe missing values. Ohe mehods fo esimaing missing values in RCBD can e found in Dempse & e al(977, Jae (978 and Muay (986. Bhad and Ahmed (0 developed a mahemaical pogamming model o esimae missing values in he design wih seveal souces of vaiaion and illusaed hei model wih RCBD wih wo missing values. Thei pape conains some flaws in he funcional consains and hus need o e coeced. Fuhemoe, he ieaive pocedue fo esimaing wo o moe missing values may ae seveal ieaions o convege and hus equie anohe mehod of achieving he same o appoximae esul in a vey sho ime and wih less compuaional effo. In his pape, we coec flaws in Bhad and Ahmed (0 and hen deived a fomula fo he esimaion of wo missing values in RCBD.. Review of Bhad and Ahmed (0 The geneal mahemaical pogamming model inoduced in Bhad and Ahmed (0 is as follows: M: Minimize {f(xi SSE} ( Sujec o: Vaiance(S i i ( x i 0 (3 47

Mahemaical Theoy and Modeling ISSN 45804 (Pape ISSN 505 (Online Vol.6, No.7, 06 www.iise.og Whee S i is he i h souce of vaiaion and he missing values. i is he vaiance of he i h souce of vaiaion wihou consideing The model M was applied o andomized complee loc designs (i.e. design having wo souces of vaiaions, loc and eamen. Thei model which consideed missing values is given as follows: M: Minimize {SSE x. j x i... i i Sujec o: (4 ( i. xi [ ( i. xi ] (5 ( j xi [ ( j xi ].. (6 x i 0 ( i j xi [ ( ij xi ] T (7 Whee, and T ae he fis souce vaiailiy, second souce vaiailiy and he oal vaiailiy.. Flaws in Bhad and Ahmed (0 and hei coecions In model M hee is no disincion eween he suscip used fo he missing values and souces of vaiaions. Also, he i h souce of vaiaion is no esiced o whee he missing value is siuaed as illusaed in he example. If an osevaion is missing in he l h level of vaiaion souce m, hen consains ( could have een expessed as follows: Va(S, m =,,..., n (8 Whee S is l h level of souce m whee a missing value occus, and souce m excluding he missing value. is he vaiance of he l h level of Wih he aove definiions, he model M should have een wien as M3: Minimize {f(xi SSE} Sujec o: Vaiance(S, m =,,..., n; fo each l. x i 0, i,,..., The esiced model M should have een fomulaed as follows: Le ij e he value in he i h level of souce and jh level of souce, whee i =,,..., ; j =,,...,. Suppose he value is missing, whee p is he level of souce and q is he level of souce, hen he vaiance of he p h level of souce, qh level of souce and he oal ae given, especively, as 48

Mahemaical Theoy and Modeling ISSN 45804 (Pape ISSN 505 (Online Vol.6, No.7, 06 www.iise.og Va(S Va(S p q Va(Toal pj pj (9 j q jq iq iq (0 i p ip ij ij ( i p j q ip jq The sum of squaes of eos (SSE should een expessed as follows: Whee SSE. q p... C ( p. = Level p oal of souce vaiaion excluding he missing value.q = Level q oal of souce vaiaion excluding he missing value.. = Gand oal excluding he missing value C = Tems independen of he missing value The coeced vesion of mahemaical pogamming model (M is heefoe:. p. Min{SSE.. q M4: } (3 Sujec o: jq pj jq pj p (4 iq iq q (5 i p ip ij ij T (6 i p j q ip jq 49

Mahemaical Theoy and Modeling ISSN 45804 (Pape ISSN 505 (Online Vol.6, No.7, 06 www.iise.og 0 (7 3. Deiving Compuaional fomula fo esimaing wo missing values Suppose wo values s and ae missing in a andomized complee loc design. We esimae he missing values y solving he unconsained opimizaion polem Mininmize {SSE = SST 0 SSB SST = f( s, } (8 whee, SSE, SST 0, SSB and SST ae he sum of squaes of eos, sum of squaes of oal, sum of squaes of locs, and sum of squaes of eamens, especively. SSE i j ij.. i i... j. j.. i, p j s,q ij s i, p. p. q. s. i.. j js,q.. i, p j s,q ij s i. i, p (. s ( p. js,q. (. (. ( s s q.. j s R s (. s ( p. ( ( q. s s. (.. s (9 Whee., p.,.s,.q and.. ae he especive oals excluding he missing values; R epesen ems no involving missing values. Diffeeniaing SSE in equaion (9 wih espec o e and especively and seing each equal o zeo we oain ( s..s.. (0 s ( p..q.. ( Now solving equaions (0 and ( simulaneously we oain s ( ( ( [(( ] (..s.. p..q.. ( 50

Mahemaical Theoy and Modeling ISSN 45804 (Pape ISSN 505 (Online Vol.6, No.7, 06 www.iise.og ( ( ( [(( ] ( p..q....s.. (3 Example We apply ou fomula o he example in Bha and Ahmed (0. The ale elow shows a andomized complee loc design wih wo missing values and 35. Bloc Teamen i. 3 4 5 6 8.5 5.7 6. 4. 3.0 3.6 9..7.9 4.4 6.9.5 68.4 + 3 5.4 6.6 5.5 0.3 35.5 89.3 + 35 4 6.5 8.6.7 5.7 6.5 8.0 98.0.j 6. 50.9 + 57.3 64.5 46.4 + 35 65.6 346.8 + + 35 In ou fomulas (equaions 0 and, = 4, = 6, =, s =, p = 3, q = 5,. = 68.4, 3. = 89.3,. = 50.9,.5 = 46.4,.. = 346.8. So, 35 ( ( (.... ( [(( ] 3..5 5(4 x 68.4 6 x 50.9 346.8 (4 x 89.3 6 x 46.4 346.8 4.3 4 ( ( ( 3..5.. ( [(( ].... 5(4 x 89.3 6 x 46.4 346.8 (4 x 68.4 6 x 50.9 346.8 8.3 4.. 5. Compaison of wo mehods The esuls oained y he wo mehods ae summaized in he ale elow: Missing Value Bhad and Ahmed (0 Fomula 4.3 4.3 35 7.0 8.3 SSE 80.7 79.6 5

Mahemaical Theoy and Modeling ISSN 45804 (Pape ISSN 505 (Online Vol.6, No.7, 06 www.iise.og 5. Conclusions As can e seen in he aove ale he esimae of missing values oained y ou fomula poduces a lowe sum of squaes of eos and hence a ee esimae han ha of Bhad and Ahmed (0. 5