Fuzzy Geometric Programming in Multivariate Stratified Sample Surveys in Presence of Non-Response with Quadratic Cost Function

Similar documents
Department of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India

VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR

Handling Missing Data on Asymmetric Distribution

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Estimation Approach to Ratio of Two Inventory Population Means in Stratified Random Sampling

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS

Estimation of Population Mean on Recent Occasion under Non-Response in h-occasion Successive Sampling

Ratio estimation using stratified ranked set sample

ch (for some fixed positive number c) reaching c

Explicit Interleavers for a Repeat Accumulate Accumulate (RAA) code construction

ESTIMATION OF A POPULATION MEAN OF A SENSITIVE VARIABLE IN STRATIFIED TWO-PHASE SAMPLING

232 Calculus and Structures

The Open Petroleum Engineering Journal

EFFICIENT REPLICATION VARIANCE ESTIMATION FOR TWO-PHASE SAMPLING

WINKLER PLATES BY THE BOUNDARY KNOT METHOD

A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

Symmetry Labeling of Molecular Energies

Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

Equilibrium and Pareto Efficiency in an exchange economy

Artificial Neural Network Model Based Estimation of Finite Population Total

Strati cation by Size Revisited

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations

Research Article New Results on Multiple Solutions for Nth-Order Fuzzy Differential Equations under Generalized Differentiability

FUNDAMENTAL ECONOMICS Vol. I - Walrasian and Non-Walrasian Microeconomics - Anjan Mukherji WALRASIAN AND NON-WALRASIAN MICROECONOMICS

NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura

LINEAR MULTISCALE AUTOREGRESSIVE MODEL FOR FORECASTING SEASONAL TIME SERIES DATA

Fast Computation of Capacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow

About the Riemann Hypothesis

Parameterized Soft Complex Fuzzy Sets

5.1 We will begin this section with the definition of a rational expression. We

GELFAND S PROOF OF WIENER S THEOREM

Dedicated to the 70th birthday of Professor Lin Qun

NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE

Numerical analysis of a free piston problem

Walrasian Equilibrium in an exchange economy

SECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES

Exponential Ratio Type Estimators In Stratified Random Sampling

Forces, centre of gravity, reactions and stability

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

RP 2.4. SEG/Houston 2005 Annual Meeting 1513

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Evaluation and Accurate Estimation from Petrophysical Parameters of a Reservoir

THE COMPLETE SOLUTION PROCEDURE FOR THE FUZZY EOQ INVENTORY MODEL WITH LINEAR AND FIXED BACK ORDER COST

Click here to see an animation of the derivative

arxiv: v3 [math.gr] 7 Dec Introduction

Periodic Orbits in the Photogravitational Restricted Problem When the Primaries Are Triaxial Rigid Bodies

LEARNING FROM MISTAKES

A = h w (1) Error Analysis Physics 141

Computational Method of Structural Reliability Based on Integration Algorithms

ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING

On Approximating Expected Log-Utility of. Mutual Fund Investments

On Local Linear Regression Estimation of Finite Population Totals in Model Based Surveys

Research Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation

Variational Localizations of the Dual Weighted Residual Estimator

lecture 26: Richardson extrapolation

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

Section 3: The Derivative Definition of the Derivative

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

ON THE CONSTRUCTION OF REGULAR A-OPTIMAL SPRING BALANCE WEIGHING DESIGNS

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

Complexity of Decoding Positive-Rate Reed-Solomon Codes

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY

Variance Estimation in Stratified Random Sampling in the Presence of Two Auxiliary Random Variables

Material for Difference Quotient

Tangent Lines-1. Tangent Lines

Financial Econometrics Prof. Massimo Guidolin

DIGRAPHS FROM POWERS MODULO p

7 Semiparametric Methods and Partially Linear Regression

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines

1 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

The Verlet Algorithm for Molecular Dynamics Simulations

arxiv: v1 [math.dg] 4 Feb 2015

INTRODUCTION AND MATHEMATICAL CONCEPTS

MANY scientific and engineering problems can be

Large deflection analysis of rhombic sandwich plates placed on elastic foundation

Poisson Equation in Sobolev Spaces

3. THE EXCHANGE ECONOMY

Quantum Numbers and Rules

Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series

Combining functions: algebraic methods

Tail Conditional Expectations for Extended Exponential Dispersion Models

School of Geomatics and Urban Information, Beijing University of Civil Engineering and Architecture, Beijing, China 2

Discriminate Modelling of Peak and Off-Peak Motorway Capacity

A Reconsideration of Matter Waves

MATH Precalculus (Revised August 2014)

DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING

Wind Turbine Micrositing: Comparison of Finite Difference Method and Computational Fluid Dynamics

INTEGRATING IMPERFECTION OF INFORMATION INTO THE PROMETHEE MULTICRITERIA DECISION AID METHODS: A GENERAL FRAMEWORK

arxiv: v1 [math.pr] 28 Dec 2018

Pre-Calculus Review Preemptive Strike

Non-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model

On One Justification on the Use of Hybrids for the Solution of First Order Initial Value Problems of Ordinary Differential Equations

ETNA Kent State University

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA EXAMINATION MODULE 5

Application of Quintic B-splines Collocation Method on Some Rosenau Type Nonlinear Higher Order Evolution Equations.

, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim

Transcription:

American Journal of Operations Researc, 04, 4, 73-88 Publised Online May 04 in SciRes. ttp://.scirp.org/ournal/aor ttp://dx.doi.org/0.436/aor.04.4307 Fuzzy Geometric Programming in Multivariate Stratified Sample Surveys in Presence of Non-Response it Quadratic Cost Function Safiulla, Moammad Faisal Kan, Irfan Ali Department of Statistics & Operations Researc, Aligar Muslim University, Aligar, India Department of Computing and Informatics Saudi Electronics University, Riyad, Saudi Arabia Email: safi.stats@gmail.com, faisalkan004@yaoo.com, irfii.ali@gmail.com Received 9 Marc 04; revised 9 April 04; accepted 6 April 04 Copyrigt 04 by autors and Scientific Researc Publising Inc. Tis ork is licensed under te Creative Commons Attribution International icense (CC BY). ttp://creativecommons.org/licenses/by/4.0/ Abstract In tis paper, te problem of non-response it significant travel costs in multivariate stratified sample surveys as been formulated of as a Multi-Obective Geometric Programming Problem (MOGPP). Te fuzzy programming approac as been described for solving te formulated MOGPP. Te formulated MOGPP as been solved it te elp of INGO Softare and te dual solution is obtained. Te optimum allocations of sample sizes of respondents and non respondents are obtained it te elp of dual solutions and primal-dual relationsip teorem. A numerical example is given to illustrate te procedure. Keyords Geometric Programming, Fuzzy Programming, Non-Response it Travel Cost, Optimum Allocations, Multivariate Stratified Sample Surveys. Introduction In sampling te precision of an estimator of te population parameters depends on te size of te sample and variability among te units of te population. If te population is eterogeneous and size of te sample depends on te cost of te survey, ten it is likely to be impossible to get a sufficiently precise estimate it te elp of simple random sampling from te entire population. In order to estimate te population mean or total it greater precision, te eterogeneous population is divided into mutually-exclusive, exaustive and non-overlap- Ho to cite tis paper: Safiulla, Kan, M.F. and Ali, I. (04) Fuzzy Geometric Programming in Multivariate Stratified Sample Surveys in Presence of Non-Response it Quadratic Cost Function. American Journal of Operations Researc, 4, 73-88. ttp://dx.doi.org/0.436/aor.04.4307

Safiulla et al. ping strata ic ill be more omogeneous tan te entire population. Te entire population is called Stratified Random Sampling. Te problem of optimum allocation in stratified random sampling for univariate population is ell knon in sampling literature; see for example Cocran [] and Sukatme et al. []. In multivariate stratified sample surveys problems te non-response can appear en te required data are not obtained. Te problem of non-response may occur due to te refusal by respondents or teir not being at ome, making te information of sample inaccessible. Te problem of non-response occurs in almost all surveys. Te extent of nonresponse depends on various factors suc as type of te target population, type of te survey and te time of survey. For te problem of non-response in stratified sampling it may be assumed tat every stratum is divided into to mutually exclusive and exaustive groups of respondents and non respondents. Hansen and Huritz [3] presented a classical non-response teory ic as first developed for surveys in ic te first attempt as made by mailing te questionnaires and a second attempt as made by personal intervies to a sub sample of te non respondents. Tey constructed te estimator for te population mean and derived te expression for its variance and also orked out te optimum sampling fraction among te non respondents. El-Badry [4] furter extended te Hansen and Huritz s tecnique by sending aves of questionnaires to te non respondent units to increase te response rate. Te generalized El-Badry s approac for different sampling design as given by Foradori [5]. Srinat [6] suggested te selection of sub samples by making several attempts. Kare [7] investigated te problem of optimum allocation in stratified sampling in presence of nonresponse for fixed cost as ell as for fixed precision of te estimate. Kan et al. [8] suggested a tecnique for te problem of determining te optimum allocation and te optimum sizes of subsamples to various strata in multivariate stratified sampling in presence of non-response ic is formulated as a Nonlinear Programming Problem (NPP). Varsney et al. [9] formulated te multivariate stratified random sampling in te presence of non-response as a Multi-obective Integer Nonlinear Programming problem and a solution procedure is developed using lexicograpic goal programming tecnique to determine te compromise allocation. Fatima and Asan [0] addressed te problem of optimum allocation in stratified sampling in te presence of non-response and formulated as an All Integer Nonlinear Programming Problem (AINPP). Varsney et al. [] ave considered te multivariate stratified population it unknon strata eigts and an optimum sampling design is proposed in te presence of non-response to estimate te unknon population means using DSS strategy and developed a solution procedure using Goal Programming tecnique and obtained an integer solution directly by te optimization softare INGO. Ragav et al. [] as discussed te various multi-obective optimization tecniques in te multivariate stratified sample surveys in case of non-response. Geometric Programming (GP) is a smoot, systematic and an effective non-linear programming metod used for solving te problems of sample surveys, management, transportation, engineering design etc. tat takes te form of convex programming. Te convex programming problems occurring in GP are generally represented by an exponential or poer function. GP as certain advantages over te oter optimization metods because it is usually muc simpler to ork it te dual tan te primal one. Te degree of difficulty (DD) plays a significant role for solving a non-linear programming problem by GP metod. Geometric Programming (GP) as been knon as an optimization tool for solving te problems in various fields from 960 s. Duffin, Peterson and Zener [3] and also Zener [4] ave discussed te basic concepts and teories of GP it application in engineering in teir books. Beigtler, C.S., and Pililps, D.T. [5], ave also publised a famous book on GP and its application. Davis and Rudolp [6] applied GP to optimal allocation of integrated samples in quality control. Amed and Carles [7] applied geometric programming to obtain te optimum allocation problems in multivariate double sampling. Oa, A.K. and Das, A.K. [8] ave taken te Multi- Obective Geometric Programming Problem being cost coefficient as continuous function it eigted mean and used te geometric programming tecnique for te solutions. Maqbool et al. [9] and Safiulla et al. [0] ave discussed te geometric programming approac to find te optimum allocations in multivariate to-stage sampling and tree-stage sample surveys respectively. In many real-orld decision-making problems of sample surveys, environmental, social, economical and tecnical areas are of multiple-obectives. It is significant to realize tat multiple obectives are often non-commensurable and in conflict it eac oter in optimization problems. Te multi-obective models it fuzzy obectives are more realistic tan deterministic of it. Te concept of fuzzy set teory as firstly given by Zade []. ater on, Bellman and Zade [] used te fuzzy set teory for te decision-making problem. Tanaka et al. [3] introduces te obective as fuzzy goal over te α-cut of a fuzzy constraint set and Zimmermann [4] gave te concept to solve multi-obective linear-programming problem. Bisal [5] and Verma [6] developed fuzzy 74

Safiulla et al. geometric programming tecnique to solve Multi-Obective Geometric Programming (MOGP) problem. Islam [7] [8] as discussed modified geometric programming problem and its applications and also anoter fuzzy geometric programming tecnique to solve MOGPP and its applications. Fuzzy matematical programming as been applied to several fields. In tis paper, e ave formulated te problem of non-response it significant travel cost ere te cost is quadratic in n in multivariate stratified sample surveys as a Multi-Obective Geometric Programming problem (MOGPP). Te fuzzy programming approac as been described for solving te formulated MOGPP and te optimum allocations of sample sizes of respondents and non respondents are obtained. A numerical example is given to illustrate te procedure.. Formulation of te Problem In stratified sampling te population of N units is first divided into non-overlapping subpopulation called strata, of sizes N, N,, N,, N it N = N and te respective sample sizes itin strata are dran = it independent simple random sampling denoted by n, n,, n,, n it n. = n = t et for te stratum: N : Stratum size. Y : Stratum mean. S : Stratum variance. Wˆ ˆ = N N : te estimated stratum eigt among respondents. Wˆ ˆ = N N : te estimated stratum eigt among non-respondents. N : te sizes of te respondents. N = N N : te sizes of non respondents groups. t n : Units are dran from te stratum. Furter let out of n, n units belong to te respondents group. n = n n : Units belong to te non respondents group. n = n = : Te total sample size. A more careful second attempt is made to obtain information on a random subsample of size r out of n non respondents for te representation from te non respondents group of te sample. r = n k; =,,, : Subsamples of sizes at te second attempt to be dran from n non-respondent group of te t stratum. Were k and k denote te sampling fraction among non respondents. Since N and N are random variables ence teir unbiased estimates are given as: Nˆ = n N n : Te unbiased estimates of te respondents group. Nˆ = n N n : Te unbiased estimate of te non respondents group. y ;,, = p: denotes te sample means of t caracteristic measured on te n respondents at te first attempt. y ;,, r = p: denotes te r sub sampled units from non respondents at te second attempt. t Using te estimator of Hansen and Huritz [3], te stratum mean Y for caracteristic in te t stratum may be estimated by n y ny( r ) y = () n t t It can be seen tat y is an unbiased estimate of te stratum mean Y of te stratum for te caracteristic it a variance. W S W S v( y ) = S () n N r n t t ere S is te stratum variance of caracteristic in te stratum; =,,, p and =,,, given as: N S = yi Y N i= 75

Safiulla et al. ere y denote te value of te i N i= i t i unit of te Y = N y : is te stratum mean of y. i S is te stratum variance of te t stratum for t caracteristic in te Nˆ = ˆ i N i= t caracteristic. t stratum among non respondents, given by: Sˆ y Y Nˆ Y ˆ = N y i= i is te stratum mean of y among non respondents. i ˆ ˆ t W = N N is stratum eigt of non respondents in stratum. If te true values of S and S are not knon tey can be estimated troug a preliminary sample or te value of some previous occasion, if available, may be used. Furtermore, te variance of y =, Wy (ignoring fpc) is given as: = ( ˆ ˆ ) W S W S WWˆ Sˆ V y = W v y = = f nr ( ) ( ) = = n = r 0 (, ) (3) ere y is an unbiased estimate of te overall population mean Y of te t caracteristic and V ( y ) is as given in Equation (). Assuming a linear cost function te total cost C of te sample survey may be given as: C = c n c n c n 0 = = = ere c 0 = te per unit cost of making te first attempt, c = c = is te per unit cost for processing te t results of all te p caracteristics on te n selected units from respondents group in te stratum in te first p attempt and c = c : te per unit cost for measuring and processing te results of all te p caracteristics = t on te r units selected from te non respondents group in te stratum in te second attempt. Also, c t and c are per unit costs of measuring te caracteristic in first and second attempts respectively. As n is not knon until te first attempt as been made, te quantity W n may be used as its expected value. Te total expected cost Ĉ of te survey may be given as: ( 0 ) 0 (4) Cˆ = c c W n c r t n t r = = = = Te problem terefore reduces to find te optimal values of sample sizes of respondents respondents r ic are expressed as: ( ) W S W S WW S Min ( f0 ) = = n = r Subect to, =,,, p ( c0 c W) n cr t0 n t r C0 = = = = n, r 0 and =,,, 3. Geometric Programming Formulation n and non- Te folloing Multi-obective Nonlinear Programming Problem (MNPP) it te cost function quadratic in n and significant travel cost is defined in Equations (6) as follos: (5) 76

Safiulla et al. Min V Subect to ( ˆ π s ) a = = n c c W n c r t n t r C n,r 0 and =,,, 0 0 0 = = = = =,,, P Similarly, te expression (6) can be expressed in te standard Primal GPP it cost function quadratic in n ere te travel cost is significant is given as follos: ( nr) Max f0, Subect fq nr, =,, p n, r 0, =,,, { } p [ ] i p ere,, 0,,,, i pi p i p i pi fq nr = di n n n r r r q= k i q p i p i or fq( nr, ) = di n r, di > 0, n > 0, q= 0,,,, k, i q [ ] = = p i : arbitrary real numbers, d i : positive and fq ( n ) : posinomials et for simplicity f0 = a and di = a = ( c c W ) C = c C = t C = t C ( ) 0 0 0 0 0 0 c c W c t t ere f nr, n r n r r 0 0 q = = C0 = C0 = C0 = C0 Te dual form of te Primal GPP ic is stated in (7) can be given as: i i k k d [ ] i q i Max v0 = i ( i) q 0 i q [ ] = i q= i q [ ] Subect to i = ( ii) i [ 0] k pii ( iii) q i q [ ] i 0, q,,, k and i =,,, mk ( iv) =,, p. Te above formulated dual GPP (8) can be solved in te folloing to steps: Step : For te Optimum value of te obective function, te obective function alays takes te form: C 0 ( x ) ( s in te last constraint s) 0 0 Coeffi. of first term Coeffi. of Second term = Coeffi. of last term k 0 0 k ( s in te first constraint s) s in te last constraint s Te Multi-Obective obective function for our problem is: k i k d i q [ ] i i [ ] i i q [ ] q i q q= i s in te first constraint s Step : Te equations tat can be used for GPP for te eigts are given belo: i in te obective function = (Normality condition ) and for eac primal variable n & i & o [ ] r aving m terms. (6) (7) (8) n and 77

Safiulla et al. m k i= ( i ) ( n n r r ) for eac term exponent on & and & in tat term (Ortogonality condition) and i 0 (Positivity condition). Te above problem (8) as been solved it te elp of steps (-) discussed in Section (3) and te corresponding solutions 0i is te unique solution to te dual constraints; it ill also maximize te obective function for te dual problem. Next, te solution of te primal problem ill be obtained using primal-dual relationsip teorem ic is given belo. 4. Primal-Dual Relationsip Teorem If 0i primal problem (7) satisfies te system of equations: is a maximizing point for dual problem (8), eac minimizing points (,,, and,,, ) f 0 ( n) [ ] 0 iv, i J 0, = i, i J, v( 0 i) ere ranges over all positive integers for ic v( 0i) > 0. Te optimal values of sample sizes of te respondents ( n ) and non-respondents ( ) it te elp of te primal-dual relationsip teorem (9). 5. Fuzzy Geometric Programming Approac [ ] n n n n r r r r for 3 4 3 4 (9) r can be calculated Te solution procedure to solve te problem (5) consists of te folloing steps: Step-: Solve te MOGPP as a single obective problem using only one obective at a time and ignoring te oters. Tese solutions are knon as ideal solution. Step-: From te results of step-, determine te corresponding values for every obective at eac solution ( ) ( ) p p derived. et n, r, n, r,, n, r,, n, r are te ideal solutions of te obective functions ( ) ( ) ( ) ( ) ( ) ( ) ( p) ( p) (, ), (, ),, (, ),, p(, ). ( ) ( ) ( ) ( ) ( p) ( p) U Max { f ( n, r ), f ( n, r ),, f (, )} 0 0 0 p n r f n r f n r f n r f n r 0 0 0 0 So = and ( ) ( ) 0 = f n, r, =,,, p. t U and be te upper and loer bonds of te obective function f0 nr,, =,,, p. Step 3: Te membersip function for te given problem can be define as: µ ( ) 0, if f0 nr, U U ( nr, ) f0 ( nr, ) f nr, =, if f ( nr, ) U, =,,, p (0) (, ) (, ), if f0 ( nr, ) 0 0 U nr nr Here U ( nr, ) is a strictly monotonic decreasing function it respect to Te membersip functions in Equation () i.e., µ f0 ( nr, ), =,,, p Terefore te general aggregation function can be defined as { p } ( nr, ) = ( f ( nr, )), ( f ( nr, )),, f ( nr, ) 0 0 0 µ µ µ µ µ D D p Te fuzzy multi-obective formulation of te problem can be defined as: ( nr) 0 0 n r n r = C0 = C0 = C0 = C0 f nr. 0, Max µ, D Subect to ( c c W ) c t t ; n, r 0 and =,,, p. () 78

Te problem to find te optimal values of (, ) Safiulla et al. n r for tis convex-fuzzy decision based on addition operator (like Tiari et. al. [9]). Terefore te problem () is reduced according to max-addition operator as: Te above problem () reduces to Te problem (3) maximizes if te function F ( n) ( 0 (, )) p p U f nr Max µ D( n, r ) = µ ( f0 ( nr, )) = = = U Subect to fq ( nr, ) ; 0 µ ( f0 ( nr, )) n, r 0 and =,,, p. ( f0 ( nr, )) p U Max µ D ( n, r ) = = U U Subect to fq ( nr, ) ; n, r 0 and =,,, p. ( f0 ( nr, )) o = U problem (3) reduces into te primal problem (4) define as: p Min Fo ( nr, ) = 0 Subect to fq ( nr, ) ; n, r 0 and =,,, Te dual form of te Primal GPP ic is stated in (6) can be given as: i i k k d i q [ ] i Max v = i ( i) q 0 i q [ ] = i q= i q [ ] Subect i = ( ii) i [ 0] k pii iii q i q [ ] i 0, q,,, k and i =,,, mk ( iv ) () (3) attain te minimum values. Terefore te Te optimal values of sample sizes of te respondents ( n ) and non-respondents ( ) it te elp of te primal-dual relationsip teorem (9). 6. Numerical (4) (5) r can be calculated In Table, te stratum sizes, stratum eigts, stratum standard deviations, measurement costs and te travel costs itin te stratum are given for to caracteristics under study in a population stratified in four strata. Te data are mainly from Kan et al. [8]. Te travelling costs t 0 and t are assumed. Te total budget available for te survey is taken as C 0 = 5000. Te relative values of te variances of te non-respondents and respondents, tat is S S is assumed to be constant and equal to 0.5 for =, and =,,3,4. Hoever, tese ratios may vary from stratum to stratum and from caracteristic to caracteristic and can be andled accordingly. For solving MOGPP by using fuzzy programming, e sall first solve te to sub-problems: 79

Safiulla et al. Table. Data for four Strata and to caracteristics. N S S c c c t t 0 0 4 487.7 8.5 0.7 0.30 3 0.5 8 65.6 763.5 0.80 0.0 3 4 0.5.5 3 08 3066.6 456.4 0.75 0.5 4 5 0.5 3 4 786 607.5 6977.7 0.7 0.8 5 6 0.5 4.5 Sub problem : On substituting te table values in sub-problem, e ave obtained te expressions given belo: 456. 3344 6. 8965 09. 559 30. 9097 Min f0 = n n n3 n 4. 000688. 75680566 3. 49547875 4. 866484 r r r3 r 4 Subect to 0.00048n 0.00068n 0.0008n3 0.0009n 4 0.0006r 0.0008r 0.00r3 0.00r4 0.000 n 0.000 n 0.000 n3 0.000 n 4 0.0004 r 0.0005 r 0.0006 r3 0.0009 r4 n 0, r 0, n, r are integers; =,,, Te dual of te above problem (6) is obtained as: 0 03 0 04 v( 0 i ) = (( 0 ) ) (( 0 ) ) (( 03 ) ) (( 04 ) ) 05 06 07 08 ((.0007 05 ) ) ((.756806 06 ) ) (( 3.49548 07 ) ) (( 4.866484 08 ) ) Max 456.3344 6.8965 09.559 30.9097 0.00048 0.00068 3 4 0.0008 5 0.0009 0.0006 3 4 5 6 7 8 9 0 0.0008 0.00 0.00 0.000 0.000 0.000 6 7 8 9 0 3 4 5 6 0.000 0.0004 0.0005 0.0006 0.0009 3 4 5 6 ^ (( 3 4 5 6 7 8 9 0 3 4 5 6 ) ( 3 4 5 6 7 8 9 0 3 4 5 6 )); ( i) Subect to 0 0 03 04 05 06 07 08 = ; ( normality condition) ( ii) 0 9 0 0 03 3 04 4 ( ortogonality condition) 05 5 ( ) 3 06 6 4 07 7 5 ( ) 0 ( iii) 08 8 6 = 0, 0, 03, 04, 05, 06, 07, 08 > 0; ( positivity condition) ( iv),, 3, 4, 5, 6, 7, 8, 9, 0,,, 3, 4, 5, 6 0 (6) (7) 80

Safiulla et al. For ortogonality condition defined in expression 7(iii) are evaluated it te elp of te payoff matrix ic is defined belo 0 0 03 04 05 06 07 08 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 8 9 0 3 4 5 6 0 9 0 0 03 3 04 4 05 5 3 06 6 4 07 7 5 08 8 6 Solving te above formulated dual problem (7) it te elp of ingo softare, e ave te corresponding dual solutions as follos: 0.306558, 0.0798, 03.07047, 04.70049, 05.04094539,.037640,.098686,.0383407, and v = 4.75039. 06 07 08 Using te primal dual-relationsip teorem (9), e ave te optimal solution of primal problem: i.e., te optimal sample sizes of respondents and non respondents are computed as follos: (, ) = f nr v 0 0i 0i In expression (6), e first keep te r constant and calculate te values of n as: 8

Safiulla et al. ( i) ( i) ( i) ( i) f0 n, r v 0 n 474 f0 n, r v 0 n 30 f n, r = v n 49 f n, r = v n 44 03 3 03 0 3 04 4 04 0 4 No, from te expression (6), e keep te n constant and calculate te values of r as: f0 ( nr, ) v( 0 i) r 65 f0 ( nr, ) v( 0 i) r 8 f nr, = v r 8 f nr, = v r 3 ( i) ( i) 03 3 03 0 3 04 4 04 0 4 Te optimal values and te obective function value are given belo: n = 474, n = 30, n3 = 49 and n4 = 44; r = 65, r = 8, r = 8 and r = 3 and te obective value of te primal problem is 4.75039. 3 4 Sub problem : On substituting te table values in sub-problem, e ave obtained te expressions given belo: 769. 353 38. 9684 99. 53584 59. 57 Min f0 = n n n3 n 4 8.796 3.357563.65893065 5.4705348 ; r r r3 r4 Subect to 0.00048n 0.00068n 0.0008n3 0.0009n4 0.0006r 0.0008r 0.00r3 0.00r4 0.000 n 0.000 n 0.000 n3 0.000 n4 0.0004 r 0.0005 r 0.0006 r 3 0.0009 r4 n 0, r 0, n, r are integers; =,,, Te dual of te above problem (8) is obtained as follos: 0 03 0 04 v( 0 i ) = (( 0 ) ) (( 0 ) ) (( 03 ) ) (( 04 ) ) 05 06 07 08 (( 8.796 / 05 ) ) (( 3.357563 / 06 ) ) ((.65893065 / 07 ) ) (( 5.4705348 /08 ) ) Max 769.353 38.9684 99.53584 59.57 0.00048 0.00068 3 4 5 6 0. 0008 0.0009 0.0006 0.0008 3 4 5 6 7 8 9 0.00 0.00 0.000 0 0.000 0.000 0.000 7 8 9 0 3 4 5 6 0.0004 0.0005 0.0006 0.0009 3 4 5 6 ^ (( 3 4 5 6 7 8 9 0 3 4 5 6 ) ( 3 4 5 6 7 8 9 0 3 4 5 6 )); ( i) Subect to 0 0 03 04 05 06 07 08 = ; ( normality condition) ( ii) 0 9 0 0 03 3 04 4 05 5 ( ) 3 06 6 4 07 7 5 ( ) 0 ( ortogonality condition) ( iii) 08 8 6 = 0, 0, 03, 04, 05, 06, 07, 08 > 0; ( positivity condition) ( iv),, 3, 4, 5, 6, 7, 8, 9, 0,,, 3, 4, 5, 6 0 (8) (9) 8

Safiulla et al. For ortogonality condition defined in expression 9(iii) are evaluated it te elp of te payoff matrix ic is defined belo: 0 0 03 04 05 06 07 08 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 8 9 0 3 4 5 6 0 9 0 0 03 3 04 4 05 5 3 06 6 4 07 7 5 08 8 6 Solving te above formulated dual problems, e ave te corresponding solution as:.854095,.87303,.35806,.94348,.04989059, 0 0 03 04 05 06 07 08.05057448,.049875,.0974803, and v = 4. 59398. Te optimal values of sample sizes of respondents and non-respondents (, ) n r can be calculated it te 83

Safiulla et al. elp of te primal-dual relationsip teorem (9) as e ave calculated in te sub-problem are given as follos: n = 587, n = 37, n3 = 63 and n4 = 46; r = 8, r = 9, r = 8 and r = 3 and te obective value of te primal problem is 4.59398. 3 4 No te pay-off matrix of te above problems is given belo: (, ) (, ) f nr f nr 0 0 ( ) ( ) ( n, r ) ( ) ( ) ( n, r ) 4.75039 4.79368 4.40394 4.59398 f Te loer and upper bond of (, ) and (, ) f nr f nr can be obtained from te pay-off matrix 0 0 4.75039 f nr, 4.40394 and 4.59398 f nr, 4.79368. 0 0 et µ ( nr, ) and µ ( nr, ) be te fuzzy membersip function of te obective function f0 (, ) ( nr ) respectively and tey are defined as: 0, ( nr), if f0, 4.75039 4.40394 f0 ( nr, ) µ ( nr, ) =, if 4.75039 f0 ( nr, ) 4.40394 0.855 0, if f0 ( nr, ) 4.40394 ( nr), if Z, 4.59398 4.79368 Z ( nr, ) µ ( nr, ) =, if 4.59398 Z( nr, ) 4.79368 0.995 0, if Z ( nr, ) 4.79368 nr and On applying te max-addition operator, te MOGPP, te standard primal problem reduces to te problem as: f0 n,r f0 n,r Maximize 43. 347 0. 855 0. 995 Subect to 0.00048n 0.00068n 0.0008n 0.0009n 0.0006r 0.0008r 0.00r 0.00r 3 4 3 4 0.000 n 0.000 n 0.000 n3 0.000 n4 0.0004 r 0.0005 r 0.0006 r3 0.0009 r4 n 0, r 0, n, r are integers; =,,, (0) In order to maximize te above problem, e ave to minimize f0 nr f0 nr, subect to te con- 0.855 0.995 straints as described belo: (, ) (, ) 84

Safiulla et al. 5859.8847 748.334 47. 67 34.3744 n n n3 n 4 Min 4.5377 8.987 3.67 48.776 r r r3 r 4 Subect to 0.00048n 0.00068n 0.0008n3 0.0009n 4 0.0006r 0.0008r 0.00r3 0.00r 4 0.000 n 0.000 n 0.000 n3 0.000 n4 0.0004 r 0.0005 r 0.0006 r3 0.0009 r4 n 0, r 0, n, r are integers; =,,, () Degree of Difficulty of te problem () is = (4 (8 ) =5. Hence te dual problem of te above final formulated problem () is given as: 0 03 0 04 v( 0 i ) = (( 0 ) ) (( 0 ) ) (( 03 ) ) (( 04 ) ) 05 06 07 08 (( 4.5377 05 ) ) (( 8.987 06 ) ) (( 3.67 07 ) ) (( 48.776 08 ) ) Max 5859.8847 748.334 47.67 34.3744 0.00048 0.00068 3 4 5 0.0008 6 0.0009 0.0006 0.0008 3 4 5 6 7 8 9 0 0.00 0.00 0.000 0.000 0.000 0.000 7 8 9 0 3 4 5 6 0.0004 0.0005 0.0006 0.0009 3 4 5 6 (( ) 3 4 5 6 7 8 9 0 3 4 5 6 ( 3 4 5 6 7 8 9 0 3 4 5 6 )); ( i) Subect to 0 0 03 04 05 06 07 08 = ; ( normality condition) ( ii) 0 9 0 0 03 3 04 4 05 5 3 06 6 4 07 7 5 ( ) 0 ( ortogonality condition) ( iii) 08 8 6 = 0, 0, 03, 04, 05, 06, 07, 08 > 0;,, 3, 4, 5, 6, 7, 8, 9, 0,,, 3, 4, 5, 6 0 ( positivity condition) ( iv) For ortogonality condition defined in expression (iii) are evaluated it te elp of te payoff matrix ic is defined belo: ^ () 85

Safiulla et al. 0 0 03 04 05 06 07 08 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 0 7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( ) 8 9 0 3 4 5 6 0 9 0 0 03 3 04 4 05 5 3 06 6 4 07 7 5 08 8 6 After solving te formulated dual problem () using lingo softare e obtain te folloing values of te dual variables ic are given as:.60739,.7708,.657664,.7033,.0469584, 0 0 03 04 05 06.043073, 07.046037, 08.0383880, and v( ) = 4.83433. Te optimal values of sample sizes of respondents and non-respondents (, ) n r can be calculated it te elp of te primal-dual relationsip teorem (9) as e ave calculated in te sub-problem are given as follos: 86

Safiulla et al. n = 537, n = 309, n = 04 and n = 44; r = 74, r = 8, r = 3 and r = 3 3 4 3 4 and te obective value of te primal problem is 4.83433. 7. Conclusion Tis paper provides an insigtful study of fuzzy programming for solving te multi-obective geometric programming problem (MOGPP). Te problem of non-response it significant travel costs ere te cost is quadratic in n in multivariate stratified sample surveys as been formulated of as a Multi-Obective Geometric Programming Problem (MOGPP). Te fuzzy programming approac as described for solving te formulated MOGPP. Te formulated MOGPP as been solved it te elp of INGO Softare [30] and te dual solution is obtained. Te optimum allocations of sample sizes of respondents and non respondents are obtained it te elp of dual solutions and primal-dual relationsip teorem. To ascertain te practical utility of te proposed metod in sample surveys problem in presence of non-response it significant travel cost ere te cost is quadratic in n, a numerical example is also given to illustrate te procedure. References [] Cocran, W.G. (977) Sampling Tecniques. 3rd Edition, Jon Wiley and Sons, Ne York. [] Sukatme, P.V., Sukatme, B.V., Sukatme, S. and Asok, C. (984) Sampling Teory of Surveys it Applications. Ioa State University Press, Ames and Indian Society of Agricultural Statistics, Ne Deli. [3] Hansen, M.H. and Huritz, W.N. (946) Te Problem of Non-Response in Sample Surveys. Journal of te American Statistical Association, 4, 57-59. ttp://dx.doi.org/0.080/06459.946.050894 [4] El-Badry, M.A. (956) A Sampling Procedure for Mailed Questionnaires. Journal of te American Statistical Association, 5, 09-7. ttp://dx.doi.org/0.080/06459.956.0503 [5] Fordori, G.T. (96) Some Non-Response Sampling Teory for To Stage Designs. Mimeograp Series No. 97, Nort Carolina State University, Raleig. [6] Srinat, K.P. (97) Multiple Sampling in Non-Response Problems. Journal of te American Statistical Association, 66, 583-586. ttp://dx.doi.org/0.080/06459.97.04830 [7] Kare, B.B. (987) Allocation in Stratified Sampling in Presence of Non-Response. Metron, 45, 3-. [8] Kan, M.G.M., Kan, E.A. and Asan, M.J. (008) Optimum Allocation in Multivariate Stratified Sampling in Presence of Non-Response. Journal of te Indian Society of Agricultural Statistics, 6, 4-48. [9] Varsney, R., Asan, M.J. and Kan, M.G.M. (0) An Optimum Multivariate Stratified Sampling Design it Non- Response: A exicograpic Goal Programming Approac. Journal of Matematical Modeling and Algoritms, 65, 9-96. [0] Fatima, U. and Asan, M.J. (0) Nonresponse in Stratified Sampling: A Matematical Programming Approac. Te Sout Pacific Journal of Natural and Applied Sciences, 9, 40-4. [] Varsney, R., Namussear and Asan, M.J. (0) An Optimum Multivariate Stratified Double Sampling Design in Presence of Non-Response. Optimization etters, 6, 993-008. [] Ragav, Y.S., Ali, I. and Bari, A. (04) Multi-Obective Nonlinear Programming Problem Approac in Multivariate Stratified Sample Surveys in Case of Non-Response. Journal of Statistical Computation and Simulation, 84, -36. ttp://dx.doi.org/0.080/00949655.0.69370 [3] Duffin, R.J., Peterson, E.. and Zener, C. (967) Geometric Programming: Teory & Applications. Jon Wiley & Sons, Ne York. [4] Zener, C. (97) Engineering Design by Geometric Programming. Jon Wiley & Sons, Ne York. [5] Beigtler, C.S. and Pilips, D.T. (976) Applied Geometric Programming. Wiley, Ne York. [6] Davis, M. and Scartz, R.S. (987) Geometric Programming for Optimal Allocation of Integrated Samples in Quality Control. Communications in Statistics Teory and Metods, 6, 335-354. ttp://dx.doi.org/0.080/03609870889568 [7] Amed, J. and Bonam, C.D. (987) Application of Geometric Programming to Optimum Allocation Problems in Multivariate Double Sampling. Applied Matematics and Computation,, 57-69. ttp://dx.doi.org/0.06/0096-3003(87)9004-5 [8] Oa, A.K. and Das, A.K. (00) Multi-Obective Geometric Programming Problem Being Cost Coefficients as Con- 87

Safiulla et al. tinuous Function it Weigted Mean. Journal of Computing,, 5-967. [9] Maqbool, S., Mir, A.H. and Mir, S.A. (0) Geometric Programming Approac to Optimum Allocation in Multivariate To-Stage Sampling Design. Electronic Journal of Applied Statistical Analysis, 4, 7-8. [0] Safiulla, Ali, I. and Bari, A. (03) Geometric Programming Approac in Tree Stage Sampling Design. International Journal of Scientific & Engineering Researc (France), 4, 9-558. [] Zade,.A. (965) Fuzzy Sets. Information and Control, 8, 338-353. ttp://dx.doi.org/0.06/s009-9958(65)904-x [] Bellman, R.E. and Zade,.A. (970) Decision-Making in a Fuzzy Environment. Management Sciences, 7, B4- B64. ttp://dx.doi.org/0.87/mnsc.7.4.b4 [3] Tanaka, H., Okuda, T. and Asai, K. (974) On Fuzzy-Matematical Programming. Journal of Cybernetics, 3, 37-46. ttp://dx.doi.org/0.080/0969773085459 [4] Zimmermann, H.J. (978) Fuzzy Programming and inear Programming it Several Obective Functions. Fuzzy Sets and Systems,, 45-55. ttp://dx.doi.org/0.06/065-04(78)9003-3 [5] Bisal, M.P. (99) Fuzzy Programming Tecnique to Solve Multi-Obective Geometric Programming Problems. Fuzzy Sets and Systems, 5, 67-7. ttp://dx.doi.org/0.06/065-04(9)90076-g [6] Verma, R.K. (990) Fuzzy Geometric Programming it Several Obective Functions. Fuzzy Sets and Systems, 35, 5-0. ttp://dx.doi.org/0.06/065-04(90)9004-z [7] Islam, S. and Roy, T.K. (005) Modified Geometric Programming Problem and Its Applications. Journal of Applied Matematics and Computing, 7, -44. [8] Islam, S. (00) Multi-Obective Geometric Programming Problem and Its Applications. Yugoslav Journal of Operations Researc, 0, 3-7. ttp://dx.doi.org/0.98/yjor003i [9] Tiari, R.N., Darman, S. and Rao, J.R. (987) Fuzzy Goal Programming An Additive Model. Fuzzy Sets and Systems, 4, 7-34. ttp://dx.doi.org/0.06/065-04(87)90-4 [30] indo Systems Inc. (03) INGO User s Guide. indo Systems Inc., Cicago. 88