Tests Concerning Equicorrelation Matrices with Grouped Normal Data
|
|
- Brice O’Connor’
- 5 years ago
- Views:
Transcription
1 Tests Concerning Equicorrelation Matrices with Grouped Normal Data Paramjit S Gill Department of Mathematics and Statistics Okanagan University College Kelowna, BC, Canada, VV V7 pgill@oucbcca Sarath G Banneheka Department of Statistics and Computer Science University of Sri Jayewardenepura Gangodawila, Nugegoda, Sri Lanka bannehek@sjpaclk Tim B Swartz Department of Statistics and Actuarial Science Simon Fraser University Burnaby, BC, Canada, V5A S6 tim@statsfuca Key Words: intraclass correlation; likelihood ratio test; maximum likelihood estimation; score test; small-sample inference ABSTRACT This paper considers three practical hypotheses involving the equicorrelation matrix for grouped normal data We obtain statistics and computing formulae for common test procedures such as the score test and the likelihood ratio test In addition, statistics and computing formulae are obtained for various small sample procedures as proposed in Skovgaard (200) The properties of the tests for each of the three hypotheses are compared using Monte Carlo simulations
2 Introduction A common data structure in applied statistics occurs when a small number of multivariate measurements are collected from a finite number of groups or classes For example, instruments known as computer-analyzed corneal topographers (CACT) are used to obtain multivariate observations on curvature points where the sample size is typically around 0 (Viana, Olkin and McMahon 993) In such situations, a parsimonious model for the covariance matrix is desired Also, in order to increase sample size, one may wish to combine data from different sources (eg instruments) Therefore, there is a need to test the equality of covariance matrices arising from different groups In this paper we consider tests of hypotheses regarding a special form of the covariance matrix called the equicorrelation or intraclass correlation matrix Such a covariance matrix is suitable when subjects are related (such as a family or a litter of animals) or when measurements are made repeatedly on the same subject More formally, suppose that we have samples of p-dimensional measurements from K groups Let n i be the number of observations from group i =, 2,, K and let n = K i= n i Let X i = x i x i2 x ip x i2 x i22 x i2p x ini x ini 2 x ini p denote the i th data matrix with x ij = (x ij, x ij2,,x ijp ) as the vector of observations on the j th sample from group i Assume x ij N p (µ i, Σ i ), i =, 2,, K, j =, 2,, n i where denotes the unit vector of dimension p and the equicorrelation or intraclass correlation matrix is given by Σ i = σ 2 i ρ i ρ i ρ i ρ i, ( p) < ρ i < ρ i ρ i 2
3 Note that for each group, the measurements share a common mean and correlation Although this is a restrictive structure, it is still appropriate in many contexts The intraclass correlation coefficient ρ i measures the degree of resemblance between the members in the i th group This covariance pattern also arises in the one-way random effects (within each group i) linear model x ijk = µ i + α ij + ε ijk where α ij is the random effect due to the j th unit from group i with E(α ij ) = 0, V ar(α ij ) = σ 2 iα, ε ijk is the random error with E(ε ijk ) = 0, V ar(ε ijk ) = σ 2 iε and Cov(α ij, ε ijk ) = 0 It is easily seen that the covariance matrix of x ij is Σ i with σ 2 i = σ 2 iα+σ 2 iε and ρ i = σ 2 iα/(σ 2 iα+σ 2 iε) (see, for example, Donner and Koval 980) We are interested in the following three practical hypotheses concerning the equicorrelation matrix Case : H 0 : ρ = = ρ K Case 2: H 02 : σ = = σ K Case 3: H 03 : ρ = = ρ K, σ = = σ K The Case problem was first addressed by Konishi and Gupta (989) using a maximum likelihood approach with combined moment estimators for the common ρ Paul and Barnwal (990) derived a score-type test based on Neyman s C(α) While estimating the intraclass correlation coefficient, Ahmed, Gupta, Khan and Nicol (200) also derived a score-type test statistic for the Case problem The Case 2 problem does not seem to have been studied For the Case 3 problem, Han (975) obtained a modified likelihood ratio test Our main concern regarding the Case, Case 2 and Case 3 problems involves data sets where the sample sizes are small For small sample sizes, asymptotic tests such as those described above may be inappropriate For example, the actual type error of a statistic may differ considerably from its nominal level In this regard, we have derived various small sample tests based on the general theory of Skovgaard (200) and we have compared these tests to the standard tests Some of these standard tests have not previously appeared in 3
4 the literature and are derived here This paper serves as a summary of the variety of tests that one might consider in connection with the Case, Case 2 and Case 3 problems Small sample tests are typically derived by considering additional terms in asymptotic expansions These corrections attempt to hasten the convergence of statistics to their asymptotic distributions The downside of such corrections is that they sometimes involve difficult calculations Well known examples of correction terms include Bartlett and Edgeworth series type corrections, saddlepoint corrections and the signed loglikelihood ratio correction Reviews of asymptotic methods in statistics are given by Barndorff-Nielsen and Cox (994), Skovgaard (200) and Reid (2003) Notably, Skovgaard (200) generalized the Barndorff-Nielsen (99) correction to the multi-parameter case; we follow Skovgaard s approach to derive small sample corrections for likelihood ratio tests involving equicorrelation matrices In section 2, we provide a general description of the modified likelihood ratio test due to Skovgaard (200) and the score test We also discuss unrestricted maximum likelihood estimation for the problem considered in the paper In sections 3, 4 and 5, Cases, 2 and 3 are discussed in detail with an emphasis on deriving asymptotic test procedures In section 6, simulation results are provided to compare the various tests under each of the three hypotheses 2 PRELIMINARIES 2 The Skovgaard Modifications Suppose that a random variable X belongs to the exponential family having a density of the form f(x; θ) = b(x) exp{θ t(x) κ(θ)} Then θ is called the canonical parameter and t is the corresponding canonical sufficient statistic Let ω = (ν, ψ ) be a (possibly nonlinear) function of θ where ψ ν = (ψ,, ψ q ) is the vector of parameters of interest and = (ν,,ν r q ) is a vector of nuisance parameters Suppose one is interested in testing the null hypothesis H 0 : ψ = = ψ q = ψ 0 (unspecified) Let ω = ( ν, ψ ) denote the maximum likelihood (ML) estimate of the full parameter vector and ω 0 =, ψ ( ν 0 ) be the 4
5 maximum likelihood estimate of the free parameters under the null hypothesis Let l( ω) and l( ω 0 ) be the corresponding maximum values of the loglikelihood Denote J νν and I νν as the observed and expected Fisher information matrices for the ν-part of the parameter vector and let J νν and Ĩνν denote the corresponding estimates obtained by replacing the unknown parameter values by their restricted ML estimates Further, let τ and Ω denote the mean and variance matrix of t; τ and Ω are the corresponding estimates when the unknown parameters are replaced by the restricted ML estimates and Ω is the estimate under unrestricted ML estimation Let w = 2(l( ω) l( ω 0 )) denote the likelihood ratio test statistics for testing H 0 Following Skovgaard (200), we define δ = {(t τ) Ω (t τ)} (q )/2 w (q )/2 (ˆθ θ) (t τ) The corrected likelihood ratio is then given by ( Ω J ) /2 νν () Ω Ĩνν w = w ( w ln δ ) 2 (2) where it is hoped that w χ 2 q is a better approximation under H 0 than the likelihood ratio test w χ 2 q Another corrected χ2 q test statistic (Skovgaard 200) is which may assume a negative value, especially when w is small w = w 2 lnδ (3) 22 The Score Test The score test is attractive for its relative simplicity of computations as we need maximum likelihood estimation only under the null hypothesis For the r-dimensional parameter vector ω = (ν, ψ ), let U = (U ν, U ψ) be the score vector where U ν = l ν, U ψ = l ψ Then under the null hypothesis H 0 : ψ = = ψ q = ψ 0, 5
6 U N r 0, T S S R, where T is a (r q) (r q) matrix with T(i, j) = E, ν i ν j H0 S is a q (r q) matrix with S(i, j) = E and R is a q q matrix with R(i, j) = E The score test statistic is then given by ψ i ν j ( 2 l ψ i ψ j ) H0 H0 ξ = U ψ(r ST S ) U ψ and under the null hypothesis, ξ is asymptotically distributed as χ 2 q Note that all the terms involved in ξ are evaluated at the ML estimates of ω under H 0 In this paper, we compare the performance of the score test with other tests using Monte Carlo simulations 23 Unrestricted Maximum Likelihood Estimation In Cases, 2 and 3, tests in this paper require maximum likelihood estimation under the alternative hypotheses We refer to this as unrestricted maximum likelihood estimation Let A be the p p matrix A = p p p p p( p) p( p) p( p) p( p) p (p ) p( p) 6
7 Define Z i = X i A, i =, 2,, K or equivalently, z ij = Ax ij, i =, 2,, K, j =, 2,, n i Then z ij N p (µ i, Σ i), where µ = i ρ i + p ( ρ i ) µ i 0 ρ i , Σ i = AΣ ia = σ 2 i 0 0 ρ i ρ i Note that z ij N (µ i, σ 2 i {ρ i + p ( ρ i )}) and z ij2,, z ijp iid N (0, σ 2 i { ρ i}), i =, 2,, K ; j =, 2,, n i Also, z ij is independent of z ij2,,z ijp The log-likelihood in terms of µ i, σ i and ρ i is l = l(µ,,µ K, σ 2,,σ 2 K, ρ,,ρ K ) = K 2 p n i ln σ 2 i K n i ln{ρ i= 2 i + p ( ρ i )} K i= 2 (p ) n i ln( ρ i ) i= K n i (z ij µ i ) 2 K 2σ 2 i {ρ i + p ( ρ i )} n i p zijk 2 2σ 2 i( ρ i ) i= j= i= j= k=2 It follows that the ML estimates ˆµ i, ˆσ 2 i and ˆρ i of µ i, σ 2 i and ρ i are given by ˆµ i = n i n i z ij j= and ˆσ 2 i = ni j=(z ij ˆµ i ) 2 ni n i p{ˆρ i + p ( ˆρ i )} + j= k=2 z2 ijk n i p( ˆρ i ) ˆρ i = p n i j=(z ij ˆµ i ) 2 (p ) n i p n i j=(z ij ˆµ i ) 2 + n i j= k=2 z2 ijk j= k=2 z2 ijk (4) 3 CASE : H 0 : ρ = = ρ K The Case problem with grouped normal data has been well studied in the literature We fill in some gaps concerning the likelihood ratio test and derive the Skovgaard (200) test For comparison purposes, we present the score test (Paul and Barnwall 990) 7
8 3 Likelihood Ratio Test Under the hypothesis H 0 : ρ = ρ 2 = = ρ K = ρ (an unspecified value), the log-likelihood is l 0 = l 0 (µ,,µ K, σ 2,,σ2 k, ρ) = 2 p K i= K n i i= j= n i ln σ 2 i 2 n ln{ρ + p ( ρ)} n(p ) ln( ρ) 2 (z ij µ i ) 2 2σ 2 i {ρ + p ( ρ)} K n i p i= j= k=2 z 2 ijk 2σ 2 i( ρ) The maximum likelihood estimates for µ i and σ 2 i are given by and σ 2 i = The score equation for ρ µ i = n i z ij = ˆµ n i i j= U(ρ) = l 0 ρ = n(p ) n(p ) + 2p{ρ + p ( ρ)} 2( ρ) ni j=(z ij µ i ) 2 ni n i p{ ρ + p ( ρ)} + j= k=2 z2 ijk (5) n i p( ρ) 2( ρ) 2 K n i p zijk 2 /σ2 i = 0 i= j= k=2 (p ) Ki= ni + j=(z ij µ i ) 2 /σ 2 i 2p{ρ + p ( ρ)} 2 does not admit an explicit solution Estimates are obtained by the scoring method ρ (m+) = ρ (m) + U(ρ(m) ), m = 0,, (6) I(ρ (m) ) where I(ρ) = E ( 2 l 0 ρ 2 ) is the Fisher information To obtain the Fisher information, we have 2 l 0 ρ = n(p ) 2 n(p ) 2 2p 2 {ρ + p ( ρ)} 2+ 2( ρ) (p )2 K ni i= j=(z ij µ i ) 2 /σ 2 i 2 p 2 {ρ + p ( ρ)} 3 Since E(z ij µ i ) 2 = σ 2 i {ρ + p ( ρ)} and E(zijk 2 ) = σ2 i ( ρ), it follows that I(ρ) = E ( 2 ) l 0 = ρ 2 n(p ) 2 n(p ) + 2p 2 {ρ + p ( ρ)} 2 2( ρ) 2 Ki= ni j= k=2 z2 ijk /σ2 i ( ρ) 3 We summarize the estimation procedure for µ,, µ K, σ 2,,σ2 K algorithm: and ρ in the following 8
9 Let m = 0 and ρ (0) = K i= ˆρ i /K where ˆρ i is given by (4) 2 Let (σ 2 i) (m) be the m th approximation to σ 2 i obtained by substituting ρ (m) for ρ in (5) 3 Calculate ρ (m+) using (6) where (σ 2 i )(m) and ρ (m) replace σ 2 i and ρ respectively 4 Let m = m + 5 If U(ρ (m) ) < 0 6 then stop Else, go to step 2 and continue In all of the simulations that we have considered, we have not experienced any difficulty with the above algorithm Let l = l( ω) = l(ˆµ,, ˆµ K, ˆσ 2,, ˆσ 2 K, ˆρ,,ˆρ K ) and l 0 = l 0 ( ω 0 ) = l 0 ( µ,, µ K, σ 2,, σ 2 K, ρ) be the maximized values of the loglikelihood Then, under the null hypothesis H 0, the likelihood ratio test statistic w = 2( l l 0 ) follows the chi-square distribution with (K ) degrees of freedom for large sample sizes 32 Modified Likelihood Ratio Tests In the following, we derive the Skovgaard (200) modifications to the likelihood ratio test statistic which is intended to provide more accurate inference in small-sample situations For our problem, the canonical parameter θ = (θ,,θ K, θ K+,,θ 2K, θ 2K+,,θ 3K ) and canonical sufficient statistic t = (t,,t K, t K+,,t 2K, t 2K+,,t 3K ) are as follows θ i = λ 2, t i = zij 2 /2, i j= θ K+i = µ n i i λ 2, t K+i = z ij, i j= θ 2K+i = n i p β 2, t 2K+i = zijk 2 /2 i j= k=2 where i =,, K, λ 2 i = σ 2 i {ρ i + p ( ρ i )} and β 2 i = σ 2 i ( ρ i) Then the terms in τ = E(t) are τ i = n i E(z 2 ij )/2 = n i(µ 2 i + λ2 i )/2, n i 9
10 τ K+i = n i E(z ij ) = n i µ i, and the elements in Ω = V ar(t) are τ 2K+i = n i (p )E(z 2 ij2)/2 = n i (p )β 2 i/2 Ω(i, i) = n i V ar(z 2 ij)/4 = n i (λ 2 i + 2µ 2 i)λ 2 i/2, Ω(K + i, K + i) = n i V ar(z ij ) = n i λ 2 i, Ω(2K + i, 2K + i) = n i (p ) V ar(z 2 ij2)/4 = n i (p )β 4 i/2, n i Ω(K + i, i) = Ω(K + i, i) = Cov where i =,,K All the other covariances are zero n i zij, 2 z ij j= j= /2 = n i µ i λ 2 i In order to obtain the Skovgaard (200) term δ in (), we need to evaluate J νν and Ĩνν for the nuisance parameter ν = (µ,,µ K, σ 2,, σ 2 K) We have J (µ i, µ i ) = 2 l n i µ 2 = i ω= ω σ 2 i { ρ + p ( ρ)}, J (µ i, σ 2 i) = 2 l ni j=(z ij µ µ i σ 2 = i ) i ω= ω σ 4 i { ρ + p ( ρ)} = 0, J (σ 2 i, σ 2 i) = 2 l σ 4 = n ni ip j=(z ij µ i ω= ω 2 σ 4 + i ) 2 ni i σ 6 i { ρ + p ( ρ)} + j= k=2 z2 ijk σ 6 = i ( ρ) In addition, J (µ i, µ j ) = J (µ i, σ 2 j ) = J (σ 2 i, σ2 j ) = 0, i j since Further, 2 l µ i µ j = 2 l σ 2 i σ 2 j = 2 l µ i σ 2 j = 0, i j n i Ĩ(µ i, µ i ) = E = µ 2 i ω= ω σ 2 i { ρ + p ( ρ)}, ( Ĩ(µ i, σ 2 2 ) i) = E = 0, µ i σ 2 i ( ω= ω Ĩ(σ 2 2 ) i, σ2 i ) = E = n ip σ 4 i ω= ω 2 σ 4, i Ĩ(µ i, µ j ) = Ĩ(µ i, σ 2 j) = Ĩ(σ2 i, σ 2 j) = 0, i j n ip 2 σ 4 i using equation (5) Note that as J νν = Ĩνν, it is not necessary to evaluate Ĩνν and J νν as they cancel out in δ 0
11 33 Score Test In order to evaluate the score test statistic ( ξ = U ψ R ST S ) Uψ we require the vector of scores U ψ = (U,,U K ) where U i = l ρ i ω= ω and the K K matrix R, the K 2K matrix S and the 2K 2K matrix T where R(i, j) = E, ρ i ρ j ( ω= ω 2 ) S(i, j) = E, ρ i µ j ( ω= ω 2 ) S(i, K + j) = E, ρ i σ 2 j ( ω= ω 2 ) T(i, j) = E, µ i µ j ( ω= ω 2 ) T(i, K + j) = T K+i j = E, µ i σ 2 j ( ω= ω 2 ) T(K + i, K + j) = E σ 2 i σ 2 j ω= ω Explicit expressions for the components of ξ are given by U i n i = 2{ρ + p ( ρ)} n i(p ) 2( ρ) (p ) ni j=(x ij x i ) 2 ni + 2pσ 2 i {ρ + p ( ρ)} j= k=2 z2 ijk, 2 2σ 2 i( ρ) 2 R = diag ( 3ni (p ) 2( ρ) n ) i(3p 2)(p ), 2 2p 2 {ρ + p ( ρ)} 2 K K S = [ ( 0 K K diag ) n i (p )ρ 2σ 2 i ( ρ){ρ+p ( ρ)} K K ], T = ( diag n i σ 2 i {ρ+p ( ρ)} ) K K 0 K K diag 0 K K ( ) n i p 2σ 4 i K K As the matrices R, S and T have a special structure, the inverses required in the score test statistic ξ can be easily computed
12 4 CASE 2: H 02 : σ = = σ K As mentioned previously, there has been considerable amount of work directed at the Case problem The Case 2 problem does not seem to have been studied It is conceivable that group correlations may differ yet group variances may be the same For example, if the same type of instruments are used across groups and the source of variability is due to measurement error, then group variances are comparable However, if additional sources of variability are at play, then H 02 may be false For the three Case 2 statistics that follow, the asymptotic distributions follow a χ 2 K distribution 4 Likelihood Ratio Test Under H 02 : σ = = σ K = σ (an unspecified value), the log-likelihood is l 0 = l 0 (µ,,µ K, σ 2, ρ,,ρ K ) = 2 np ln σ2 K n i ln{ρ 2 i + p ( ρ i )} K 2 (p ) n i ln( ρ i ) n i K σ 2 i= j= i= (z ij µ i ) 2 2{ρ i + p ( ρ i )} σ 2 K n i p i= j= k=2 i= z 2 ijk 2( ρ i ) The maximum likelihood estimates for µ,, µ K and σ 2 are given by and The equations σ 2 = np { K n i i= j= l 0 n i (p ) = ρ i 2p{ρ i + p ( ρ i )} + n i(p ) 2( ρ i ) µ i = n i z ij = ˆµ n i i j= (z ij µ i ) 2 { ρ i + p ( ρ i )} + p k=2 zijk 2 } (7) ( ρ i ) (p ) ni + j=(z ij µ i ) 2 ni 2pσ 2 {ρ i + p ( ρ i )} j= k=2 z2 ijk = 0 2 2σ 2 ( ρ i ) 2 do not lead to an explicit solution for ρ i, i =,, K To obtain the ML estimate ρ i of ρ i under H 02, our first instinct was to use the scoring method However, we experienced some convergence problems Instead, the method of bisection was successfully used on equation (8) Because of the interdependence between σ 2 and ρ i, the following algorithm was used for estimation 2 (8)
13 Let m = 0 and ρ (0) i = ˆρ i where ˆρ i is given by (4), i =,, K 2 Let (σ 2 ) (m) be the m th approximation to σ 2 obtained by substituting ρ (m) i (7) for ρ i in 3 Calculate ρ (m+) i using bisection on (8) where (σ 2 ) (m) replaces σ 2, i =,, K 4 Let m = m + 5 If l 0 ρ i < 0 6 for i =,,K then stop Else, go to step 2 and continue The likelihood ratio test statistic is given by w = 2( l l 0 ) where l is the maximized unrestricted loglikelihood and l 0 is the maximized restricted loglikelihood 42 Modified Likelihood Ratio Tests Everything except the terms Ĩνν and J νν in the δ formula () remains the same as in Case The matrices J νν and Ĩνν for the nuisance parameter ν = (µ,,µ K, ρ,,ρ K ) are as follows For i =,,K, J νν (i, i) = 2 l n i µ 2 = i ω= ω σ 2 { ρ i + p ( ρ i )} ni J νν (i, K + i) = 2 l (p ) = j=(z ij µ i ) µ i ρ i ω= ω p σ 2 { ρ i + p ( ρ i )} = 0 2 J νν (K + i, K + i) = 2 l n i (p ) 2 ρ 2 = i ω= ω 2p 2 { ρ i p ( ρ i )} n i(p ) 2 2( ρ i ) 2 + (p )2 n i j=(z ij µ i ) 2 ni p 2 σ 2 { ρ i + p ( ρ i )} + frac zijk σ 2 2 ( ρ 3 i ) 3 j= k=2 All other terms are zero since 2 l = 2 l = µ i µ j ρ i ρ j 2 l µ i ρ j = 0, i j For i =,,K, n i Ĩ νν (i, i) = E = µ 2 i ω= ω σ 2 { ρ i + p ( ρ i )} Ĩ νν (i, K + i) = E = 0 µ i ρ i ( ω= ω 2 ) n i (p ) 2 Ĩ νν (K + i, K + i) = E = ω= ω 2p 2 { ρ i + p ( ρ i )} + n i(p ) 2 2( ρ i ) 2 ρ 2 i 3
14 and all other terms are zero 43 Score Test The score vector U ψ and the matrices R, S and T are given by the following expressions where i, j =,,K U i = l σ 2 i = n ip ω= ω 2 σ 2 + ni j=(z ij µ i ) 2 ni 2 σ 4 { ρ i + p ( ρ i )} + j= k=2 z2 ijk 2 σ 4, ( ρ i ) R(i, i) = E = n ip ω= ω 2 σ 4, σ 4 i R(i, j) = E = 0, i j, σ 2 i σ 2 j ω= ω S(i, i) = E = 0, σ 2 i µ i ω= ω S(i, K + i) = E ( 2 l σ 2 i ρ i ) ω= ω = (p )n i 2p σ 2 { ρ i + p ( ρ i )} (p )n i 2 σ 2 ( ρ i ), S(i, j) = E = 0, i j, σ 2 i µ j ω= ω S(i, K + j) = E = 0, i j, σ 2 i ρ j ω= ω T(i, j) = E = ν i ν j Ĩνν(i, j) ω= ω where Ĩνν is given in Subsection 42 The score statistic for testing H 02 : σ = = σ K is then given by ξ = U ψ (R ST S ) U ψ 4
15 5 CASE 3: H 03 : σ = = σ K, ρ = = ρ K Since H 0 (Case) and H 02 (Case2) are common hypotheses of interest with respect to the grouped normal data problem, it follows that the composite hypothesis H 03 is also of interest Although there have been a number of test procedures proposed for H 03 (Han 975), we make use of the methodologies developed for Cases and 2 to derive the Skovgaard (200) modifications The likelihood ratio test and the score test are presented for comparison purposes For the three Case 3 statistics that follow, the asymptotic distributions follow a χ 2 2K 2 distribution 5 Likelihood Ratio Test Under H 03 : σ = = σ K = σ and ρ = = ρ K = ρ, the log-likelihood is l 0 = l 0 (µ,,µ K, σ 2, ρ) = 2 np ln σ2 2 n ln{ρ + p ( ρ)} n(p ) ln( ρ) 2 K n i (z 2σ 2 {ρ + p ij µ ( ρ)} i ) 2 K n i 2σ 2 ( ρ) The ML estimates are given by µ i = n i z ij = ˆµ n i, i σ 2 = np j= K n i i= j= i= j= { (zij µ i ) 2 ρ + p ( ρ) + p k=2 zijk 2 }, ( ρ) ρ = p K ni i= j=(z ij µ i ) 2 (p ) K ni i= p K ni i= j=(z ij µ i ) 2 + K ni i= j= k=2 z2 ijk 52 Modified Likelihood Ratio Tests p zijk 2 i= j= k=2 j= k=2 z2 ijk For the Skovgaard (200) modifications, the matrices J νν and Ĩνν for the nuisance parameter ν = (µ,,µ K ) are as follows J νν (i, i) = 2 l n µ 2 = i ω= ω σ 2 { ρ + p ( ρ)}, J νν (i, j) = 2 l = 0, i j, µ i µ j ω= ω 5
16 Ĩ νν (i, i) = E = ω= ω µ 2 i Ĩ νν (i, j) = E ( 2 l µ i µ j ) n σ 2 { ρ + p ( ρ)}, ω= ω = 0, i j Note that Ĩνν = J νν providing a simplification in the δ formula () 53 Score Test The score vector U ψ = (U,,U 2K ) and the matrices R and S are given by U i = l σ 2 i = n ip ω= ω 2 σ 2 + ni j=(z ij µ i ) 2 ni 2 σ 4 { ρ + p ( ρ)} + j= k=2 z2 ijk 2 σ 4, ( ρ) U K+i = l n i (p ) = ρ i ω= ω 2p{ ρ + p ( ρ)} +n i(p ) ) ni +(p j=(z ij µ i ) 2 ni 2( ρ) 2p σ 2 { ρ + p ( ρ)} j= k=2 z2 ijk 2 2 σ 2, ( ρ) 2 R(i, i) = E = n ip σ 4 i ω= ω 2 σ 4, (p )n i R(i, K + i) = R(K + i, i) = E = σ 2 i ρ i ω= ω 2p σ 2 { ρ + p ( ρ)} (p )n i 2 σ 2 ( ρ), n(p ) 2 n(p ) R(K + i, K + i) = E = + ω= ω 2p 2 { ρ + p ( ρ)} 2 2( ρ) 2 ρ 2 i and all other terms are zero since ( 2 ) ( 2 ) E = E = E = 0, i j σ 2 i σ 2 j ω= ω ρ i ρ j ω= ω σ 2 i ρ j ω= ω Further, since S = E = 0 ψ i ν j ω= ω ( 2 ) E = E = 0, σ 2 i µ i ω= ω ρ i µ i ( ω= ω 2 ) ( 2 ) E = E = 0, i j σ 2 i µ j ω= ω ρ i µ j ω= ω Then, the score statistic for testing H 03 reduces to ξ = U ψ (R ST S ) U ψ = U ψ R U ψ 6
17 6 SIMULATION RESULTS AND DISCUSSION We present the results of some Monte Carlo simulations to investigate the adequacy of the test statistics discussed in this paper The following tables show the Type I error probabilities of the likelihood ratio statistic (w), the modified likelihood ratio statistics (w and w ) and the score statistic (ξ) based on N =00,000 simulations In each case, we used the nominal value 005 and note that the results are comparable for the nominal value 00 In addition to investigating the Type I error rate, we are also interested in the adequacy of the χ 2 approximations over the entire range of the test statistics To do so, we calculate the Cramér-von Mises statistic (Anderson and Darling 954) W 2 = /(2N) + {Z (i) (2i )/(2N)} 2 where Z (i) is the i th order statistic corresponding to Z = F(Y ), F is the cumulative distribution function (cdf) of the relevant χ 2 distribution and Y is the test statistic obtained when generating from the equicorrelation model in question The Cramér-von Mises statistic is one in the wide class of discrepancy measures given by the Cramér-von Mises family Q = n (F n (y) F(y)) 2 Ψ(y)dF(y) where F and F n = F n (Y,,Y n ) are the cdf and the empirical distribution function (edf) respectively, and Ψ is a suitable function which gives weights to the squared difference (F n (y) F(y)) 2 Larger values of the statistic Q (and in our case W 2 ), provide greater evidence of the lack of fit between the data Y,,Y n and the proposed F Theory and practical issues associated with the Cramér-von Mises family of statistics are described in detail in Stephens (986) The simulation results are presented in Tables, 2 and 3 for K = 5 groups where n = = n K Before discussing the particular simulation results, we make some general comments It appears that when the sample size n i reaches 20, all four statistics are adequate as the Type I error rates are close to the nominal values Therefore, our discussion is focused on small sample problems and we choose n = = n K = 5 Also, except in a few instances, the likelihood ratio statistic w is inferior to the w, w and ξ Finally, there seems to be little difference in the results when the number of variables changes from p = 5 to p = 0 7
18 In Table, we present the simulation results with respect to the Case scenario The w and ξ statistics perform best with respect to the Type I error rate However, our recommendation is the modified likelihood ratio statistic w based on its superior goodness of fit In Table 2, we present the simulation results with respect to the Case 2 scenario We observe that the performance of all of the test statistics depends on the correlation ρ i Notably, w, w and ξ perform worse at the extreme values of ρ i Again, our recommendation is w although its performance is not as good as in Case In Table 3, we present the simulation results with respect to the Case 3 scenario In Case 3, the score test ξ has a slight edge over the other tests However, we note that the score test is anti-conservative whereas the two modified likelihood ratio tests are conservative ACKNOWLEDGMENTS Gill and Swartz were partially supported by grants from the Natural Sciences and Engineering Research Council of Canada The work was initiated during Banneheka s sabbatical visit to Simon Fraser University in 2003 Authors thank a referee for valuable comments References [] Ahmed, S E; Gupta, A K; Khan, S M; Nicol, C J Simultaneous estimation of several intraclass correlation coefficients Annals of the Institute of Statistical Mathematics 200, 53, [2] Anderson, T W; Darling, D A Asymptotic theory of certain goodness-of-fit criteria based on stochastic processes Annals of Mathematical Statistics, 954, 23, [3] Barndorff-Nielsen, O E Modified signed loglikelihood ratio Biometrika, 99, 78, [4] Barndorff-Nielsen, O E; Cox, D R Inference and Asymptotics Chapman and Hall, London, 994 8
19 [5] Donner, A; Koval, J J The estimation of intraclass correlation in the analysis of family data Biometrics, 980, 36, 9-25 [6] Han, C -P Testing the equality of covariance matrices under intraclass correlation models Annals of the Institute of Statistical Mathematics, 975, 27, [7] Konishi, S; Gupta, A K Testing the equality of several intraclass correlation coefficients Journal of Statistical Planning and Inference, 989, 2, [8] Paul, S R; Barnwal, R K Maximum likelihood estimation and a C(α) test for a common intraclass correlation The Statistician, 990, 39, 9-24 [9] Reid, N Asymptotics and the theory of inference Annals of Statistics, 2003, 3, [0] Skovgaard, I M Likelihood asymptotics Scandinavian Journal of Statistics, 200, 28, 3-32 [] Stephens, M A Tests based on EDF statistics Goodness of Fit Techniques eds R B D Agostino and M A Stephens, Eds; Marcel Dekker, New York, 986 [2] Viana, M A G; Olkin, I; McMahon, T T Multivariate assessment of computeranalyzed corneal topographers Journal of the Optical Society of America Series A, 993, 0,
20 Table Case (H 0 : ρ = = ρ K where σ 2 = = σ2 K = 0 and K = 5) Type I error rates Cramer-von Mises statistic p = 5 ρ i w w w ξ w w w ξ (n i = 5) (n i = 20) p = 0 (n i = 5) (n i = 20)
21 Table 2 Case 2 (H 02 : σ = = σ K where K = 5) Type I error rates Cramer-von Mises statistic p = 5 ρ i w w w ξ w w w ξ (n i = 5) (n i = 20) p = 0 (n i = 5) (n i = 20)
22 Table 3 Case 3 (H 03 : σ = = σ K, ρ = = ρ K where K = 5) Type I error rates Cramer-von Mises statistic p = 5 ρ i w w w ξ w w w ξ (n i = 5) (n i = 20) p = 0 (n i = 5) (n i = 20)
Testing Some Covariance Structures under a Growth Curve Model in High Dimension
Department of Mathematics Testing Some Covariance Structures under a Growth Curve Model in High Dimension Muni S. Srivastava and Martin Singull LiTH-MAT-R--2015/03--SE Department of Mathematics Linköping
More informationResearch Article Inference for the Sharpe Ratio Using a Likelihood-Based Approach
Journal of Probability and Statistics Volume 202 Article ID 87856 24 pages doi:0.55/202/87856 Research Article Inference for the Sharpe Ratio Using a Likelihood-Based Approach Ying Liu Marie Rekkas 2 and
More informationASSESSING A VECTOR PARAMETER
SUMMARY ASSESSING A VECTOR PARAMETER By D.A.S. Fraser and N. Reid Department of Statistics, University of Toronto St. George Street, Toronto, Canada M5S 3G3 dfraser@utstat.toronto.edu Some key words. Ancillary;
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationApproximate Inference for the Multinomial Logit Model
Approximate Inference for the Multinomial Logit Model M.Rekkas Abstract Higher order asymptotic theory is used to derive p-values that achieve superior accuracy compared to the p-values obtained from traditional
More informationCanonical Correlation Analysis of Longitudinal Data
Biometrics Section JSM 2008 Canonical Correlation Analysis of Longitudinal Data Jayesh Srivastava Dayanand N Naik Abstract Studying the relationship between two sets of variables is an important multivariate
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 informationChapter 1 Likelihood-Based Inference and Finite-Sample Corrections: A Brief Overview
Chapter 1 Likelihood-Based Inference and Finite-Sample Corrections: A Brief Overview Abstract This chapter introduces the likelihood function and estimation by maximum likelihood. Some important properties
More informationNotes on the Multivariate Normal and Related Topics
Version: July 10, 2013 Notes on the Multivariate Normal and Related Topics Let me refresh your memory about the distinctions between population and sample; parameters and statistics; population distributions
More informationLikelihood and p-value functions in the composite likelihood context
Likelihood and p-value functions in the composite likelihood context D.A.S. Fraser and N. Reid Department of Statistical Sciences University of Toronto November 19, 2016 Abstract The need for combining
More informationStatement: With my signature I confirm that the solutions are the product of my own work. Name: Signature:.
MATHEMATICAL STATISTICS Homework assignment Instructions Please turn in the homework with this cover page. You do not need to edit the solutions. Just make sure the handwriting is legible. You may discuss
More informationConditional Inference by Estimation of a Marginal Distribution
Conditional Inference by Estimation of a Marginal Distribution Thomas J. DiCiccio and G. Alastair Young 1 Introduction Conditional inference has been, since the seminal work of Fisher (1934), a fundamental
More informationTesting Goodness Of Fit Of The Geometric Distribution: An Application To Human Fecundability Data
Journal of Modern Applied Statistical Methods Volume 4 Issue Article 8 --5 Testing Goodness Of Fit Of The Geometric Distribution: An Application To Human Fecundability Data Sudhir R. Paul University of
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 informationA comparison study of the nonparametric tests based on the empirical distributions
통계연구 (2015), 제 20 권제 3 호, 1-12 A comparison study of the nonparametric tests based on the empirical distributions Hyo-Il Park 1) Abstract In this study, we propose a nonparametric test based on the empirical
More informationDepartment of Statistics
Research Report Department of Statistics Research Report Department of Statistics No. 05: Testing in multivariate normal models with block circular covariance structures Yuli Liang Dietrich von Rosen Tatjana
More information3. (a) (8 points) There is more than one way to correctly express the null hypothesis in matrix form. One way to state the null hypothesis is
Stat 501 Solutions and Comments on Exam 1 Spring 005-4 0-4 1. (a) (5 points) Y ~ N, -1-4 34 (b) (5 points) X (X,X ) = (5,8) ~ N ( 11.5, 0.9375 ) 3 1 (c) (10 points, for each part) (i), (ii), and (v) are
More informationBootstrap and Parametric Inference: Successes and Challenges
Bootstrap and Parametric Inference: Successes and Challenges G. Alastair Young Department of Mathematics Imperial College London Newton Institute, January 2008 Overview Overview Review key aspects of frequentist
More informationRecall the Basics of Hypothesis Testing
Recall the Basics of Hypothesis Testing The level of significance α, (size of test) is defined as the probability of X falling in w (rejecting H 0 ) when H 0 is true: P(X w H 0 ) = α. H 0 TRUE H 1 TRUE
More informationMATH5745 Multivariate Methods Lecture 07
MATH5745 Multivariate Methods Lecture 07 Tests of hypothesis on covariance matrix March 16, 2018 MATH5745 Multivariate Methods Lecture 07 March 16, 2018 1 / 39 Test on covariance matrices: Introduction
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 informationTesting for a unit root in an ar(1) model using three and four moment approximations: symmetric distributions
Hong Kong Baptist University HKBU Institutional Repository Department of Economics Journal Articles Department of Economics 1998 Testing for a unit root in an ar(1) model using three and four moment approximations:
More informationThe formal relationship between analytic and bootstrap approaches to parametric inference
The formal relationship between analytic and bootstrap approaches to parametric inference T.J. DiCiccio Cornell University, Ithaca, NY 14853, U.S.A. T.A. Kuffner Washington University in St. Louis, St.
More informationAdaptive Designs: Why, How and When?
Adaptive Designs: Why, How and When? Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj ISBS Conference Shanghai, July 2008 1 Adaptive designs:
More informationUNIVERSITÄT POTSDAM Institut für Mathematik
UNIVERSITÄT POTSDAM Institut für Mathematik Testing the Acceleration Function in Life Time Models Hannelore Liero Matthias Liero Mathematische Statistik und Wahrscheinlichkeitstheorie Universität Potsdam
More informationHypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations
Hypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations Takeshi Emura and Hisayuki Tsukuma Abstract For testing the regression parameter in multivariate
More informationAnalysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems
Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Jeremy S. Conner and Dale E. Seborg Department of Chemical Engineering University of California, Santa Barbara, CA
More informationAccurate directional inference for vector parameters
Accurate directional inference for vector parameters Nancy Reid February 26, 2016 with Don Fraser, Nicola Sartori, Anthony Davison Nancy Reid Accurate directional inference for vector parameters York University
More informationTesting equality of two mean vectors with unequal sample sizes for populations with correlation
Testing equality of two mean vectors with unequal sample sizes for populations with correlation Aya Shinozaki Naoya Okamoto 2 and Takashi Seo Department of Mathematical Information Science Tokyo University
More informationSample size determination for logistic regression: A simulation study
Sample size determination for logistic regression: A simulation study Stephen Bush School of Mathematical Sciences, University of Technology Sydney, PO Box 123 Broadway NSW 2007, Australia Abstract This
More informationNonparametric Tests for Multi-parameter M-estimators
Nonparametric Tests for Multi-parameter M-estimators John Robinson School of Mathematics and Statistics University of Sydney The talk is based on joint work with John Kolassa. It follows from work over
More informationSample size calculations for logistic and Poisson regression models
Biometrika (2), 88, 4, pp. 93 99 2 Biometrika Trust Printed in Great Britain Sample size calculations for logistic and Poisson regression models BY GWOWEN SHIEH Department of Management Science, National
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 informationAccurate directional inference for vector parameters
Accurate directional inference for vector parameters Nancy Reid October 28, 2016 with Don Fraser, Nicola Sartori, Anthony Davison Parametric models and likelihood model f (y; θ), θ R p data y = (y 1,...,
More informationAn Approximate Test for Homogeneity of Correlated Correlation Coefficients
Quality & Quantity 37: 99 110, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. 99 Research Note An Approximate Test for Homogeneity of Correlated Correlation Coefficients TRIVELLORE
More informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
More informationCONVERTING OBSERVED LIKELIHOOD FUNCTIONS TO TAIL PROBABILITIES. D.A.S. Fraser Mathematics Department York University North York, Ontario M3J 1P3
CONVERTING OBSERVED LIKELIHOOD FUNCTIONS TO TAIL PROBABILITIES D.A.S. Fraser Mathematics Department York University North York, Ontario M3J 1P3 N. Reid Department of Statistics University of Toronto Toronto,
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationA TEST OF FIT FOR THE GENERALIZED PARETO DISTRIBUTION BASED ON TRANSFORMS
A TEST OF FIT FOR THE GENERALIZED PARETO DISTRIBUTION BASED ON TRANSFORMS Dimitrios Konstantinides, Simos G. Meintanis Department of Statistics and Acturial Science, University of the Aegean, Karlovassi,
More informationTesting Equality of Two Intercepts for the Parallel Regression Model with Non-sample Prior Information
Testing Equality of Two Intercepts for the Parallel Regression Model with Non-sample Prior Information Budi Pratikno 1 and Shahjahan Khan 2 1 Department of Mathematics and Natural Science Jenderal Soedirman
More informationDEFNITIVE TESTING OF AN INTEREST PARAMETER: USING PARAMETER CONTINUITY
Journal of Statistical Research 200x, Vol. xx, No. xx, pp. xx-xx ISSN 0256-422 X DEFNITIVE TESTING OF AN INTEREST PARAMETER: USING PARAMETER CONTINUITY D. A. S. FRASER Department of Statistical Sciences,
More informationA copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2
based on the conditional probability integral transform Daniel Berg 1 Henrik Bakken 2 1 Norwegian Computing Center (NR) & University of Oslo (UiO) 2 Norwegian University of Science and Technology (NTNU)
More informationObtaining Critical Values for Test of Markov Regime Switching
University of California, Santa Barbara From the SelectedWorks of Douglas G. Steigerwald November 1, 01 Obtaining Critical Values for Test of Markov Regime Switching Douglas G Steigerwald, University of
More informationStandardizing the Empirical Distribution Function Yields the Chi-Square Statistic
Standardizing the Empirical Distribution Function Yields the Chi-Square Statistic Andrew Barron Department of Statistics, Yale University Mirta Benšić Department of Mathematics, University of Osijek and
More informationAn Improved Specification Test for AR(1) versus MA(1) Disturbances in Linear Regression Models
An Improved Specification Test for AR(1) versus MA(1) Disturbances in Linear Regression Models Pierre Nguimkeu Georgia State University Abstract This paper proposes an improved likelihood-based method
More informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationThird-order inference for autocorrelation in nonlinear regression models
Third-order inference for autocorrelation in nonlinear regression models P. E. Nguimkeu M. Rekkas Abstract We propose third-order likelihood-based methods to derive highly accurate p-value approximations
More informationTESTING FOR NORMALITY IN THE LINEAR REGRESSION MODEL: AN EMPIRICAL LIKELIHOOD RATIO TEST
Econometrics Working Paper EWP0402 ISSN 1485-6441 Department of Economics TESTING FOR NORMALITY IN THE LINEAR REGRESSION MODEL: AN EMPIRICAL LIKELIHOOD RATIO TEST Lauren Bin Dong & David E. A. Giles Department
More information2.6.3 Generalized likelihood ratio tests
26 HYPOTHESIS TESTING 113 263 Generalized likelihood ratio tests When a UMP test does not exist, we usually use a generalized likelihood ratio test to verify H 0 : θ Θ against H 1 : θ Θ\Θ It can be used
More informationLikelihood Inference in the Presence of Nuisance Parameters
PHYSTAT2003, SLAC, September 8-11, 2003 1 Likelihood Inference in the Presence of Nuance Parameters N. Reid, D.A.S. Fraser Department of Stattics, University of Toronto, Toronto Canada M5S 3G3 We describe
More informationPRINCIPLES OF STATISTICAL INFERENCE
Advanced Series on Statistical Science & Applied Probability PRINCIPLES OF STATISTICAL INFERENCE from a Neo-Fisherian Perspective Luigi Pace Department of Statistics University ofudine, Italy Alessandra
More informationLikelihood Inference in the Presence of Nuisance Parameters
Likelihood Inference in the Presence of Nuance Parameters N Reid, DAS Fraser Department of Stattics, University of Toronto, Toronto Canada M5S 3G3 We describe some recent approaches to likelihood based
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationAdjusted Empirical Likelihood for Long-memory Time Series Models
Adjusted Empirical Likelihood for Long-memory Time Series Models arxiv:1604.06170v1 [stat.me] 21 Apr 2016 Ramadha D. Piyadi Gamage, Wei Ning and Arjun K. Gupta Department of Mathematics and Statistics
More informationExact goodness-of-fit tests for censored data
Exact goodness-of-fit tests for censored data Aurea Grané Statistics Department. Universidad Carlos III de Madrid. Abstract The statistic introduced in Fortiana and Grané (23, Journal of the Royal Statistical
More informationCentral Limit Theorem ( 5.3)
Central Limit Theorem ( 5.3) Let X 1, X 2,... be a sequence of independent random variables, each having n mean µ and variance σ 2. Then the distribution of the partial sum S n = X i i=1 becomes approximately
More informationDSGE Methods. Estimation of DSGE models: GMM and Indirect Inference. Willi Mutschler, M.Sc.
DSGE Methods Estimation of DSGE models: GMM and Indirect Inference Willi Mutschler, M.Sc. Institute of Econometrics and Economic Statistics University of Münster willi.mutschler@wiwi.uni-muenster.de Summer
More informationA COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky
A COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky Empirical likelihood with right censored data were studied by Thomas and Grunkmier (1975), Li (1995),
More informationAsymptotic Statistics-VI. Changliang Zou
Asymptotic Statistics-VI Changliang Zou Kolmogorov-Smirnov distance Example (Kolmogorov-Smirnov confidence intervals) We know given α (0, 1), there is a well-defined d = d α,n such that, for any continuous
More informationHypothesis testing in multilevel models with block circular covariance structures
1/ 25 Hypothesis testing in multilevel models with block circular covariance structures Yuli Liang 1, Dietrich von Rosen 2,3 and Tatjana von Rosen 1 1 Department of Statistics, Stockholm University 2 Department
More informationMultivariate Regression
Multivariate Regression The so-called supervised learning problem is the following: we want to approximate the random variable Y with an appropriate function of the random variables X 1,..., X p with the
More informationMultivariate Statistics
Multivariate Statistics Chapter 2: Multivariate distributions and inference Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2016/2017 Master in Mathematical
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationModified Simes Critical Values Under Positive Dependence
Modified Simes Critical Values Under Positive Dependence Gengqian Cai, Sanat K. Sarkar Clinical Pharmacology Statistics & Programming, BDS, GlaxoSmithKline Statistics Department, Temple University, Philadelphia
More informationCORRELATION ANALYSIS OF ORDERED OBSERVATIONS FROM A BLOCK-EQUICORRELATED MULTIVARIATE NORMAL DISTRIBUTION
CORRELATION ANALYSIS OF ORDERED OBSERVATIONS FROM A BLOCK-EQUICORRELATED MULTIVARIATE NORMAL DISTRIBUTION MARLOS VIANA AND INGRAM OLKIN Abstract. Inferential procedures are developed for correlations between
More informationApplication of Parametric Homogeneity of Variances Tests under Violation of Classical Assumption
Application of Parametric Homogeneity of Variances Tests under Violation of Classical Assumption Alisa A. Gorbunova and Boris Yu. Lemeshko Novosibirsk State Technical University Department of Applied Mathematics,
More informationLecture 3. Inference about multivariate normal distribution
Lecture 3. Inference about multivariate normal distribution 3.1 Point and Interval Estimation Let X 1,..., X n be i.i.d. N p (µ, Σ). We are interested in evaluation of the maximum likelihood estimates
More informationStatistics and econometrics
1 / 36 Slides for the course Statistics and econometrics Part 10: Asymptotic hypothesis testing European University Institute Andrea Ichino September 8, 2014 2 / 36 Outline Why do we need large sample
More informationStatistical Inference with Monotone Incomplete Multivariate Normal Data
Statistical Inference with Monotone Incomplete Multivariate Normal Data p. 1/4 Statistical Inference with Monotone Incomplete Multivariate Normal Data This talk is based on joint work with my wonderful
More informationStatistical Inference: Estimation and Confidence Intervals Hypothesis Testing
Statistical Inference: Estimation and Confidence Intervals Hypothesis Testing 1 In most statistics problems, we assume that the data have been generated from some unknown probability distribution. We desire
More informationAN EMPIRICAL LIKELIHOOD RATIO TEST FOR NORMALITY
Econometrics Working Paper EWP0401 ISSN 1485-6441 Department of Economics AN EMPIRICAL LIKELIHOOD RATIO TEST FOR NORMALITY Lauren Bin Dong & David E. A. Giles Department of Economics, University of Victoria
More informationA Multiple Testing Approach to the Regularisation of Large Sample Correlation Matrices
A Multiple Testing Approach to the Regularisation of Large Sample Correlation Matrices Natalia Bailey 1 M. Hashem Pesaran 2 L. Vanessa Smith 3 1 Department of Econometrics & Business Statistics, Monash
More informationTheory and Methods of Statistical Inference. PART I Frequentist theory and methods
PhD School in Statistics cycle XXVI, 2011 Theory and Methods of Statistical Inference PART I Frequentist theory and methods (A. Salvan, N. Sartori, L. Pace) Syllabus Some prerequisites: Empirical distribution
More informationFlat and multimodal likelihoods and model lack of fit in curved exponential families
Mathematical Statistics Stockholm University Flat and multimodal likelihoods and model lack of fit in curved exponential families Rolf Sundberg Research Report 2009:1 ISSN 1650-0377 Postal address: Mathematical
More informationBARTLETT IDENTITIES AND LARGE DEVIATIONS IN LIKELIHOOD THEORY 1. By Per Aslak Mykland University of Chicago
The Annals of Statistics 1999, Vol. 27, No. 3, 1105 1117 BARTLETT IDENTITIES AND LARGE DEVIATIONS IN LIKELIHOOD THEORY 1 By Per Aslak Mykland University of Chicago The connection between large and small
More informationExact goodness-of-fit tests for censored data
Ann Inst Stat Math ) 64:87 3 DOI.7/s463--356-y Exact goodness-of-fit tests for censored data Aurea Grané Received: February / Revised: 5 November / Published online: 7 April The Institute of Statistical
More informationOn testing the equality of mean vectors in high dimension
ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 17, Number 1, June 2013 Available online at www.math.ut.ee/acta/ On testing the equality of mean vectors in high dimension Muni S.
More informationIntegrated likelihoods in models with stratum nuisance parameters
Riccardo De Bin, Nicola Sartori, Thomas A. Severini Integrated likelihoods in models with stratum nuisance parameters Technical Report Number 157, 2014 Department of Statistics University of Munich http://www.stat.uni-muenchen.de
More informationImproved Inference for First Order Autocorrelation using Likelihood Analysis
Improved Inference for First Order Autocorrelation using Likelihood Analysis M. Rekkas Y. Sun A. Wong Abstract Testing for first-order autocorrelation in small samples using the standard asymptotic test
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 informationTheory and Methods of Statistical Inference
PhD School in Statistics cycle XXIX, 2014 Theory and Methods of Statistical Inference Instructors: B. Liseo, L. Pace, A. Salvan (course coordinator), N. Sartori, A. Tancredi, L. Ventura Syllabus Some prerequisites:
More informationAnalysis of variance, multivariate (MANOVA)
Analysis of variance, multivariate (MANOVA) Abstract: A designed experiment is set up in which the system studied is under the control of an investigator. The individuals, the treatments, the variables
More informationModeling the scale parameter ϕ A note on modeling correlation of binary responses Using marginal odds ratios to model association for binary responses
Outline Marginal model Examples of marginal model GEE1 Augmented GEE GEE1.5 GEE2 Modeling the scale parameter ϕ A note on modeling correlation of binary responses Using marginal odds ratios to model association
More informationMaster s Written Examination
Master s Written Examination Option: Statistics and Probability Spring 05 Full points may be obtained for correct answers to eight questions Each numbered question (which may have several parts) is worth
More informationHypothesis testing:power, test statistic CMS:
Hypothesis testing:power, test statistic The more sensitive the test, the better it can discriminate between the null and the alternative hypothesis, quantitatively, maximal power In order to achieve this
More informationMISCELLANEOUS TOPICS RELATED TO LIKELIHOOD. Copyright c 2012 (Iowa State University) Statistics / 30
MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD Copyright c 2012 (Iowa State University) Statistics 511 1 / 30 INFORMATION CRITERIA Akaike s Information criterion is given by AIC = 2l(ˆθ) + 2k, where l(ˆθ)
More information,..., θ(2),..., θ(n)
Likelihoods for Multivariate Binary Data Log-Linear Model We have 2 n 1 distinct probabilities, but we wish to consider formulations that allow more parsimonious descriptions as a function of covariates.
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 informationGeneralized Quasi-likelihood versus Hierarchical Likelihood Inferences in Generalized Linear Mixed Models for Count Data
Sankhyā : The Indian Journal of Statistics 2009, Volume 71-B, Part 1, pp. 55-78 c 2009, Indian Statistical Institute Generalized Quasi-likelihood versus Hierarchical Likelihood Inferences in Generalized
More informationGARCH Models Estimation and Inference
Università di Pavia GARCH Models Estimation and Inference Eduardo Rossi Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationIrr. Statistical Methods in Experimental Physics. 2nd Edition. Frederick James. World Scientific. CERN, Switzerland
Frederick James CERN, Switzerland Statistical Methods in Experimental Physics 2nd Edition r i Irr 1- r ri Ibn World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI CHENNAI CONTENTS
More informationEmpirical Likelihood Tests for High-dimensional Data
Empirical Likelihood Tests for High-dimensional Data Department of Statistics and Actuarial Science University of Waterloo, Canada ICSA - Canada Chapter 2013 Symposium Toronto, August 2-3, 2013 Based on
More informationTheory and Methods of Statistical Inference. PART I Frequentist likelihood methods
PhD School in Statistics XXV cycle, 2010 Theory and Methods of Statistical Inference PART I Frequentist likelihood methods (A. Salvan, N. Sartori, L. Pace) Syllabus Some prerequisites: Empirical distribution
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
More informationGoodness-of-fit tests for randomly censored Weibull distributions with estimated parameters
Communications for Statistical Applications and Methods 2017, Vol. 24, No. 5, 519 531 https://doi.org/10.5351/csam.2017.24.5.519 Print ISSN 2287-7843 / Online ISSN 2383-4757 Goodness-of-fit tests for randomly
More informationModified Kolmogorov-Smirnov Test of Goodness of Fit. Catalonia-BarcelonaTECH, Spain
152/304 CoDaWork 2017 Abbadia San Salvatore (IT) Modified Kolmogorov-Smirnov Test of Goodness of Fit G.S. Monti 1, G. Mateu-Figueras 2, M. I. Ortego 3, V. Pawlowsky-Glahn 2 and J. J. Egozcue 3 1 Department
More informationApplication of Variance Homogeneity Tests Under Violation of Normality Assumption
Application of Variance Homogeneity Tests Under Violation of Normality Assumption Alisa A. Gorbunova, Boris Yu. Lemeshko Novosibirsk State Technical University Novosibirsk, Russia e-mail: gorbunova.alisa@gmail.com
More informationTest of fit for a Laplace distribution against heavier tailed alternatives
DEPARTMENT OF STATISTICS AND ACTUARIAL SCIENCE University of Waterloo, 200 University Avenue West Waterloo, Ontario, Canada, N2L 3G1 519-888-4567, ext. 00000 Fax: 519-746-1875 www.stats.uwaterloo.ca UNIVERSITY
More informationModification and Improvement of Empirical Likelihood for Missing Response Problem
UW Biostatistics Working Paper Series 12-30-2010 Modification and Improvement of Empirical Likelihood for Missing Response Problem Kwun Chuen Gary Chan University of Washington - Seattle Campus, kcgchan@u.washington.edu
More informationMinimum Hellinger Distance Estimation in a. Semiparametric Mixture Model
Minimum Hellinger Distance Estimation in a Semiparametric Mixture Model Sijia Xiang 1, Weixin Yao 1, and Jingjing Wu 2 1 Department of Statistics, Kansas State University, Manhattan, Kansas, USA 66506-0802.
More information