Modified Systematic Sampling in the Presence of Linear Trend
|
|
- Wesley Mitchell
- 6 years ago
- Views:
Transcription
1 Modified Systeatic Sapling in the Presence of Linear Trend Zaheen Khan, and Javid Shabbir Keywords: Abstract A new systeatic sapling design called Modified Systeatic Sapling (MSS, proposed by ] is ore general than Linear Systeatic Sapling (LSS and Circular Systeatic Sapling (CSS In the present paper, this schee is further extended for populations having a linear trend Expressions for ean and variance of saple ean are obtained for the population having perfect linear trend aong population values Expression for the average variance is also obtained for super population odel Further, efficiency of MSS with respect to CSS is obtained for different saple size LSS, CSS, MSS, linear trend, super population 000 AMS Classification: 6DO5 Introduction In survey sapling, Linear Systeatic Sapling (LSS is a coonly used design Generally, it is useful when population size N is a ultiple of saple size n, ie N = nk ( k is the sapling interval Thus, we have k saples each of size n However, LSS is not beneficial when population size N is not a ultiple of the saple size n, ie N nk Because in this case, LSS cannot provide a constant saple size n, thus, estiate of population ean (total is biased Therefore, Circular Systeatic Sapling (CSS was introduced by Lahiri in 95 (cited in, p39] Contrary to LSS, CSS is not advantageous when population size N is a ultiple of the saple size n, ie N = nk as in this case, CSS produces n replicates of k saples Further, in CSS, the nuber of saples also rapidly increase to N as copared to k saples of LSS To iprove the efficiency of systeatic sapling, researchers proposed several odifications in the selection procedure The considerable work is done by 4], 6] and 7] In the recent years, 8] proposed Diagonal Systeatic Sapling (DSS under the condition that n k as a copetitor of LSS Later, the condition n k for DSS have relaxed by 9] The generalization of DSS is suggested by 0] Soe odification in LSS are proposed by ], in which odd and even saple sizes are dealt separately Further odification on LSS is also proposed by ] Diagonal Circular Systeatic Sapling (DCSS proposed by 5] is an extension of DSS to Departent of Matheatics and Statistics, Federal Urdu University of Arts, Science and Technology, Islaabad, Pakistan Eail: zkurdu@gailco, Corresponding Auther Departent of Statistics, Quaid-i-Aza University, Islaabad 4530, Pakistan, Eail: javidshabbir@gailco
2 the circular version of systeatic sapling A note on DCSS has been proposed by 3] However, soe of these schees are applicable when N = nk while other can be used only when N nk A new systeatic sapling design called Modified Systeatic Sapling (MSS proposed by ], which is applicable in both situations, whether N = nk or N nk According to this design, first copute least coon ultiple of N and n, ie L, then find k,, s and k, k = L n, = L N, s = N k and k = k /] or k = N/n] is rounded off to integer Consequently, s = n, which eans that there are sets and each set contains s units in a saple Thus, in MSS the j th unit of the i th set of a saple of n units can be written as: ( y (r ij = r+(i k+(j k hn if hn < r+(i k+(j k (h+n for h = 0,, ; i =,, 3,, and j =,, 3,, s This sapling schee reduces to LSS if L = N or N = nk and CSS if L = N n, the detail is given below If N = nk then L = N, k = k, = and s = n Thus, Equation ( reduces to ( y (r j = r + (j k, j =,, 3,, n which is LSS Siilarly, if L = N n, then k = N, = n and s = So, Equation ( can be written as (3 y (r i = r + (i k hn if hn < r + (i k (h + N for h = 0,, ; i =,, 3,, n Which is CSS To study the characteristics of MSS, we use an alternative ethod by partitioning the total nuber of saples into different sets of siilar saples To develop an alternative ethod, let us assue that k can be written as k = qk + r, q and r are quotient and reainder respectively Further, we assue that w = if ( q and, w = ( q if ( q > In both cases, there are two types of partitioning, ie between saples and within saples(see detail in Subsections and When w = In this case partitioning between saples and within saples are given in the Subsections and Partitioning between saples In this case, k possible saples are ainly partitioned into two groups The first group consists of initial k ( k saples and second group contains last ( k saples However, in the second group, there are ( subgroups, each attains k saples If a rando nuber r is selected fro the first k units of a population, there is a possibility that it is selected fro the first group, ie k ( k or it is selected fro the ( subgroups of the second group, ie k ( uk < r k ( u k such that u =,,, (, integer u is selected corresponding to a rando nuber r
3 Saple article to present hjs class 3 Partitioning within saples Furtherore, in the first group, all s units of all sets in each saple are labeled as r + (i k + (j k, such that i =,,, and j =,,, s; while in the second group, all s units of the first ( u sets are labeled as r + (i k + (j k such that i =,,, ( u and j =,,, s and in each of the last u sets, first (s units are labeled as r + (i k + (j k such that i = ( u +,, ; j =,,, (s and last unit is labeled as r+(i k+(j k N such that i = ( u+,, and j = s When w = (n q > In this case partitioning between saples and within saples are given in the Subsections and Partitioning between saples In this case k saples are ainly partitioned into two groups, the first group consists of the nuber of saples in which r k (w k + r The econd group contains the nuber of saples in which r > k (w k + r The first group is further partitioned into ( w (w + subgroups in which, there are r nuber of saples in each of the first and the last subgroups, and k saples in each of the iddle ( w (w subgroups In each subgroup of the first group, corresponding to a rando nuber r, an integer u is picked in such a way that u = (w if r k ( u k, u = w, w +, w +,, ( w if k ( uk < r k ( u k and u = ( w+ if k ( uk < r k ( uk + r The second group consists of last (w k r = (k r + (w k saples, which is the cobination of the first (k r saples and the last (w sets of k saples These (w sets of k saples further partitioned in such a way that the first r of every k saples fors the first subgroup and the last (k r saples of every k saples together with the first (k r saples of this group fors the second subgroup However, when w =, then we have only (k r saples in the second group Partitioning within saples In each saple of the first group, all s units of the first ( u sets are labeled as r+(i k+(j k such that i =,,, ( u and j =,,, s, and in each of the last u sets, the first (s units are labeled as r + (i k + (j k such that i = ( u +,,, j =,,, (s and the last unit is labeled as r + (i k + (j k N such that i = ( u +,, and j = s In each saple of the first subgroup of the second group, all s units of the first (w x sets are labeled as r + (i k + (j k such that i =,,, (w x and j =,,, s; the units of iddle ( w + sets are labeled in such a way that, the first (s units of each set are labeled as r + (i k + (j k such that i (w x +,, ( x +, j =,,, (s and the last unit of each set is labeled as r + (i k + (j k N such that i (w x +,, ( x + and j = s; the units of the last (x sets are labeled in such a way that, the first (s units are labeled as r+(i k+(j k such that i ( x+,,, j =,,, (s and the last two units in each set is labeled as r+(i k+(j k N such that i ( x +,, and j = (s, s However, when s =, the units
4 4 in these (x sets are labeled as r + (i k + (j k N The possible values of x are, 3,, (w Note: If w =, then this set of saples does not exist In the second subgroup of the second group, all s units of the the first (w x sets are labeled as r + (i k + (j k such that i =,,, (w x and j =,,, s; The units of iddle ( w sets are labeled in such a way that the first (s units of each set are labeled as r + (i k + (j k such that i (w x+,, ( x and j =,,, (s, the last unit of each set is labeled as r + (i k + (j k N such that i (w x+,, ( x and j = s, the units of the last (x sets are labeled in such a way that, the first (s units are labeled as r +(i k +(j k such that i ( x+,,, j =,,, (s and the last two units in each set is labeled as r + (i k + (j k N such that i ( x +,, and j = (s, s However, when s =, the units in these (x sets are labeled as r + (i k + (j k N, the possible values of x are,, (w Mean and variance of MSS for population having linear trend The following linear odel of hypothetical population is to be considered as ( Y t = α + β t, t =,, 3,, N α and β are the intercept and slope of the odel respectively Mean of MSS The saple ean for both cases, ie w = and w > are given below (see detial in Appendix A Case (i when w = ( ȳ MSS = α + β Case (ii when w > r + (s k + ( k ], if r k ( k ] r + (s k + ( k u k, u =,,, ( if k ( uk < r k ( u k (3 ȳ MSS = α + β ] r + (s k + ( k u N, n u = (w if r k ( u k u = w, w +,, ( w if N ( uk < r k ( u k, u = ( w + if k ( k < r k ( uk ] + r r + (s k + ( k ( w + x N, n x =,, (w if k (w xk < r k (w xk ] + r r + (s k + ( k ( w + x N, n x =,, 3,, (w if k (w xk + r < r k (w xk + k
5 Saple article to present hjs class 5 If w =, the Equation (3 will reduce to (4 ȳ MSS = α + β ] r + (s k + ( k u N n u = (w if r k ( u k, u = w, w +,, ( w if N ( uk < r k ( u k, u = ( w + if k ( k < r k ( uk ] + r r + (s k + ( k ( w + x N, n x =,, 3,, (w if k (w xk + r < r k (w xk + k Unbiasedness of ȳ MSS (see detial in Appendix A The saple ean (ȳ MSS is an unbiased estiator of population ean (Ȳ N + (5 E(ȳ MSS = α + β = Ȳ Variance of ȳ MSS (see detial in Appendix A (i when w = (6 V (ȳ MSS = b (k + ( k(k k ] Note: In this case, if N = nk then L = N, so =, thus, V (ȳ MSS = b (k This is a variance of linear systeatic sapling (ii when w > (7 (ȳ MSS = b (k + ( k(k k +4w(w k 3k (3 w + k] 3 Yates corrected estiator Yates corrected estiator of population ean for MSS is derived below 3 Yates corrected estiator for MSS The corrected estiator ȳmss c of population ean using MSS is given by (3 ȳmss c = n λ l Y r + ] n l= Y rl + λ l Y rn, λ l and λ l are selected so that saple ean coincides with the population ean in the presence of linear trend for all choices of r,,, k Alternatively Equation (3 can be written as (3 ȳ c MSS = ȳ MSS + a l (r (Y r Y rn, a l (r = λ l n = λ l n
6 6 Under the odel given in (, the population ean is (33 Ȳ = α + β N + As entioned earlier, that there are two cases, ie (i w = and (ii w > First, we consider the Case (i 3 Case (i: when w = In this situation, a rando start r is selected fro k units such that r k ( k or r > k ( k If r k ( k, then l = and the last value of each saple is labeled r + ( k + (s k Thus, (3 becoes (34 ȳ c MSS = ȳ MSS + a (r ( Y r Y r+( k+(s k Under the linear odel given in (, we have ȳ MSS = α+β r + (s k + ( k ], Y r = α + β r and Y r+( k+(s k = α + β r + ( k + (s k Putting these values in (34, we have (35 ȳmss c = α + β r + ] ( k + (s k a(r ( k + (s k Coparing the coefficients of α and β in (33 and (35 and solving for a (r, we have r + ( k k a (r = ( k + (s k Putting a (r in (34, we have (36 ȳmss c = ȳ MSS + r + ( k k ( Yr Y r+( k+(s k ( k + (s k If r > k ( k, then l = and the last value of each saple is labeled r + ( k + (s k N Thus, (3 becoes (37 ȳ c MSS = ȳ MSS + a (r ( Y r Y r+( k+(s k N Under the linear odel (, we have ȳ MSS = α+β r + (s k ] + ( k u k, u =,,, ( ifk ( uk < r k ( u k, Y r = α+β r and Y r+(s k+( k N = α + β r + (s k + ( k N Putting these values in (37, we have (38 ȳmss c = α + β r + ] k (s k + ( k u a(r( k k Coparing the coefficients of α and β in (33 and (38 and solving for a (r, we have r (k + + ( k uk / (39 a (r =, ( k k u =,,, (, which are picked corresponding to a rando nuber r such that k ( uk < r k ( u k
7 Saple article to present hjs class 7 Putting a (r in (37, we have r ȳmss c (k + + ( k uk / = ȳ MSS + (30 ( k k ( Yr Y r+( k+(s k N 3 Case (ii: when w > As entioned earlier in Section, when s =, MSS becoes CSS (see ] Therefore, we focus the MSS for s > It is also entioned in Subsection, all k saples are partitioned into two groups The first group contains the saples r k (w k + r and the second group consist of the saples in which r > k (w k + r The corrected saple ean for each saple in the first group is siilar to the corrected saple ean found in Subsection 3, r > k ( k, because the pattern of saples in both situations is siilar Further, the weights assigned to the first and the last units of each saple in this group will be siilar to the weights given in (39, ie ȳ c MSS = ȳ MSS + a (r ( Y r Y r+( k+(s k N, r (k + + ( k uk /n a (r =, ( k k u = (w corresponding to a rando nuber r such that r k ( u k, u = w, w+,, ( w if k ( uk < r k ( u k and u = ( w + if k ( uk < r k ( uk + r Thus, (3 y c MSS = ȳ MSS + r (k++( k uk /n ( k k ( Y r Y r+( k+(s k N The second group having saples in which r > k (w k + r, and the first subgroup consists of the nuber of saples in which k (w xk < r k (w xk + r such that x =,, (w The Yates corrected estiator with l = 3 in (3, for the saples of the first subgroup can be written as (3 ȳ c MSS = ȳ MSS + a 3 (r ( Y r Y r+( k+(s k N Under a linear odel ( ȳ MSS = α + βr + (s k + ( k ( w + x k ], x =,, (w if k (w xk < r k (w xk + r, Y r = α + β r and Y r+(s k+( k N = α + βr + (s k + ( k N Putting these values in (3, we have (33 ȳmss r c = α + β + (s k + ( k ( w + x k ] +a 3 (r( k k Coparing the coefficients of α and β given in (33 and (33 and solving for a 3 (r, we have r (k + + ( k k ( w + x / (34 a 3 (r = ( k k
8 8 Putting a 3 (r in the corrected estiator given in (3, we have ȳmss c = ȳ (35 MSS + r (k ++( k k ( w+x / ( k k ( Y r Y r+( k+(s k N, x =,, (w, which are picked corresponding to a rando nuber r such that k (w xk < r k (w xk + r Siilarly, the second subgroup consists of the nuber of saples in which k (w xk + r < r k (w xk + k such that x =,,, (w The Yates corrected estiator with l = 4 in (3, for saples of this subgroup, can be written as (36 ȳ c MSS = ȳ MSS + a 4 (r ( Y r Y r+( k+(s k N Under the linear odel (, ȳ MSS = α+βr + (s k +( k ( w + x k ], x =,,, (w if k (w xk < r k (w xk + r, Y r = α + βr and Y r+(s k+( k N = α + β r + (s k + ( k N Putting these values in (36, we have (37 ȳmss r c = α + β + (s k + ( k ( w + x k ] +a 4 (r( k k Coparing the coefficients of α and β given in (33 and (37 and solving for a 4 (r, we have r (k + + ( k k ( w + x/ (38 a 4 (r = ( k k Putting a 4 (r in the corrected estiator given in (36, we have y (39 MSS c = ȳ MSS + r (k++( k k ( w+x/ ( k k ( Y r Y r+( k+(s k N, x =,,, (w, which are picked corresponding to a rando nuber r such that k (w xk < r k (w xk + r 4 Average variance In real life application, we hardly found such population exhibiting perfect linear trend Therefore, it is necessary to study the average variance of the corrected estiator under MSS using following super population odel (4 Y t = α + β t + e t, E(e t = 0, V (e t = E(e t = σ t g, Cov(e t, e v = 0, t v =,, 3,, N and g is the predeterined constant The average variance of ȳ (r MSS under odified systeatic sapling for population odeled by Y t = α + β t + e t is given by
9 Saple article to present hjs class 9 Case (i when w = (see detial in Appendix B (4 E V (ȳ (r MSS = σ r= χ (u, r k ( k k k ( u k + u= / N +k t= tg r=k χ ( uk+ (u, r N, χ (u, r = δ + (rrg + θ s i= j= r + (i k + (j k g +δ (rr + ( k + (s k g, χ (u, r = δ + u (rrg + θ s n j= (r + (i k + (j k g i= ( s + ( i= u+ j= (r + (i k + (j k g + r + (i k g +(s k N + δ (r (r + ( k + (s k N g, δ + l (r = a l(ra l (r + ( n N, δ l (r = a l(ra l (r ( n ( N and θ = n n N, such that l =, Case (ii When w > (see detial in Appendix B (43 ( E V ȳ (r MSS = σ + w u=w + w+ w k ( u k k u=w k ( u k r= χ (u, r r=k ( uk+ χ (u, r k ( u k u= w+ k (w xk+r + w x= r=k ( uk+ χ (u, r r=k (w xk+ χ 3 (x, r k (w xk+k + w x= r=k χ (w xk+r + 4 (x, r / ] N + k t= tg N, χ (u, r = δ + u (rrg + θ s n j= (r + (i k + (j k g + i= u+ ( s i= j= (r + (i k + (j k g +(r + (i k + (s k N g + δ (r (r + ( k + (s k N g, χ 3 (x, r = δ 3 + w x s (rrg + θ j= r + (i k + (j k g i= ( s + x+ i=w x+ j= r + (i k + (j k g +r + (i k + (s k N g ( s + i= x+ j= r + ( k + (j k g + s j=s r + (i k + (j k N g +δ3 (rr + ( k + (s k N g,
10 0 χ 4 (x, r = δ + 4 (rrg + θ w x i= ( s s j= r + (i k + (j k g + x i=w x+ j= r + (i k + (j k g +r + (i k + (s k N g ( s + i= x+ j= r + ( k + (j k g + s j=s r + (i k + (j k N g +δ4 (rr + ( k + (s k N g, δ + l (r = a l(ra l (r + ( n N, δ l (r = a l(ra l (r ( n ( N and θ = n n N, for l =, 3, 4 5 Epirical study Due to the coplex nature of the derived expressions, the average variances of MSS and CSS cannot be theoretically copared Therefore, in this paper, a coputer based efficiency coparison of MSS and CSS is ade nuerically under super population odel (4 The nuerical coparison has been ade for N =, N = 30, N = 50 and N = 78 As entioned earliar, if L = N then MSS reduces to LSS and if L = (N n then MSS becoes CSS Therefore, the choice of a saple size considered in this paper is based on the fact that N < L < (N n The relative efficiency of MSS over CSS is presented in Table under g = 0,,, 3 This table includes 40 different cobinations of N and n each at g = 0,, and 3 which are to be considered for efficiency coparison, and it is observed that CSS is not applicable for 4 cobinations Thus, we have 36 4 = 44 results of efficiency coparison and found that MSS is ore efficient than CSS in 35 cases Further, it is to be noted, whenever N n = ( n +, the efficiency of MSS over CSS is draatically increased 5 Natural Population We use the following natural population for efficiency coparison The results are given in Table Population : Source:, page8] Table reflacts that MSS is ore efficient than CSS 6 Conclusion Modified Systeatic Sapling (MSS is a ore general schee than LSS and CSS Because, when least coon ultiple of N and n is equal to lower extree, ie L = N, MSS coincides with LSS If it is equal to upper extree, ie L = (N n, then MSS coincides with CSS However, when L lies between these two extree values, ie N < L < (N n, MSS is advantageous over CSS In this case, the nuber of saples is considerably reduced in MSS as copared to CSS, ie iniu reduction is half of the saples Contrary to the CSS, the explicit expressions for ean and variance of ean are derived for population having perfect linear trend aong the population values Further, nuerical coparison is carried out in this paper clearly favors the use of MSS over CSS for population odeled by a super population odel with linear trend as well as for natural population
11 Saple article to present hjs class Table Rercent Relative Efficiency (PRE of MSS over CSS under linear trend N n g = 0 g = g = g = 3 N n g = 0 g = g = g = The sybol (- indicates that CSS is not possible Table Percent Relative Efficiency (PRE of MSS over CSS for Population N = 80 Variance Eff = V (CSS V (MSS 00 n MSS CSS Here, V (MSS = Variance of odified systeatic sapling and V (CSS = Variance of circular systeatic sapling Acknowledgeent: We would like to thank the two anonyous Referees whose coents helped to iprove the presentation of this paper
12 References ] Murthy, MN Sapling Theory and Method, Statistical Publishing Society, Calcutta, 967 ] Khan, Z Shabbir, J and Gupta, S A new sapling design for systeatic sapling, Counications in Statistics - Theory and Methods,4, , 03 3] Khan, Z Gupta, S and Shabbir, J A Note on Diagonal Circular Systeatic Sapling, Journal of Statistical Theory and Practice, 8, , 04 4] Madow, WG On the theory of systeatic sapling III, Annals of Matheatics and Statistics, 4, 0-06, 953 5] Sapath, S and Varalakshi, V Diagonal circular systeatic sapling, Model Assisted Statistics and Applications, 3 (4, , 008 6] Sethi, V K On optiu pairing of units, Sankhya B, 7, 35-30, 965 7] Singh, D Jindal, KK and Grag, JNOn odified systeatic sapling, Bioetrika, 55, 4-546, 968 8] Subraani, J Diagonal systeatic sapling schee for finite populations, Journal of the Indian Society of Agriculture Statistics, 53 (, 87-95, 000 9] Subraani, J Further results on diagonal systeatic sapling schee for finite populations, Journal of the Indian Society of Agricultural Statistics, 63 (3, 77-8, 009 0] Subraani, JGeneralization of diagonal systeatic sapling schee for finite populations, Model Assisted Statistics and Applications, 5, 7-8, 00 ] Subraani, J A Modification on linear systeatic sapling, Model Assisted Statistics and Applications, 8 (3, 5-7, 03a ] Subraani, J A further odification on linear systeatic sapling for finite populations, Journal of Statistical Theory and Practice, 7, , 03b
13 Saple article to present hjs class 3 Appendix A Mean and variance of MSS for population having linear trend The following linear odel of hypothetical population is to be considered (A Y t = α + β t, t =,, 3,, N, α and β are the intercept and slope of the odel respectively A Mean of MSS The saple ean for both cases, ie w = and w > are separately discussed below: Case (i when w = If r (k ( k, the ean, ȳ MSS can be written as ȳ MSS = α + β s w+ i= s r + (i k + (j k j= After siplification, we have ȳ MSS = α + β r + ] (s k + ( k If k ( uk < r k ( u k for u =,,,, then u ȳ MSS = α + β s s i= j= r + (i k + (j k + s i= u+ j= r + (i k + (j k ] + r + (i k + (s k N After siplification, we have ȳ MSS = α + β r + (s k + ( k u k ] Thus ȳ MSS is a piecewise function of r, ie r + (s k + ( k ] if r k ( k ] (A ȳ MSS = α + β r + (s k + ( k u k u =,,, ( if k ( uk < r k ( u k Case (ii when w > If r k (w k + r, then r ust belongs to any one of the three subgroups which have been discussed in Section Therefore, corresponding to a rando nuber r, an integer u is picked in such a way that u = (w if r k ( u k; u = w, w +, w +,, ( w if k ( uk < r k ( u k and u = ( w + if k ( uk < r k ( uk + r
14 4 (A3 For each subgroup, ȳ MSS can be written as u ȳ MSS = α + β s s i= j= r + (i k + (j k + s i= u+ j= r + (i k + (j k ] + r + (i k + (s k N After few steps, we have ȳ MSS = α + β r + (s k + ( k u k ] If w >, then it is also possible that k (w xk < r k (w xk+r, such that x =, 3,, w So, w x ȳ MSS = α + β s s i= j= r + (i k + (j k + x+ s i=w x+ j= r + (i k + (j k + r + (i k + (s k N + i= x+ s j= r + (i k + (j k + s j=s r + (i k + (j k N When s =, then Equation (A3 can be expressed as w x ȳ MSS = α + β s s i= j= r + (i k + (j k + x+ s i=w x+ j= r + (i k + (j k + r + (i k + (s k N + s ] i= x+ j=s r + (i k + (j k N Also, when s =, then Equation (A3 can be expressed as w x ȳ MSS = α + β s s i= j= r + (i k + (j k + x+ i=w x+ r + (i k + (s k N + i= x+ r + (i k + (s k N ] After siplifying of Equation (A3 for each case, ie s =, s = and s >, we have ȳ MSS = α + β r + (s k + ( k ( w + x k ] If k (w xk + r < r k (w xk + k, then w x ȳ MSS = α + β s s i= j= r + (i k + (j k + x s i=w x+ j= r + (i k + (j k (A4 + r + (i k + (s k N + i= x+ s j= r + (i k + (j k ] + s j=s r + (i k + (j k N
15 Saple article to present hjs class 5 When s =, then Equation (A4 can be expressed as w x ȳ MSS = α + β s s i= j= r + (i k + (j k + x s i=w x+ j= r + (i k + (j k + r + (i k + (s k N + i= x+ s j=s r + (i k + (j k N When s =, then Equation (A4 can be expressed as, w x ȳ MSS = α + β s s i= j= r + (i k + (j k + x i=w x+ r + (i k + (s k N + ] i= x+ r + (i k + (s k N After siplification of Equation (A4 for each case, ie s =, s = and s >, we have ȳ MSS = α + β r + (s k + ( k ( w + x k ] Thus, ean of MSS for above odel of hypothetical population with rando start r is given by: r + (s k + ( k u ] N n u = (w if r k ( u k, u = w, w +,, ( w if N ( uk < r k ( u k, u = ( w + if (A5 ȳ MSS = α + β k ( uk < r k ( uk + r r + (s k + ( k ( w + x ] N n, x =,, (w if k (w xk < r k (w xk + r r + (s k + ( k ( w + x ] N n, x =,, 3,, (w if k (w xk + r < r k (w xk + k If w =, then Equation (A5 reduces to r + (s k ] + ( k u N n u = (w if r k ( u k, u = w, w +,, ( w if N ( uk < r k ( u k, ȳ MSS = α+β u = ( w + if k ( uk < r k ( uk + r r + (s k ] + ( k ( w + x N n, x =,, 3,, (w if k (w xk + r < r k (w xk + k A Unbiasedness of saple ean ȳ MSS We have two cases: ]
16 6 Case (i when w = : Taking the expected value of (A, we have (k ( k E(ȳ MSS = k r= r + (s k + ( k ] + u= k ( u k r=(k ( uk+ r + (s k + ( k u k ] ] As sk = N, then (k ( k E(ȳ MSS = k r= α + β r + (N k + ( k ] + k ( u k u= r=(k ( uk+ α + β r + (N k + ( k ]] u k After a little algebra, we have N + E(ȳ MSS = α + β = Ȳ, which shows that ȳ MSS is an unbiased estiator of Ȳ Case (ii when w > : If w >, we take the expected value of (A5, we have w E(ȳ MSS = k ( u k k u=w r= α + β r + ((s k ] +( k u k + w k ( u k u=w r=k ( uk+ α + β ] +( k u k + w+ k ( uk+r u= w+ r=k ( uk+ α + β ] +( k u k + w k (w xk+r x= r=k (w xk+ α + β ] +( k ( w +xk + w k (w xk+k x= r=k (w xk+r + α + β ]] +( k ( w+xk r + ((s k r + ((s k r + ((s k r + ((s k But if w =, we take the expected value of (A5, given by w E(ȳ MSS = k ( u k k u=w r= α + β r + ((s k ] +( k u k + w u=w k ( u k r=k ( uk+ ] +( k u k + w+ u= w+ ] +( k u k + w x= α + β k ( uk+r r=k ( uk+ k (w xk+k r=k (w xk+r + +( k ( w+xk ]] r + ((s k α + β α + β r + ((s k r + ((s k
17 Saple article to present hjs class 7 After few steps, we have E(ȳ MSS = α + β k r k + wk + wr ] +kw k w kw + N + E(ȳ MSS = α + β (N + k (w ] +r (w wk(w + k(w ] N + (w (A6 E(ȳ MSS = α + β + (r k wk + k When w = ( q in (A6, we have ] N + ( q E(ȳ MSS = α + β + (r k ( qk + k, ] N + ( q E(ȳ MSS = α + β + (qk + r k k + k, N + ( q E(ȳ MSS = α + β + (qk + r k, N + E(ȳ MSS = α + β = Ȳ the above equation shows that ȳ MSS is unbiased estiator of Ȳ as k = qk+r Note: Putting w = in (A6, we also have N + E(ȳ MSS = α + β = Ȳ A The variance of ȳ MSS V (ȳ MSS = E(ȳ MSS Ȳ = k k r= (ȳ r(mss Ȳ (i when w = k ( k V (ȳ MSS = k ȳr(mss Ȳ + u= r= k ( u k r=k ( uk+ ȳr(mss Ȳ ] k ( k V (ȳ MSS = k r= α + β r + ((s k + ( k ] α + β N+ + k ( u k u= r=k ( uk+ α + βr + ((s k + ( k u k ] ] α + β N+ k ( k V (ȳ MSS = k r= β r + (( k k ] β + k ( u k u= r=k ( uk+ β r + (( k k ] ] u k β
18 8 After siplification, we have (A7 V (ȳ MSS = b (k + ( k(k k ], Note: If N = nk, then L = N, so =, thus V (ȳ MSS = b (k This is a variance of linear systeatic sapling (ii when w > If w >, then V (ȳ MSS will be expressed as: V (ȳ MSS = k w k ( u k u=w r= k ( u k r=k ( uk+ k ( uk+r u= w+ r=k ( uk+ k (w xk+r r=(k (w xk+ + w u=w + w+ + w x= + w x= ȳr(mss Ȳ ȳr(mss Ȳ ȳr(mss Ȳ ȳr(mss Ȳ k (w xk+k r=(k (w xk+r + ȳr(mss Ȳ ] If w =, then the ter w k (w xk+r x= r=(k (w xk+ ȳr(mss Ȳ will be oitted fro V (ȳ MSS w V (ȳ MSS = ( k ( u k k u=w r= α + β r + (s k +( k u k ] α + β N+ + w ( k ( u k u=w r=k ( uk+ α + β r + (s k +( k u k ] α + β N+ + w+ ( k ( uk+r u= w+ r=k ( uk+ α + β r + (s k +( k u k ] α + β N+ + w k (w xk+r x= r=(k (w xk+ α + β r + (s k +( k ( w + x k ] α + β N+ + w k (w xk+k x= r=(k (w xk+r + α + β r + (s k +( k ( w + x k ] ] α + β N+,
19 (A8 w V (ȳ MSS = k u k β + w u=w β Saple article to present hjs class 9 u=w ] u k + w+ u k k ( u k r= β r + (( k k k ( u k r=k ( uk+ u= w+ ] β k (w xk+r r=k (w xk+ + w x= ( w + x k + w x= β r + (( k k k ( uk+r r=k ( uk+ β r + (( k k β α + β ] k (w xk+k r=k (w xk+r + ( w + x k β ] r + (( k k α + β r + (( k k If w =, then Equation (A8 reduces to w V (ȳ MSS = k u k β + w u=w β u=w ] u k + w+ u k k ( u k r= β r + (( k k k ( u k r=k ] ( uk+ β u= w+ ] β k (w xk+k r=k (w xk+r + + w x= ( w + x k After siplification, we have r + (( k k k ( uk+r r=k ( uk+ β r + (( k k β ] ] α + β r + (( k k V (ȳ MSS = b (k + ( k(k k +4(w 3k( q w (k qk + (k k ] +k w 3k (3 w + k The ter (w 3k( q w (k qk + (k k] will be vanished in both situations, when w = or w = ( q So, we are left with (A9 V (ȳ MSS = b (k + ( k(k k +4w(w k 3k (3 w + k]
20 0 Appendix B Average variance In real life application, we hardly found such population exhibiting perfect linear trend, therefore, it is necessary to study the average variance of the corrected estiator under MSS using following super population odel (B Y t = α + βt + e t, E(e t = 0, V (e t = E(e t = σ t g, Cov(e t, e v = 0, t v =,, 3,, Nand g is a predeterine constant Under the above super population odel (B, the average variance expression of MSS is given below: Case (i when w = Consider that l th su of squares (SSl are given by (B SSl = (B3 ȳ (r MSS Ȳ ] = ȳmss + a l (r (Y r Y rn Ȳ ], l = if r k ( k and l = if r > k ( k When r k ( k, the expressions of ȳ MSS, Ȳ, Y rand Y rn under the odel Y t = α + β t + e t, can be expressed as: ȳ MSS = α + β r + (s k + ( k ] + s i= s j= e r+(i k+(j k, Ȳ = α + β N+ + N N t= e t, Y r = α + βr + e r α + β r + ( k + (s k + e r+( k+(s k and Y rn = Substituting these expressions in (B, we have ] SS = ȳ (r MSS Ȳ = s n i= j= e r+(i k+(j k +na (r ( e r e r+( k+(s k N N t= t] e Siilarly, if r > k ( k the expressions of ȳ MSS, Ȳ, Y r and Y rn under the super population odel Y t = α + β t + e t, can be written as ȳ MSS = α + β r + (s k ] + ( k u k u + s s i= j= e r+(i k+(j k + s ] i= u+ j= e r+(i k+(j k + e r+(i k+(s k N, Ȳ = α + β N+ + N N t= e t, Y r = α + βr + e r α + β r + ( k + (s k N + e r+( k+(s k N and Y rn = Thus, SS = ȳ (r MSS Ȳ ] = n + i= u+ u i= ( s s j= e r+(i k+(j k j= e r+(i k+(j k + e r+(i k+(s k N + na (r(e r e r+( k+(s k N N N t= t] e The average variance of the corrected saple ean can be written as: E V (ȳ (r MSS k ( k = k + u= r= E (SS k ( u k r=k E (SS ( uk+
21 Saple article to present hjs class under the assuptions of super population odel s E (SS = n i= j= E(e r+(i k+(j k + n a (r E(e r +E(e r+( k+(s k + na (r E(e r E(e r+( k+(s k + N N t= E(e t s nn i= j= E(e r+(i k+(j k ] +na (r E(e r E(e r+( k+(s k, E (SS = σ n i= s j= (r + (i k + (j k g +n a (r r g + (r + ( k + (s k g +na (r r g (r + ( k + (s k g s nn i= j= (r + (i k + (j k g +na (r r g (r + ( k + (s k g + N N t= tg ] (B4 E (SS = σ a (r ( a (r + ( n N r g ( + n n s N i= j= r + (i k + (j k g +a (r ( a (r ( n N r + ( k + (s k g + N N t= tg Siilarly, E (SS = u n i=( s + i= u+ + i= u+ +E(e r+(i k+(s k N s j= E(e r+(i k+(j k j= E(e r+(i k+(j k +E(e r+(i k+(s k N + n a (r E(e r +E(e r+( k+(s k N + na (r (E(e r E(e r+( k+(s k N u s nn i= j= E(e r+(i k+(j k ( s j= E(e r+(i k+(j k + na (r (E(e r E(e r+( k+(s k N + ] N N t= E(e t,
22 (B5 E (SS = u n i=( s s j= (r + (i k + (j k g + i= u+ j= (r + (i k + (j k g + (r + (i k + (s k N g + n a (r r g + (r + ( k + (s k N g +na (r r g (r + ( k + (s k N g u s j= (r + (i k + (j k g nn i= ( s + i= u+ j= (r + (i k + (j k g + (r + (i k + (s k N g + na (r r g (r + ( k + (s k N g + ] N N t= tg, E (SS = a (r ( a (r + ( n N r g ( n u s N n j= (r + (i k + (j k g + n i= ( s + i= u+ j= (r + (i k + (j k g + (r + (i k + (s k N g +a (r ( a (r ( n N ( r + ( k + (s k N g + N N t= tg ] Equations (B4 and (B5 can be written as: and E (SS = σ δ + (rrg + θ s i= j= r + (i k + (j k g +δ (rr + ( k + (s k g + N N t= tg E (SS = σ δ + u (rrg + θ s n i= j= (r + (i k + (j k g + ( s i= u+ j= (r + (i k + (j k g + (r + (i k + (s k N g +δ (r (r + ( k + (s k N g + N N t= ], tg δ + l (r = a l(ra l (r + ( n N, δ l (r = a l(ra l (r ( n ( N and θ = n n N, such that l =, Also (B6 E (SS = σ χ (u, r + N N t= tg and (B7 E (SS = σ χ (u, r + N N t= tg,
23 Saple article to present hjs class 3 χ (u, r = δ + (rrg + θ s i= j= r + (i k + (j k g +δ (rr + ( k + (s k g, and χ (u, r = δ + u (rrg + θ s n j= (r + (i k + (j k g i= ( s + i= u+ + (r + (i k + (s k N g j= (r + (i k + (j k g +δ (r (r + ( k + (s k N g Substituting the values of E (SS and E (SS in (B3, we have E V (ȳ (r MSS k ( k = σ k r= χ (u, r (B8 + k ( u k u= r=k χ ( uk+ (u, r / N +k t= tg N Case (ii when w > We can write ] (B9 SSl = ȳ (r MSS Ȳ = ȳmss + a l (r (Y r Y rn Ȳ ], l = if r k (w k + r, l = 3 if k (w xk < r k (w xk + r such that x =,, ( and l = 4 if k (w xk+r < r k (w xk+k such that x =,,, Furtherore, when r k (w k + r, we realize whether r k ( u k such that u = w, k ( uk < r k ( u k such that u = w, w +,, ( w or k ( uk < r k ( uk + r such that u = ( w + However, for each of these subgroups E (SS will be used Thus, the average variance of the corrected saple ean can be expressed as (B0 E V ( ȳ (r MSS ] w = k ( u k N u=w k ( u k + n w u=w + w+ r= E SS] r=k ( uk+ E SS] k ( uk+r u= w+ k (w xk+r + w x= + w x= r=k ( uk+ E SS] r=k E SS3] (w xk+ k (w xk+k r=k E SS4] (w xk+r + ] The E (SS is already obtained in case of w =, ie (B E (SS = χ (u, r + N N t= tg Now consider (B E (SS3 = E ȳ MSS + a 3 (r (Y r Y rn Y ]
24 4 Under the super population odel, we have ȳ MSS = α + β r + (s k + ( k ( w + x k w x + s n i= j= e r+(i k+(j k + x+ s i=w x+ j= e r+(i k+(j k + e r+(i k+(s k N + i= x+ s j= e r+(i k+(j k + s j=s e r+(i k+(s k N ], Ȳ = α + β N+ + N N t= e t, Y r = α + βr + e r β r + ( k + (s k + e r+( k+(s k and Y rn = α + Substituting these expressions in (B, we have E (SS3 = E n w x + x+ i=w x+ i= s s s j= e r+(i k+(j k j= e r+(i k+(j k + e r+(i k+(s k N + i= x+ j= e r+(i k+(j k + s + na (r(e r e r+( k+(s k N N N t= t] e j=s e r+(i k+(j k N Applying the assuption of super population odel, we have E (SS3 = w x n i= s s j= E(e r+(i k+(j k + x+ i=w x+ +E(e r+(i k+(s k N s j= E(e r+(i k+(j k + i= x+ j= E(e r+(i k+(j k + s j=s E(e r+(i k+(j k N +n a (r E(e r + E(e r+( k+(s k N +na (r E(e r E(e r+( k+(s k N nn w x i= + x+ i=w x+ +E(e r+(i k+(s k N s j= E(e r+(i k+(j k s j= E(e r+(i k+(j k s + i= x+ j= E(e r+(i k+(j k + s j=s E(e r+(i k+(j k N +na (r E(e r E(e r+( k+(s k N + N N t= E(e t ]
25 Saple article to present hjs class 5 E (SS3 = w x n i= s s j= (r + (i k + (j k g + x+ i=w x+ j= (r + (i k + (j k g + (r + (i k + (j k N g s + i= x+ j= (r + (i k + (j k g + s j=s (r + (i k + (j k N g +n a (r r g + (r + (i k + (j k N g +na (r r g (r + (i k + (j k N g w x s j= (r + (i k + (j k g nn i= s + x+ i=w x+ j= (r + (i k + (j k g + (r + (i k + (j k N g s + i= x+ j= (r + (i k + (j k g + s j=s (r + (i k + (j k N g +na (r r g (r + (i k + (j k N g + N N t= tg ], E (SS3 = a 3 (r ( a 3 (r + ( n N k (w xk+r r=k (w xk+ + w x= +(i k + (j k g ( s r g n ( n w x s N i= j= r + x+ i=w x+ j= r + (i k + (j k g +r + (i k + (s k N g ( s + i= x+ j= r + ( k + (j k g + s j=s r + (i k + (j k N g +a 3 (r ( a 3 (r ( n N r + ( k + (s k N g + N N t= tg, E (SS3 = δ + 3 (rrg + θ w x i= ( s s j= r + (i k + (j k g + x+ i=w x+ j= r + (i k + (j k g +r + (i k + (s k N g ( s + i= x+ j= r + ( k + (j k g + s j=s r + (i k + (j k N g +δ3 (rr + ( k + (s k N g + N N t= tg, δ 3 + (r = a 3(ra 3 (r + ( n N and δ 3 (r = a 3 (ra 3 (r ( n N
26 6 Also (B3 E (SS3 = χ 3 (x, r + N N t= tg, χ 3 (x, r = δ 3 + w x s (rrg + θ j= r + (i k + (j k g Siilarly, i= ( s + x+ i=w x+ j= r + (i k + (j k g +r + (i k + (s k N g ( s + i= x+ j= r + ( k + (j k g + s j=s r + (i k + (j k N g +δ3 (rr + ( k + (s k N g (B4 E (SS4 = χ 4 (x, r + N N t= tg, χ 4 (x, r = δ 4 + w x s (rrg + θ j= r + (i k + (j k g i= ( s + x i=w x+ j= r + (i k + (j k g +r + (i k + (s k N g ( s + i= x+ j= r + ( k + (j k g + s j=s r + (i k + (j k N g +δ4 (rr + ( k + (s k N g Putting E (SSl for l =, 3, 4 in (B0, we have ( E V ȳ (r w MSS = σ k ( u k k u=w r= χ (u, r + w k ( u k u=w r=k χ ( uk+ (u, r + w+ k ( u k u= w+ r=k χ ( uk+ (u, r + w k (w xk+r x= r=k χ (w xk+ 3 (x, r + w k (w xk+k x= r=k χ (w xk+r + 4 (x, r / N +k t= tg N ]
MODIFIED SYSTEMATIC SAMPLING WITH MULTIPLE RANDOM STARTS
RESTAT Statistical Journal olume 6, Number, April 08, 87 MODIFIED SYSTEMATIC SAMPLING WITH MULTIPLE RANDOM STARTS Authors: Sat Gupta Department of Mathematics and Statistics, University of North Carolina,
More informationAN EFFICIENT CLASS OF CHAIN ESTIMATORS OF POPULATION VARIANCE UNDER SUB-SAMPLING SCHEME
J. Japan Statist. Soc. Vol. 35 No. 005 73 86 AN EFFICIENT CLASS OF CHAIN ESTIMATORS OF POPULATION VARIANCE UNDER SUB-SAMPLING SCHEME H. S. Jhajj*, M. K. Shara* and Lovleen Kuar Grover** For estiating the
More informationCurious Bounds for Floor Function Sums
1 47 6 11 Journal of Integer Sequences, Vol. 1 (018), Article 18.1.8 Curious Bounds for Floor Function Sus Thotsaporn Thanatipanonda and Elaine Wong 1 Science Division Mahidol University International
More informationEstimation of the Mean of the Exponential Distribution Using Maximum Ranked Set Sampling with Unequal Samples
Open Journal of Statistics, 4, 4, 64-649 Published Online Septeber 4 in SciRes http//wwwscirporg/ournal/os http//ddoiorg/436/os4486 Estiation of the Mean of the Eponential Distribution Using Maiu Ranked
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationEstimation of the Population Mean Based on Extremes Ranked Set Sampling
Aerican Journal of Matheatics Statistics 05, 5(: 3-3 DOI: 0.593/j.ajs.05050.05 Estiation of the Population Mean Based on Extrees Ranked Set Sapling B. S. Biradar,*, Santosha C. D. Departent of Studies
More informationTesting equality of variances for multiple univariate normal populations
University of Wollongong Research Online Centre for Statistical & Survey Methodology Working Paper Series Faculty of Engineering and Inforation Sciences 0 esting equality of variances for ultiple univariate
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationPolygonal Designs: Existence and Construction
Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationIn this chapter, we consider several graph-theoretic and probabilistic models
THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions
More informationEstimating Parameters for a Gaussian pdf
Pattern Recognition and achine Learning Jaes L. Crowley ENSIAG 3 IS First Seester 00/0 Lesson 5 7 Noveber 00 Contents Estiating Paraeters for a Gaussian pdf Notation... The Pattern Recognition Proble...3
More informationMoments of the product and ratio of two correlated chi-square variables
Stat Papers 009 50:581 59 DOI 10.1007/s0036-007-0105-0 REGULAR ARTICLE Moents of the product and ratio of two correlated chi-square variables Anwar H. Joarder Received: June 006 / Revised: 8 October 007
More informationThe Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters
journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn
More informationTEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES
TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES S. E. Ahed, R. J. Tokins and A. I. Volodin Departent of Matheatics and Statistics University of Regina Regina,
More informationBIVARIATE NONCENTRAL DISTRIBUTIONS: AN APPROACH VIA THE COMPOUNDING METHOD
South African Statist J 06 50, 03 03 BIVARIATE NONCENTRAL DISTRIBUTIONS: AN APPROACH VIA THE COMPOUNDING METHOD Johan Ferreira Departent of Statistics, Faculty of Natural and Agricultural Sciences, University
More informationA Simple Regression Problem
A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where
More information1 Proof of learning bounds
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #4 Scribe: Akshay Mittal February 13, 2013 1 Proof of learning bounds For intuition of the following theore, suppose there exists a
More informationBiostatistics Department Technical Report
Biostatistics Departent Technical Report BST006-00 Estiation of Prevalence by Pool Screening With Equal Sized Pools and a egative Binoial Sapling Model Charles R. Katholi, Ph.D. Eeritus Professor Departent
More informationBipartite subgraphs and the smallest eigenvalue
Bipartite subgraphs and the sallest eigenvalue Noga Alon Benny Sudaov Abstract Two results dealing with the relation between the sallest eigenvalue of a graph and its bipartite subgraphs are obtained.
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More informationAn Introduction to Meta-Analysis
An Introduction to Meta-Analysis Douglas G. Bonett University of California, Santa Cruz How to cite this work: Bonett, D.G. (2016) An Introduction to Meta-analysis. Retrieved fro http://people.ucsc.edu/~dgbonett/eta.htl
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More informationSampling How Big a Sample?
C. G. G. Aitken, 1 Ph.D. Sapling How Big a Saple? REFERENCE: Aitken CGG. Sapling how big a saple? J Forensic Sci 1999;44(4):750 760. ABSTRACT: It is thought that, in a consignent of discrete units, a certain
More informationNOTES AND CORRESPONDENCE. Two Extra Components in the Brier Score Decomposition
752 W E A T H E R A N D F O R E C A S T I N G VOLUME 23 NOTES AND CORRESPONDENCE Two Extra Coponents in the Brier Score Decoposition D. B. STEPHENSON School of Engineering, Coputing, and Matheatics, University
More informationare equal to zero, where, q = p 1. For each gene j, the pairwise null and alternative hypotheses are,
Page of 8 Suppleentary Materials: A ultiple testing procedure for ulti-diensional pairwise coparisons with application to gene expression studies Anjana Grandhi, Wenge Guo, Shyaal D. Peddada S Notations
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
More informationTwo New Unbiased Point Estimates Of A Population Variance
Journal of Modern Applied Statistical Methods Volue 5 Issue Article 7 5--006 Two New Unbiased Point Estiates Of A Population Variance Matthew E. Ela The University of Alabaa, ela@baa.ua.edu Follow this
More informationIN modern society that various systems have become more
Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto
More informationA new type of lower bound for the largest eigenvalue of a symmetric matrix
Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box
More informationGeneral Properties of Radiation Detectors Supplements
Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent
More informationUpper bound on false alarm rate for landmine detection and classification using syntactic pattern recognition
Upper bound on false alar rate for landine detection and classification using syntactic pattern recognition Ahed O. Nasif, Brian L. Mark, Kenneth J. Hintz, and Nathalia Peixoto Dept. of Electrical and
More informationBootstrapping Dependent Data
Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly
More informationNote on generating all subsets of a finite set with disjoint unions
Note on generating all subsets of a finite set with disjoint unions David Ellis e-ail: dce27@ca.ac.uk Subitted: Dec 2, 2008; Accepted: May 12, 2009; Published: May 20, 2009 Matheatics Subject Classification:
More informationResearch Article Rapidly-Converging Series Representations of a Mutual-Information Integral
International Scholarly Research Network ISRN Counications and Networking Volue 11, Article ID 5465, 6 pages doi:1.54/11/5465 Research Article Rapidly-Converging Series Representations of a Mutual-Inforation
More informationC na (1) a=l. c = CO + Clm + CZ TWO-STAGE SAMPLE DESIGN WITH SMALL CLUSTERS. 1. Introduction
TWO-STGE SMPLE DESIGN WITH SMLL CLUSTERS Robert G. Clark and David G. Steel School of Matheatics and pplied Statistics, University of Wollongong, NSW 5 ustralia. (robert.clark@abs.gov.au) Key Words: saple
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationpaper prepared for the 1996 PTRC Conference, September 2-6, Brunel University, UK ON THE CALIBRATION OF THE GRAVITY MODEL
paper prepared for the 1996 PTRC Conference, Septeber 2-6, Brunel University, UK ON THE CALIBRATION OF THE GRAVITY MODEL Nanne J. van der Zijpp 1 Transportation and Traffic Engineering Section Delft University
More informationMODIFICATION OF AN ANALYTICAL MODEL FOR CONTAINER LOADING PROBLEMS
MODIFICATIO OF A AALYTICAL MODEL FOR COTAIER LOADIG PROBLEMS Reception date: DEC.99 otification to authors: 04 MAR. 2001 Cevriye GECER Departent of Industrial Engineering, University of Gazi 06570 Maltepe,
More informationResearch Article Some Formulae of Products of the Apostol-Bernoulli and Apostol-Euler Polynomials
Discrete Dynaics in Nature and Society Volue 202, Article ID 927953, pages doi:055/202/927953 Research Article Soe Forulae of Products of the Apostol-Bernoulli and Apostol-Euler Polynoials Yuan He and
More informationarxiv: v1 [cs.ds] 3 Feb 2014
arxiv:40.043v [cs.ds] 3 Feb 04 A Bound on the Expected Optiality of Rando Feasible Solutions to Cobinatorial Optiization Probles Evan A. Sultani The Johns Hopins University APL evan@sultani.co http://www.sultani.co/
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationBest Linear Unbiased and Invariant Reconstructors for the Past Records
BULLETIN of the MALAYSIAN MATHEMATICAL SCIENCES SOCIETY http:/athusy/bulletin Bull Malays Math Sci Soc (2) 37(4) (2014), 1017 1028 Best Linear Unbiased and Invariant Reconstructors for the Past Records
More informationa a a a a a a m a b a b
Algebra / Trig Final Exa Study Guide (Fall Seester) Moncada/Dunphy Inforation About the Final Exa The final exa is cuulative, covering Appendix A (A.1-A.5) and Chapter 1. All probles will be ultiple choice
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationBoosting with log-loss
Boosting with log-loss Marco Cusuano-Towner Septeber 2, 202 The proble Suppose we have data exaples {x i, y i ) i =... } for a two-class proble with y i {, }. Let F x) be the predictor function with the
More informationCombining Classifiers
Cobining Classifiers Generic ethods of generating and cobining ultiple classifiers Bagging Boosting References: Duda, Hart & Stork, pg 475-480. Hastie, Tibsharini, Friedan, pg 246-256 and Chapter 10. http://www.boosting.org/
More informationPOD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION EQUATION IN NEURO-MUSCULAR SYSTEM
6th European Conference on Coputational Mechanics (ECCM 6) 7th European Conference on Coputational Fluid Dynaics (ECFD 7) 1115 June 2018, Glasgow, UK POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationDEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS
ISSN 1440-771X AUSTRALIA DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS An Iproved Method for Bandwidth Selection When Estiating ROC Curves Peter G Hall and Rob J Hyndan Working Paper 11/00 An iproved
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE
DRAFT Proceedings of the ASME 014 International Mechanical Engineering Congress & Exposition IMECE014 Noveber 14-0, 014, Montreal, Quebec, Canada IMECE014-36371 ANALYTICAL INVESTIGATION AND PARAMETRIC
More informationFast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials
Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter
More informationWhen Short Runs Beat Long Runs
When Short Runs Beat Long Runs Sean Luke George Mason University http://www.cs.gu.edu/ sean/ Abstract What will yield the best results: doing one run n generations long or doing runs n/ generations long
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More informationEXPLICIT CONGRUENCES FOR EULER POLYNOMIALS
EXPLICIT CONGRUENCES FOR EULER POLYNOMIALS Zhi-Wei Sun Departent of Matheatics, Nanjing University Nanjing 10093, People s Republic of China zwsun@nju.edu.cn Abstract In this paper we establish soe explicit
More informationAN OPTIMAL SHRINKAGE FACTOR IN PREDICTION OF ORDERED RANDOM EFFECTS
Statistica Sinica 6 016, 1709-178 doi:http://dx.doi.org/10.5705/ss.0014.0034 AN OPTIMAL SHRINKAGE FACTOR IN PREDICTION OF ORDERED RANDOM EFFECTS Nilabja Guha 1, Anindya Roy, Yaakov Malinovsky and Gauri
More informationLectures 8 & 9: The Z-transform.
Lectures 8 & 9: The Z-transfor. 1. Definitions. The Z-transfor is defined as a function series (a series in which each ter is a function of one or ore variables: Z[] where is a C valued function f : N
More informationMulticollision Attacks on Some Generalized Sequential Hash Functions
Multicollision Attacks on Soe Generalized Sequential Hash Functions M. Nandi David R. Cheriton School of Coputer Science University of Waterloo Waterloo, Ontario N2L 3G1, Canada 2nandi@uwaterloo.ca D.
More informationSoft Computing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis
Soft Coputing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis Beverly Rivera 1,2, Irbis Gallegos 1, and Vladik Kreinovich 2 1 Regional Cyber and Energy Security Center RCES
More informationUniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval
Unifor Approxiation and Bernstein Polynoials with Coefficients in the Unit Interval Weiang Qian and Marc D. Riedel Electrical and Coputer Engineering, University of Minnesota 200 Union St. S.E. Minneapolis,
More informationBest Procedures For Sample-Free Item Analysis
Best Procedures For Saple-Free Ite Analysis Benjain D. Wright University of Chicago Graha A. Douglas University of Western Australia Wright s (1969) widely used "unconditional" procedure for Rasch saple-free
More information1 Identical Parallel Machines
FB3: Matheatik/Inforatik Dr. Syaantak Das Winter 2017/18 Optiizing under Uncertainty Lecture Notes 3: Scheduling to Miniize Makespan In any standard scheduling proble, we are given a set of jobs J = {j
More informationImpact of Imperfect Channel State Information on ARQ Schemes over Rayleigh Fading Channels
This full text paper was peer reviewed at the direction of IEEE Counications Society subject atter experts for publication in the IEEE ICC 9 proceedings Ipact of Iperfect Channel State Inforation on ARQ
More informationFairness via priority scheduling
Fairness via priority scheduling Veeraruna Kavitha, N Heachandra and Debayan Das IEOR, IIT Bobay, Mubai, 400076, India vavitha,nh,debayan}@iitbacin Abstract In the context of ulti-agent resource allocation
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More informationA MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION
A eshsize boosting algorith in kernel density estiation A MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION C.C. Ishiekwene, S.M. Ogbonwan and J.E. Osewenkhae Departent of Matheatics, University
More informationON SEQUENCES OF NUMBERS IN GENERALIZED ARITHMETIC AND GEOMETRIC PROGRESSIONS
Palestine Journal of Matheatics Vol 4) 05), 70 76 Palestine Polytechnic University-PPU 05 ON SEQUENCES OF NUMBERS IN GENERALIZED ARITHMETIC AND GEOMETRIC PROGRESSIONS Julius Fergy T Rabago Counicated by
More informationMA304 Differential Geometry
MA304 Differential Geoetry Hoework 4 solutions Spring 018 6% of the final ark 1. The paraeterised curve αt = t cosh t for t R is called the catenary. Find the curvature of αt. Solution. Fro hoework question
More informationA Comparison of Three Different Procedures for Estimating Variance Components
International Journal of Statistics and Applications 2016, 6(4): 241-252 DOI: 10.5923/j.statistics.20160604.05 A Coparison of Three Different Procedures for Estiating Variance Coponents Zakaria A. Abdel
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More information3.8 Three Types of Convergence
3.8 Three Types of Convergence 3.8 Three Types of Convergence 93 Suppose that we are given a sequence functions {f k } k N on a set X and another function f on X. What does it ean for f k to converge to
More informationEnzyme kinetics: A note on negative reaction constants in Lineweaver-Burk plots
Enzye kinetics: A note on negative reaction constants in Lineweaver-Burk plots Sharistha Dhatt# and Kaal Bhattacharyya* Departent of Cheistry, University of Calcutta, Kolkata 700 009, India #pcsdhatt@gail.co
More informationCompression and Predictive Distributions for Large Alphabet i.i.d and Markov models
2014 IEEE International Syposiu on Inforation Theory Copression and Predictive Distributions for Large Alphabet i.i.d and Markov odels Xiao Yang Departent of Statistics Yale University New Haven, CT, 06511
More informationLost-Sales Problems with Stochastic Lead Times: Convexity Results for Base-Stock Policies
OPERATIONS RESEARCH Vol. 52, No. 5, Septeber October 2004, pp. 795 803 issn 0030-364X eissn 1526-5463 04 5205 0795 infors doi 10.1287/opre.1040.0130 2004 INFORMS TECHNICAL NOTE Lost-Sales Probles with
More informationOn Constant Power Water-filling
On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives
More informationFast Structural Similarity Search of Noncoding RNAs Based on Matched Filtering of Stem Patterns
Fast Structural Siilarity Search of Noncoding RNs Based on Matched Filtering of Ste Patterns Byung-Jun Yoon Dept. of Electrical Engineering alifornia Institute of Technology Pasadena, 91125, S Eail: bjyoon@caltech.edu
More informationLeast squares fitting with elliptic paraboloids
MATHEMATICAL COMMUNICATIONS 409 Math. Coun. 18(013), 409 415 Least squares fitting with elliptic paraboloids Heluth Späth 1, 1 Departent of Matheatics, University of Oldenburg, Postfach 503, D-6 111 Oldenburg,
More informationTesting the lag length of vector autoregressive models: A power comparison between portmanteau and Lagrange multiplier tests
Working Papers 2017-03 Testing the lag length of vector autoregressive odels: A power coparison between portanteau and Lagrange ultiplier tests Raja Ben Hajria National Engineering School, University of
More informationEstimation for the Parameters of the Exponentiated Exponential Distribution Using a Median Ranked Set Sampling
Journal of Modern Applied Statistical Metods Volue 14 Issue 1 Article 19 5-1-015 Estiation for te Paraeters of te Exponentiated Exponential Distribution Using a Median Ranked Set Sapling Monjed H. Sau
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationComparing Probabilistic Forecasting Systems with the Brier Score
1076 W E A T H E R A N D F O R E C A S T I N G VOLUME 22 Coparing Probabilistic Forecasting Systes with the Brier Score CHRISTOPHER A. T. FERRO School of Engineering, Coputing and Matheatics, University
More informationAlgorithms for parallel processor scheduling with distinct due windows and unit-time jobs
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 57, No. 3, 2009 Algoriths for parallel processor scheduling with distinct due windows and unit-tie obs A. JANIAK 1, W.A. JANIAK 2, and
More informationInference on Downton s Bivariate Exponential Distribution Based on Moving Extreme Ranked Set Sampling
AUSTRIAN JOURNAL OF STATISTICS Volue 42 203), Nuber 3, 6 79 Inference on Downton s Bivariate Exponential Distribution Based on Moving Extree Ranked Set Sapling Ahed Ali Hanandeh and Mohaad Fraiwan Al-Saleh
More informationWorld's largest Science, Technology & Medicine Open Access book publisher
PUBLISHED BY World's largest Science, Technology & Medicine Open Access book publisher 2750+ OPEN ACCESS BOOKS 95,000+ INTERNATIONAL AUTHORS AND EDITORS 88+ MILLION DOWNLOADS BOOKS DELIVERED TO 5 COUNTRIES
More informationDESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *
IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS
More informationDepartment of Electronic and Optical Engineering, Ordnance Engineering College, Shijiazhuang, , China
6th International Conference on Machinery, Materials, Environent, Biotechnology and Coputer (MMEBC 06) Solving Multi-Sensor Multi-Target Assignent Proble Based on Copositive Cobat Efficiency and QPSO Algorith
More informationE. Alpaydın AERFAISS
E. Alpaydın AERFAISS 00 Introduction Questions: Is the error rate of y classifier less than %? Is k-nn ore accurate than MLP? Does having PCA before iprove accuracy? Which kernel leads to highest accuracy
More informationCHAPTER ONE. Physics and the Life Sciences
Solution anual for Physics for the Life Sciences 2nd Edition by Allang Link download full: http://testbankair.co/download/solution-anual-forphysics-for-the-life-sciences-2nd-edition-by-allang/ CHAPTER
More informationCompleteness of Bethe s states for generalized XXZ model, II.
Copleteness of Bethe s states for generalized XXZ odel, II Anatol N Kirillov and Nadejda A Liskova Steklov Matheatical Institute, Fontanka 27, StPetersburg, 90, Russia Abstract For any rational nuber 2
More informationSymbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm
Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat
More informationA general forulation of the cross-nested logit odel Michel Bierlaire, Dpt of Matheatics, EPFL, Lausanne Phone: Fax:
A general forulation of the cross-nested logit odel Michel Bierlaire, EPFL Conference paper STRC 2001 Session: Choices A general forulation of the cross-nested logit odel Michel Bierlaire, Dpt of Matheatics,
More informationRandomized Recovery for Boolean Compressed Sensing
Randoized Recovery for Boolean Copressed Sensing Mitra Fatei and Martin Vetterli Laboratory of Audiovisual Counication École Polytechnique Fédéral de Lausanne (EPFL) Eail: {itra.fatei, artin.vetterli}@epfl.ch
More informationClassical and Bayesian Inference for an Extension of the Exponential Distribution under Progressive Type-II Censored Data with Binomial Removals
J. Stat. Appl. Pro. Lett. 1, No. 3, 75-86 (2014) 75 Journal of Statistics Applications & Probability Letters An International Journal http://dx.doi.org/10.12785/jsapl/010304 Classical and Bayesian Inference
More informationMeasures of average are called measures of central tendency and include the mean, median, mode, and midrange.
CHAPTER 3 Data Description Objectives Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance,
More informationAn improved self-adaptive harmony search algorithm for joint replenishment problems
An iproved self-adaptive harony search algorith for joint replenishent probles Lin Wang School of Manageent, Huazhong University of Science & Technology zhoulearner@gail.co Xiaojian Zhou School of Manageent,
More informationarxiv: v3 [cs.ds] 22 Mar 2016
A Shifting Bloo Filter Fraewor for Set Queries arxiv:1510.03019v3 [cs.ds] Mar 01 ABSTRACT Tong Yang Peing University, China yangtongeail@gail.co Yuanun Zhong Nanjing University, China un@sail.nju.edu.cn
More information