ON ALMOST LINEARITY OF LOW-DIMENSIONAL PROJECTIONS FROM HIGH-DIMENSIONAL DATA

This paper studies the shapes of low dimensional projections from high dimensional data. After standardization, let x be a p-dimensional random variable with mean zero and identity covariance. For a projection beta'x, \\beta\\ = 1, find another direction b so that the regression curve of b'x against beta'x is as nonlinear as possible. We show that when the dimension of x is large, for most directions beta even the most nonlinear regression is still nearly linear. Our method depends on the construction of a pair of p-dimensional random variables, w1, w2, called the rotational twin, and its density function with respect to the standard normal density. With this, we are able to obtain closed form expressions for measuring deviation from normality and deviation from linearity in a suitable sense of average. As an interesting by-product, from a given set of data we can find simple unbiased estimates of E(f(beta'x)(t)/phi1(t)-1)2 and E[(\\E(x\beta, beta'x = t)\\2-t2)f(beta'x)2(t)/phi1(2)t)], where phi1 is the standard normal density, f(beta'x) is the density for beta'x and the ''E'' is taken with respect to the uniformly distributed beta. This is achieved without any smoothing and without resorting to any laborious projection procedures such as grand tours. Our result is related to the work of Diaconis and Freedman. The impact of our result on several fronts of data analysis is discussed. For example, it helps establish the validity of regression analysis when the link function of the regression model may be grossly wrong. A further generalization, which replaces beta'x by B'x with B = (beta1,...,beta(k)) for k randomly selected orthonormal vectors (beta(i), i = 1,...,k), helps broaden the scope of application of sliced inverse regression (SIR).