K. Krishnamoorthy a & Jie Peng b a Department of Mathematics, University of Louisiana at
|
|
- Ruth Stephens
- 6 years ago
- Views:
Transcription
1 This article was downloaded by: [ ] On: 4 March 014, At: 05:08 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: Registered office: Mortimer House, Mortimer Street, London W1T 3JH, UK Journal of Statistical Computation and Simulation Publication details, including instructions for authors and subscription information: Approximate one-sided tolerance limits in random effects model and in some mixed models and comparisons K. Krishnamoorthy a & Jie Peng b a Department of Mathematics, University of Louisiana at Lafayette, Lafayette, LA 70504, USA b Department of Finance, Economics and Decision Science, St. Ambrose University, Davenport, IA 5803, USA Published online: 17 Mar 014. To cite this article: K. Krishnamoorthy & Jie Peng 014): Approximate one-sided tolerance limits in random effects model and in some mixed models and comparisons, Journal of Statistical Computation and Simulation, DOI: / To link to this article: PLEASE SCROLL DOWN FOR ARTICLE Taylor & Francis makes every effort to ensure the accuracy of all the information the Content ) contained in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and should be independently verified with primary sources of information. Taylor and Francis shall not be liable for any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of the Content. This article may be used for research, teaching, and private study purposes. Any substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to anyone is expressly forbidden. Terms &
2 Conditions of access and use can be found at
3 Journal of Statistical Computation and Simulation, Approximate one-sided tolerance limits in random effects model and in some mixed models and comparisons K. Krishnamoorthy a and Jie Peng b a Department of Mathematics, University of Louisiana at Lafayette, Lafayette, LA 70504, USA; b Department of Finance, Economics and Decision Science, St. Ambrose University, Davenport, IA 5803, USA Received 10 July 013; accepted 1 January 014) An approximate closed-form one-sided tolerance limit TL) in a general mixed model is proposed. One-sided TLs for the distribution of observable random variable and for the distribution of unobservable random variable in one-way random model are obtained as special cases from the one for the general mixed model. Applications to a two-way nested random model are also given. The merits of the TLs are evaluated using Monte Carlo simulation and compared with the existing ones. Our comparison studies indicate that the approximate TLs are quite satisfactory for all parameter and sample size configurations, and better than the existing ones in some cases. Approximate confidence intervals for exceedance probabilities in one-way random effects model are also proposed. The procedures are illustrated using three examples. Keywords: exceedance probability; generalized variable; Graybill Wang approximation; MOVER; Satterthwaite approximation; nested model; survival probability; variance components 1. Introduction Tolerance intervals TIs) are widely used in many areas of sciences including engineering, quality assessment, pharmaceutical industry. Apart from these areas of sciences, TIs are also used to assess the exposure level in a workplace or pollution level in an environment.[1,] The problems of constructing one-sided tolerance limits TLs) and two-sided TIs in one-way random models and in some general mixed models have been investigated by several authors for the past three decades. For example, see Mee and Owen,[3] Beckman and Tietjen,[4] Bhaumik and Kulkarni,[5] Vangel,[6] Liao and Iyer,[7] Krishnamoorthy and Mathew,[8] and Liao et al.,[9] and the references therein. So far all available methods are approximate, and no single approach is satisfactory for all parameter sample size configurations. Krishnamoorthy and Mathew [8] have proposed one-sided TLs based on the generalized variable approach, and an approximate one based on the noncentral t percentiles for one-way random model with balanced or unbalanced data. More specifically, they proposed approximate TLs for the distribution of observable random variable and for the distribution of unobservable true value. Monte Carlo evaluation and comparison studies by these authors indicated that none of the methods, including Vangel s [6] methods, is satisfactory for all situations. The approximate Corresponding author. krishna@louisiana.edu 014 Taylor & Francis
4 K. Krishnamoorthy and J. Peng TLs based on the noncentral t percentiles are simple to compute and are satisfactory when the intraclass correlation ρ.5, which maybe a realistic condition in many applications. In some cases, the generalized TLs [8,9] and Vangel s [6] TLs could be very conservative, yielding TIs that are unnecessarily wide. Recently, Hoffman [10] has proposed an approximate closed-form TLs for one-way random models which are quite conservative; a bootstrap-based adjustment was suggested to improve the TLs, and the improved ones are comparable with the generalized TLs by Krishnamoorthy and Mathew.[8] An excellent review of available procedures for constructing one-sided as well as two-sided TIs in one-way random models and in general mixed models, and applications of TIs, see Chapters 4 6 of the book by Krishnamoorthy and Mathew.[11] In this article, we shall provide simple closed-form approximate one-sided TLs in a general mixed model using the method of variance estimate recovery MOVER) by Zou and Donner [1] and Zou et al.[13,14] In the context of Bland-Altman limits of agreement, Zou [15] has used the MOVER to develop confidence intervals CIs) for quantiles in a one-way random model. Noting that a one-sided TL is a one-sided confidence limit CL) for an appropriate quantile, Zou s result can be used to find approximate one-sided TLs for the distribution of observable random variable in a one-way random model. We here generalize Zou s result by obtaining one-sided TLs in a general mixed model, and show applications in finding TLs for the unobservable random variable in one-way random models, and for finding TLs for observable as well as for unobservable random variables in a two-way nested model with random effects. Furthermore, we shall compare the approximate TLs with the popular ones proposed in the literature, and the generalized TLs proposed by Fonseca et al. [16] for a two-way nested model. A related problem of estimating the exceedance probability survival probability) in a one-way random model is also addressed and an approximate solution is provided. The rest of the article is organized as follows. In the following section, we describe the MOVER approach for finding CLs for a linear combination of parameters, and for finding TLs in a general mixed model. The TLs in one-way random model with balanced or unbalanced data are described, and their coverage properties are studied in Section 3. Also, a closed-form solutions to the related problem of finding CLs for the survival probability or exceedance probability) are given. Approximate TLs in a two-way nested model with random effects are given in Section 4, and they are compared with the generalized TLs given in Fonseca et al.[16] The results are illustrated using three examples in Section 5, and some concluding remarks are given in Section 6.. MOVER TLs in general mixed model The MOVER CI by Zou and Donner,[1] Zou et al. [13,14] for a linear combination of parameters θ 1,..., θ g is described as follows. Let ˆθ i be an unbiased estimate of θ i, i = 1,..., k. Assume that ˆθ 1,..., ˆθ g are independent. Further, let l i, u i ) denote the 1 α CI for θ i, i = 1,..., k. The 1 α MOVER CI L, U) for k c iθ i can be expressed as L = g c i ˆθ i g ci ˆθ i li ) with li = { l i if c i > 0, u i if c i < 0 1) and U = g c i ˆθ i + g ci ˆθ i ui ) with ui = { u i if c i > 0, l i if c i < 0. )
5 Journal of Statistical Computation and Simulation 3 It should be noted that Graybill and Wang [17] first obtained the above CI for a linear combinations variance components. Zou and co-authors gave a different argument so as to the above CI is valid for any parameters. Finding a p,1 α) upper TL in one-way random model or in a general mixed model simplifies i ) distribution, where a i s are to finding a 1 α upper CL for the p quantile of a Nμ, g a iσ known constants, and σ i s are linear combinations of variance components. Noting that the p g quantile of the Nμ, g a iσi ) distribution is given by μ + z p a iσi, an upper TL can be g readily obtained using 1 α upper CLs for μ and a iσi. Let U μ and U v denote 1 α upper CLs for μ and g a iσi, respectively. The p,1 α) upper TL for the Nμ, g a iσi ) distribution can be obtained from Equation )as g ˆμ + z p a i σ i ) 1/ + + U μ ˆμ) + zp U v g a i ˆσ i ) 1/ + 1/, 3) where x + = max{x,0}, and ˆσ i s are independent unbiased estimates of σi s. The p,1 α) MOVER lower TL is similarly obtained from Equation 1) as g ˆμ z p a i ˆσ i ) 1/ + + L μ ˆμ) + zp U v g a i ˆσ i ) 1/ + 1/, 4) where L μ isa1 α lower CL for μ. The required upper CL U v for g a i ˆσ i )1/ can be obtained using the MOVER method. We also note that U μ ˆμ) = L μ ˆμ) for all the cases considered in the sequel, as a result, the expressions within curly parentheses of Equations 3) and 4) will be the same. One-sided TLs for various problems involving random models with balanced or unbalanced data can be easily constructed in a straightforward manner using the above general form. In fact, for each problem, we simply need to identify the coefficients a i and the individual CLs for μ, and σ i, i = 1,..., g. 3. One-way random models Let Y ij denote the jth observation corresponding to the ith level, assumed to follow the one-way random model Y ij = μ + τ i + e ij, j = 1,,..., n i, i = 1,,..., k, 5) where μ is an unknown general mean, τ i s represent random effects, and e ij s represent error terms. It is assumed that τ i s and e ij s are all mutually independent with τ i N0, στ ) and e ij N0, σe ). Thus, the observable random variable Y ij Nμ, στ + σ e ), and the unobservable true value associated with the ith level μ + τ i Nμ, στ ). We shall consider the problems of constructing one-sided TLs for the distribution Nμ, στ + σ e ), and for the distribution Nμ, σ τ ).
6 4 K. Krishnamoorthy and J. Peng Define ñ = 1 k SS e = k n 1 i, Ȳ i. = 1 n i n i j=1 Y ij, Ȳ = 1 k k Ȳ i., SSȳ = k Ȳ i. Ȳ), and k n i Y ij Ȳ i. ). 6) j=1 Then, Ȳ i. N μ, στ + σ e ), n i Ȳ N μ, σ τ +ñσ e ), and SS e σe k χ N k with N = k n i. It can be verified that ESSȳ) = k 1)στ +ñσ e ). Furthermore, SSȳ σ τ +ñσ e χ k 1. 7) When n 1 = =n k, the above distributional result is exact, and Ȳ, SSȳ and SS e are mutually independent. For the unbalanced case, the distributional result in Equation 7) is approximate,[18] Ȳ and SS e are independently distributed, and SSȳ and SS e are independently distributed. However, Ȳ and SSȳ are not independent. As noted in Krishnamoorthy and Mathew,[8] and as will be seen in the sequel, this dependence does not pose any serious problems to obtain a satisfactory TLs for the unbalanced case One-sided TLs for Nμ, σ τ + σ e ) MOVER TLs To construct a TL for the Nμ, στ + σ e ) distribution, let σ 1 = σ τ +ñσ e and σ = σ e, so that σ τ = σ 1 ñσ. 8) In these notations, we see that Y ij Nμ, σ τ + σ e ) = Nμ, a 1σ 1 + a σ ) with a 1, a ) = 1, 1 ñ), and ˆμ = Ȳ Nμ, σ1 /k). To obtain a 1 α upper CL for a 1σ1 + a σ, note that ˆσ 1 = SS ȳ k 1 σ 1 χ k 1 k 1 and ˆσ = SS e N k σ χ N k N k. Furthermore, ˆσ 1 and ˆσ are independent with E ˆσ 1 ) = σ 1 and E ˆσ ) = σ. Using these distributional results, we obtain a 1 α upper CL for μ as U μ = Ȳ + t k 1;1 α ˆσ 1 k. 9) The above CL is exact when the data are balanced, and is an approximate when the data are unbalanced. The upper CL for a 1 σ1 + a σ can be obtained as ) ) 1/ U v = a 1 ˆσ 1 + a ˆσ ) + a1 ˆσ 1 4 k 1 χk 1;α 1 + a ˆσ 4 N k χn k;α 1. 10)
7 Journal of Statistical Computation and Simulation 5 Using Equations 9) and 10) in Equation 3), a p,1 α) upper TL for Nμ, στ + σ e ) can be obtained as { Ȳ + z p a 1 ˆσ 1 + a ˆσ )1/ + tk 1;1 α ˆσ 1 k + z p a 1 ˆσ 1 + a ˆσ ) 1/ } U v ), 11) where U v is given in Equation 10), and a 1 = 1 and a = 1 ñ. Furthermore, it follows from Equation 4) that the p,1 α) MOVER lower TL for the Nμ, στ + σ e ) is given by { Ȳ z p a 1 ˆσ 1 + a ˆσ )1/ + tk 1;1 α ˆσ 1 k + z p a 1 ˆσ 1 + a ˆσ ) 1/ } U v ). 1) TLs based on noncentral t percentiles The approximate upper TL given in Krishnamoorthy and Mathew [8] is expressed as SSȳ Ȳ + t k 1;1 α δ) kk 1), with δ = z p k + kk 1)1 ñ) N k ) SS 1/ e F k 1,N k;α, 13) SSȳ where F m,n;q denotes the qth quantile of an F distribution with dfs m and n. We shall refer to the above TL as the NCT TL Vangel s TLs Vangel [6] proposed approximate TLs for the balanced case, and they can be described as follows. Define Ȳ i. = 1 n n Y ij, j=1 Ȳ = 1 kn k n Y ij, SS τ =n j=1 k Ȳ i. Ȳ.. ) and SS e = k n Y ij Ȳ i. ). Further, define MS τ = SS τ /k 1), MS e = SS e /kn 1), ˆσ τ = 1/n)MS τ MS e ) and ˆσ e = MS e. Vangel s lower TL is expressed as Ȳ F ˆσ τ +ˆσ e, and the corresponding upper TL is given by Ȳ + F ˆσ τ +ˆσ e, where the F is the tolerance factor given by j=1 F V = a 1 + a W + a 3 W + a 4 W 3, 14) and W =[1 + n 1)/MS τ /MS e )] 1/. The coefficients a 1, a, a 3 and a 4 are to be determined numerically, and they are tabulated in Vangel s paper for p,1 α) = 0.90, 0.95) and 0.99, 0.95), and some selected values of k, n). 3.. One-sided TLs for Nμ, σ τ ) MOVER TLs To construct a p,1 α)upper TL for Nμ, στ ), we note from Equation 8) that σ τ = a 1σ1 + a σ with a 1 = 1 and a = ñ. To find an upper CL for μ + z p σ τ = μ + z p a 1 σ1 + a σ, we can
8 6 K. Krishnamoorthy and J. Peng use the upper CL for μ as defined in Equation 9), and we need to find only the upper CL for a 1 σ1 + a σ with a 1 = 1 and a = ñ. The 1 α upper CL for a 1 σ1 + a σ is expressed as ) ) 1/ U v = a 1 ˆσ 1 + a ˆσ ) + + a1 ˆσ 1 4 a a ˆσ 4 N a 1. 15) χ a 1;α χ N a;1 α In terms of U μ in Equation 9) and U v in Equation 15), the p,1 α) upper TL for Nμ, στ ) is given by { Ȳ + z p a1 ˆσ 1 + a ˆσ ) 1/ t + + k 1;1 α ˆσ 1 k + z p [a 1 ˆσ 1 + a ˆσ ) +] 1/ ) 1/ } U v ), 16) where a 1 = 1 and a = ñ. The p,1 α) lower TL is given by Equation 16) with the plus sign after Ȳ replaced by the minus sign Noncentral t TLs The p,1 α) upper TL for a Nμ, στ ) distribution based on the noncentral t approximation [8] is given by SSȳ Ȳ + t k 1;1 α δ) kk 1), 17) where δ = z p { k 3.3. CLs for a survival probability ñkk 1) N k } SS 1/ e F k 1,N k;α. SSȳ Consider the one-way random model in Equation 5). Let θ = PY > T), where T is a specified value, and Y Nμ, στ + σ e ). In industrial hygiene, T is the occupational exposure limit OEL), and the probability is referred to as the exceedance probability. In the manufacturing industry, this probability is referred to as the unilateral conformance proportion, which is an important index for assessing the efficiency of a manufacturing process; see Lee and Liao.[19]A1 α lower CL for θ is the value of p for which the p,1 α) lower TL is equal to T. That is, the value of p that satisfies the equation Ȳ z p a 1 ˆσ 1 + a ˆσ t k 1;1 α ˆσ 1 k + z p a 1 ˆσ 1 + a ˆσ ) 1/ U v ) = T. It is easy to check that z p satisfies the above equation and it is given by + A = b b ac), 18) a where a = a 1 ˆσ 1 + a ˆσ a 1 ˆσ 1 + a ˆσ U v ), b = a 1 ˆσ 1 + a ˆσ )1/ Ȳ T) and c = Ȳ T) tk 1;1 α ˆσ 1 /k, and U v is given in Equation 10). The approximate lower CL for θ is given by A), where denotes the cumulative distribution function of the standard normal random variable.
9 Journal of Statistical Computation and Simulation 7 There are practical situations where one is interested in finding an upper CL for θ. For instance, in exposure data analysis, it is desired to find upper CL for the probability that the exposure level of a random worker exceeds the OEL. A 1 α upper CL for θ is obtained similarly by equating p,1 α) upper TL to T, and solving the resulting equation for p. This yields upper CL for θ as A ), with A = b b ac)/a, where a, b, and c are as defined in the preceding paragraph Coverage studies To determine the performance of the new TLs described in the preceding sections, and to compare with the noncentral t TLs, we estimated the coverage probabilities and expectations of both TLs by Monte Carlo method with 100,000 runs. Estimates are obtained for various values of k, n and the intraclass correlation ρ = στ /σ τ + σ e ). In our simulation study, without loss of generality, we set μ = 0 and στ + σ e = 1, so that μ + z p σ τ + σe = z p. The results in Table 1 correspond to.90,.95) upper TLs based on balanced data. Furthermore, we note that the sample size and parameter configurations are chosen as in Table of Krishnamoorthy and Mathew.[8] From Table 1, we see that the MOVER TLs are conservative for small k, n and ρ<.5. Even though the coverage probabilities of the MOVER TLs are larger than.95, their expected values are comparable, and shorter than the corresponding NCT and Vangel s TLs for some cases. For example, for k, n, ρ) = 3, 4,.5). In situations where the coverage probabilities of all TLs Table 1. Monte carlo estimates of the coverage probabilities and expected values for the p,1 α) NCT, MOVER and Vangel s upper TLs for Nμ, στ + σ e ) based on balanced data with μ =0,ρ = σ τ /σ τ + σ e ), σ τ + σ e =1, p, γ)= 0.90, 0.95); the expected values are given in parentheses. k, n) ρ NCT MOVER Vangel k, n) NCT MOVER Vangel k, n) NCT MOVER Vangel 3,4) ,) ,5) ) 3.11) 3.0) 1.84).00) 1.9) 1.85).00) 1.94) ) 3.64) 3.65) 1.93).03) 1.99).19).4).5) ) 4.31) 4.44).08).10).11).63).60).66) ) 4.71) 4.91).16).14).18).88).81).88) ) 5.19) 5.30).8).).7) 3.16) 3.05) 3.15) 3,10) ,10) ,4) ).6).5) 1.41) 1.54) 1.53) 1.8) 1.97) 1.9) ) 3.14) 3.13) 1.63) 1.70) 1.71).07).14).13) ) 4.09) 4.17) 1.91) 1.9) 1.95).41).38).44) ) 4.60) 4.75).07).04).09).58).5).63) ) 5.1) 5.35).5).0).6).84).74).81) 3, 10 3 ) ,10 3 ) ,0) ) 1.34).1) 1.8) 1.31) 1.30) 1.31) 1.38) 1.40) ).76) 3.16) 1.55) 1.60) 1.64) 1.4) 1.47) 1.50) ) 3.96) 4.19) 1.88) 1.88) 1.91) 1.56) 1.58) 1.61) ) 4.49) 4.7).04).0) ) 1.65) 1.68) ) 5.17) 5.34).5).0).6) 1.74) 1.73) 1.75)
10 8 K. Krishnamoorthy and J. Peng Table. Monte carlo estimates of coverage probabilities and expected value of the.90,.95) upper TL for Nμ, στ + σ e ) based on unbalanced data with μ =0,ρ = στ /σ τ + σ e ), σ τ + σ e =1; the expected values are given in parentheses. ρ n 1 n n 3 NCT MOVER n 1 n n 3 NCT MOVER n 1 n n 3 n 4 NCT MOVER ) 3.17).56).6).63).77) ) 3.68) 3.36) 3.34).9).96) ) 4.37) 4.6) 4.17) 3.30) 3.5) ) 4.78) 4.78) 4.66) 3.5) 3.4) ) 5.3) 5.41) 5.5) 3.77) 3.63) ).78) 1.79) 1.89) 1.47) 1.58) ) 3.45).96).95).19).0) ) 4.7) 4.09) 4.01).89).84) ) 4.69) 4.68) 4.56) 3.3) 3.) ) 5.4) 5.9) 5.13) 3.75) 3.61) a b c d ρ NCT MOVER NCT MOVER NCT MOVER NCT MOVER ) 1.57) 1.67) 1.80) 1.47) 1.56) 1.35) 1.44) ) 1.69) 1.78) 1.86) 1.64) 1.69) 1.56) 1.61) ) 1.86) 1.94) 1.95) 1.86) 1.86) 1.81) 1.8) ) 1.96).04).0) 1.99) 1.97) 1.96) 1.95) ).09).13).08).14).09).13).09) Notes: a, n = 3, 15, 30, 14,, 3, 13,, 8, 6, 9, 11);b,n = 3, 4, 3, 4,, 3, 3,,,,, ) c, n = 13, 40, 7, 14,, 30, 3,, 1,, 1, 4);d,n = 0, 30, 40, 30, 0, 10, 40, 0, 10, 5, 3, 40). are close to the nominal level, the MOVER TLs are expected to be smaller than the other TLs; see for the cases k, n, ρ) = 3, 4,.7) and 3, 10 3,.5). For moderate-to-large samples, the MOVER TLs are superior to the NCT TLs in terms of coverage probabilities and precision. We also observe that MOVER TL is quite comparable with the Vangel TL in most cases in terms of coverage probabilities and precision. However, the MOVER TLs are simple to compute and does not require any table values, whereas the Vangel TLs require table values the coefficients a i s in Equation 14)) and are available only for the balanced case. The estimated coverage probabilities reported in Table are for the unbalanced data. We once again see that the NCT TLs could be very liberal for smaller values of ρ where as the MOVER TLs are conservative. Nevertheless, the expected values of the MOVER TLs are comparable or smaller than those of the NCT TLs. For example, when ρ, n 1, n, n 3 ) =., 4, 6, 7), the coverage probability of the MOVER TL is.969 with the expected value 3.45, and those of the NCT TL are.94 and 3.45, respectively. For ρ.5, which could be the case in many applications, both TLs are satisfactory. Even in these cases, the expected values of the MOVER TLs are slightly smaller than the corresponding NCT TLs.
11 Journal of Statistical Computation and Simulation 9 Table 3. Monte carlo estimates of the coverage probabilities and expected values for the.90,.95) NCT and MOVER upper TLs for Nμ, στ ) based on balanced data with μ = 0, ρ = σ τ /σ τ + σ e ), σ τ + σ e = 1; the expected values are given in parentheses. ρ k, n) NCT MOVER k, n) NCT MOVER k, n) NCT MOVER ,) ,) ,4) ) 3.66) 1.34) 1.34) 0.68) 0.68) ) 4.0) 1.58) 1.51) 1.05) 1.00) ) 4.55) 1.88) 1.79) 1.43) 1.39) ) 4.86).05) 1.97) 1.64) 1.60) ) 5.).5).19) 1.86) 1.83) ,5) ,5) ,10) ).34) 0.8) 0.8) 0.43) 0.43) ) 3.14) 1.7) 1.1) 0.93) 0.90) ) 4.09) 1.7) 1.66) 1.38) 1.35) ) 4.57) 1.96) 1.90) 1.61) 1.58) ) 5.19).4).17) 1.85) 1.8) ,10) ,10) ,0) ) 0.48) 0.58) 0.58) 0.6) 0.6) ) 0.97) 1.15) 1.11) 0.83) 0.81) ) 1.44) 1.67) 1.6) 1.7) 1.5) ) 1.68) 1.94) 1.89) 1.48) 1.46) ) 1.94).3).17) 1.7) 1.70) Estimates of coverage probabilities and expected values of upper TLs for Nμ, στ ) are given in Table 3 balanced case) and Table 4 unbalanced case). These two table values indicate that both TLs are conservative for small values of ρ, and control the coverage probability satisfactorily for ρ.5. We see in general that the expected values of the MOVER TLs are slightly smaller than those of NCT TLs. Overall, we see that the MOVER TLs are not only simple to compute, but also they are very satisfactory. In many cases, we noticed that the MOVER TLs are better than those based on other methods. 4. A two-way nested model Consider the design of experiment where there are two factors A and B with the levels of B nested within the levels of A. Suppose there are a levels of A, and b levels of B are nested within each level of A, with n observations obtained on each combination of levels. Let Y ijl denote the lth observation corresponding to the jth level of B nested within the ith level of A. Thus, we have Y ijl = μ + τ i + β ji) + e ijl, i = 1,..., a; j = 1,..., b; l = 1,..., n, 19) where μ is an overall mean effect, τ i is the main effect due to the ith level of A, β ji) is the effect due to the jth level of B nested within the ith level of A, and e ijl s are the error terms with
12 10 K. Krishnamoorthy and J. Peng Table 4. Monte carlo estimates of the confidence level and expected value of the.90,.95) NCT and MOVER upper TLs for Nμ, στ ) based on unbalanced data with μ = 0.0, σ e = 1.0, Q = μ + z pσ τ, ρ = στ /σ τ + σ e ); the expected values are given in parentheses. ρ n 1 n n 3 NCT MOVER n 1 n n 3 n 4 n 5 NCT MOVER ) 3.5).7).7) ) 3.71).50).43) ) 4.36).78).67) ) 4.73).97).84) ) 5.1) 3.18) 3.06) ).93) 1.1) 1.0) ) 3.5) 1.79) 1.71) ) 4.4).44).34) ) 4.69).76).65) ) 5.15) 3.13) 3.01) e ijl N0, σe ). If the levels of A and B are randomly selected, then we have a random effects model, and in addition to the above distributional assumptions, we also assume τ i N0, στ ) and β ji) N0, σβ ), and τ i, β ji) and e ijl are all independent. The required statistics to find the MOVER TLs are defined as follows: SS τ = bn Ȳ ij. = 1 n n l=1 Y ijl, Ȳ i.. = 1 bn a Ȳ i.. Ȳ... ), SS β =n a b j=1 n l=1 Y ijl, Ȳ... = 1 abn a b Ȳ ij. Ȳ i.. ), and SS e = j=1 b j=1 a n Y ijl, l=1 b j=1 n Y ijl Ȳ ij. ). 0) Let σ1 = bnσ τ + nσ β + σ e, σ = nσ β + σ e and σ 3 = σ e so that σ τ = σ 1 σ )/bn), σ β = σ σ 3 )/n and σ e = σ 3. Under the model assumptions, we have l=1 ˆμ = Ȳ... N σ ) 1 μ,, abn SS τ σ 1 χ a 1, SS β σ χ ab 1), and SS e σ 3 χ abn 1), 1) where all the random variables are mutually independent. It follows from the above distributional results that ˆσ 1 = SS τ a 1, ˆσ = SS β ab 1), and ˆσ 3 = SS e abn 1) are independent unbiased estimators of σ 1, σ and σ 3, respectively.
13 4.1. TLs for Nμ, σ τ + σ β + σ e ) MOVER TLs Journal of Statistical Computation and Simulation 11 We shall find p,1 α)one-sided upper TL for the distribution Nμ, στ + σ β + σ e ), equivalently, 3 a1 α upper CL for μ + z p a iσi where a 1 = 1/bn), a = 1 1/b)/n and a 3 = 1 1/n). The upper CL for μ is U μ = Ȳ... + t a 1 ˆσ 1 abn, ) and the MOVER upper CL for 3 a iσi is given by U v = 3 a i ˆσ i + 3 ai ˆσ i 4 m i χm i ;α 1), 3) where m 1 = a 1, m = ab 1) and m 3 = abn 1). Finally, the p,1 α) MOVER upper TL for the Nμ, 3 a iσi ) distribution is expressed as 3 ) 1/ Ȳ... + z p a i ˆσ i + t a 1;1 α ˆσ 1 abn + z 3 p a i ˆσ i U v 1/. 4) A p,1 α) lower TL is obtained from Equation 4) with the plus sign after Ȳ... replaced by the minus sign Generalized TLs The p,1 α) upper TL based on the generalized variable approach [16] is expressed as follows. Let ȳ..., ss τ, ss β, ss e ) be an observed value of Ȳ..., SS τ, SS β, SS e ). In terms of these notations, the generalized pivotal quantity is given by ȳ... + Z ssτ χa 1 abn + z [ p ss τ bn χa 1 + b 1) ss β χ ab 1) + bn 1) ss e χ abn 1) ] 1/, 5) where Z N0, 1), and all the random variables in Equation 5) are mutually independent. The 1001 α) percentile is the desired upper TL for the Nμ, στ + σ β + σ e ) distribution. For agivenȳ..., ss τ, ss β, ss e ), the percentiles of Equation 5) can be estimated using Monte Carlo simulation. 4.. TLs for Nμ, σ τ + σ β ) MOVER TLs An upper TL for the distribution Nμ, στ + σ β ) can be readily obtained from Equation 4) as follows. Note that σβ + σ τ = a 1σ1 + a σ + a 3σ3, where a 1 = bn) 1, a = b 1)/nb), and
14 1 K. Krishnamoorthy and J. Peng a 3 = 1/n. The MOVER upper CL for 3 a iσi is given by U v = 3 a i ˆσ i + a 1 ˆσ 1 4 m 1 χm 1 ;α 1) + a ˆσ 4 m χm ;α 1) + a 3 ˆσ 4 3 m 3 χm 3 ;1 α 1), 6) where m i s are as defined in Equation 4). The MOVER upper TL for the distribution Nμ, σ τ + σ β ) is given by Equation 4) with U v defined in Equation 6) Generalized TLs The generalized upper TL due to Fonseca et al. [16] is given as follows. Let ȳ..., ss τ, ss β, ss e ) be an observed value of Ȳ..., SS τ, SS β, SS e ). In terms of these notations, the generalized pivotal quantity is given by ȳ... + Z ssτ χa 1 abn + z [ p ss τ bn χa 1 + b 1) ss β χ ab 1) b ss e χ abn 1) ] 1/ +, 7) where Z N0, 1), and all the random variables in Equation 7) are mutually independent. The 1001 α) percentile is the desired upper TL for the Nμ, στ + σ β ) distribution. For a given ȳ..., ss τ, ss β, ss e ), the percentiles of Equation 7) can be estimated using Monte Carlo simulation Coverage studies In order to judge the performance of the MOVER TLs for the two-way nested model, and to compare them with the generalized TLs proposed in Fonseca et al.,[16] we estimated the coverage probabilities of the.90,.95) TLs using Monte carlo simulation consisting of 100,000 runs. For simulation studies, we used the same sample size and parameter configurations given in Fonseca et al.[16] Specifically, the simulations were carried out fixing n = 0, and the quantities a and b were chosen to be 10 or 0. In Table 5, we reported the estimated coverage probabilities of Table 5. Coverage probabilities of.90,.95) upper TLs for Nμ, στ + σ β + σ e ) in a two-way nested model. ρ a = 5, b = 5 GV MOVER a = 5, b = 0 GV MOVER a = 0, b = 5 GV MOVER a = 0, b = 0 GV MOVER Notes: σβ = σ e = 1; ρ = σ τ /σ τ + σ β + σ e ); n = 0.
15 Journal of Statistical Computation and Simulation 13 Table 6. Coverage probabilities of.90,.95) upper TLs for Nμ, στ + σ β ) in a two-way nested model. ρ a = 5, b = 5 GV MOVER a = 5, b = 0 GV MOVER a = 0, b = 5 GV MOVER a = 0, b = 0 GV MOVER Notes: σβ = σ e = 1; ρ = σ τ /σ τ + σ β + σ e ); n = 0. MOVER TLs for Nμ, στ + σ β + σ e ) along with those for the generalized TLs reported in [16, Table ]) in Table 5. The estimated coverage probabilities clearly indicated that the MOVER TLs are less conservative than the generalized TLs. For the case of a = b = n = 0, the coverage probabilities of the MOVER TLs are very close to the nominal level.95. The estimated coverage of.90,.95) TLs for Nμ, στ + σ β ) are given in Table 6. The coverage probabilities of the generalized TLs are taken from Table of Fonseca et al.[16] We once again note that the MOVER TLs are less conservative than the generalized TLs, and in some cases the coverage probabilities are slightly less than.95, but not less than.94. Overall, we see that the MOVER TLs are less conservative than the generalized TLs, and they do not require simulation. So the MOVER TLs are preferable to the generalized TLs in terms of accuracy and simplicity. 5. Examples Example 1 This example is taken from Vangel,[6] and is concerned with tensile strength measurements made on k = 5 batches of composite materials. Each batch consists of n = 5 specimens. The data are given in Vangel [6, Table 4] and they are reproduced here in Table 7. The summary statistics to compute TLs are Ȳ 1. = 375.6, Ȳ. = 370.4, Ȳ 3. = 401.6, Ȳ 4. = 394.0, Ȳ 5. = 400., Ȳ = , SSȳ = 83.67, and SS e = Other quantities to compute the.90,.95) MOVER lower TL are a 1 = 1, a =.8, ˆσ 1 = 08.17, ˆσ = 78.9, t k 1;.95 = and U v = 136.0, where U v is computed using Equation 9). Using these values along with z.9 = 1.8 in Equation 1), we found the desired lower TL as This limit along with the ones [8] based on the generalized variable approach, noncentral t approximation, and Vangel s [6] approach are given in Table 7. Note that the MOVER TL is also the largest among the four. Since we are computing a lower TL, the larger the TL, the better. Suppose it is desired to assess the percentage of tensile strength measurements that exceed 350. The required quantities see Equation 18)) are a = , b = , c = 18.7, and
16 14 K. Krishnamoorthy and J. Peng Table 7. The tensile strength data and lower TLs. Batch Observations Methods.90,.95) lower TLs GV approach Noncentral t Vangel s MOVER Table 8. The moisture content data. Storage condition Moisture content %) 1 7.3, 8.3, 7.6, 8.4, , 7.4, , , 9.5, A = The 95% lower confidence for the exceedance probability is.9581) =.831. This means that at least 83% of tensile strengths exceed 350 with confidence 95%. Example This example is taken from Ostle and Mensing [0, p. 96] and is on a study of the effect of storage conditions on the moisture content of white pine lumber. Five different storage conditions k = 5) were studied with a varying number of sample boards stored under each condition. The example is also considered by Bhaumik and Kulkarni,[5] and the data are reproduced here in Table 8. These authors assumed the one-way random model for the purpose of estimating a 0.90, 0.95) upper TL for moisture content. We see from Table 8 that n 1 = 5, n = 3, n 3 =, n 4 = 3 and n 5 = 1. The summary statistics are Ȳ = 7.6, SSȳ = 3.80, SS e = 7.17, ñ =.4733, ˆσ 1 =.9500, and ˆσ = To compute the MOVER upper TL, we found the.95 upper CL for μ as and the.95 upper CL for a 1 σ 1 + a σ as 5.84, where a 1 = 1 and a = 1 ñ. Using these values along with z.9 = 1.8, we calculated the MOVER upper TL in Equation 16) as The NCT upper TL and the one based on the generalized variable approach are reported in Krishnamoorthy and Mathew [8] as and 11.1, respectively. We also evaluated 0.90, 0.95) upper TL for the distribution of the true moisture content, that is, the distribution of μ + τ i. the MOVER upper TL is computed as Krishnamoorthy and Mathew s [8] calculation showed that the NCT TL and the one based on the generalized variable approach are very close, given by We once again see that the MOVER TL is the smallest among these three TLs. Example 3 To show the applications of TLs in a two-way nested random model, we shall use the example of a breeding experiment given in Fonseca et al.[16] In this experiment, the breeding values of five sires in rising pigs were evaluated. Each sire was mated to two randomly selected dams, and the average weight gain of two pigs from each litter was recorded. Note that a = 5, b = and n = 5. Other summary statistics are ȳ 1.. =.67, ȳ.. =.53, ȳ 3.. =.63, ȳ 4.. =.47, ȳ 5.. =.57, ss β =.56, ss e =.39 ss τ =.05.
17 Journal of Statistical Computation and Simulation 15 Assuming that both factors are random, Fonseca et al. obtained the.90,.95) generalized upper TL for Nμ, στ + σ β + σ e ) as To find the.90,.95) MOVER upper TL, we calculated U v in Equation 3) as.1516, and Ȳ... + z p 3 a i ˆσ i )1/ as Using these numbers, and other statistics in Equation 4), we obtained the upper TL as , which is the same as generalized upper TL Fonseca et al. [16] have calculated the.90,.95) generalized upper TL for Nμ, στ + σ β ) as.99. To find the.90,.95) MOVER upper TL, we calculated U v in Equation 6) as.1074 and Ȳ... + z p 3 a i ˆσ i )1/ as.71. Finally, substituting these numbers in Equation 4), we obtained the upper TL as.99, which is the same as the generalized upper TL. 6. Concluding remarks In this article, we have developed simple closed-form one-sided TLs in one-way random models with balanced or unbalanced data. Although the proposed MOVER method is a large sample method, our simulation studies clearly indicate that the method works satisfactorily even for small samples. We also notice that the TLs based on the MOVER method are not only simple to compute, but also comparable with those based on other methods and produce better TLs in some cases, especially when sample sizes are moderate to large. Overall, we observe in our simulation study that the MOVER TLs are seldom liberal, and so they can be safely used in applications. It is evident that the performance of the MOVER TLs depends on the individual CIs for the overall mean effect μ and the variance components. Regarding the CI for μ, we have chosen the usual one; however, there are choices for a linear combination of variance components. For example, one could use the generalized CIs proposed in Liao et al. [9] or the ones based on the Satterthwaite, [1] approximation to find a CI for a i σi. However, as noted in Graybill and Wang,[17] the MOVER CIs are better than the ones based the Satterthwaite approximation, and the MOVER method is applicable even when some a i s are negative. Our preliminary coverage studies of TLs based on the Satterthwaite CIs of a i σi ) for Nμ, σ τ + σ e ) indicated that such TLs are liberal for ρ.. The generalized CIs proposed in Liao et al. [9] are numerically involved, and require simulation to obtain them. Furthermore, as noted in the illustrative examples in Krishnamoorthy and Mathew [11, Chapters 4 and 5] and Krishnamoorthy and Lian,[] coverage properties of the MOVER CIs and those of the generalized CIs for a i σ i the MOVER CIs for a i σ i should be similar. As are easy to compute, we recommend the TLs based on them. References [1] Lyles RH, Kupper LL, Rappaport SM. Assessing regulatory compliance of occupational exposures via the balanced one-way random effects ANOVA model. J Agric Biol Environ Stat. 1997a;: [] Lyles RH, Kupper LL, Rappaport SM. A lognormal distribution based exposure assessment method for unbalanced data. Ann Occup Hyg. 1997b;41: [3] Mee RW, Owen DB. Improved factors for one-sided tolerance limits for balanced one-way ANOVA random model. J Amer Statist Assoc. 1983;78: [4] Beckman RJ, Tietjen GL. Two-sided tolerance limits for balanced random-effects ANOVA models Correction: 1998, 40, p. 69). Technometrics. 1989;31: [5] Bhaumik DK, Kulkarni PM. One-sided tolerance limits for unbalanced one-way ANOVA random effects model. Comm Statist Theory Methods. 1991;0: [6] Vangel MG. New methods for one-sided tolerance limits for a one-way balanced random-effects ANOVA model. Technometrics. 199;34: [7] Liao C-T, Iyer HK. A tolerance interval for the normal distribution with several variance components. Statist Sinica. 004;14:17 9. [8] Krishnamoorthy K, Mathew T. One-sided tolerance limits in balanced and unbalanced one-way random models based on generalized confidence limits. Technometrics. 004;46:44 5.
18 16 K. Krishnamoorthy and J. Peng [9] Liao C-T, Lin TY, Iyer HK. One- and two-sided tolerance intervals for general balanced mixed models and unbalanced one-way random models. Technometrics. 005;47: [10] Hoffman D. One-sided tolerance limits for balanced and unbalanced random effects models. Technometrics. 010;5: [11] Krishnamoorthy K, Mathew T. Statistical tolerance regions: theory, applications and computation. Hoboken, NJ: Wiley; 009. [1] Zou GY, Donner A. Construction of confidence limits about effect measures: a general approach. Stat Med. 008;7: [13] Zou GY, Taleban J, Huo CY. Simple confidence intervals for lognormal means and their differences with environmental applications. Environmetrics. 009a;0: [14] Zou GY, Huo CY, Taleban J. Confidence interval estimation for lognormal data with application to health economics. Comput Statist Data Anal. 009b;53: [15] Zou GY. Confidence interval estimation for the Bland-Altman limits of agreement with multiple observations per individual. Stat Methods Med Res. 011;: doi: / [16] Fonseca M, Mathew T, Mexia JT, Zmyslony R. Tolerance intervals in a two-way nested model with mixed or random effects. Statistics. 007;41: [17] Graybill FA, Wang CM. Confidence intervals on nonnegative linear combinations of variances. J Amer Statist Assoc. 1980;75: [18] Thomas JD, Hultquist RA. Interval estimation for the unbalanced case of the one-way random effects model. Ann Statist. 1978;6: [19] Lee H-I, Liao C-T. Interval estimation for conformance proportions in a normal variance components models. J Qual Technol. 01;44: [0] Ostle B, Mensing RW. Statistics in research: basic concepts and techniques for research workers. 3rd ed. Ames, IA: Iowa State University Press; [1] Satterthwaite FE. An approximate distribution of estimates of variance components. Biom Bull. 1946;: [] Krishnamoorthy K, Lian X. Closed-form approximate tolerance intervals for some general linear models and comparison studies. J Stat Comput Simul. 01;8:
Assessing occupational exposure via the one-way random effects model with unbalanced data
Assessing occupational exposure via the one-way random effects model with unbalanced data K. Krishnamoorthy 1 and Huizhen Guo Department of Mathematics University of Louisiana at Lafayette Lafayette, LA
More informationA Simple Approximate Procedure for Constructing Binomial and Poisson Tolerance Intervals
This article was downloaded by: [Kalimuthu Krishnamoorthy] On: 11 February 01, At: 08:40 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 107954 Registered office:
More informationPublished online: 17 May 2012.
This article was downloaded by: [Central University of Rajasthan] On: 03 December 014, At: 3: Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 107954 Registered
More informationUse and Abuse of Regression
This article was downloaded by: [130.132.123.28] On: 16 May 2015, At: 01:35 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationGilles Bourgeois a, Richard A. Cunjak a, Daniel Caissie a & Nassir El-Jabi b a Science Brunch, Department of Fisheries and Oceans, Box
This article was downloaded by: [Fisheries and Oceans Canada] On: 07 May 2014, At: 07:15 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationModified Normal-based Approximation to the Percentiles of Linear Combination of Independent Random Variables with Applications
Communications in Statistics - Simulation and Computation ISSN: 0361-0918 (Print) 1532-4141 (Online) Journal homepage: http://www.tandfonline.com/loi/lssp20 Modified Normal-based Approximation to the Percentiles
More informationOnline publication date: 12 January 2010
This article was downloaded by: [Zhang, Lanju] On: 13 January 2010 Access details: Access Details: [subscription number 918543200] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationUniversity, Wuhan, China c College of Physical Science and Technology, Central China Normal. University, Wuhan, China Published online: 25 Apr 2014.
This article was downloaded by: [0.9.78.106] On: 0 April 01, At: 16:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 10795 Registered office: Mortimer House,
More informationOnline publication date: 01 March 2010 PLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [2007-2008-2009 Pohang University of Science and Technology (POSTECH)] On: 2 March 2010 Access details: Access Details: [subscription number 907486221] Publisher Taylor
More informationTesting Goodness-of-Fit for Exponential Distribution Based on Cumulative Residual Entropy
This article was downloaded by: [Ferdowsi University] On: 16 April 212, At: 4:53 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationImproved Confidence Intervals for the Ratio of Coefficients of Variation of Two Lognormal Distributions
Journal of Statistical Theory and Applications, Vol. 16, No. 3 (September 017) 345 353 Improved Confidence Intervals for the Ratio of Coefficients of Variation of Two Lognormal Distributions Md Sazib Hasan
More informationFull terms and conditions of use:
This article was downloaded by:[rollins, Derrick] [Rollins, Derrick] On: 26 March 2007 Access Details: [subscription number 770393152] Publisher: Taylor & Francis Informa Ltd Registered in England and
More informationPark, Pennsylvania, USA. Full terms and conditions of use:
This article was downloaded by: [Nam Nguyen] On: 11 August 2012, At: 09:14 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationDistribution Theory. Comparison Between Two Quantiles: The Normal and Exponential Cases
Communications in Statistics Simulation and Computation, 34: 43 5, 005 Copyright Taylor & Francis, Inc. ISSN: 0361-0918 print/153-4141 online DOI: 10.1081/SAC-00055639 Distribution Theory Comparison Between
More informationGeorge L. Fischer a, Thomas R. Moore b c & Robert W. Boyd b a Department of Physics and The Institute of Optics,
This article was downloaded by: [University of Rochester] On: 28 May 2015, At: 13:34 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationPublished online: 05 Oct 2006.
This article was downloaded by: [Dalhousie University] On: 07 October 2013, At: 17:45 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationOpen problems. Christian Berg a a Department of Mathematical Sciences, University of. Copenhagen, Copenhagen, Denmark Published online: 07 Nov 2014.
This article was downloaded by: [Copenhagen University Library] On: 4 November 24, At: :7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 72954 Registered office:
More informationThe American Statistician Publication details, including instructions for authors and subscription information:
This article was downloaded by: [National Chiao Tung University 國立交通大學 ] On: 27 April 2014, At: 23:13 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationModified Large Sample Confidence Intervals for Poisson Distributions: Ratio, Weighted Average and Product of Means
Modified Large Sample Confidence Intervals for Poisson Distributions: Ratio, Weighted Average and Product of Means K. KRISHNAMOORTHY a, JIE PENG b AND DAN ZHANG a a Department of Mathematics, University
More informationNEW APPROXIMATE INFERENTIAL METHODS FOR THE RELIABILITY PARAMETER IN A STRESS-STRENGTH MODEL: THE NORMAL CASE
Communications in Statistics-Theory and Methods 33 (4) 1715-1731 NEW APPROXIMATE INFERENTIAL METODS FOR TE RELIABILITY PARAMETER IN A STRESS-STRENGT MODEL: TE NORMAL CASE uizhen Guo and K. Krishnamoorthy
More informationUniversity, Tempe, Arizona, USA b Department of Mathematics and Statistics, University of New. Mexico, Albuquerque, New Mexico, USA
This article was downloaded by: [University of New Mexico] On: 27 September 2012, At: 22:13 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationCommunications in Algebra Publication details, including instructions for authors and subscription information:
This article was downloaded by: [Professor Alireza Abdollahi] On: 04 January 2013, At: 19:35 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationInference on reliability in two-parameter exponential stress strength model
Metrika DOI 10.1007/s00184-006-0074-7 Inference on reliability in two-parameter exponential stress strength model K. Krishnamoorthy Shubhabrata Mukherjee Huizhen Guo Received: 19 January 2005 Springer-Verlag
More informationOnline publication date: 30 March 2011
This article was downloaded by: [Beijing University of Technology] On: 10 June 2011 Access details: Access Details: [subscription number 932491352] Publisher Taylor & Francis Informa Ltd Registered in
More informationOnline publication date: 22 March 2010
This article was downloaded by: [South Dakota State University] On: 25 March 2010 Access details: Access Details: [subscription number 919556249] Publisher Taylor & Francis Informa Ltd Registered in England
More informationGuangzhou, P.R. China
This article was downloaded by:[luo, Jiaowan] On: 2 November 2007 Access Details: [subscription number 783643717] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number:
More informationCOMPARISON OF FIVE TESTS FOR THE COMMON MEAN OF SEVERAL MULTIVARIATE NORMAL POPULATIONS
Communications in Statistics - Simulation and Computation 33 (2004) 431-446 COMPARISON OF FIVE TESTS FOR THE COMMON MEAN OF SEVERAL MULTIVARIATE NORMAL POPULATIONS K. Krishnamoorthy and Yong Lu Department
More informationDissipation Function in Hyperbolic Thermoelasticity
This article was downloaded by: [University of Illinois at Urbana-Champaign] On: 18 April 2013, At: 12:23 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationCCSM: Cross correlogram spectral matching F. Van Der Meer & W. Bakker Published online: 25 Nov 2010.
This article was downloaded by: [Universiteit Twente] On: 23 January 2015, At: 06:04 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationAnkara, Turkey Published online: 20 Sep 2013.
This article was downloaded by: [Bilkent University] On: 26 December 2013, At: 12:33 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationCharacterizations of Student's t-distribution via regressions of order statistics George P. Yanev a ; M. Ahsanullah b a
This article was downloaded by: [Yanev, George On: 12 February 2011 Access details: Access Details: [subscription number 933399554 Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationDerivation of SPDEs for Correlated Random Walk Transport Models in One and Two Dimensions
This article was downloaded by: [Texas Technology University] On: 23 April 2013, At: 07:52 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationNacional de La Pampa, Santa Rosa, La Pampa, Argentina b Instituto de Matemática Aplicada San Luis, Consejo Nacional de Investigaciones Científicas
This article was downloaded by: [Sonia Acinas] On: 28 June 2015, At: 17:05 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationThe Homogeneous Markov System (HMS) as an Elastic Medium. The Three-Dimensional Case
This article was downloaded by: [J.-O. Maaita] On: June 03, At: 3:50 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 07954 Registered office: Mortimer House,
More informationTo cite this article: Edward E. Roskam & Jules Ellis (1992) Reaction to Other Commentaries, Multivariate Behavioral Research, 27:2,
This article was downloaded by: [Memorial University of Newfoundland] On: 29 January 2015, At: 12:02 Publisher: Routledge Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationOF SCIENCE AND TECHNOLOGY, TAEJON, KOREA
This article was downloaded by:[kaist Korea Advanced Inst Science & Technology] On: 24 March 2008 Access Details: [subscription number 731671394] Publisher: Taylor & Francis Informa Ltd Registered in England
More informationPrecise Large Deviations for Sums of Negatively Dependent Random Variables with Common Long-Tailed Distributions
This article was downloaded by: [University of Aegean] On: 19 May 2013, At: 11:54 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationThe Fourier transform of the unit step function B. L. Burrows a ; D. J. Colwell a a
This article was downloaded by: [National Taiwan University (Archive)] On: 10 May 2011 Access details: Access Details: [subscription number 905688746] Publisher Taylor & Francis Informa Ltd Registered
More informationDiatom Research Publication details, including instructions for authors and subscription information:
This article was downloaded by: [Saúl Blanco] On: 26 May 2012, At: 09:38 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer House,
More informationErciyes University, Kayseri, Turkey
This article was downloaded by:[bochkarev, N.] On: 7 December 27 Access Details: [subscription number 746126554] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number:
More informationDresden, Dresden, Germany Published online: 09 Jan 2009.
This article was downloaded by: [SLUB Dresden] On: 11 December 2013, At: 04:59 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationVersion of record first published: 01 Sep 2006.
This article was downloaded by: [University of Miami] On: 27 November 2012, At: 08:47 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationG. S. Denisov a, G. V. Gusakova b & A. L. Smolyansky b a Institute of Physics, Leningrad State University, Leningrad, B-
This article was downloaded by: [Institutional Subscription Access] On: 25 October 2011, At: 01:35 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationMSE Performance and Minimax Regret Significance Points for a HPT Estimator when each Individual Regression Coefficient is Estimated
This article was downloaded by: [Kobe University] On: 8 April 03, At: 8:50 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 07954 Registered office: Mortimer House,
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Uniwersytet Slaski] On: 14 October 2008 Access details: Access Details: [subscription number 903467288] Publisher Taylor & Francis Informa Ltd Registered in England and
More informationFull terms and conditions of use:
This article was downloaded by:[smu Cul Sci] [Smu Cul Sci] On: 28 March 2007 Access Details: [subscription number 768506175] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered
More informationConfidence Intervals for Process Capability Indices Using Bootstrap Calibration and Satterthwaite s Approximation Method
CHAPTER 4 Confidence Intervals for Process Capability Indices Using Bootstrap Calibration and Satterthwaite s Approximation Method 4.1 Introduction In chapters and 3, we addressed the problem of comparing
More informationGeometric View of Measurement Errors
This article was downloaded by: [University of Virginia, Charlottesville], [D. E. Ramirez] On: 20 May 205, At: 06:4 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number:
More informationApproximate and Fiducial Confidence Intervals for the Difference Between Two Binomial Proportions
Approximate and Fiducial Confidence Intervals for the Difference Between Two Binomial Proportions K. Krishnamoorthy 1 and Dan Zhang University of Louisiana at Lafayette, Lafayette, LA 70504, USA SUMMARY
More informationAcyclic, Cyclic and Polycyclic P n
This article was downloaded by: [German National Licence 2007] On: 15 December 2010 Access details: Access Details: [subscription number 777306419] Publisher Taylor & Francis Informa Ltd Registered in
More informationTong University, Shanghai , China Published online: 27 May 2014.
This article was downloaded by: [Shanghai Jiaotong University] On: 29 July 2014, At: 01:51 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationFB 4, University of Osnabrück, Osnabrück
This article was downloaded by: [German National Licence 2007] On: 6 August 2010 Access details: Access Details: [subscription number 777306420] Publisher Taylor & Francis Informa Ltd Registered in England
More informationTolerance limits for a ratio of normal random variables
Tolerance limits for a ratio of normal random variables Lanju Zhang 1, Thomas Mathew 2, Harry Yang 1, K. Krishnamoorthy 3 and Iksung Cho 1 1 Department of Biostatistics MedImmune, Inc. One MedImmune Way,
More informationMarko A. A. Boon a, John H. J. Einmahl b & Ian W. McKeague c a Department of Mathematics and Computer Science, Eindhoven
This article was downloaded by: [Columbia University] On: 8 March 013, At: 1: Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 107954 Registered office: Mortimer
More informationUniversity of Thessaloniki, Thessaloniki, Greece
This article was downloaded by:[bochkarev, N.] On: 14 December 2007 Access Details: [subscription number 746126554] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number:
More informationFrance. Published online: 23 Apr 2007.
This article was downloaded by: [the Bodleian Libraries of the University of Oxford] On: 29 August 213, At: 6:39 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number:
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationPublished online: 10 Apr 2012.
This article was downloaded by: Columbia University] On: 23 March 215, At: 12:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationGeometrical optics and blackbody radiation Pablo BenÍTez ab ; Roland Winston a ;Juan C. Miñano b a
This article was downloaded by: [University of California, Merced] On: 6 May 2010 Access details: Access Details: [subscription number 918975015] ublisher Taylor & Francis Informa Ltd Registered in England
More informationOn Selecting Tests for Equality of Two Normal Mean Vectors
MULTIVARIATE BEHAVIORAL RESEARCH, 41(4), 533 548 Copyright 006, Lawrence Erlbaum Associates, Inc. On Selecting Tests for Equality of Two Normal Mean Vectors K. Krishnamoorthy and Yanping Xia Department
More informationBy Godase, Shirke, Kashid. Published: 26 April 2017
Electronic Journal of Applied Statistical Analysis EJASA, Electron. J. App. Stat. Anal. http://siba-ese.unisalento.it/index.php/ejasa/index e-issn: 2070-5948 DOI: 10.1285/i20705948v10n1p29 Tolerance intervals
More informationON COMBINING CORRELATED ESTIMATORS OF THE COMMON MEAN OF A MULTIVARIATE NORMAL DISTRIBUTION
ON COMBINING CORRELATED ESTIMATORS OF THE COMMON MEAN OF A MULTIVARIATE NORMAL DISTRIBUTION K. KRISHNAMOORTHY 1 and YONG LU Department of Mathematics, University of Louisiana at Lafayette Lafayette, LA
More informationPLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:
This article was downloaded by: [Stanford University] On: 20 July 2010 Access details: Access Details: [subscription number 917395611] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationA Note on Coverage Probability of Confidence Interval for the Difference between Two Normal Variances
pplied Mathematical Sciences, Vol 6, 01, no 67, 3313-330 Note on Coverage Probability of Confidence Interval for the Difference between Two Normal Variances Sa-aat Niwitpong Department of pplied Statistics,
More informationA note on adaptation in garch models Gloria González-Rivera a a
This article was downloaded by: [CDL Journals Account] On: 3 February 2011 Access details: Access Details: [subscription number 922973516] Publisher Taylor & Francis Informa Ltd Registered in England and
More informationLecture 2: Basic Concepts and Simple Comparative Experiments Montgomery: Chapter 2
Lecture 2: Basic Concepts and Simple Comparative Experiments Montgomery: Chapter 2 Fall, 2013 Page 1 Random Variable and Probability Distribution Discrete random variable Y : Finite possible values {y
More informationDiscussion on Change-Points: From Sequential Detection to Biology and Back by David Siegmund
This article was downloaded by: [Michael Baron] On: 2 February 213, At: 21:2 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Sun Yat-Sen University] On: 1 October 2008 Access details: Access Details: [subscription number 781195905] Publisher Taylor & Francis Informa Ltd Registered in England
More informationMean. Pranab K. Mitra and Bimal K. Sinha. Department of Mathematics and Statistics, University Of Maryland, Baltimore County
A Generalized p-value Approach to Inference on Common Mean Pranab K. Mitra and Bimal K. Sinha Department of Mathematics and Statistics, University Of Maryland, Baltimore County 1000 Hilltop Circle, Baltimore,
More informationA. H. Abuzaid a, A. G. Hussin b & I. B. Mohamed a a Institute of Mathematical Sciences, University of Malaya, 50603,
This article was downloaded by: [University of Malaya] On: 09 June 2013, At: 19:32 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationInferences on a Normal Covariance Matrix and Generalized Variance with Monotone Missing Data
Journal of Multivariate Analysis 78, 6282 (2001) doi:10.1006jmva.2000.1939, available online at http:www.idealibrary.com on Inferences on a Normal Covariance Matrix and Generalized Variance with Monotone
More informationAccepted author version posted online: 28 Nov 2012.
This article was downloaded by: [University of Minnesota Libraries, Twin Cities] On: 15 May 2013, At: 12:23 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationPublished online: 27 Aug 2014.
This article was downloaded by: [49.50.78.27] On: 26 August 2015, At: 22:51 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: 5 Howick
More informationAndrás István Fazekas a b & Éva V. Nagy c a Hungarian Power Companies Ltd., Budapest, Hungary. Available online: 29 Jun 2011
This article was downloaded by: [András István Fazekas] On: 25 July 2011, At: 23:49 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [University of Santiago de Compostela] On: 6 June 2009 Access details: Access Details: [subscription number 908626806] Publisher Taylor & Francis Informa Ltd Registered
More informationThe use of the Normalized Difference Water Index (NDWI) in the delineation of open water features
This article was downloaded by: [University of Thessaly] On: 07 May 2015, At: 06:17 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Los Alamos National Laboratory] On: 21 July 2009 Access details: Access Details: [subscription number 908033413] Publisher Taylor & Francis Informa Ltd Registered in England
More informationTore Henriksen a & Geir Ulfstein b a Faculty of Law, University of Tromsø, Tromsø, Norway. Available online: 18 Feb 2011
This article was downloaded by: [Bibliotheek van het Vredespaleis] On: 03 May 2012, At: 03:44 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationA process capability index for discrete processes
Journal of Statistical Computation and Simulation Vol. 75, No. 3, March 2005, 175 187 A process capability index for discrete processes MICHAEL PERAKIS and EVDOKIA XEKALAKI* Department of Statistics, Athens
More informationPLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:
This article was downloaded by: On: 2 January 211 Access details: Access Details: Free Access Publisher Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered
More informationOnline publication date: 29 April 2011
This article was downloaded by: [Volodin, Andrei] On: 30 April 2011 Access details: Access Details: [subscription number 937059003] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationBootstrap Procedures for Testing Homogeneity Hypotheses
Journal of Statistical Theory and Applications Volume 11, Number 2, 2012, pp. 183-195 ISSN 1538-7887 Bootstrap Procedures for Testing Homogeneity Hypotheses Bimal Sinha 1, Arvind Shah 2, Dihua Xu 1, Jianxin
More informationChiao Tung University, Taiwan b Department of Industrial Management, National Taiwan
This article was downloaded by: [National Chiao Tung University 國立交通大學 ] On: 6 April 14, At: :59 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 17954 Registered
More informationT. Runka a, M. Kozielski a, M. Drozdowski a & L. Szczepańska b a Institute of Physics, Poznan University of
This article was downloaded by: [Politechnika Poznanska] On: 04 October 2013, t: 04:02 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationTaiwan Published online: 11 Sep To link to this article:
This article was downloaded by: [National Chiao Tung University 國立交通大學 ] On: 8 April 04, At: 6:09 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 07954 Registered
More informationC.K. Li a a Department of Mathematics and
This article was downloaded by: [Brandon University GST] On: 31 August 2013, At: 07:46 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationTOLERANCE INTERVALS FOR DISCRETE DISTRIBUTIONS IN EXPONENTIAL FAMILIES
Statistica Sinica 19 (2009), 905-923 TOLERANCE INTERVALS FOR DISCRETE DISTRIBUTIONS IN EXPONENTIAL FAMILIES Tianwen Tony Cai and Hsiuying Wang University of Pennsylvania and National Chiao Tung University
More informationMelbourne, Victoria, 3010, Australia. To link to this article:
This article was downloaded by: [The University Of Melbourne Libraries] On: 18 October 2012, At: 21:54 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationPLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:
This article was downloaded by: [Ferdowsi University of Mashhad] On: 7 June 2010 Access details: Access Details: [subscription number 912974449] Publisher Taylor & Francis Informa Ltd Registered in England
More informationPublished online: 04 Oct 2006.
This article was downloaded by: [Max-Planck Institut Fur Polymerforschung] On: 04 August 2015, At: 02:21 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationOsnabrück, Germany. Nanotechnology, Münster, Germany
This article was downloaded by:[universitatsbibliothek Munster] On: 12 December 2007 Access Details: [subscription number 769143588] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales
More informationHorizonte, MG, Brazil b Departamento de F sica e Matem tica, Universidade Federal Rural de. Pernambuco, Recife, PE, Brazil
This article was downloaded by:[cruz, Frederico R. B.] On: 23 May 2008 Access Details: [subscription number 793232361] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered
More information35-959, Rzeszów, Poland b Institute of Computer Science, Jagiellonian University,
This article was downloaded by: [Polska Akademia Nauk Instytut Matematyczny] On: 07 March 03, At: 03:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 0795 Registered
More informationSimultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions
Statistics Preprints Statistics 10-2014 Simultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions Yimeng Xie Virginia Tech Yili Hong Virginia Tech Luis A. Escobar Louisiana
More informationResearch Article A Nonparametric Two-Sample Wald Test of Equality of Variances
Advances in Decision Sciences Volume 211, Article ID 74858, 8 pages doi:1.1155/211/74858 Research Article A Nonparametric Two-Sample Wald Test of Equality of Variances David Allingham 1 andj.c.w.rayner
More informationP. C. Wason a & P. N. Johnson-Laird a a Psycholinguistics Research Unit and Department of
This article was downloaded by: [Princeton University] On: 24 February 2013, At: 11:52 Publisher: Routledge Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationA nonparametric two-sample wald test of equality of variances
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 211 A nonparametric two-sample wald test of equality of variances David
More informationQuan Wang a, Yuanxin Xi a, Zhiliang Wan a, Yanqing Lu a, Yongyuan Zhu a, Yanfeng Chen a & Naiben Ming a a National Laboratory of Solid State
This article was downloaded by: [Nanjing University] On: 07 April 2012, At: 19:39 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationA Strongly Convergent Method for Nonsmooth Convex Minimization in Hilbert Spaces
This article was downloaded by: [IMPA Inst de Matematica Pura & Aplicada] On: 11 November 2011, At: 05:10 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More information