Estimation of the large covariance matrix with two-step monotone missing data
|
|
- Florence Allen
- 5 years ago
- Views:
Transcription
1 Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo University of Science 2 Deartment of Health Sciences, Oita University of Nursing and Health Sciences 3 Deartment of Mathematics, KTH Royal Institute of Technology We suggest a shrinkage based technique for estimating covariance matrix in the high-dimensional normal model with missing data. Our asymtotic framework allows the dimensionality grow to infinity together with the samle size, N, and extends the methodology of Ledoit and Wolf 2004 to the case of 2-ste monotone missing data. Two new shrinkage tye estimators are derived and their dominance roerties over the Ledoit and Wolf 2004 estimator are shown under the exected quadratic loss. We erform a simulation study and conclude that the roosed estimators are successful for a range of missing data scenarios. Key Words Monotone missing data; High-dimensional estimation; Shrinkage tye estimation. Mathematics Subject Classification 62H2; 62F2.
2 Introduction Estimation of the covariance matrix Σ in high dimensions has been imortant for many statistical rocedures such as rincial comonent analysis, linear regression and discriminant analysis. Some of these rocedures require estimation of the covariance structure and eigenvalues of Σ while the others need the inverse. It has long been known that standard estimators such as maximum likelihood or least squares develoed for large samle sizes N and relatively small number of variables do not rovide a satisfactory solution in high dimensions. Therefore modification of these estimators is of crucial imortance for high-dimensional roblems. During recent years, many methods for imroving the estimators of Σ and its inverse in the context of high dimensional data, i.e. if the number of variables is much larger than N, have been develoed; see e.g., Kubokawa and Srivastava 2008 and Ledoit and Wolf However, all of these methods have been focused on the case where all data is observed. Since in ractice data sets often contain missing atterns, it is imortant to develo aroaches for high dimensional covariance estimation in resence of missing data. In this study, we focus on a articular tye of statistical missing data, a monotone missing scheme with missing values atterns occurring comletely at random. In other words, we assume that the underlying missing-data mechanism is ignorable under arametric estimation. This tye of missing data roblem has been extensively studied in the statistical literature; see e.g. Anderson 957, Bhargava 975. The samle scheme treated in them is called a k-ste monotone missing data. In articular, in Anderson and Olkin 985, the MLE of Σ for the case of 2-ste monotone missing data was derived exlicitly. Further, Kanda and Fujikoshi 998 studied asymtotic roerties of MLE of Σ constructed by 2-ste monotone missing data and extend these results to the general case of k-ste monotone missing data. However, in high dimensional situations, the MLE of Σ is known to erform oorly even for the case of comlete data, since it is usually imossible to collect 2
3 enough observations to fulfill n condition. In this aer, we resent a solution to this roblem by extending the idea of Ledoit and Wolf 2004 to the case of 2-ste monotone missing data and suggest a new estimator of the covariance matrix that is adjusted to high dimensionality in combination with missing data. The rest of the aer is organized as follows. In Section 2, we set u 2-ste monotone missing data scheme and derive some auxiliary unbiased estimators, which then will be used in Section 3 where we resent our new shrinkage estimators. Section 4 reorts numerical studies on the risk erformance of the suggested estimators and comarison with results by Ledoit and Wolf We conclude in Section 5 and resent the roof of key results in Aendix A.2, A.3 and A.4. 2 Auxiliary unbiased estimators and their asymtotic roerty under the 2-ste monotone missing data To secify the roblem considered here, let x be the -dimensional random vector distributed as N µ, Σ, where µ is mean vector and Σ is covariance matrix. Assume that x can be decomosed as x, x 2, where x and x 2 are and 2 -dimensional vectors, resectively. Suose now that we have N indeendent observations on the full set of variables, x, and N 2 indeendent observations on x. We introduce some notation for the samle means and covariance matrices. Let x denote the samle, mean of x based on the N observations, and x = x, x 2 x i : i, i =, 2. Let further x 2 denote the -dimensional samle mean vector for based on the N 2 observations. Throughout this aer, we use the letter j only as running suffix for samle observations. With above notations, the 2-ste monotone missing data scheme, described in Shutoh et al. 20 can be resented by x, x 2,..., x N, x 2, x 2 2,..., x 2 N 2, where x j j =,..., N denotes the -dimensional samle vector and x 2 j j = 3
4 ,..., N 2 denotes the -dimensional samle vector. The grahical reresentation of this scheme is given in Figure.. Please insert Figure. around here. Further the samle covariance matrices based on the N and N 2 observations are exressed as S = N x j x x j x, S 2 = n n 2 j= resectively, where n i = N i, i =, 2. N 2 j= x 2 j x 2 x 2 Let the artitions of µ, Σ and S corresonding to ones of x be µ Σ Σ µ =, Σ = 2, S S S 2 = µ 2 Σ 2 Σ 22 S 2 S. 22, j x 2, resectively. Let ˆµ and ˆΣ denote the MLE of µ and Σ, resectively, and assume that ˆµ and ˆΣ can be artitioned in the same way as µ and Σ. From Anderson and Olkin 985, we can exress the MLE s ˆµ and ˆΣ as follows where ˆµ = N N x + N 2 x 2, ˆµ 2 = x 2 ˆΣ ˆΣ 2 x ˆµ, ˆΣ = W + W 2, ˆΣ2 = N ˆΣ W W 2, ˆΣ 22 = N W 22 + ˆΣ 2 ˆΣ ˆΣ 2, N = N + N 2, W = n S, W 2 = n 2 S 2 + N N 2 N W W 2 W = W 2 W 22 x x 2 x x 2,, W 22 = W 22 W 2 W W 2. In Lemma 2. stated below, we roose two auxiliary unbiased estimators of Σ under 2-ste monotone missing data. Lemma 2.. Let ˆΣ and ˆΣ 2 be defined as ˆΣ, ˆΣ,2 ˆΣ2, ˆΣ2,2 ˆΣ =, ˆΣ2 =, ˆΣ,2 ˆΣ,22 ˆΣ 2,2 ˆΣ2,22 4
5 where ˆΣ, = W + W 2, ˆΣ,2 = N ˆΣ W W 2, ˆΣ,2 = ˆΣ,2, ˆΣ,22 = b W 22 + ˆΣ ˆΣ 2 ˆΣ2, ˆΣ 2, = W + W 2, ˆΣ2,2 = N ˆΣ 2,2 = ˆΣ 2,2, ˆΣ2,22 = N W 22, N W 2, and Then b = N N N 2. N N 2 E[ˆΣ ] = E[ˆΣ 2 ] = Σ. Proof : See Aendix A.2. Due to the unbiasedness roerty, both estimators, ˆΣ and ˆΣ 2 are exected to erform better than S which is based on the comlete data art only. The following lemma rovides exact results on the relative erformance accuracy of these three estimators in terms of exected quadratic loss. Lemma 2.2. Let β 0 = a 2 + a 2, n n where a i = tr Σ i / and Σ i denotes i-th ower of the matrix Σ for i =, 2. Then, i E[ S Σ 2 F ] = β 0, ii E[ ˆΣ Σ 2 F ] = β 0 m, iii E[ ˆΣ 2 Σ 2 F ] = β 0 m 2, where 2 F denotes the normalized Frobenious norm, m = N { 2 a 2 tr Σ 2 22 tr Σ 22 2 tr Σ tr Σ 22 + a 2 k k 2 2k 3 nn + tr Σ 22 Σ 2 Σ Σ 2 + tr Σ 22 tr Σ 2 Σ } Σ 2, m 2 = N 2 a 2 + a 2 nn, 5
6 and k = k 2 = k 3 = 3r + 4r 2 2r 3 r 2 4 r + 6r 2 r r r 2 n n 2 r n r 3 }, 2. n { r 4r + 7r 2 4r 3 2rr rr r3 r 2 4rr 2 2 5r 2 + 6r 5 + 8r2 r 2 2r r 2 2r + 7 n n 2 + 2r { 2 r 2 3 r r } n r 3n, 2.2 n 3 r n, 2.3 r = n, r 2 = n n. Proof : See Aendix A.3. It follows from that k i = + On for all i, assuming that is fixed and n is relatively large. Then the following relationshi between m and m 2 holds m 2 = m + N 2 trσ2 Σ Σ Σ 2Σ 2 nn + tr Σ tr Σ 2 Σ Σ o, n + tr Σ 2Σ Σ 2 2 which imlies that ˆΣ is more accurate estimator than ˆΣ 2 in a large samle case. However, ˆΣ will be worse than S since k, k 2 and k 3 grow when /N. On the other hand, since ˆΣ 2 does not deend on the relationshi between and N, its accuracy is better than S. By this reason, it can be said that ˆΣ 2 is a flexibly accurate estimator. 3 The shrinkage estimator for high-dimensional data with 2-ste monotone missing values On the comlete samles, Ledoit and Wolf 2004 rovide a shrinkage estimator for Σ given by ˆΣ f = ˆβ f ˆδ f ˆµ f I + ˆα f ˆδ f S, 6
7 where ˆα f = â 2 â 2, ˆβf = n â 2, ˆδ f = â 2 + â 2 n, ˆµ f = â. In this study, we extend this aroach to the case of 2-ste monotone missing data using auxiliary estimators ˆΣ and ˆΣ 2 derived in Section 2. Ledoit and Wolf tye estimator is derived by the linear combination Σ = ρ I + ρ 2 ˆΣb of the identity matrix I and the unbiased estimator ˆΣ b = ˆΣ or ˆΣ 2 whose exected quadratic loss is minimum. Observed that ˆΣ = ˆΣ f if ˆΣ b = S. Following the technique of Ledoit and Wolf 2004, consider now the otimization roblem min E[ Σ Σ 2 F ] s.t. Σ = ρ I + ρ 2 ˆΣb, 3. ρ, ρ 2 where the coefficients ρ and ρ 2 are non-random. The solution of 3. is given by Σ = β δ µi + α δ ˆΣ b, where µ = tr Σ/, α = Σ µi 2 F, β = E[ ˆΣ b Σ 2 F ], δ = E[ ˆΣ b µi 2 F ]. With this notations we obtain µ = a, α = a 2 a The suggested estimator based on ˆΣ A crucial ste of our estimation is to get the exlicit exressions of E[ ˆΣ Σ 2 F ] and E[ ˆΣ µi 2 F ] in terms of a i and a 22 i consider the following conditions which are functions of tr Σ 22. We C0 : n, n 2,,, C : n c 0 [0,, 0 < lim a i = lim tr Σ i / <, i =,..., 4. Using Lemma A., under the conditions C0 and C, n i c i 0,, i =, 2. E[ ˆΣ Σ 2 F ] = β + o, E[ ˆΣ µi 2 F ] = δ + o, 7
8 where β = a 2 m n, δ = a 2 + a 2 m n, m = m a, a 22 = N { 2 a 2 2k 3 a 22 a + 2k 3 k 2 a 22 2 nn and a 22 = tr Σ 22 /. Thus, the asymtotically otimal linear combination of the identity matrix I and ˆΣ is given by Σ = β δ µi + α δ ˆΣ. } Observe however that Σ deends on unknown arameters. Therefore, the next ste is to relace these arameters their consistent estimators, and show that the asymtotic roerties of the resulting estimator are unchanged. Under the conditions C0 and C, the consistent estimators of α, β and δ are given by ˆα = â 2 â 2, ˆβ = n â 2 ˆm, ˆδ = /n â 2 + â 2 ˆm, 3.2 where ˆm = m â, â 22, and â, â 2 and â 22 are consistent estimators of the corresonding arameters â = â 22 tr S, â 2 = n 2 n n + 2 tr S 2 tr S 2, n = n tr S 22 n. 3.3 Now, by combining 3.2 and 3.3, we define the our estimator as ˆΣ = ˆβ ˆδ ˆµI + ˆαˆδ ˆΣ. One can see that ˆΣ has the same asymtotic roerties as Σ : ˆΣ Σ 2 F P 0 and E[ ˆΣ Σ 2 F ] E[ Σ Σ 2 F ] under conditions C0 and C. Consider further the relationshi between ˆΣ and ˆΣ f by Ledoit and Wolf 2004 when using comlete data art only. The following roosition rovides conditions for the estimator ˆΣ f to dominate ˆΣ f under the exected quadratic loss. 8
9 Proosition 3.. Under the conditions C and n, n 2,, with /n 0, /n i c i 0,, i =, 2 and n i /n + n 2 c i 0,, i = 3, 4 E[ ˆΣ f Σ 2 F ] E[ ˆΣ Σ 2 F ] c c 4 α δ 0 δ a 2 a Proof : See Aendix A The suggested estimator based on ˆΣ 2 The dominance roerty of ˆΣ over ˆΣ f follows from the condition /n 0, which is very strong. Furthermore, if > n, ˆΣ is not alicable due to singularity of ˆΣ. In this subsection, we get another condition for the estimator of Σ to dominate ˆΣ f in the case of > n. Using the technique of subsection 3. and by considering C3 : n, n 2,,, n c 0 [0,, n i /n + n 2 c i 0,, i = 3, 4, n i c i 0,, i =, 2, we obtain E[ ˆΣ 2 Σ 2 F ] = β 2 + o, E[ ˆΣ 2 µi 2 F ] = δ 2 + o, where β 2 = n a 2 m 2, δ 2 = a 2 + a 2 m n 2, m 2 = m 2a = N 2 a 2 nn and a = tr Σ /. Observe that C3 resents more flexible asymtotic framework than C0. Now, the asymtotically otimal linear combination of the identity matrix I and ˆΣ 2 is given by Σ 2 = β 2 δ 2 µi + α δ 2 ˆΣ2. Of this, we roose our estimator of the form ˆΣ 2 = ˆβ 2 ˆδ 2 ˆµI + ˆαˆδ2 ˆΣ2, 9
10 where ˆα = â 2 â 2, ˆβ2 = n â 2 m 2â, ˆδ 2 = /n â 2 + â 2 m 2â and â = tr S /. By the same arguments as in subsection 3., Σ 2 and ˆΣ 2 have the same asymtotic roerties because ˆΣ 2 Σ 2 2 F conditions C and C3. P 0 and E[ ˆΣ 2 Σ 2 F ] E[ Σ 2 Σ 2 F ] under The dominance roerty of ˆΣ 2 over ˆΣ f follows from the following roosition. Proosition 3.2. Under the conditions C and C3, E[ ˆΣ f Σ 2 F ] E[ ˆΣ 2 Σ 2 F ] c c 4 α δ 0 δ 2 a 2. Proof : The roof of this roosition which we omit, is straightforward alication of the technique derived in Aendix A.4. 4 Numerical results Using Monte Carlo simulations, we comare erformance of our new estimators ˆΣ and ˆΣ 2 to the shrinkage estimator by Ledoit and Wolf The simulation exeriment is lanned as follows. The data are generated from the normal model with zero mean vector and covariance matrix Σ = diagσ, σ 2,..., σ ρ i j diagσ, σ 2,..., σ, where σ i is a random variable following continuous uniform distribution on the interval 0, 5. For all simulations we assume N = N 2 for convenience. We set = 00, 50, 200, and for each one consider three different values of ρ = 0.2, 0.5, 0.8 to vary the covariance arameters. The first setting reresents the case when missing art of the data asymtotically exceeds the comlete art. We fix = 0, the size of comlete art, and assume that N = N 2 = 50. This imlies that 2, the dimension of the missing art grows with. 0
11 In the second setting, we vary the roortion of the data that are missing. To ut it more recisely, for each = 00, 50, 200, we set = [q ], for q = 0.5, 0.7, 0.9, and Σ = diagσ, σ 2,..., σ 0.8 i j diagσ, σ 2,..., σ, where σ i is a random variable having continuous uniform distribution on the interval 0, 5. Here, [a] denotes the integer art of a. Within the first setting, three following estimators are evaluated: ˆΣ f, ˆΣ, ˆΣ 2, whereas in the second setting we only focus on ˆΣ f and ˆΣ 2. To evaluate erformance accuracy we use the loss function based on normalized Frobenius norm, LΣ, = Σ 2 F, and define the risk function as the exected loss, RΣ, = E[LΣ, ], when estimating Σ by. We asses the relative risk rate as follows RR = RΣ, /RΣ, S, where S is used reference estimator for all comarisons. For all the settings, the value of RR is comuted by 00,000 relications. Tables reort the risk rates for all three estimators in the for the first setting, and figures lot the risk rates for ˆΣ f and ˆΣ 2 as functions of for the second setting. Overview through the simulation results reveals that our new estimators rovide significant imrovements over ˆΣ f for both tyes of missing scenarios. Although RRs of ˆΣ and ˆΣ 2 are very close to each other for all ρ, ˆΣ 2 dominates ˆΣ in the case when missing art grows. As exected, ˆΣ 2 yields substantial imrovement over ˆΣ f even in the case of fixed missing art. It could be observed that both estimators are sensitive to the size of ρ 0, which is naturally exected in a AR model. In conclusion, both suggested estimators, ˆΣ and ˆΣ 2 are shown to have dominance roerties over ˆΣ f in high-dimensional normal model with 2-ste monotone missing data. Please insert Tables and Figures around here.
12 5 Conclusion This aer rovided the estimators for Σ in the case of 2-ste monotone missing data with high-dimensional setting by extending Ledoit and Wolf s 2004 estimator. The result rovided a solution to the roblem such that the existing MLE on the basis of 2-ste monotone missing data is not also alicable in the case of high-dimensional setting. Further, it could be observed that our estimators were suerior to Ledoit and Wolf s 2004 estimator in view of the relative risk rate for the erformed simulations. In the future, we may also derive the theoretical results of the rocedures with our estimators which are alicable to the high-dimensional and missing data. Aendix To rove roositions 2. and 2.2 we need two additional lemmas stated below. For these roofs we need the some reliminary results formulated in Lemma A. and Lemma A.2. A. Preliminary results Lemma 5.. Suose that A follows Wishart distribution W Σ, n, where A admits the similar artition as that of Σ in Section 2, and let A 22 = A 22 A 2 A A 2. Then i A 22 W Σ 22, n, and A 22 is indeendent of A and A 2, ii A 2 A N 2 V eca Σ Σ 2, Σ 22 A Proof. For the roof, see Muirhead 982 and Siotani et al Now we state Lemma A.2. 2
13 Lemma 5.2. For ˆΣ defined in Section 2, it holds that E[tr ˆΣ 2 ] = n + n tr Σ 2 + n tr Σ 2 N 2 nn {tr Σ 2 + tr Σ 2 k tr Σ 2 22 k 2 tr Σ k 3 tr Σ 22 tr Σ + tr Σ 22 Σ 2 Σ Σ 2 + tr Σ 22 tr Σ 2 Σ Σ 2 }, where k i, i =, 2, 3 is defined in Lemma 2.2. Proof. Using the following exression E[tr ˆΣ 2 ] = E[trˆΣ 2, + ˆΣ,2 ˆΣ,2 ] + E[trˆΣ 2,22 + ˆΣ,2 ˆΣ,2 ] A. and the results of Lemma A., we obtain E[trˆΣ 2, + ˆΣ,2 ˆΣ,2 ] = n + n trσ2 + Σ 2 Σ 2 + n tr Σ tr Σ N 2 + nn tr Σ 22 tr Σ, A.2 For the second term in A. we get E[trˆΣ 2,22 + ˆΣ,2 ˆΣ,2 ] = n 2 E [tr ˆΨ 2 W 2 ˆΨ 2 ] + b 2 E + b n E [tr ˆΨ 2 W ˆΨ2 W 22 ] + b n E [tr W 22 ˆΨ 2 W ˆΨ2 ] ] [tr W n 2 E [tr ˆΨ 2 W ˆΨ2 ˆΨ2 W ˆΨ2 ], A.3 3
14 where ˆΨ 2 = W 2 W. From Lemma A., E[tr ˆΨ 2 W 2 ˆΨ 2 ] = nn + tr Σ 2 Σ 2 + ntr Σ tr Σ 2 Σ Σ 2 + tr Σ 22 tr Σ + E[tr W 2 22 ] = n n tr Σ n tr Σ 22 2, E[tr ˆΨ 2 W ˆΨ2 W 22 ] = n + E[tr ˆΨ 2 W ˆΨ2 2 ] = +nn tr Σ 2 Σ { N 2 nn 2 n tr Σ tr Σ 22, N 2 n } Σ 2 Σ 22 tr Σ 2 22 n + c 2N 2 + N2 2 2 } tr Σ c 2 N 2 + N N2 2 { N 2 A.4, A.5 n + c 22N 2 + N2 2 2 } +c 2 n N 2 + N2 2 + N 2 2 +{tr Σ tr Σ 22 2 } + 2 n + N 2 + N N 2 + N 2 n +2n + N 2 n +n 2 + n trσ 2 Σ Σ 2 2 tr Σ 22 Σ 2 Σ Σ 2 tr Σ 22 tr Σ 2 Σ Σ 2 +n tr Σ 2 Σ Σ 2 2. A.6 Substituting A.4-A.6 into A.3, the roof of Lemma A. is comlete. 4
15 A.2 Proof of Lemma 2. Using the moment roerties of the Wishart distribution and alying Lemma A., we get E[ˆΣ, ] = [ E N = ] W + W 2 N {N Σ + N 2 Σ } = Σ, A.7 [ ] E[ˆΣ,2 ] = N E W + W 2 W W 2 [ ] = N E W 2 + W 2 W W 2 = N {N Σ 2 + N 2 Σ 2 } = Σ 2, A.8 E[ˆΣ,22 ] = N N N 2 E[W N N 2 22 ] + [ ] N E W 2 W W + W 2 W W 2.A.9 Then, for the first term of A.9 we have N = N N 2 N N 2 Σ 22. N N 2 N N 2 For the second term of A.9 we obtain [ N E W 2 W = Σ 22 W + W 2 W N N 2 N N 2 E[W 22 ] ] W 2 Σ 22. A.0 A. Substituting A.0 and A. into A.9, we have E[ˆΣ,22 ] = Σ 22. A.2 Combining A.7, A.8 and A.2, we get E[ˆΣ ] = Σ. The unbiasedness of ˆΣ 2 is easily shown by calculating the first moment of the Wishart matrices W and W 2. 5
16 A.3 Proof of Lemma 2.2 The result in i is established in see, e.g., Ledoit and Wolf To show ii we first exress it as E[ ˆΣ Σ 2 F ] = E[tr ˆΣ 2 ]/ a 2. By using Lemma A.2, we can calculate E[tr ˆΣ 2 ], which gives the statement ii. It should be noted that i is shown assuming that N 2 = 0. Finaly, we give the roof of iii, by first exressing E[ ˆΣ 2 Σ 2 F ] = E[tr ˆΣ 2 2]/ a 2 A.3 and using the decomosition of ˆΣ 2 as n n S + N 2 n S2 S 2 S 2 S. 22 Then E [tr ˆΣ 2 2] = = { E [ n tr n S + N 2 ]} [tr S 2 n S2 + E 22 { n + tr Σ 2 + tr Σ 2 n n 2 ] } + 2 E[tr S 2 S 2 ] N 2 nn { tr Σ 2 + tr Σ 2}. A.4 Substituting A.4 into A.3, we obtain iii. It should be noted that i is shown assuming that N 2 = 0. A.4 Proof of Proosition 3. Under the conditions C and assuming that n, n 2,, with /n 0 and /n i c i 0,, i =, 2, E[ ˆΣ Σ 2 F ] αβ δ and E[ ˆΣ f Σ 2 F ] αβ 0 δ 0, A.5 6
17 where δ 0 = α+β 0. By using the continuous maing theorem and A.5, we obtain E[ ˆΣ f Σ 2 F ] E[ ˆΣ Σ 2 F ] αβ 0 δ 0 αβ δ. A.6 On the other hand, under the same assumtions as above, αβ 0 δ 0 αβ δ = α δ 0 δ β 0 β = α δ 0 δ m N 2α nn δ 0 δ a 2 a A.7 Combining A.6 and A.7, Proosition 3. follows. Acknowledgments The authors would like to extend their sincere gratitude to the referee who gave invaluable comments and suggestions, which have greatly enhanced this aer. Any remaining errors are the authors resonsibility. This research was in art suorted by grants from the Forum for Asian Studies , Stockholm University, Sweden. M.Hyodo s research was in art suorted by Grant-in-Aid for JSPS Fellow T. Seo s research was in art suorted by Grant-in-Aid for Scientific Research C T. Pavlenko s research was in art suorted by grants from Swedish Research Council References Anderson, T. W Maximum likelihood estimates for multivariate normal distribution when some observations are missing. Journal of the American Statistical Association, 52, Anderson, T.W. and Olkin, I Maximum likelihood estimation of the arameters of multivariate normal distribution. Linear Algebra and Its Alications, 70, Bhargava, R. P Some one-samle hyothesis testing roblems when there is monotone samle from a multivariate normal oulation. Annals of the Institute of Statistical Mathematics, 27,
18 Fujikoshi, Y. and Seo, T Asymtotic aroximations for EPMC s of the linear and the quadratic discriminant function when the samle size and the dimension are large. Random Oerators and Stochastic Equations, 6, Jinadasa, K. G. and Tracy, D. S Maximum likelihood estimation for multivariate normal distribution with monotone samle. Communications in Statistics A, Theory and Methods, 2, Kanda, T. and Fujikoshi, Y Some basic roerties of the MLE s for a multivariate normal distribution with monotone missing data. American Journal of Mathematical and Management Sciences, 8, Kubokawa, T. and Srivastava, M. S Estimation of the recision matrix of a singular Wishart distribution and its alication in high-dimensional data. Journal of Multivariate Analysis, 99, Little, R. J. A. and Rubin, D. R Statistical Analysis with Missing Data, Wiley, New York. Ledoit, O. and Wolf, M A well-conditioned estimator for large dimensional covariance matrices. Journal of Multivariate Analysis, 88, Muirhead, R. J Asects of multivariate statistical theory, Wiley, New York. Shutoh, N., Hyodo, M. and Seo, T. 20. An asymtotic aroximation for EPMC in linear discriminant analysis based on two-ste monotone missing samles. Journal of Multivariate Analysis, 02, Siotani, M., Hayakawa, T. and Fujikoshi, Y Modern Multivariate Statistical Analysis, A Graduate Course and Handbook. Ohio:American Science Press. 8
19 Table 4.: The risk rates when ρ = 0.2 = 00 = 50 = 200 RRˆΣ f RRˆΣ RRˆΣ Table 4.2: The risk rates when ρ = 0.5 = 00 = 50 = 200 RRˆΣ f RRˆΣ RRˆΣ Table 4.3: The risk rates when ρ = 0.8 = 00 = 50 = 200 RRˆΣ f RRˆΣ RRˆΣ
20 N N 2 missing data 2 Figure.: 2-ste monotone missing data RR lotted as a function of the number of variables for = ˆΣ fsolid and = ˆΣ 2dashed. Fig.4.. q = 0.5. Fig.4.2. q = 0.7. Fig.4.3. q =
ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE
J Jaan Statist Soc Vol 34 No 2004 9 26 ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE Yasunori Fujikoshi*, Tetsuto Himeno
More informationA Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression
Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi
More informationarxiv: v2 [stat.me] 3 Nov 2014
onarametric Stein-tye Shrinkage Covariance Matrix Estimators in High-Dimensional Settings Anestis Touloumis Cancer Research UK Cambridge Institute University of Cambridge Cambridge CB2 0RE, U.K. Anestis.Touloumis@cruk.cam.ac.uk
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationA New Asymmetric Interaction Ridge (AIR) Regression Method
A New Asymmetric Interaction Ridge (AIR) Regression Method by Kristofer Månsson, Ghazi Shukur, and Pär Sölander The Swedish Retail Institute, HUI Research, Stockholm, Sweden. Deartment of Economics and
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationUsing the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process
Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund
More informationChapter 3. GMM: Selected Topics
Chater 3. GMM: Selected oics Contents Otimal Instruments. he issue of interest..............................2 Otimal Instruments under the i:i:d: assumtion..............2. he basic result............................2.2
More informationHotelling s Two- Sample T 2
Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationIntroduction to Probability and Statistics
Introduction to Probability and Statistics Chater 8 Ammar M. Sarhan, asarhan@mathstat.dal.ca Deartment of Mathematics and Statistics, Dalhousie University Fall Semester 28 Chater 8 Tests of Hyotheses Based
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationOn split sample and randomized confidence intervals for binomial proportions
On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have
More informationExpected probabilities of misclassification in linear discriminant analysis based on 2-Step monotone missing samples
Expected probabilities of misclassification in linear discriminant analysis based on 2-Step monotone missing samples Nobumichi Shutoh, Masashi Hyodo and Takashi Seo 2 Department of Mathematics, Graduate
More informationAsymptotically Optimal Simulation Allocation under Dependent Sampling
Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationTests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test)
Chater 225 Tests for Two Proortions in a Stratified Design (Cochran/Mantel-Haenszel Test) Introduction In a stratified design, the subects are selected from two or more strata which are formed from imortant
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationGeneral Linear Model Introduction, Classes of Linear models and Estimation
Stat 740 General Linear Model Introduction, Classes of Linear models and Estimation An aim of scientific enquiry: To describe or to discover relationshis among events (variables) in the controlled (laboratory)
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More informationSOME TESTS CONCERNING THE COVARIANCE MATRIX IN HIGH DIMENSIONAL DATA
J. Jaan Statist. Soc. Vol. 35 No. 005 5 7 SOME TESTS CONCERNING THE COVARIANCE MATRIX IN HIGH DIMENSIONAL DATA Muni S. Srivastava* In this aer, tests are develoed for testing certain hyotheses on the covariance
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationCOMMUNICATION BETWEEN SHAREHOLDERS 1
COMMUNICATION BTWN SHARHOLDRS 1 A B. O A : A D Lemma B.1. U to µ Z r 2 σ2 Z + σ2 X 2r ω 2 an additive constant that does not deend on a or θ, the agents ayoffs can be written as: 2r rθa ω2 + θ µ Y rcov
More informationLower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data
Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment
More informationEstimating function analysis for a class of Tweedie regression models
Title Estimating function analysis for a class of Tweedie regression models Author Wagner Hugo Bonat Deartamento de Estatística - DEST, Laboratório de Estatística e Geoinformação - LEG, Universidade Federal
More informationŽ. Ž. Ž. 2 QUADRATIC AND INVERSE REGRESSIONS FOR WISHART DISTRIBUTIONS 1
The Annals of Statistics 1998, Vol. 6, No., 573595 QUADRATIC AND INVERSE REGRESSIONS FOR WISHART DISTRIBUTIONS 1 BY GERARD LETAC AND HELENE ` MASSAM Universite Paul Sabatier and York University If U and
More informationSupplementary Materials for Robust Estimation of the False Discovery Rate
Sulementary Materials for Robust Estimation of the False Discovery Rate Stan Pounds and Cheng Cheng This sulemental contains roofs regarding theoretical roerties of the roosed method (Section S1), rovides
More informationImproved Bounds on Bell Numbers and on Moments of Sums of Random Variables
Imroved Bounds on Bell Numbers and on Moments of Sums of Random Variables Daniel Berend Tamir Tassa Abstract We rovide bounds for moments of sums of sequences of indeendent random variables. Concentrating
More informationMULTIVARIATE SHEWHART QUALITY CONTROL FOR STANDARD DEVIATION
MULTIVARIATE SHEWHART QUALITY CONTROL FOR STANDARD DEVIATION M. Jabbari Nooghabi, Deartment of Statistics, Faculty of Mathematical Sciences, Ferdowsi University of Mashhad, Mashhad-Iran. and H. Jabbari
More informationTesting Block-Diagonal Covariance Structure for High-Dimensional Data
Testing Block-Diagonal Covariance Structure for High-Dimensional Data MASASHI HYODO 1, NOBUMICHI SHUTOH 2, TAKAHIRO NISHIYAMA 3, AND TATJANA PAVLENKO 4 1 Department of Mathematical Sciences, Graduate School
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationSign Tests for Weak Principal Directions
Sign Tests for Weak Princial Directions DAVY PAINDAVEINE, JULIEN REMY and THOMAS VERDEBOUT arxiv:1812.09367v1 [math.st] 21 Dec 2018 Université libre de Bruxelles ECARES and Déartement de Mathématique Avenue
More informationA multiple testing approach to the regularisation of large sample correlation matrices
A multile testing aroach to the regularisation of large samle correlation matrices Natalia Bailey Queen Mary, University of London M. Hashem Pesaran University of Southern California, USA, and rinity College,
More informationVIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES
Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The
More informationUnobservable Selection and Coefficient Stability: Theory and Evidence
Unobservable Selection and Coefficient Stability: Theory and Evidence Emily Oster Brown University and NBER August 9, 016 Abstract A common aroach to evaluating robustness to omitted variable bias is to
More informationCovariance Matrix Estimation for Reinforcement Learning
Covariance Matrix Estimation for Reinforcement Learning Tomer Lancewicki Deartment of Electrical Engineering and Comuter Science University of Tennessee Knoxville, TN 37996 tlancewi@utk.edu Itamar Arel
More informationAsymptotic non-null distributions of test statistics for redundancy in high-dimensional canonical correlation analysis
Asymtotic non-null distributions of test statistics for redundancy in high-dimensional canonical correlation analysis Ryoya Oda, Hirokazu Yanagihara, and Yasunori Fujikoshi, Deartment of Mathematics, Graduate
More informationYixi Shi. Jose Blanchet. IEOR Department Columbia University New York, NY 10027, USA. IEOR Department Columbia University New York, NY 10027, USA
Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelsach, K. P. White, and M. Fu, eds. EFFICIENT RARE EVENT SIMULATION FOR HEAVY-TAILED SYSTEMS VIA CROSS ENTROPY Jose
More informationDeriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.
Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationAN OPTIMAL CONTROL CHART FOR NON-NORMAL PROCESSES
AN OPTIMAL CONTROL CHART FOR NON-NORMAL PROCESSES Emmanuel Duclos, Maurice Pillet To cite this version: Emmanuel Duclos, Maurice Pillet. AN OPTIMAL CONTROL CHART FOR NON-NORMAL PRO- CESSES. st IFAC Worsho
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationarxiv: v2 [stat.me] 27 Apr 2018
Nonasymtotic estimation and suort recovery for high dimensional sarse covariance matrices arxiv:17050679v [statme] 7 Ar 018 1 Introduction Adam B Kashlak and Linglong Kong Deartment of Mathematical and
More informationHigh-Dimensional Asymptotic Behaviors of Differences between the Log-Determinants of Two Wishart Matrices
TR-No. 4-0, Hiroshima Statistical Research Grou, 9 High-Dimensional Asymtotic Behaviors of Differences between the Log-Determinants of Two Wishart Matrices Hirokazu Yanagihara,2,3, Yusuke Hashiyama and
More informationESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION
STATISTICA, anno LXVIII, n., 008 ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION 1. INTRODUCTION The Inverse Gaussian distribution was first introduced
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationThe spectrum of kernel random matrices
The sectrum of kernel random matrices Noureddine El Karoui Deartment of Statistics, University of California, Berkeley Abstract We lace ourselves in the setting of high-dimensional statistical inference,
More informationAn Improved Generalized Estimation Procedure of Current Population Mean in Two-Occasion Successive Sampling
Journal of Modern Alied Statistical Methods Volume 15 Issue Article 14 11-1-016 An Imroved Generalized Estimation Procedure of Current Poulation Mean in Two-Occasion Successive Samling G. N. Singh Indian
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationOn the sphericity test with large-dimensional observations. Citation Electronic Journal of Statistics, 2013, v. 7, p
Title On the shericity test with large-dimensional observations Authors Wang, Q; Yao, J Citation Electronic Journal of Statistics, 203, v. 7,. 264 292 Issued Date 203 URL htt://hdl.handle.net/0722/89449
More informationON UNIFORM BOUNDEDNESS OF DYADIC AVERAGING OPERATORS IN SPACES OF HARDY-SOBOLEV TYPE. 1. Introduction
ON UNIFORM BOUNDEDNESS OF DYADIC AVERAGING OPERATORS IN SPACES OF HARDY-SOBOLEV TYPE GUSTAVO GARRIGÓS ANDREAS SEEGER TINO ULLRICH Abstract We give an alternative roof and a wavelet analog of recent results
More informationOn the asymptotic sizes of subset Anderson-Rubin and Lagrange multiplier tests in linear instrumental variables regression
On the asymtotic sizes of subset Anderson-Rubin and Lagrange multilier tests in linear instrumental variables regression Patrik Guggenberger Frank Kleibergeny Sohocles Mavroeidisz Linchun Chen\ June 22
More informationApplicable Analysis and Discrete Mathematics available online at HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS
Alicable Analysis and Discrete Mathematics available online at htt://efmath.etf.rs Al. Anal. Discrete Math. 4 (010), 3 44. doi:10.98/aadm1000009m HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS Zerzaihi
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationMaxisets for μ-thresholding rules
Test 008 7: 33 349 DOI 0.007/s749-006-0035-5 ORIGINAL PAPER Maxisets for μ-thresholding rules Florent Autin Received: 3 January 005 / Acceted: 8 June 006 / Published online: March 007 Sociedad de Estadística
More informationElements of Asymptotic Theory. James L. Powell Department of Economics University of California, Berkeley
Elements of Asymtotic Theory James L. Powell Deartment of Economics University of California, Berkeley Objectives of Asymtotic Theory While exact results are available for, say, the distribution of the
More informationThe power performance of fixed-t panel unit root tests allowing for structural breaks in their deterministic components
ATHES UIVERSITY OF ECOOMICS AD BUSIESS DEPARTMET OF ECOOMICS WORKIG PAPER SERIES 23-203 The ower erformance of fixed-t anel unit root tests allowing for structural breaks in their deterministic comonents
More informationMULTIVARIATE STATISTICAL PROCESS OF HOTELLING S T CONTROL CHARTS PROCEDURES WITH INDUSTRIAL APPLICATION
Journal of Statistics: Advances in heory and Alications Volume 8, Number, 07, Pages -44 Available at htt://scientificadvances.co.in DOI: htt://dx.doi.org/0.864/jsata_700868 MULIVARIAE SAISICAL PROCESS
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationAsymptotic Analysis of the Squared Estimation Error in Misspecified Factor Models
Asymtotic Analysis of the Squared Estimation Error in Missecified Factor Models Alexei Onatski July 7, 2013 Abstract In this aer, we obtain asymtotic aroximations to the squared error of the least squares
More informationUsing a Computational Intelligence Hybrid Approach to Recognize the Faults of Variance Shifts for a Manufacturing Process
Journal of Industrial and Intelligent Information Vol. 4, No. 2, March 26 Using a Comutational Intelligence Hybrid Aroach to Recognize the Faults of Variance hifts for a Manufacturing Process Yuehjen E.
More informationOn the Maximum and Minimum of Multivariate Pareto Random Variables
ISSN 1392 124X (rint) ISSN 2335 884X (online) INFORMATION TECHNOLOGY AND CONTROL 213 Vol.42 No.3 On the Maximum Minimum of Multivariate Pareto Rom Variables Jungsoo Woo Deartment of Statistics Yeungnam
More informationProximal methods for the latent group lasso penalty
Proximal methods for the latent grou lasso enalty The MIT Faculty has made this article oenly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher Villa,
More informationSpectral Analysis by Stationary Time Series Modeling
Chater 6 Sectral Analysis by Stationary Time Series Modeling Choosing a arametric model among all the existing models is by itself a difficult roblem. Generally, this is a riori information about the signal
More informationMinimax Design of Nonnegative Finite Impulse Response Filters
Minimax Design of Nonnegative Finite Imulse Resonse Filters Xiaoing Lai, Anke Xue Institute of Information and Control Hangzhou Dianzi University Hangzhou, 3118 China e-mail: laix@hdu.edu.cn; akxue@hdu.edu.cn
More informationPositivity, local smoothing and Harnack inequalities for very fast diffusion equations
Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Dedicated to Luis Caffarelli for his ucoming 60 th birthday Matteo Bonforte a, b and Juan Luis Vázquez a, c Abstract
More informationLocation of solutions for quasi-linear elliptic equations with general gradient dependence
Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations
More informationSome Measures of Agreement Between Close Partitions
Some Measures of Agreement Between Close Partitions Genane Youness and Gilbert Saorta CEDRIC CNAM, BP - 466, Beirut, Lebanon, genane99@hotmail.com Chaire de Statistique Aliquée- CEDRIC, CNAM, 9 rue Saint
More informationA Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split
A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com
More informationBEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH
BEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH DORIN BUCUR, ALESSANDRO GIACOMINI, AND PAOLA TREBESCHI Abstract For Ω R N oen bounded and with a Lischitz boundary, and
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationAn Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices
Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem
More informationApproximation of the Euclidean Distance by Chamfer Distances
Acta Cybernetica 0 (0 399 47. Aroximation of the Euclidean Distance by Chamfer Distances András Hajdu, Lajos Hajdu, and Robert Tijdeman Abstract Chamfer distances lay an imortant role in the theory of
More information2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized
BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,
More informationElements of Asymptotic Theory. James L. Powell Department of Economics University of California, Berkeley
Elements of Asymtotic Theory James L. Powell Deartment of Economics University of California, Berkeley Objectives of Asymtotic Theory While exact results are available for, say, the distribution of the
More informationCONVOLVED SUBSAMPLING ESTIMATION WITH APPLICATIONS TO BLOCK BOOTSTRAP
Submitted to the Annals of Statistics arxiv: arxiv:1706.07237 CONVOLVED SUBSAMPLING ESTIMATION WITH APPLICATIONS TO BLOCK BOOTSTRAP By Johannes Tewes, Dimitris N. Politis and Daniel J. Nordman Ruhr-Universität
More informationOn-Line Appendix. Matching on the Estimated Propensity Score (Abadie and Imbens, 2015)
On-Line Aendix Matching on the Estimated Proensity Score Abadie and Imbens, 205 Alberto Abadie and Guido W. Imbens Current version: August 0, 205 The first art of this aendix contains additional roofs.
More informationAdaptive Estimation of the Regression Discontinuity Model
Adative Estimation of the Regression Discontinuity Model Yixiao Sun Deartment of Economics Univeristy of California, San Diego La Jolla, CA 9293-58 Feburary 25 Email: yisun@ucsd.edu; Tel: 858-534-4692
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationRatio Estimators in Simple Random Sampling Using Information on Auxiliary Attribute
ajesh Singh, ankaj Chauhan, Nirmala Sawan School of Statistics, DAVV, Indore (M.., India Florentin Smarandache Universit of New Mexico, USA atio Estimators in Simle andom Samling Using Information on Auxiliar
More informationarxiv: v4 [math.nt] 11 Oct 2017
POPULAR DIFFERENCES AND GENERALIZED SIDON SETS WENQIANG XU arxiv:1706.05969v4 [math.nt] 11 Oct 2017 Abstract. For a subset A [N], we define the reresentation function r A A(d := #{(a,a A A : d = a a }
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analysis of Variance and Design of Exeriment-I MODULE II LECTURE -4 GENERAL LINEAR HPOTHESIS AND ANALSIS OF VARIANCE Dr. Shalabh Deartment of Mathematics and Statistics Indian Institute of Technology Kanur
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More informationPublished: 14 October 2013
Electronic Journal of Alied Statistical Analysis EJASA, Electron. J. A. Stat. Anal. htt://siba-ese.unisalento.it/index.h/ejasa/index e-issn: 27-5948 DOI: 1.1285/i275948v6n213 Estimation of Parameters of
More informationRobust Beamforming via Matrix Completion
Robust Beamforming via Matrix Comletion Shunqiao Sun and Athina P. Petroulu Deartment of Electrical and Comuter Engineering Rutgers, the State University of Ne Jersey Piscataay, Ne Jersey 8854-858 Email:
More informationShadow Computing: An Energy-Aware Fault Tolerant Computing Model
Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
More informationSome results of convex programming complexity
2012c12 $ Ê Æ Æ 116ò 14Ï Dec., 2012 Oerations Research Transactions Vol.16 No.4 Some results of convex rogramming comlexity LOU Ye 1,2 GAO Yuetian 1 Abstract Recently a number of aers were written that
More informationMorten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014
Morten Frydenberg Section for Biostatistics Version :Friday, 05 Setember 204 All models are aroximations! The best model does not exist! Comlicated models needs a lot of data. lower your ambitions or get
More informationUnderstanding and Using Availability
Understanding and Using Availability Jorge Luis Romeu, Ph.D. ASQ CQE/CRE, & Senior Member Email: romeu@cortland.edu htt://myrofile.cos.com/romeu ASQ/RD Webinar Series Noviembre 5, J. L. Romeu - Consultant
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationsubstantial literature on emirical likelihood indicating that it is widely viewed as a desirable and natural aroach to statistical inference in a vari
Condence tubes for multile quantile lots via emirical likelihood John H.J. Einmahl Eindhoven University of Technology Ian W. McKeague Florida State University May 7, 998 Abstract The nonarametric emirical
More information