arxiv: v2 [stat.ml] 3 Sep 2015

Transcription

1 Nonparametric Independence Testing for Small Sample Sizes arxiv:46.9v [stat.ml] 3 Sep 5 Aaditya Ramdas Dept. of Statistics and Machine Learning Dept. Carnegie Mellon University aramdas@cs.cmu.edu September 4, 5 Abstract Leila Wehbe Machine Learning Dept. Carnegie Mellon University lwehbe@cs.cmu.edu This paper deals with the problem of nonparametric independence testing, a fundamental decision-theoretic problem that asks if two arbitrary (possibly multivariate) random variables X, Y are independent or not, a question that comes up in many fields like causality and neuroscience. While quantities like correlation of X, Y only test for (univariate) linear independence, natural alternatives like mutual information of X, Y are hard to estimate due to a serious curse of dimensionality. A recent approach, avoiding both issues, estimates norms of an operator in Reproducing Kernel Hilbert Spaces (RKHSs). Our main contribution is strong empirical evidence that by employing shrunk operators when the sample size is small, one can attain an improvement in power at low false positive rates. We analyze the effects of Stein shrinkage on a popular test statistic called HSIC (Hilbert-Schmidt Independence Criterion). Our observations provide insights into two recently proposed shrinkage estimators, SCOSE and FCOSE - we prove that SCOSE is (essentially) the optimal linear shrinkage method for estimating the true operator; however, the non-linearly shrunk FCOSE usually achieves greater improvements in test power. This work is important for more powerful nonparametric detection of subtle nonlinear dependencies for small samples. Introduction The problem of nonparametric independence testing deals with ascertaining if two random variables are independent or not, making no parametric assumptions about their underlying distributions. Formally, given n samples (x i, y i ) for i {,..., n} where x i R p, y i R q, that are drawn from a joint distribution P XY supported on X Y R p+q, we want to decide between the null and alternate hypotheses H : P XY = P X P Y vs. H : P XY P X P Y where P X, P Y are the marginals of P XY w.r.t. X, Y. A test is a function from the data to {, }. Tests aim to have high power (probability of detecting dependence, when it exists) at a prespecified allowable type- error rate α (probability of detecting dependence when there isn t any). Independence testing is often a precursor to further analysis. Consider for instance conditional independence testing for inferring causality, say by the PC algorithm Spirtes et al. (), whose first step is (unconditional) independence testing. It is also useful for scientific discovery like in Both authors have equally contributed to this paper.

2 neuroscience, to see if a stimulus X (say an image) is independent of the brain activity Y (say fmri) in a relevant part of the brain. Since detecting nonlinear correlations is much easier than estimating a nonparametric regression function (of Y onto X), it can be done at smaller sample sizes, with further samples collected for estimation only if an effect is detected by the hypothesis test. For such situations, correlation only tests for univariate linear independence, while other statistics like mutual information that do characterize multivariate independence are hard to estimate from data, suffering from a serious curse of dimensionality. A recent popular approach for this problem (and a related two-sample testing problem) involve the use of quantities defined in reproducing kernel Hilbert spaces (RKHSs) - see Gretton et al. (6), Harchaoui et al. (7), Gretton, Herbrich, Smola, Bousquet & Schölkopf (5), Gretton, Bousquet, Smola & Schölkopf (5). This paper will concern itself with increasing the statistical power at small samples of a popular kernel statistic called HSIC, by using shrunk empirical estimators of the unknown population quantity (introduced below).. Hilbert Schmidt Independence Criterion Due to limited space, familiarity with RKHS terminology is assumed - see Scholkopf & Smola () for an introduction. Let k : X X R and l : Y Y R be two positive-definite reproducing kernels that correspond to RKHSs H k and H l respectively with inner-products, k and, l. Let k, l arise from (implicit) feature maps φ : X H k and ψ : Y H l. In other words, φ, ψ are not functions, but mappings to the Hilbert space. i.e. φ(x) H k, ψ(y) H l respectively. These functions, when evaluated at points in the original spaces, must satisfy φ(x)(x ) = φ(x), φ(x ) k = k(x, x ) and ψ(y)(y ) = ψ(y), ψ(y ) l = l(y, y ). The mean embedding of P X and P Y are defined as µ X := E x PX φ(x) H k and µ Y := E y PY ψ(y) H l whose empirical estimates are µ X := n n i= φ(x i) and µ Y := n n i= ψ(y i). Finally, the cross-covariance operator of X, Y is defined as Σ XY := E (x,y) PXY (φ(x) µ X ) (ψ(y) µ Y ) where is an outer-product. For unfamiliar readers, if we used the linear kernel k(x, x ) = x T x and l(y, y ) = y T y, then the cross-covariance operator is just the cross-covariance matrix. The plug-in empirical estimator of Σ XY is S XY := n n (φ(x i ) µ X ) (ψ(y i ) µ Y ) i= For conciseness, define φ(x i ) = φ(x i ) µ X, ψ(y i ) = ψ(y i ) µ Y, k(x, x ) = φ(x), φ(x ) k and l(y, y ) = ψ(y), ψ(y ) l. The test statistic Hilbert-Schmidt Independence Criterion (HSIC) defined in Gretton, Bousquet, Smola & Schölkopf (5) is the squared Hilbert-Schmidt norm of S XY, and can be calculated using centered kernel matrices K, L, where K ij = k(x i, x j ), L ij = l(y i, y j ), as HSIC := S XY HS = n tr( K L) () For unfamiliar readers, if we used the linear kernel, this just corresponds to the Frobenius norm of the cross-covariance matrix. The most important property is: when the kernels k, l are characteristic, then the corresponding population statistic Σ XY HS is zero iff X, Y are independent Gretton, Bousquet, Smola & Schölkopf (5). This gives rise to a natural test - calculate S XY HS and reject the null if it is large. ( ) x x γ and Laplace Examples of characteristic kernels include Gaussian k(x, x ) = exp ( k(x, x ) = exp x x γ ), for any bandwidth γ, while the aforementioned linear kernel is not

3 characteristic the corresponding HSIC tests only linear relationships, and a zero cross-covariance matrix characterizes independence only for multivariate Gaussian distributions. Working with the infinite dimensional operator with characteristic kernels, allows us to identify any general nonlinear dependence (in the limit) between any pair of distributions, not just Gaussians.. Independence Testing using HSIC A permutation-based test is described in Gretton, Bousquet, Smola & Schölkopf (5), and proceeds in the following manner. From the given data, calculate the test statistic T := S XY HS. Keeping the order of x,..., x n fixed, randomly permute y,..., y n a large number of times, and recompute the permuted HSIC each time. This destroyed any dependence between x, y simulating a draw from the product of marginals, making the empirical distribution of the permuted HSICs behave like the null distribution of the test statistic (distribution of HSIC when H is true). For a pre-specified type- error α, calculate threshold t α in the right tail of the null distribution. Reject H if T > t α. This test was proved to be consistent against any fixed alternative, meaning for any fixed type- error α, the power goes to as n. Empirically, the power can be calculated using simulations by repeating the above permutation test many times for a fixed P XY (for which dependence holds), and reporting the empirical probability of rejecting the null (detecting the dependence). Note that the power depends on P XY (unknown to the user of the test)..3 Shrunk Estimators of S XY Even though S XY is an unbiased estimator of Σ XY, it typically has high variance at low sample sizes. The idea of Stein shrinkage Stein (956) is to trade-off bias and variance, first introduced in the context of Gaussian mean estimation. This strategy of introducing some bias and decreasing the variance to get different estimators of Σ XY was followed by Muandet et al. (4) who define a linear shrinkage estimator of S XY called SCOSE (Simple Covariance Shrinkage Estimator) and a nonlinear shrinkage estimator called FCOSE (Flexible Covariance Shrinkage Estimator). When we refer to shrunk estimators, we implicitly mean SCOSE and FCOSE. We will describe these briefly in Section..4 Contributions Our first contribution is the following :. We provide evidence that employing shrunk estimators of Σ XY, instead of S XY, to calculate the aforementioned test statistic, can increase the power of the associated independence test at low false positive rates, when the sample size is small (there is higher variance in estimating infinitedimensional operators). Our second contribution is to analyze the effect of shrinkage on the test statistic, to provide some practical insight.. The effect of shrinkage on the test-statistic is very similar to soft-thresholding (see Section 4), shrinking very small statistics to zero, and shrinking other values nearly (but not) linearly, and nearly (but not) monotonically. Our last contribution is an insight on the two estimators considered in this paper, SCOSE and FCOSE. 3. We prove that SCOSE is (essentially, up to lower order terms) the optimal/oracle linear shrinkage estimator with respect to quadratic risk (see Section 5). However, we observe that FCOSE typically achieves higher power than SCOSE. This indicates that it may be useful to search for the optimal estimator in a larger class than linearly shrunk estimators, and also that quadratic loss may not be the right loss function for the purposes of test power. 3

4 The rest of this paper is organized as follows. Section introduces SCOSE, FCOSE and their corresponding shrunk test statistics. Section 3 presents illuminating experiments that bring out the statistically significant improvement in power over HSIC. Section 4 conducts a deeper investigation into the effect of shrinkage and proves the oracle optimality of SCOSE under quadratic risk. Shrunk Estimators and Test Statistics Let HS(H k, H l ) represent the set of Hilbert-Schmidt operators from H k to H l. We first note that S XY can be written as the solution to the following optimization problem. n S XY := min φ(x i ) Z HS(H k,h l ) n ψ(y i ) Z HS Using this idea Muandet et al. (4) suggest the following two shrunk/regularized estimators. From SCOSE to HSIC S This is derived in Muandet et al. (4) by solving n min φ(x i ) Z HS(H k,h l ) n ψ(y i ) Z + λ Z HS HS i= and the optimal solution (called SCOSE) is S S XY := i= ( λ ) S XY + λ where λ (and hence the shrinkage intensity) is estimated by leave-one-out cross-validation (LOOCV), in closed form as ( ) λ ρ S CV := + λ CV n K i= ii Lii n K ] n i,j= ij Lij [ n = (n ) n K n i,j= ij Lij + n K n i= ii Lii Observing the expression for λ CV in Muandet et al. (4), the denominator can be negative (for example, with the Gaussian kernel for small bandwidths, resulting in a kernel matrix close to the identity). This can cause λ CV to be negative, and ρ S to be (unintentionally) outside the range [, ]. Though not discussed in Muandet et al. (4), we shall follow the convention ( that when ) ρ S <, we shall use ρ S = and if ρ S >, we use ρ S =. Indeed, one can show that λ +λ S XY dominates λ +λ ( ) + S XY where (x) + = max{x, }. In Section 4, we prove that SXY S is (essentially) the optimal/oracle linear shrinkage estimator with respect to quadratic risk. We can now calculate the corresponding shrunk statistic HSIC S = SXY S HS = n n i= K ii Lii HSIC (n )HSIC + n K HSIC () n i= ii Lii While the above expression looks daunting, one thing to note is that the amount that HSIC is shrunk (i.e. the multiplicative factor) depends on the value of HSIC. As we shall see in section 4, small HSIC values get shrunk to zero, but as can be seen above, the shrinkage of HSIC is nonmonotonic. n + 4

5 From FCOSE to HSIC F The Flexible Covariance Shrinkage Estimator is derived by relying on the Representer theorem, see Scholkopf & Smola (), to instead minimize n n φ(x i ) ψ(y n β i i ) n φ(x i ) ψ(y i ) + λ β i= i= over all β R n, and the optimal solution (called FCOSE) is S F XY := n i= β λ i n φ(x i ) ψ(y i ) HS where β λ = ( K L + λi) K L where denotes elementwise (Hadamard) product, is the vector [,,..., ] T, and as before the best λ is determined by LOOCV. The procedure to evaluate the optimal λ efficiently is described by Muandet et al. (4) - a single eigenvalue decomposition of K L costing O(n 3 ) can be done, following which evaluating LOOCV is only O(n ) per λ, see Muandet et al. (4), section 3. for more details. As before, after picking the λ by LOOCV, we can derive the corresponding shrunk test statistic as HSIC F = S S XY HS = n tr(m(m + λi) M(M + λi) M) where M = K L. Note here that the shrinkage is not linear, and the effect on HSIC cannot be seen immediately. Similar to SCOSE, we shall see in section 4, small HSIC values get shrunk to zero (LOOCV chooses a large λ). 3 Linear Shrinkage and Quadratic Risk In this section, we prove that SCOSE is (essentially) optimal within a particular class of estimators. Such oracle arguments also exist elsewhere in the literature, like Ledoit & Wolf (4), so we provide only a brief proof outline. Proposition. The oracle (with respect to quadratic risk) linear shrinkage estimator and intensity is defined as S, ρ := and is given by S := ( ρ )S XY where arg min Z Σ XY HS Z HS,Z=( ρ)s XY, ρ ρ := E S XY Σ XY HS E S XY Proof. Define α = Σ XY HS, β = E S XY Σ XY HS, δ = E S XY. Since E[S XY ] = Σ XY, it is easy to verify that α + β = δ. Substituting and expanding the objective, we get: E Z Σ XY HS = E ρs XY + (S XY Σ XY ) HS = ρ δ + β ρ(δ α ) = ρ α + ( ρ) β Differentiating and equating to zero, gives ρ = β δ. 5

6 This ρ appears in terms of quantities that depend on the unknown underlying distribution (hence the term oracle estimator). We use plugin estimates b, d for β, δ. Let d = S XY HS = n K n i,j= ij Lij = HSIC. Since β is the variance of S XY, let b be the sample variance of S XY, i.e. b = n n n k= φ(x i ) ψ xi S XY = n into S and simplifying, we see that HSIC := S HS is ( HSIC n n i= = K ii Lii HSIC nhsic [ n n i= K ii Lii n n i,j= K ij Lij ]. Plugging these ) HSIC (3) Comparing Eq.(3) with Eq.() shows that SCOSE is essentially S, up to a factor in the denominator which is of the same order as the bias of the HSIC empirical estimator (see Theorem in Gretton, Bousquet, Smola & Schölkopf (5)). In other words, SCOSE just corresponds to using a slightly different estimator for δ than the simple plugin d, which varies on the same order as the bias δ Ed. Hence SCOSE, as estimated via regularization and LOOCV, is (essentially) the optimal linear shrinkage estimator under quadratic risk. To the best of our knowledge, this is the first such characterization of optimality of an estimator achieved through leave-one-out cross-validation. We are only able to prove this because one can explicitly calculate both the oracle linear shrinkage intensity ρ as well as the optimal λ CV (as mentioned in Section ). This raises a natural open question can we find other situations where the LOOCV estimator is optimal with respect to some risk measure? (perhaps when explicit calculations are not possible, like ridge regression). 4 Experiments In this section, we run three kinds of experiments: a) to verify that SCOSE has better quadratic risk than FCOSE and original sample estimator, b) detailed synthetic experiments to verify that shrinkage does improve power, across interesting regimes of α = {.,.5,.}, and c) real data obtained from MNIST, to show that we shrinkage detect dependence at much lower samples than the original data size. 4. Quadratic Risk Figure shows that SCOSE is indeed much better than both S XY and FCOSE with respect to quadratic risk. Here, we calculate E Z Σ XY HS for the distribution given in dataset (A) for Z {S XY, SXY S, SF XY }. The expectation is calculated by repeating the experiment times. Each time Z is calculated according to N {, 5, } samples and Σ XY is approximated by the empirical cross-covariance matrix on 5, samples. The four panels use four different kernels which are linear, polynomial, Laplace and Gaussian from top to bottom. The shrunk estimators are always better than the unshrunk, with a larger difference between SCOSE and FCOSE for finitedimensional feature spaces (top two). In infinite-dimensional feature spaces (bottom two), SCOSE and FCOSE are much better than the unshrunk estimator but very similar to each other. The differences between all estimators decreases with increasing n, since the sample cross-covariance operator itself becomes very accurate. 4. Synthetic Data We perform synthetic experiments in a wide variety of settings to demonstrate that the shrunk test statistics achieve higher power than HSIC in a variety of settings. We follow the schema provided in HSIC and HSIC HSIC/n C/n both converge to population HSIC at same rate determined by the dominant term (HSIC). 6

7 quadratic risk quadratic risk.5. 5 N.. 5 N quadratic risk quadratic risk.. 5 N.. 5 N Figure : All panels show quadratic risk E X Σ XY HS for X {S XY, SXY S, SF XY }. Dataset (A) was used in all four panels, but the kernels were varied - from top to bottom is the linear, quadratic, Gaussian and Laplace kernel. the introduction for independence testing and calculating power. We only consider difficult distributions with nonlinear dependence between X, Y, on which linear methods like correlation are shown to fail to detect dependence (some of them were used in previous papers on independence testing like Gretton et al. (7) and Chwialkowski & Gretton (4)). For all experiments, α {.,.5,.} is chosen as the type- error (for choosing the threshold level of the null distribution s right tail). For every setting of parameters of each experiment, power is calculated as the percentage of rejection over repetitions (independent trials), with permutations per repetition (permutation testing to find the null distribution threshold at level α). We use the Gaussian kernel where the bandwidth is chosen by the common median heuristic Scholkopf & Smola (). Table is a representative sample from what we saw on other examples - either large, small or no improvement in power was seen but almost never a worsening of power. The improvements in power may not always be huge, but they are statistically significant - it is difficult to detect such non-linear dependencies at low sample sizes, so any increase in power can be important in scientific applications. Remark. A more appropriate way than using error bars to assess significance is by the Wilcoxon rank sum test, omitted for lack of space, though it yields more favorable results. 4.3 Real Data We use two real datasets - the first is a good example where shrinkage helps a lot, but in the second it does not help (we show it on purpose). Like the synthetic datasets, for most real datasets it either helps or does not hurt (being very rarely worse; see remark in the discussion section). The first is the Eckerle dataset Eckerle (979) from the NIST Statistical Reference Datasets (NIST StRD) for Nonlinear Regression, data from a NIST study of circular interference transmittance (n=35, Y is transmittance, X is wavelength). A plot of the data in Figure reveals a nonlinear relationship between X, Y (though the correlation is.35 with p-value.84). We subsample the data to see how often we can detect a relationship at %, %, 3% of the original data size, when the false positive level is always controlled at.5. The second is the Aircraft dataset Bowman & Azzalini (4) (n=79, X is log(speed), Y is log(span)). Once again, correlation is low, with a p-value of over.8, and we subsample the data to 5%, %, % of the original data size. 7

8 α =. α =.5 α =. HSIC HSIC S HSIC F HSIC HSIC S HSIC F HSIC HSIC S HSIC F ±.. ±.. ±.. ±.. ±.. ± Table : The first column shows scatterplots of X vs Y (all having dependence between X, Y ). There are 3 sets of 5 columns each - for α =.,.5,. (controlled by running permutations). In eachs set, the first three columns show the power of HSIC, HSIC S, HSIC F (with standard deviation over repetitions below). The fourth column shows when HSIC S is significantly better than HSIC, and the fifth column when HSIC F has significantly higher power than HSIC. A blank means the powers are not significantly better or worse. In the first dataset (A) (top 4) we show how the power varies with increasing n (becomes easier). In the second dataset (B) (second 4) we show how the power varies with rotation (goes from near-independence to clear dependence). In the third dataset (C) (third 4), we demonstrate a case where shrinkage does not help much, which is a circle with a hole. In the last dataset (D) (last 4), we demonstrate a case where HSIC S does as well as HSIC F. We tried many more datasets, these are a few representative samples. 8

9 transmitance log(span) wavelength log(speed) HSIC. HSIC S HSIC F % % 3% HSIC. HSIC S HSIC F 5% % % Figure : Top Row: The left figure shows a plot of wavelength against transmittance. The right figure shows the power of HSIC, HSIC S, HSIC F when the data are subsampled to %, %, 3% (error bars over repetitions). Bottom Row: The left figure shows a plot of log(wingspan) vs log(airspeed). The right figure shows the power of HSIC, HSIC S, HSIC F when the data are subsampled to 5%, %, % (error bars over repetitions). 5 Discussion Why might shrinkage improve power? Let us examine the net effect of using shrunk estimators on the value of HSIC, i.e. let us compare HSIC S and HSIC F to HSIC by computing these over all the repetitions of the permutation testing procedure described in the introduction. In Fig. 3, both estimators are visually similar in transforming the actual test statistic. Perhaps the more interesting phenomenon is that Fig. 3 is reminiscent of the graph of a soft-thresholding operator ST t (x) = max{, x t}. Intuitively, if the unshrunk HSIC value is small, the shrinkage methods deem it to be noise and it is shrunk to zero. Looking at the X-axis scaling of the top and bottom row, the size of the region that gets shrunk to zero decreases with n - as expected, shrinkage has less effect when S XY has low variance). The shrinkage being non-monotone (more so for n = than n = 5 in Figure 3) is key to achieving an improvement in power. Using the intuition from the above figure, we can finally piece together why shrinkage may yield benefits. A rejection of H occurs when the test statistic stands out in the right tail of its null distribution. Typically, when the alternative is true (this is when rejecting the null improves power) the unshrunk test statistics calculated from the permuted samples is smaller than the unshrunk HSIC calculated on the original sample. However, the effect of shrinking the small statistics towards zero, and setting the smallest ones to zero, is that the unpermuted test statistic under the alternative 9

10 HSIC S.5 HSIC F.5.5 HSIC.5 HSIC HSIC S. HSIC F.. HSIC. HSIC Figure 3: The top row corresponds to n =, and the bottom row has n = 5. The left plots compare HSIC S to HSIC, and the right plots compare HSIC F to HSIC. Each cross mark corresponds to the shrunk and unshrunk HSIC calculated during a single permutation of a permutation test. distribution stands out more in the right tail of the null. In other words, relative to the unshrunk null distribution and the unshrunk test statistic, the tail of the null distribution is shrunk more towards zero than the unpermuted test statistic, causing the latter to have a higher quantile in the right tail of the former (relative to the quantile before shrinkage). Let us verify this experimentally. In Fig.4 we plot for each of the datasets in Table, the average ratio of unpermuted statistic T to the 95th percentile of the permuted statistics, for T {HSIC, HSIC S, HSIC F }. Recall that for dataset (C), we didn t see much of an improvement in power, but for (A),(B),(D) it is clear from Fig. 4 that the unpermuted statistic is shrunk less than its null distribution s 95th quantile. Remark. In our experiments, real and synthetic, shrinkage usually improves (and almost never worsens) power in false-positive regimes that we usually care about. Will shrinkage always improve power? Possibly not. Even though shrunk the shrunk S XY dominates S XY for estimation error, it may not be the case that shrunk HSIC always dominates unshrunk HSIC for test power (i.e. the latter may not be inadmissible). However, just as no single classifier always outperforms another, it is still beneficial to add techniques like shrinkage, that seem to consistently yield benefits in practice, to the practitioner s array of tools. 6 Conclusion We presented evidence for an important phenomenon - using biased but lower variance shrunk estimators of cross-covariance operators can often significantly improve test power of HSIC at small sample sizes. This observation (that shrinkage can improve power) has rarely been made in the statistics and machine learning testing literature. We think the reason is that most test statistics for independence testing cannot be immediately expressed as the norm of an empirical operator, making it less obvious how to apply shrinkage to improve their power at low sample sizes. We also showed the optimality (among linear shrinkage estimators) of SCOSE, but observe that the nonlinear shrinkage of FCOSE usually yields higher power. To the best of our knowledge, there seems to be no current literature showing that the choice made by leave-one-out cross-validation (SCOSE) explicitly leads to an estimator that is optimal in some sense (among linear shrinkage estimators). This may be because it is often not possible to explicitly calculate the form of the

11 unpermuted HSIC 95%HSIC permuted unpermuted HSIC 95%HSIC permuted estimator estimator unpermuted HSIC 95%HSIC permuted 3 unpermuted HSIC 95%HSIC permuted estimator estimator Figure 4: All panels show the ratio of the unpermuted HSIC to the 95th percentile of the null distribution based on HSICs calculated from the permuted data. (see Table ) The top row has datasets (C) with radius., (B) with angle 3 π/3, and the bottom row has (D) with N = 5, (A) with N = 4. These observations were qualitatively the same in all other synthetic data parameter settings, and also for other percentiles than 95th, and since the figures look identical in spirit, they were omitted due to lack of space. LOOCV estimator, nor the explicit form of the best linear shrinkage estimator, as can both be done in this simple setting. Since even the best possible linear shrinkage estimator (as represented by SCOSE) is usually worse than FCOSE, this result indicates that in order to improve upon FCOSE, it will be necessary to further study the class of non-linear shrinkage estimators for our infinite dimensional operators, as done for finite dimensional covariance matrices in Ledoit & Wolf () and other papers by the same authors. We ended with a brief investigation into the effect of shrinkage on HSIC and why shrinkage may intuitively improve power. We think that our work will be important for more powerful nonparametric detection of subtle nonlinear dependencies at low sample sizes, a common problem in scientific applications. Acknowledgments We would like to thank Arthur Gretton and Jessica Chemali for useful feedback on an earlier draft of the paper, Larry Wasserman for pointing out some useful references, and Krikamol Muandet for sharing his code. AR was supported in part by ONR MURI grant N495. References Bowman, A. W. & Azzalini, A. (4), R package sm: nonparametric smoothing methods (version.-5.4), University of Glasgow, UK and Università di Padova, Italia. Chwialkowski, K. & Gretton, A. (4), A kernel independence test for random processes, in Proceedings of The 3st International Conference on Machine Learning, pp

12 Eckerle, K. (979), Circular interference transmittance study, National Institute of Standards and Technology (NIST), US Department of Commerce, USA. Gretton, A., Borgwardt, K. M., Rasch, M., Schölkopf, B. & Smola, A. J. (6), A kernel method for the two-sample-problem, in Neural Information Processing Systems, pp Gretton, A., Bousquet, O., Smola, A. & Schölkopf, B. (5), Measuring statistical dependence with Hilbert-Schmidt norms, in Algorithmic Learning Theory, Springer, pp Gretton, A., Fukumizu, K., Teo, C. H., Song, L., Schölkopf, B. & Smola, A. J. (7), A kernel statistical test of independence, Neural Information Processing Systems. Gretton, A., Herbrich, R., Smola, A., Bousquet, O. & Schölkopf, B. (5), Kernel methods for measuring independence, The Journal of Machine Learning Research 6, Harchaoui, Z., Bach, F. & Moulines, E. (7), Testing for homogeneity with kernel fisher discriminant analysis, Arxiv preprint Ledoit, O. & Wolf, M. (4), A well-conditioned estimator for large-dimensional covariance matrices, Journal of Multivariate Analysis 88(), Ledoit, O. & Wolf, M. (), Nonlinear shrinkage estimation of large-dimensional covariance matrices, Institute for Empirical Research in Economics University of Zurich Working Paper (55). Muandet, K., Fukumizu, K., Sriperumbudur, B., Gretton, A. & Schoelkopf, B. (4), Kernel mean estimation and stein effect, in Proceedings of The 3st International Conference on Machine Learning, pp. 8. Scholkopf, B. & Smola, A. (), Learning with kernels, MIT press Cambridge. Spirtes, P., Glymour, C. N. & Scheines, R. (), Causation, prediction, and search, Vol. 8, MIT press. Stein, C. (956), Inadmissibility of the usual estimator for the mean of a multivariate normal distribution, Proceedings of the Third Berkeley symposium on mathematical statistics and probability (399), 97 6.