Ratio Estimator in Adaptive Cluster Sampling without Replacement of ...

0 downloads 0 Views 2MB Size Report
May 6, 2014 - Adaptive cluster sampling, proposed by Thompson [1], is an efficient method for sampling rare and hidden clustered populations. In adaptiveΒ ...
Hindawi Publishing Corporation Journal of Probability and Statistics Volume 2014, Article ID 726398, 6 pages http://dx.doi.org/10.1155/2014/726398

Research Article Ratio Estimator in Adaptive Cluster Sampling without Replacement of Networks Nipaporn Chutiman and Monchaya Chiangpradit Department of Mathematics, Faculty of Science, Mahasarakham University, Maha Sarakham 44150, Thailand Correspondence should be addressed to Nipaporn Chutiman; [email protected] Received 6 January 2014; Revised 6 May 2014; Accepted 6 May 2014; Published 27 May 2014 Academic Editor: Shein-chung Chow Copyright Β© 2014 N. Chutiman and M. Chiangpradit. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. In this paper, we study the estimators of the population total in adaptive cluster sampling by using the information of the auxiliary variable. The numerical examples showed that the ratio estimator in adaptive cluster sampling without replacement of networks is more efficient than the ratio estimators in adaptive cluster sampling without replacement of units.

1. Introduction Adaptive cluster sampling, proposed by Thompson [1], is an efficient method for sampling rare and hidden clustered populations. In adaptive cluster sampling, an initial sample is selected by simple random sampling without replacement of units (Thompson [1]). The unbiased estimators for adaptive cluster sampling with an initial sample taken by simple random sampling without replacement of units are HansenHurwitz and Horvitz-Thompson estimator. Although the initial sample is selected by simple random sampling without replacement of units, some networks in the sample may be selected more than once. An adaptive cluster sampling design exists with selection without replacement of networks (Salehi and Seber [2]). The unbiased estimator for adaptive cluster sampling without replacement of networks is Des Raj estimator. The estimators mentioned above are provided by the variable of interest (𝑦) only. Sometimes other variables are related to the variable of interest. We can obtain additional information for estimating the population total. The use of an auxiliary variable is a common method to improve the precision of estimates of a population total. Dryver and Chao [3] proposed ratio estimators based on Hansen-Hurwitz and Horvitz-Thompson estimator in adaptive cluster sampling without replacement of units. In this study, we will study

the estimator of population total in adaptive cluster sampling without replacement of networks using an auxiliary variable. Some comparisons are made using a simulation.

2. Adaptive Cluster Sampling without Replacement of Units As in the finite population sampling situation, the population consists of 𝑁 units (𝑒1 , 𝑒2 , . . . , 𝑒𝑁) indexed by their labels (1, 2, . . . , 𝑁). Unit 𝑖 is associated as a variable of interest 𝑦𝑖 . The vector of the population total 𝑦-values may be written as y = {𝑦1 , 𝑦2 , . . . , 𝑦𝑁}. The estimator of the population total πœπ‘¦ = βˆ‘π‘ 𝑖=1 𝑦𝑖 is studied. For adaptive cluster sampling, an initial sample of units is selected by simple random sampling without replacement. The unbiased estimator of the population total of the variable of interest is formed by using Hansen-Hurwitz estimator and Horvitz-Thompson estimator. 2.1. The Hansen-Hurwitz Estimator. The Hansen-Hurwitz estimator is based on draw-by-draw probabilities that a unit’s network is intersected by the initial sample. Let 𝑛 denote the initial sample size and ] denote the final sample size. Let πœ“π‘– denote the network that includes unit 𝑖 and let π‘šπ‘– be the number of units in that network.

2

Journal of Probability and Statistics

The Hansen-Hurwitz estimator of the population total of the variable of interest can be written as (Thompson and Seber [4]) (Μ‚ πœπ‘¦ )HH =

𝑁 𝑛 βˆ‘(𝑀 ) , 𝑛 𝑖=1 𝑦 𝑖

(1)

where (𝑀𝑦 )𝑖 is the average of the 𝑦-values in the network that includes unit 𝑖 of the initial sample; that is, (𝑀𝑦 )𝑖 =

2.3. Ratio Estimators in Adaptive Cluster Sampling. Dryver and Chao [3] proposed the ratio estimator in adaptive cluster sampling based on Hansen-Horvitz’s ratio estimator as

1 βˆ‘π‘¦ . π‘šπ‘– π‘—βˆˆπœ“π‘– 𝑗

(2)

(Μ‚ πœπ‘¦ )HH R =

𝑉(Μ‚ πœπ‘¦ )HH

𝑉(Μ‚ πœπ‘¦ )HH R = (3)

Μ‚ πœπ‘¦ ) = 𝑁 (𝑁 βˆ’ 𝑛) βˆ‘((𝑀𝑦 ) βˆ’ 𝑉(Μ‚ HH 𝑖 𝑛 (𝑛 βˆ’ 1) 𝑖=1

(Μ‚ πœπ‘¦ )HH 𝑁

2

).

(4)

2 𝑁 (𝑁 βˆ’ 𝑛) 𝑁 βˆ‘((𝑀𝑦 )𝑖 βˆ’ 𝑅(𝑀π‘₯ )𝑖 ) . 𝑛 (𝑁 βˆ’ 1) 𝑖=1

𝑁 (𝑁 βˆ’ 𝑛) 𝑛 Μ‚ πœπ‘¦ ) Μ‚ HH (𝑀π‘₯ ) )2 . 𝑉(Μ‚ = βˆ‘((𝑀𝑦 )𝑖 βˆ’ 𝑅 𝑖 HH R 𝑛 (𝑛 βˆ’ 1) 𝑖=1

πœ… π‘¦π‘˜. πΌπ‘˜ 𝑦 = βˆ‘ π‘˜. , π›Όπ‘˜ 𝛼 π‘˜=1 π‘˜=1 π‘˜

(Μ‚ πœπ‘¦ )HT R =

(Μ‚ πœπ‘¦ )HT (Μ‚ 𝜏π‘₯ )HT

(5)

where π‘¦π‘˜. is the sum of 𝑦-values in network π‘˜; that is, π‘¦π‘˜. = βˆ‘π‘–βˆˆπœ“π‘˜ 𝑦𝑖 . 𝐾 is the total number of distinct networks in the population and πœ… is the total number of distinct networks in the sample. πΌπ‘˜ is an indicator variable that takes a value of 1 with probability π›Όπ‘˜ if the initial sample intersects network π‘˜ and 0 otherwise. The inclusion probability of network π‘˜ is 𝑁 π‘˜ π›Όπ‘˜ = 1 βˆ’ [( π‘βˆ’π‘š 𝑛 ) / ( 𝑛 )], and the probability that networks 𝑗 and π‘˜ are both intersected is π‘βˆ’π‘šπ‘— βˆ’π‘šπ‘˜ π‘βˆ’π‘šπ‘˜ 𝑗 )] [( π‘βˆ’π‘š 𝑛 )+( 𝑛 )βˆ’( 𝑛

(𝑁 𝑛 )

}.

(6)

Μ‚ HT = 𝑅 (Μ‚ 𝜏π‘₯ )HT

𝑉(Μ‚ πœπ‘¦ )HT = βˆ‘ βˆ‘

𝛼𝑗 π›Όπ‘˜

𝑗=1 π‘˜=1

(12)

(Μ‚ πœπ‘¦ )HT (Μ‚ 𝜏π‘₯ )HT

, (13)

πœ…

π‘₯ = βˆ‘ ( π‘˜. ) . π›Όπ‘˜ π‘˜=1

𝐾 𝐾

𝑉(Μ‚ πœπ‘¦ )HT R = βˆ‘ βˆ‘

𝑒𝑗. π‘’π‘˜. (π›Όπ‘—π‘˜ βˆ’ 𝛼𝑗 π›Όπ‘˜ ) 𝛼𝑗 π›Όπ‘˜

𝑗=1 π‘˜=1

,

πœ…

Μ‚ πœπ‘¦ ) = βˆ‘ βˆ‘ 𝑉(Μ‚ HT

𝑗=1 π‘˜=1

𝑦𝑗. π‘¦π‘˜. π›Όπ‘—π‘˜

(

π›Όπ‘—π‘˜ 𝛼𝑗 π›Όπ‘˜

(14)

where (15)

The unbiased estimator of 𝑉(Μ‚ πœπ‘¦ )HT R is ,

(7)

πœ…

πœ…

Μ‚ πœπ‘¦ ) 𝑉(Μ‚ = βˆ‘βˆ‘ HT R

𝑗=1 π‘˜=1

and the unbiased estimator of 𝑉(Μ‚ πœπ‘¦ )HT is πœ…

Μ‚ HT 𝜏π‘₯ , 𝜏π‘₯ = 𝑅

𝑒𝑗. = 𝑦𝑗. βˆ’ 𝑅π‘₯𝑗. . 𝑦𝑗. π‘¦π‘˜. (π›Όπ‘—π‘˜ βˆ’ 𝛼𝑗 π›Όπ‘˜ )

(11)

The approximated variance of (Μ‚ πœπ‘¦ )HT R is

The variance of (Μ‚ πœπ‘¦ )HT is 𝐾

(10)

where

𝐾

(Μ‚ πœπ‘¦ )HT = βˆ‘

𝐾

(9)

The ratio estimator in adaptive cluster sampling based on Horvitz-Thompson’s ratio estimator is

2.2. The Horvitz-Thompson Estimator. The Horvitz-Thompson estimator is based on probabilities of the initial sample intersecting networks. The unbiased estimator of the population total of the variable of interest can be written as (Thompson and Seber [4])

π›Όπ‘—π‘˜ = 1 βˆ’ {

Μ‚ HH 𝜏π‘₯ , 𝜏π‘₯ = 𝑅

The unbiased estimator of 𝑉(Μ‚ πœπ‘¦ )HH R is

The unbiased estimator of 𝑉(Μ‚ πœπ‘¦ )HH is 𝑛

(Μ‚ 𝜏π‘₯ )HH

Μ‚ HH = (Μ‚ πœπ‘¦ )HH /(Μ‚ 𝜏π‘₯ )HH , (Μ‚ 𝜏π‘₯ )HH = (𝑁/𝑛) βˆ‘π‘›π‘–=1 (𝑀π‘₯ )𝑖 , where 𝑅 and (𝑀π‘₯ )𝑖 is the average of the auxiliary variable π‘₯ in the network that includes unit 𝑖 of the initial sample; that is, (𝑀π‘₯ )𝑖 = (1/π‘šπ‘– ) βˆ‘π‘—βˆˆπœ“π‘– π‘₯𝑗 . The approximated variance of (Μ‚ πœπ‘¦ )HH R is

The variance of (Μ‚ πœπ‘¦ )HH is πœπ‘¦ 2 𝑁 (𝑁 βˆ’ 𝑛) 𝑁 = βˆ‘((𝑀𝑦 )𝑖 βˆ’ ) . 𝑛 (𝑁 βˆ’ 1) 𝑖=1 𝑁

(Μ‚ πœπ‘¦ )HH

σΈ€  𝑒𝑗.σΈ€  π‘’π‘˜.

π›Όπ‘—π‘˜

(

π›Όπ‘—π‘˜ 𝛼𝑗 π›Όπ‘˜

βˆ’ 1) ,

(16)

where βˆ’ 1) .

(8)

Μ‚ HT π‘₯𝑗. . 𝑒𝑗.σΈ€  = 𝑦𝑗. βˆ’ 𝑅

(17)

Journal of Probability and Statistics

3

3. Adaptive Cluster Sampling without Replacement of Networks Although the initial sample is selected by simple random sampling without replacement of units, some networks in the sample may be selected more than once. An adaptive cluster sampling design exists with selection without replacement of networks (Salehi and Seber [2]). Salehi and Seber [2] proposed a new sampling design, that is, adaptive cluster sampling with networks selected without replacement. Their procedure is as follows. The first initial sample unit is selected by simple random sampling from the population. Then the first network is created and this network is removed from the population. Next, the second initial unit is selected by simple random sampling without replacement of the remaining units; then a second network is created. And finally, the process is continued until all of 𝑛 networks are selected. Let 𝑛 denote the initial sample size. Let πœ“π‘– denote the network that includes unit 𝑖 and let π‘šπ‘– be the number of units in that network. Further, let 𝑝𝑖 be the first-draw probabilities for network which includes unit 𝑖. Thus 𝑝𝑖 = π‘šπ‘– /𝑁, where π‘šπ‘– is the number of units in the network which includes unit 𝑖. So 𝑝𝑖 /(1 βˆ’ βˆ‘π‘–βˆ’1 𝑗=1 𝑝𝑗 ) is the conditional 𝑖th draw probability for the 𝑖th network which includes unit 𝑖 in the sample given the first 𝑖 βˆ’ 1 networks selection. By using a modified Des Raj estimator (Salehi and Seber [2]), an unbiased estimator for the population total of the variable of interest is (Μ‚ πœπ‘¦ )Raj =

1 𝑛 βˆ‘(𝑧 ) , 𝑛 𝑖=1 𝑦 𝑖

4. Proposed Ratio Estimator in Adaptive Cluster Sampling without Replacement of Networks Some other variables are related to the variable of interest. We can obtain additional information for estimating the population total. The use of an auxiliary variable is a common method to improve the precision of estimates of a population total. In this section, proposed ratio estimator in adaptive cluster sampling without replacement of networks will be described. The proposed ratio estimator based on Des Raj’s estimator is (Μ‚ πœπ‘¦ )Raj R =

(𝑧𝑦 )𝑖 = βˆ‘π‘¦π‘—. + 𝑗=1

(1 βˆ’ βˆ‘π‘–=1 𝑗=1 𝑝𝑗 ) 𝑦𝑖. 𝑝𝑖

Μ‚ Raj βˆ’ 𝑅 = 𝑅

(Μ‚ 𝜏π‘₯ )Raj

Μ‚ Raj βˆ’ 𝑅 β‰ˆ 𝑅

(22)

βˆ’π‘…=

(Μ‚ πœπ‘¦ )Raj βˆ’ 𝑅(Μ‚ 𝜏π‘₯ )Raj (Μ‚ 𝜏π‘₯ )Raj

,

(Μ‚ πœπ‘¦ )Raj βˆ’ 𝑅(Μ‚ 𝜏π‘₯ )Raj 𝜏π‘₯

(23)

,

2

Μ‚ Raj ) = 𝐸[𝑅 Μ‚ Raj βˆ’ 𝐸 (𝑅 Μ‚ Raj )] β‰ˆ 𝐸[𝑅 Μ‚ Raj βˆ’ 𝑅] 𝑉 (𝑅

2

2

(24)

] = 𝐸[ 𝜏π‘₯ [ ] 2 1 = 2 𝐸[(Μ‚ πœπ‘¦ )Raj βˆ’ 𝑅(Μ‚ 𝜏π‘₯ )Raj ] . 𝜏π‘₯

(19) Let

1 𝑛 βˆ‘π‘‰(𝑧𝑦 )𝑖 , 𝑛2 𝑖=1

(20)

𝑑𝑖 = (𝑧𝑦 )𝑖 βˆ’ 𝑅(𝑧π‘₯ )𝑖 , 𝑛

𝑛

(25) 1 1 1 𝑛 πœπ‘¦ )Raj βˆ’ 𝑅(Μ‚ 𝜏π‘₯ )Raj . βˆ‘π‘‘π‘– = βˆ‘(𝑧𝑦 )𝑖 βˆ’ 𝑅 βˆ‘(𝑧π‘₯ )𝑖 = (Μ‚ 𝑛 𝑖=1 𝑛 𝑖=1 𝑛 𝑖=1 So, 2

𝑛 Μ‚ Raj ) = 1 𝐸( 1 βˆ‘π‘‘π‘– ) , 𝑉 (𝑅 𝜏π‘₯2 𝑛 𝑖=1

and the unbiased estimator of 𝑉(Μ‚ πœπ‘¦ )Raj is Μ‚ πœπ‘¦ ) = 𝑉(Μ‚ Raj

(Μ‚ πœπ‘¦ )Raj

(Μ‚ πœπ‘¦ )Raj βˆ’ 𝑅(Μ‚ 𝜏π‘₯ )Raj

From Raj [5] we find that (𝑧𝑦 )𝑖 are mutually uncorrelated so that the variance of (Μ‚ πœπ‘¦ )Raj (Salehi and Seber [2]) is 𝑉(Μ‚ πœπ‘¦ )Raj =

Μ‚ Raj 𝜏π‘₯ , 𝜏π‘₯ = 𝑅

where 𝑅 = πœπ‘¦ /𝜏π‘₯ is the population ratio. Μ‚ Raj ) If 𝑛 is reasonably large, (Μ‚ 𝜏π‘₯ )Raj is closed to 𝜏π‘₯ and 𝐸(𝑅 is approximately equal to 𝑅; therefore,

(18)

.

(Μ‚ 𝜏π‘₯ )Raj

Μ‚ Raj = (Μ‚ πœπ‘¦ )Raj /(Μ‚ 𝜏π‘₯ )Raj and (Μ‚ 𝜏π‘₯ )Raj is the modified Des where 𝑅 Raj estimator for 𝜏π‘₯ . To obtain the variance of (Μ‚ πœπ‘¦ )Raj R , we write

where, for 𝑖 = 1, (𝑧𝑦 )𝑖 = 𝑦1. /𝑝1 and for 𝑖 = 2, 3, . . . , 𝑛 let π‘–βˆ’1

(Μ‚ πœπ‘¦ )Raj

2

2

𝑛 1 πœπ‘¦ )Raj ) . βˆ‘ (𝑧𝑖𝑦 βˆ’ (Μ‚ 𝑛 (𝑛 βˆ’ 1) 𝑖=1

(21)

The variance of the modified Des Raj estimator, 𝑉(Μ‚ πœπ‘¦ )Raj , is less than that of the Hansen-Hurwitz estimator, 𝑉(Μ‚ πœπ‘¦ )HH (Salehi and Seber [2]).

1 𝑛 1 𝑛 1 𝑛 𝑉 [ βˆ‘π‘‘π‘– ] = 𝐸( βˆ‘π‘‘π‘– ) βˆ’ [𝐸 ( βˆ‘π‘‘π‘– )] 𝑛 𝑖=1 𝑛 𝑖=1 𝑛 𝑖=1 𝑛

2

1 = 𝐸( βˆ‘π‘‘π‘– ) , 𝑛 𝑖=1 𝑛 Μ‚ Raj ) = 1 Γ— 𝑉 [ 1 βˆ‘π‘‘π‘– ] . 𝑉 (𝑅 𝜏π‘₯2 𝑛 𝑖=1

2

(26)

4

Journal of Probability and Statistics

Table 1: The estimated MSE of the Hansen-Horvitz estimators and Horvitz-Thompson estimators for the population total of blue-winged teal. 𝑛

𝐸(])

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HH

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HT

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HH R

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HT R

5 8 10 15 20

17.09 23.44 26.14 31.27 34.71

2,359,716,325 1,369,069,022 1,076,080,659 575,418,109 348,333,275

1,855,257,461 814,459,237 527,029,260 163,370,814 42,718,623

444,212,079 155,321,999 79,499,876 11,895,009 1,512,773

444,211,995 155,321,710 79,499,398 11,894,258 1,511,986

Table 2: The estimated MSE of the modified Des Raj estimators for the population total of blue-winged teal. 𝑛

𝐸(])

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ Raj

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ Raj R

5 8 10 15 20

19.90 26.66 29.84 35.74 39.32

1,897,377,719 931,732,897 654,226,056 305,945,821 171,618,612

22,982,048 15,329,005 11,162,171 3,298,687 624,795

0 0 0 0 0

0 0 0 0 0

3 0 0 0 0

5 24 0 0 0

0 14 2 0 0

0 0 3 0 0

0 0 2 0 2

0 10 0 0 0

0 103 13639 14 0

0 0 1 122 177

Figure 1: Blue-winged teal data (𝑦-values).

Therefore, Μ‚ Raj 𝜏π‘₯ ) = 𝜏2 𝑉 (𝑅 Μ‚ Raj ) . 𝑉 ((Μ‚ πœπ‘¦ )Raj R ) = 𝑉 (𝑅 π‘₯

(27)

The approximated variance of (Μ‚ πœπ‘¦ )Raj R is 𝑉(Μ‚ πœπ‘¦ )Raj R =

𝑛

1 1 [ βˆ‘π‘‰ (𝑑𝑖 )] , 𝑛2 𝑛2 𝑖=1

(28)

and the estimator of 𝑉(Μ‚ πœπ‘¦ )Raj R is 𝑛 1 1 Μ‚ πœπ‘¦ ) Μ‚ Raj (𝑧π‘₯ ) )2 ] . (29) 𝑉(Μ‚ = [ ((𝑧𝑦 )𝑖 βˆ’ 𝑅 βˆ‘ 𝑖 2 Raj R 𝑛 𝑛 (𝑛 βˆ’ 1) 𝑖=1

5. Examples 5.1. The Blue-Winged Teal Example. As an example, the bluewinged teal data (Smith et al. [6]) were used in this study. The population region of 5,000 km2 in an area of Central Florida was partitioned into fifty units of 100 km2 . These data were set to be 𝑦-values. The total of 𝑦-values, πœπ‘¦ , is 14,121 (see Figure 1). The π‘₯-values were simulation results from the model π‘₯𝑖 = 4𝑦𝑖 βˆ’ πœ€π‘– , where πœ€π‘– ∼ 𝑁(0, 𝑦𝑖 ) (see Figure 2). For each iteration, an initial sample of units is selected by simple random sampling without replacement. The 𝑦-values are obtained for keeping the sample network. In each of the

0 0 0 0 0

0 0 0 0 0

13 0 0 0 0

19 93 0 0 0

0 59 9 0 0

0 0 10 0 0

0 0 8 0 7

0 37 0 0 0

0 419 45621 59 0

0 0 6 493 691

Figure 2: Simulated data (π‘₯-values) from the model π‘₯𝑖 = 4𝑦𝑖 βˆ’ πœ€π‘– .

sample networks, the π‘₯-values are obtained. The condition for added units in the sample is defined by 𝐢 = {𝑦 : 𝑦 > 0}. For each estimator, 20,000 iterations were performed to obtain an accuracy estimate. Initial SRS sizes were varied with 𝑛 = 5, 8, 10, 15, and 20. The estimated variance of the estimate total is Μ‚ (Μ‚ 𝑀𝑆𝐸 πœπ‘¦ ) =

2 1 20,000 πœπ‘¦ )𝑖 βˆ’ πœπ‘¦ ) , βˆ‘ ((Μ‚ 20, 000 𝑖=1

(30)

where (Μ‚ πœπ‘¦ )𝑖 is the value for the relevant estimator for sample 𝑖. For the blue-winged teal example, the estimated MSE of estimators under the adaptive cluster sampling without replacement of unit were calculated and listed in Table 1. The estimated MSE of estimators under adaptive cluster sampling without replacement of networks were shown in Table 2. 5.2. Simulation Example. As a second illustration, we have used the simulation π‘₯-values and 𝑦-values from Chutiman and Kumphon [7] were studied. The populations were shown in Figures 3 and 4. For each iteration, an initial sample of units is selected by simple random sampling without replacement. The 𝑦-values are obtained for keeping the sample network. In each sample network, the π‘₯-values are obtained. The condition for added units in the sample is defined by 𝐢 = {𝑦 : 𝑦 > 0}. For each estimator, 20,000 iterations were performed to obtain an accuracy estimate. Initial SRS sizes were varied 𝑛 = 5, 10, 15, 20, 30, 40, and 50. For the simulation example, the estimated MSE of estimators under the adaptive cluster sampling without replacement of unit were calculated and listed in Table 3. Table 4 shows the estimated MSE of estimators under adaptive cluster sampling without replacement of networks. From the examples, we see that the estimated MSE of adaptive cluster sampling without replacement of networks is smaller than the estimated MSE of adaptive cluster sampling without replacement of units; that is, Μ‚ πœπ‘¦ ) ≀ 𝑀𝑆𝐸(Μ‚ Μ‚ πœπ‘¦ ) . Μ‚ πœπ‘¦ ) ≀ 𝑀𝑆𝐸(Μ‚ 𝑀𝑆𝐸(Μ‚ Raj HT HH

(31)

Journal of Probability and Statistics

5

Table 3: The estimated MSE of the Hansen-Horvitz estimator and Horvitz-Thompson estimator for the population total of the variable of interest. 𝑛

𝐸(])

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HH

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HT

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HH R

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ HT R

5 10 15 20 30 40 50

28.999 50.596 68.027 82.526 104.183 120.132 133.058

407,050.393 197,774.349 130,815.296 95,153.494 61,714.960 45,383.260 33,434.958

398,429.963 187,902.766 120,108.018 84,773.793 52,527.501 36,945.984 25,939.191

113,906.753 55,691.349 28,580.328 15,778.559 7,242.755 5,263.785 4,451.101

113,871.977 55,554.041 28,341.920 15,465.085 6,890.563 4,962.297 4,167.722

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 7 1 1 0 0 0 0 0 0 0

0 0 0 1 0 0 0 0 0 0 65 4 7 0 0 0 0 0 0 0

0 0 0 22 0 0 0 0 0 5 0 5 0 0 0 0 0 0 0 0

0 0 0 5 2 0 0 0 1 4 4 0 0 0 0 0 0 0 0 0

0 0 0 4 0 0 0 0 0 3 5 7 0 0 0 0 0 0 0 0

0 0 2 5 4 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0

0 5 22 0 8 0 21 4 5 5 9 3 0 0 0 0 0 0 0 0

0 0 5 0 0 0 0 0 7 8 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 5 0 4 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 1 7 7 5 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 7 1 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 33 0 0 0 27 0 0 7 6 7 1 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6 3 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 3 1 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 27 0 0 21 0 29 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Figure 3: The 𝑦-values.

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 18 2 3 0 0 0 0 0 0 0

0 0 0 11 0 0 0 0 0 2 0 2 0 0 0 0 0 0 0 0

0 0 0 2 1 0 0 0 0 2 2 0 0 0 0 0 0 0 0 0

0 0 0 2 0 0 0 0 0 1 2 3 0 0 0 0 0 0 0 0

0 0 1 1 2 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0

0 2 11 0 4 0 16 2 2 2 4 1 0 0 0 0 0 0 0 0

0 0 3 0 0 0 0 0 2 3 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 2 0 2 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 3 3 2 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 3 0 0 2 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 12 2 0 0 1 2 1 0 0 0 0 0 0 0 27 0

Figure 4: Theπ‘₯-values.

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 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 12 0 0

0 0 0 0 15 2 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 12 0 0

6

Journal of Probability and Statistics

Table 4: The estimated MSE of the modified Des Raj estimators for the population total of the variable of interest. 𝑛

𝐸(])

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ Raj

Μ‚ πœπ‘¦ ) 𝑀𝑆𝐸(Μ‚ Raj R

5 10 15 20 30 40 50

31.590 56.182 75.244 91.807 116.093 133.820 148.388

395,296.216 181,264.543 116,781.850 83,718.824 51,340.406 35,017.241 25,772.744

8,146.702 7,641.379 7,033.291 6,417.127 5,359.412 4,606.211 4,041.559

For the population total in adaptive cluster sampling by using the information of the auxiliary variable, we see that the estimated MSE of the ratio estimator in adaptive cluster sampling without replacement of networks is smaller than the estimated MSE of the ratio estimator in adaptive cluster sampling without replacement of units; that is, Μ‚ πœπ‘¦ ) Μ‚ πœπ‘¦ ) Μ‚ πœπ‘¦ ) ≀ 𝑀𝑆𝐸(Μ‚ ≀ 𝑀𝑆𝐸(Μ‚ . 𝑀𝑆𝐸(Μ‚ Raj R HT R HH R

(32)

6. Conclusions We discuss using the auxiliary information, that is, the ratio estimators in adaptive cluster sampling without replacement of units and without replacement of networks. For the adaptive cluster sampling without replacement of units, HorvitzThompson’s ratio estimator is better than Hansen-Hurwitz’s ratio estimator (Dryver and Chao [3]). Although the initial units in adaptive cluster sampling are selected by simple random sampling without replacement, networks may be selected more than once. So the adaptive cluster sampling without replacement of units is equivalent to the adaptive cluster sampling without replacement of networks. Therefore, the modified Des Raj ratio estimator is better than both Hansen-Hurwitz’s ratio estimator and Horvitz-Thompson’s ratio estimator.

Conflict of Interests The authors declare that there is no conflict of interests regarding the publishing of this paper.

Acknowledgments The authors would like to profoundly thank Paveen Chutiman for his programming advice. They would also like to express their special thanks to Assistant Professor Siriluck Jemjitpornchai and Prasert Vangsantitrakul for their valuable comments and suggestions on the material of this paper and for their help in correcting English as well.

References [1] S. K. Thompson, β€œAdaptive cluster sampling,” Journal of the American Statistical Association, vol. 85, no. 412, pp. 1050–1059, 1990.

[2] M. M. Salehi and G. A. F. Seber, β€œAdaptive cluster sampling with networks selected without replacement,” Biometrika, vol. 84, no. 1, pp. 209–219, 1997. [3] A. L. Dryver and C. T. Chao, β€œRatio estimators in adaptive cluster sampling,” Environmetrics, vol. 18, no. 6, pp. 607–620, 2007. [4] S. K. Thompson and G. A. F. Seber, Adaptive Sampling, Wiley Series in Probability and Statistics: Probability and Statistics, John Wiley & Sons, New York, NY, USA, 1996. [5] D. Raj, β€œSome estimators in sampling with varying probabilities without replacement,” Journal of the American Statistical Association, vol. 51, pp. 269–284, 1956. [6] D. R. Smith, M. J. Conroy, and D. H. Brakhage, β€œEfficiency of adaptive cluster sampling for estimating density of wintering waterfowl,” Biometrics, vol. 51, no. 2, pp. 777–788, 1995. [7] N. Chutiman and B. Kumphon, β€œRatio estimator using two auxiliary variables for adaptive cluster sampling,” Journal of the Thai Statistical Association, vol. 6, no. 2, pp. 241–256, 2008.

Advances in

Operations Research Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Advances in

Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Applied Mathematics

Algebra

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Probability and Statistics Volume 2014

The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Differential Equations Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Submit your manuscripts at http://www.hindawi.com International Journal of

Advances in

Combinatorics Hindawi Publishing Corporation http://www.hindawi.com

Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of Mathematics and Mathematical Sciences

Mathematical Problems in Engineering

Journal of

Mathematics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Discrete Mathematics

Journal of

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Discrete Dynamics in Nature and Society

Journal of

Function Spaces Hindawi Publishing Corporation http://www.hindawi.com

Abstract and Applied Analysis

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Journal of

Stochastic Analysis

Optimization

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014