MEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION

Similar documents
Chapter 6 Balanced Incomplete Block Design (BIBD)

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

A New Method of Estimation of Size-Biased Generalized Logarithmic Series Distribution

Estimation of the Correlation Coefficient for a Bivariate Normal Distribution with Missing Data

Topic 5. Mean separation: Multiple comparisons [ST&D Ch.8, except 8.3]

Alternative Tests for the Poisson Distribution

Chapter 3: Theory of Modular Arithmetic 38

The Substring Search Problem

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},

APPLICATION OF MAC IN THE FREQUENCY DOMAIN

LET a random variable x follows the two - parameter

3.1 Random variables

Multiple Criteria Secretary Problem: A New Approach

Chem 453/544 Fall /08/03. Exam #1 Solutions

On the Poisson Approximation to the Negative Hypergeometric Distribution

Auchmuty High School Mathematics Department Advanced Higher Notes Teacher Version

Encapsulation theory: the transformation equations of absolute information hiding.

Determining solar characteristics using planetary data

This is a very simple sampling mode, and this article propose an algorithm about how to recover x from y in this condition.

Surveillance Points in High Dimensional Spaces

New problems in universal algebraic geometry illustrated by boolean equations

titrrvers:rtt t>1 NO~~H CAROLINA

Encapsulation theory: radial encapsulation. Edmund Kirwan *

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012

CSCE 478/878 Lecture 4: Experimental Design and Analysis. Stephen Scott. 3 Building a tree on the training set Introduction. Outline.

Safety variations in steel designed using Eurocode 3

Basic Bridge Circuits

Solution to HW 3, Ma 1a Fall 2016

Psychometric Methods: Theory into Practice Larry R. Price

MEASURING CHINESE RISK AVERSION

Berkeley Math Circle AIME Preparation March 5, 2013

Question 1: The dipole

NOTE. Some New Bounds for Cover-Free Families

Some technical details on confidence. intervals for LIFT measures in data. mining

A scaling-up methodology for co-rotating twin-screw extruders

Bayesian Analysis of Topp-Leone Distribution under Different Loss Functions and Different Priors

Fractional Zero Forcing via Three-color Forcing Games

A STUDY OF HAMMING CODES AS ERROR CORRECTING CODES

Lecture 7 Topic 5: Multiple Comparisons (means separation)

6 PROBABILITY GENERATING FUNCTIONS

Random Variables and Probability Distribution Random Variable

TESTING THE VALIDITY OF THE EXPONENTIAL MODEL BASED ON TYPE II CENSORED DATA USING TRANSFORMED SAMPLE DATA

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

Empirical Prediction of Fitting Densities in Industrial Workrooms for Ray Tracing. 1 Introduction. 2 Ray Tracing using DRAYCUB

Lecture 28: Convergence of Random Variables and Related Theorems

Fresnel Diffraction. monchromatic light source

Lecture 10. Vertical coordinates General vertical coordinate

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

No. 6 Cell image ecognition with adial hamonic Fouie moments 611 ecognition expeiment is pefomed. Results show that as good image featues, RHFMs ae su

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

radians). Figure 2.1 Figure 2.2 (a) quadrant I angle (b) quadrant II angle is in standard position Terminal side Terminal side Terminal side

THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

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

Measure Estimates of Nodal Sets of Polyharmonic Functions

7.2. Coulomb s Law. The Electric Force

Failure Probability of 2-within-Consecutive-(2, 2)-out-of-(n, m): F System for Special Values of m

Aalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández

Classical Worm algorithms (WA)

GENLOG Multinomial Loglinear and Logit Models

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

JENSEN S INEQUALITY FOR DISTRIBUTIONS POSSESSING HIGHER MOMENTS, WITH APPLICATION TO SHARP BOUNDS FOR LAPLACE-STIELTJES TRANSFORMS

4/18/2005. Statistical Learning Theory

Magnetometer Calibration Algorithm Based on Analytic Geometry Transform Yongjian Yang, Xiaolong Xiao1,Wu Liao

On Polynomials Construction

ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS. D.A. Mojdeh and B. Samadi

Diffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.

Research Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function

3.6 Applied Optimization

1. Review of Probability.

A generalization of the Bernstein polynomials

A DETAILED DESCRIPTION OF THE DISCREPANCY IN FORMULAS FOR THE STANDARD ERROR OF THE DIFFERENCE BETWEEN A RAW AND PARTIAL CORRELATION: A TYPOGRAPHICAL

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

6 Matrix Concentration Bounds

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events

Hammerstein Model Identification Based On Instrumental Variable and Least Square Methods

16 Modeling a Language by a Markov Process

An Overview of Limited Information Goodness-of-Fit Testing in Multidimensional Contingency Tables

International Journal of Mathematical Archive-3(12), 2012, Available online through ISSN

The geometric construction of Ewald sphere and Bragg condition:

Turbulent Prandtl number effect on passive scalar advection

2 K. ENTACHER seies called Es classes, Koobov[4] developed the theoy of good lattice points. Recently, in a seies of papes, Lache et al. [6, 7, 8, 9]

Nuclear Medicine Physics 02 Oct. 2007

On the integration of the equations of hydrodynamics

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

Inverse Square Law and Polarization

Journal of Inequalities in Pure and Applied Mathematics

Educational Testing Service

arxiv: v1 [math.co] 4 May 2017

PHYS 705: Classical Mechanics. Small Oscillations

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

FORMULA FOR SPIRAL CLOUD-RAIN BANDS OF A TROPICAL CYCLONE. Boris S. Yurchak* University of Maryland, Baltimore County, Greenbelt, Maryland =, (2)

INTRODUCTION. 2. Vectors in Physics 1

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

Identification of the degradation of railway ballast under a concrete sleeper

Truncated Squarers with Constant and Variable Correction

arxiv: v1 [math.co] 1 Apr 2011

CROSSTABS. Notation. Marginal and Cell Statistics

Transcription:

MEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION Walte T. Fedee 337 Waen Hall, Biometics Unit Conell Univesity Ithaca, NY 4853 and Tey P. Speed Division of Mathematics & Statistics, CSIRO G.P.O. Box 965 Canbea, ACT 260, Austalia ABSTRACT. In evaluating goodness of a class of designs, eseaches have used a measue of design efficiency poposed by F. Yates in the thities. This measue consides only intablock infomation and does not make use of the infomation contained in the inteblock vaiance. The measues of efficiency poposed hee ae dependent upon the atio of the inteblock and intablock components of vaiance, i.e., aata 2 y. The efficiency of one block design to a second may not emain in$aiant with espect to this atio. Incomplete block designs which wee inefficient unde the intablock measue, now become quite efficient fo some atios of y. Likewise, the indications ae that inteblock infomation should always be ecoveed when analyzing data fom expeiments aanged in an incomplete block design. I. INTRODUCTION. In the mid-thities Yates (e.g. in 937) intoduced an efficiency facto fo patially confounded factoials and fo incomplete block designs. The facto is computed as the atio of the aveage vaiance of a diffeence between two adjusted means (o fo factoial effects) to the vaiance of a diffeence of two means fom an othogonal design such as a completely andomized o andomized complete block design assuming no change in the eo vaiance a 2 fo the two designs. It is common pactice in statistical liteatufe to pesent this efficiency facto fo designs and to discuss optimality of classes of designs in tems of the Yates efficiency facto, which only makes use of the intablock eo vaiance. No use is made of the infomation contained in the inteblock vaiance obtained fom the incomplete blocks (eliminating teatment effects) sum of squaes. A moe pope efficiency facto should make use of the infomation contained in both the intablock and the inteblock vaiances. Some measues accomplishing this ae pesented in this pape. BU-85-M in the Technical Repot Seies of the Biometics Unit, Conell Univesity, Ithaca, NY 4853.

II. BALANCED INCOMPLETE BLOCK DESIGN. A classical balanced incomplete block design (BIBD) consists of v teatments aanged in b incomplete blocks of size k, k<v, with epetitions of each teatment, and with each and evey pai of teatments occuing togethe in an incomplete block A times. The standad elations ae bk = v, A "" (k-)/(v-), and e = (-k )/(-v ) = v(k-)/k(v-). The facto e is the Yates efficiency facto. The usual esponse model assumed fo a BlED is () whee Yhij is the esponse fo the jth teatment in the ith incomplete block in the hth complete block, h=,...,; is,... b/; j,..,v; nhij is one if the jth teatment occus in the hith incomplete block and zeo othewise; ~ is a geneal mean effect; Ph is the hth complete block effect; ahi is the hith andom incomplete block effect distibuted with mean zeo and common vaiance ~a ~j is the jth teatment effect, and thij ae andom eo effects which wee distibuted with mean zeo and vaiance ~~. An analysis of vaiance is given in Table. TABLE. Analysis of vaiance fo a esolvable BIBD. Souce of vaiation Degee of feedom Expected value of mean squae Total Coection fo mean Teatment (ignoing incomplete block effects) Within teatments Blocks (eliminating teatment effects) bk = v v- v(-) b- Complete blocks - Incomplete blocks (elim. t.) b- Intablock eo v-v-b+ ~2 E ~2 t bk-v + ~2 b- a - k~2 + a ()"2 t Expected mean squae fo Ph= 0, i.e., no complete block effects. Intablock infomation o = /~~ Inteblock infomation o' /(0"~ + k~a)

Fo intablock contasts the vaiance of a diffeence between two estimated teatment effects is 2o~/e, when e (l-/k)/(-/v). Fo inteblock contasts the vaiance of a diffeence between two teatments effects is 2(a~ + kaa> (l-e) Fo the combined e.stimato of a diffeence between two teatment effects the vaiance is 2 {7 - e + a2 + ka2 3 2a 2 { + k }.--- ( 3) + k e ' whee = aalo~. Since the intablock contast vaiance is of the fom 2a~/e, it would be lofical to have the combined estimato vaiance in the same fo, i.e., 2ai/el, whee A second measue not involving e is * + k e k(l-e) el... = - + k + k (2) ( 4) * + k e2.. - + (k+h k + +/y ( 5) 2 2 The latte measue of efficiency depends only upon = a~/oe and k; note that (4) is also a function of k and since e k/(k+l) fo v=k2 and k+. A compaison of the two measues is \iven in*table 2 fo k 3,7, and. Thee is little to choose between e and ez and it is suggested that ef be used as a measue of efficiency athe than e~. Note that as appoaches zeo the efficiency fo all k appoaches unity. When appoaches infinity, the Yates efficiency facto e is appoached fo all k. Fo small k, e is elatively low indicating an inefficient design. Howeve, ef indicates that designs with small k ae quite efficient if is /4 to /6, say. Fom these esults, it is suggested that inteblock infomation should always be ecoveed and that inefficiency of incomplete block design is not a poblem.

TABLE 2. Intablock - inteblock efficiencies fo vaious values of y fo k = 3, 7, and. k 3 7 y ef e~ ef e~ e* e~ 0 /32.98.97.98.98.98.98 /6.96.95.96.96.97.96 /4.89.88.92.92.94.94 /2.85.83.90.90.93.93.8.80.89.89.92.92 2.79.78.88.88.92.92 4. 77.76.88.88.92.92 OG.75.75.875.875.97.97 - + ke y + k y e* e~ = + k y + (k+)y Fo a andomized complete block design, the vaiance of a diffeence between two aithmetic means is 2 20' v - k e { + v - wheeas the vaiance of a diffeence between two adjusted teatment means fom a BIBD is 2a 2 e * (7) } ( 6)

Now (6) ~ (7), thei diffeence being v- k + k y + y v- + ke y v( v-k)(k-) y2 (v-)[v- + v(k-) ] (8) (8) is zeo when 0 and/o v = k. Equation (8) could be anothe measue of intablock - inteblock efficiency. Pehaps a moe appopiate measue would be a atio athe a diffeence to obtain (v-k) e~ /ef ( + ) v- The measue e~ would confom moe to the definition of efficiency oiginally pesented by Yates but would include both intablock and inteblock infomation. (9) III. OTHER BLOCK DESIGNS. A class of genealized N-ay designs wee discussed by Shafiq and Fedee (979, 983). Fo these designs the esponse model equation () is eplaced by Yghij V +Ph+ ahi + ~j + Eghij (0) whee g = O,..,nhij and when nhij = 0 thee is no esponse Yghij The othe symbols ae as defined in (). The above authos genealized the Yates efficiency facto fo this class of designs and they also poved that the Fishe inequality v~b holds fo this geneal class of balanced block designs. The efficiency facto ef fo the genealized balanced block design is e* = - c - A. (k + -) 2 whee c = ~ t nhij nhij is the numbe of times teatment j occus in block hi, and = E E nhij is the numbe of eplicates fo teatment j. h i Fo the class of esolvable incomplete block designs known as lattices, the aveage vaiance of a diffeence between two adjusted teatment means is 2 k==+l k+ --} o ()

... 20'2 k-+ e ( ) ( + ) (2) k+ - + (l+ky)- As befoe, we may take e*.., k+ ( k-+ - - + (l+ky)- + ). (3) Note that is the numbe of confounding aangements and ~2 simple (double) lattice, 3 fo the tiple lattice, etc. fo the An intablock - inteblock measue of efficiency like e~ would be e* = 3 [ + ( - (l+ky)- k+ - + (l+ky)- k [ +-y] k+ )] (4) Note that v=k2 and is the numbe of geometical factoial effects confounded. Using the aveage vaiance of a diffeence between two adjusted teatment effects, we could constuct ef and e~ fo any class of incomplete block designs. The ideas in this pape may be used to constuct efficiency factos simila to ef and e~ fo cubic lattices, fo lattice squaes, and fo othe designs. IV. LITERATURE CITED. Shafiq, M. and W. T. Fedee (979). "Genealized N-ay balanced block designs.".b.iomet:.ika 66, pp. 5-23. Shafiq, M. and W. T. Fedee ( 983). "Geneal binay patially balanced designs." Ind.ian J. Ag.ic. St:at:.ist:., XXXV(2). Yates, F. (937). "The design and analysis of factoial expeiments." Impe.ial Bu. Soil Sci., J'ech. Comm. 35, pp. -95.