Title: RECOVERING THE LOWEST LAYER OF DEEP NET-WORKS WITH HIGH THRESHOLD ACTIVATIONS Anonymous

Abstract: Giving provable guarantees for learning neural networks is a core challenge of machine learning theory. Most prior work gives parameter recovery guarantees for one hidden layer networks, however, the networks used in practice have multiple non-linear layers. In this work, we show how we can strengthen such results to deeper networks -we address the problem of uncovering the lowest layer in a deep neural network under the assumption that the lowest layer uses a high threshold before applying the activation, the upper network can be modeled as a well-behaved polynomial and the input distribution is gaussian.

1. INTRODUCTION:

Understanding the landscape of learning neural networks has been a major challege in machine learning. Various works gives parameter recovery guarantees for simple one-hidden-layer networks where the hidden layer applies a non-linear activation u after transforming the input x by a matrix W, and the upper layer is the weighted sum operator: thus f (x) = a i u(w T i x). However, the networks used in practice have multiple non-linear layers and it is not clear how to extend these known techniques to deeper networks.

We consider a multilayer neural network with the first layer activation u and the layers above represented by an unknown polynomial P such that it has non-zero non-linear components. More precisely, the function f computed by the neural network is as follows:

f W (x) = P (u(w T 1 x), u(w T 2 x), . . . , u(w T d x)) for P (X 1 , . . . , X d ) =

r∈Z d + c r • j X rj j .

We assume that the input x is generated from the standard Gaussian distribution and there is an underlying true network (parameterized by some unknown W * )1 from which the labels are generated.

In this work we strengthen previous results for one hidden layer networks to a larger class of functions representing the transform made by the upper layer functions if the lowest layer uses a high threshold (high bias term) before applying the activation: u(a -t) instead of u(a). Intuitively, a high threshold is looking for a high correlation of the input a with a direction w * i . Thus even if the function f is applying a complex transform after the first layer, the identity of these high threshold directions may be preserved in the training data generated using f . Learning with linear terms in P . Suppose P has a linear component then we show that increasing the threshold t in the lowest layer is equivalent to amplifying the coefficients of the linear part. Instead of dealing with the polynomial P it turns out that we can roughly think of it as P (µX 1 , ..., µX d ) where µ decreases exponentially in t (µ ≈ e -t 2

). As µ decreases it has the effect of diminishing the non-linear terms more strongly so that relatively the linear terms stand out. Taking advantage of this effect we manage to show that if t exceeds a certain threshold the non linear terms drop in value enough so that the directions w i can be learned by relatively simple methods. We show that we can get close to the w i applying a simple variant of PCA. While an application of PCA can be thought of as finding principal directions as the local maxima of max ||z||=1 E[f (x)(z T x) 2 ],

we instead perform max E[f (x)H2(z T x)2 ]=1 E[f (x)H 4 (z T x) 4 ]] 2 . If W * has a constant condition number then the local maxima can be used to recover directions that are transforms of w i . Theorem 1 (informal version of Claim 2, Theorem 11). If t > c √ log d for large enough constant c > 0 and P has linear terms with absolute value of coefficients at least 1/poly(d) and all coefficients at most O(1), we can recover the weight vector w i within error 1/poly(d) in time poly(d).

These approximations of w i obtained collectively can be further refined by looking at directions along which there is a high gradient in f ; for monotone functions we show how in this way we can recover w i exactly (or within any desired precision. Theorem 2. (informal version of Theorem 5) Under the conditions of the previous theorem, for monotone P , there exists a procedure to refine the angle to precision in time poly(1/ , d) starting from an estimate that is 1/poly(d) close.

The above mentioned theorems hold for u being sign and ReLU. 3When P is monotone and u is the sign function, learning W is equivalent to learning a union of half spaces. We learn W * by learning sign of P which is exactly the union of halfspaces w T i x = t. Thus our algorithm can also be viewed as a polynomial time algorithm for learning a union of large number of half spaces that are far from the origin -to our knowledge this is the first polynomial time algorithm for this problem but with this extra requirement (see earlier work Vempala (2010) for an exponential time algorithm). Refer to Appendix B.6 for more details.

Such linear components in P may easily be present: consider for example the case where P (X) = u(v T X -b) where u is say the sigmoid or the logloss function. The taylor series of such functions has a linear component -note that since the linear term in the taylor expansion of u(x) has coefficient u (0), for expansion of u(x-b) it will be u (-b) which is Θ(e -b ) in the case of sigmoid. In fact one may even have a tower (deep network) or such sigmoid/logloss layers and the linear components will still be present -unless they are made to cancel out precisely; however, the coefficients will drop exponentially in the depth of the networks and the threshold b.

Sample complexity with low thresholds and no explicit linear terms. Even if the threshold is not large or P is not monotone, we show that W * can be learned with a polynomial sample complexity (although possibly exponential time complexity) by finding directions that maximize the gradient of f . Theorem 3 (informal version of Corollary 1). If u is the sign function and w i 's are orthogonal then in poly(1/ , d) samples one can determine W * within precision if the coefficient of the linear terms in P (µ(X 1 + 1), µ(X 2 + 1), µ(X 3 + 1), . . .) is least 1/poly(d)

Learning without explicit linear terms. We further provide evidence that P may not even need to have the linear terms -under some restricted cases (section 4), we show how such linear terms may implicitly arise even though they may be entirely apparently absent. For instance consider the case when P = X i X j that does not have any linear terms. Under certain additional assumptions we show that one can recover w i as long as the polynomial P (µ(X 1 + 1), µ(X 2 + 1), µ(X 3 + 1), ..) (where µ is e -t has linear terms components larger than the coefficients of the other terms). Note that this transform when applied to P automatically introduces linear terms. Note that as the threshold increases applying this transform on P has the effect of gathering linear components from all the different monomials in P and penalizing the higher degree monomials. We show that if W * is a sparse binary matrix then we can recover W * when activation u(a) = e ρa under certain assumptions about the structure of P . When we assume the coefficients are positive then these results extend for binary low l 1 -norm vectors without any threshold. Lastly, we show that for even activations (∀a, u(a) = u(-a)) under orthogonal weights, we can recover the weights with no threshold.

Learning with high thresholds at deeper layers. We also point out how such high threshold layers could potentially facilitate learning at any depth, not just at the lowest layer. If there is any cut in the network that takes inputs X 1 , . . . , X d and if the upper layers operations can be modelled by a polynomial P , then assuming the inputs X i have some degree of independence we could use this to modularly learn the lower and upper parts of the network separately (Appendix E) Related Work. Various works have attempted to understand the learnability of simple neural networks. Despite known hardness results Goel et al. (2016) ; Brutzkus & Globerson (2017) , there has been an array of positive results under various distributional assumptions on the input and the underlying noise in the label. Most of these works have focused on analyzing one hidden layer neural networks. A line of research has focused on understanding the dynamics of gradient descent on these networks for recovering the underlying parameters under gaussian input distribution Du et al. (2017b; a) ; Li & Yuan (2017) ; Zhong et al. (2017a) ; Zhang et al. (2017) ; Zhong et al. (2017b) . Another line of research borrows ideas from kernel methods and polynomial approximations to approximate the neural network by a linear function in a high dimensional space and subsequently learning the same Zhang et al. (2015) ; Goel et al. (2016) ; Goel & Klivans (2017b; a) . Tensor decomposition methods Anandkumar & Ge (2016); Janzamin et al. (2015) have also been applied to learning these simple architectures.

The complexity of recovering arises from the highly non-convex nature of the loss function to be optimized. The main result we extend in this work is by Ge et al. (2017) . They learn the neural network by designing a loss function that allows a "well-behaved" landscape for optimization avoiding the complexity. However, much like most other results, it is unclear how to extend to deeper networks. The only known result for networks with more than one hidden layer is by Goel & Klivans (2017b) . Combining kernel methods with isotonic regression, they show that they can provably learn networks with sigmoids in the first hidden layer and a single unit in the second hidden layer in polynomial time. We however model the above layer as a multivariate polynomial allowing for larger representation. Another work Arora et al. (2014) deals with learning a deep generative network when several random examples are generated in an unsupervised setting. By looking at correlations between input coordinates they are able to recover the network layer by layer. We use some of their ideas in section 4 when W is a sparse binary matrix.

Notation. We denote vectors and matrices in bold face. || • || p denotes the l p -norm of a vector. || • || without subscript implies the l 2 -norm. For matrices || • || denotes the spectral norm and || • || F denotes the forbenius norm. N (0, Σ) denotes the multivariate gausssian distribution with mean 0 and covariance Σ. For a scalar x we will use φ(x) to denote the p.d.f. of the univariate standard normal distribution with mean zero and variance 1 .For a vector x we will use φ(x) to denote the p.d.f. of the multivariate standard normal distribution with mean zero and variance 1 in each direction. Φ denotes the c.d.f. of the standard gausssian distribution. Also define Φ c = 1 -Φ. Let h i denote the ith normalized Hermite polynomial Wikipedia contributors (2018). For a function f , let fi denote the ith coefficient in the hermite expansion of f , that is, fi = E g∼N (0,1) [f (g)h i (g)].

For a given function f computed by the neural network, we assume that the training samples (x, y) are such that x ∈ R n is distributed according to N (0, 1) and label has no noise, that is, y = f (x).

Note: Most proofs are deferred to the Appendix due to lack of space.

2. APPROXIMATE RECOVERY WITH LINEAR TERM:

In this section we consider the case when P has a positive linear component and we wish to recover the parameters of true parameters W * . The algorithm has two-steps: 1) uses existing one-hidden layer learning algorithm (SGD on carefully designed loss Ge et al. (2017) ) to recover an approximate solution , 2) refine the approximate solution by performing local search (for monotone P ). The intuition behind the first step is that high thresholds enable P to in expectation be approximately close to a one-hidden-layer network which allows us to transfer algorithms with approximate guarantees. Secondly, with the approximate solutions as starting points, we can evaluate the closeness of the estimate of each weight vector to the true weight vector using simple correlations. The intuition of this step is to correlate with a function that is large only in the direction of the true weight vectors. This equips us with a way to design a local search based algorithm to refine the estimate to small error.

For simplicity in this section we will work with P where the highest degree in any X i is 1. The degree of the overall polynomial can still be n. See Appendix B.8 for the extension to general P . More formally, Assumption 1 (Structure of network). We assume that P has the following structure

P (X 1 , . . . , X k ) = c 0 + i∈[d] c i X i + S⊆[d]:|S|>1 c S j∈S X j such that c i = Θ(1) 4 for all i ∈ [d] and for all S ⊆ [d] such that |S| > 1, |c S | ≤ O(1). W * has constant condition number. Thus f (x) = c 0 + i∈[d] c i u((w * i ) T x) + S⊆[d]:|S|>1 c S j∈S u((w * j ) T x). Denote f lin (x) = c 0 + i∈[d] c i u((w * i ) T x

) to be the linear part of f . Next we will upper bound expected value of u(x): for "high-threshold" ReLU, that is, u t (a) = max(0, a -t), E g∼N (0,σ 2 ) [u t (g)] is bounded by a function ρ(t, σ) ≈ e -t 2 2σ 2 (see Lemma 10). We also get a lower bound on |û 4 | in terms of ρ(t, σ)5 This enables us to make the following assumption. Assumption 2. Activation function u is a positive high threshold activation with threshold t, that is, the bias term is t. E g∼N (0,σ 2 ) [u t (g)] ≤ ρ(t, σ) where ρ is a positive decreasing function of t. Also,

|û k | = t Θ(1) ρ(t, 1) for k = 2, 4. Assumption 3 (Value of t). t is large enough such that ρ(t, ||W * ||) ≈ d -η and ρ(t, 1) ≈ d -pη

with for large enough constant η > 0 and p ∈ (0, 1].

For example, for high threshold ReLU, ρ(t, 1) = e -t 2 /2 and µ = ρ(t,

||W * ||) = e -t 2 /2||W * || 2

, thus t = √ 2η log d for large enough d suffices to get the above assumption (κ(W * ) is a constant).

These high-threshold activation are useful for learning as in expectation, they ensure that f is close to f lin since the product terms have low expected value. This is made clear by the following lemmas: Lemma 1. For |S| > 1, under Assumption 2 we have,

E   j∈S u t ((w * j ) T x)   ≤ ρ(t, 1) (κ(W * )ρ(t, ||W * ||)) |S|-1 . So if µ := κ(W * )ρ(t, ||W * ||), then E[ j∈S X j [x]] ≤ ρ(t, 1)µ |S|-1

Lemma 2. Let ∆(x) = f (x) -f lin (x). Under Assumptions 1, 2 and 3, if t is such that dρ(t, ||W * ||) ≤ c for some small enough constant c > 0 we have,

E[|∆(x)|] ≤ O d 3 ρ(t, 1)ρ(t, ||W * ||) = O d -(1+p)η+3 .

Note: We should point out that f (x) and f lin (x) are very different point wise; they are just close in expectation under the distribution of x. In fact, if d is some constant then even the difference in expectation is some small constant.

This closeness suggests that algorithms for recovering under the labels from f lin can be used to recover with labels from f approximately.

Learning One Layer Neural Networks using Landscape Design. Ge et al. (2017) proposed an algorithm for learning one-hidden-layer networks. Intuitively, the approach of Ge et al. ( 2017) is to design a well behaved loss function based on correlations to recover the underlying weight vectors. They show that the local minima of the following optimization corresponds to some transform of each of the w * i -thus it can be used to recover a transform of w * i , one at a time. max

z:E[f lin (x)H2(z T x)]=û2 sgn(û 4 )E[f lin (x)H 4 (z T x)]

which they optimize using the Lagrangian formulation (viewed as a minimization):

min z G lin (z) := -sgn(û 4 )E[f lin (x)H 4 (z T x)] + λ(E[f lin (x)H 2 (z T x)] -û2 ) 2 where H 2 (z T x) = ||z|| 2 h 2 z T x ||z|| = (z T x) 2 √ 2 -||z|| 2 √ 2 and H 4 (z T x) = ||z|| 4 h 4 z T x ||z|| = √ 6 (z T x) 4 12 -||z|| 2 (z T x) 2 2 + ||z|| 4 4

(see Appendix A.1 for more details). Using properties of Hermite polynomials, we have

E[f lin (x)H 2 (z T x)] = û2 i c i (z T w * i ) 2 and similarly E[f lin (x)H 4 (z T x)] = û4 i (z T w * i ) 4 . Thus G lin (z) = -|û 4 | i c i (z T w * i ) 4 + λû 2 2 i c i (z T w * i ) 2 -1 2 .

Using results from Ge et al. (2017) , it can be shown that the approximate local minima of this problem are close to columns of (TW * ) -1 where T is a diagonal matrix with

T ii = √ c i . Definition 1 (( , τ )-local minimum/maximum). z is an ( , τ )-local minimum of F if ||∇F (z)|| ≤ and λ min (∇ 2 F (z)) ≤ τ .

Claim 1 (Ge et al. (2017) ). An ( , τ )-local minima of the Lagrangian formulation z with Ge et al. (2017) do not mention û2 but it is necessary in the non-orthogonal weight vectors case for the correct reduction. Since for us, this value can be small, we mention the dependence.Note that these are not exactly the directions w * i that we need, one way to think about is that we can get the correct directions by estimating all columns and then inverting.

≤ O τ 3 /|û 4 | is such that for an index i |z T w i | = 1 ± O( /λû 2 2 ) ± O(dτ /|û 4 |) and ∀j = i, |v T w j | = O( τ /|û 4 |) where w i are columns of (TW * ) -1 .

One-hidden-layer to Deep Neural Network. Consider the loss with f instead of f lin :

min z : G(z) = -sgn(û 4 )E[f (x)H 4 (z T x)] + λ(E[f (x)H 2 (z T x)] -û2 ) 2

We previously showed that f is close to f lin in expectation due to the high threshold property. This also implies that G lin and G are close and so are the gradients and (eignevalues of) hessians of the same. This closeness implies that the landscape properties of one approximately transfers to the other function. More formally, Theorem 4. Let Z be an ( , τ )-local minimum of function A. If ||∇(B -A)(Z)|| ≤ ρ and ||∇ 2 (B -A)(Z)|| ≤ γ then Z is an ( + ρ, τ + γ)-local minimum of function B and vice-versa.

We will now apply above lemma on our G lin (z) and G(z). Claim 2. For λ = Θ(|û 4 |/û 2 2 ) ≈ d η , an ( , τ )-approximate local minima of G (for small enough

, τ ≤ d -2η ) is an (O(log d)d -(1+p)η+3 , O(log d)d -(1+p)η+3 )-approximate local minima of G lin . This implies z is such that for an index i, |z T w i | = 1 ± O(1)d -2/3pη+3 and ∀j = i, |z T w j | = O(1)d -1/3pη+3/2

where w i are columns of (TW * ) -1 (ignoring log d factors).

Note: For ReLU, setting t = √ C log d for large enough C > 0 we can get closeness 1/poly(d) to the columns of (TW * ) -1 . Refer Appendix B.7 for details for sigmoid.

The paper Ge et al. (2017) also provides an alternate optimization that when minimized simultaneously recovers the entire matrix W * instead of having to learn columns of (TW * ) -1 separately. We show how applying our methods can also be applied to that optimization in Appendix B.4 to recover W * by optimizing a single objective.

2.1. APPROXIMATE TO ARBITRARILY CLOSE FOR MONOTONE P:

Assuming P is monotone, we can show that the approximate solution from the previous analysis can be refined to arbitrarily closeness using a random search method followed by approximately finding the angle of our current estimate to the true direction.

The idea at a high level is to correlate with δ (z T x -t) where δ is the Dirac delta function. It turns out that the correlation is maximized when z is equal to one of the w i . Correlation with δ (z T x-t) is checking how fast the correlation of f with δ(z T x-t) is changing as you change t. To understand this look at the case when our activation u is the sign function then note that correlation of u t (w T x -t) with δ (w T x -t) is very high as its correlation with δ(w T x -t ) is 0 when t < t and significant when t > t. So as we change t' slightly from tto t + there is a sudden increase. If z and w differ then it can be shown that correlation of u t (w T x -t) with δ (z T x -t) essentially depends on cot(α) where α is the angle between w and z (for a quick intuition note that one can prove that E[u t (w T x)δ (z T x)] = c cot(α). See Lemma 16 in Appendix). In the next section we will show how the same ideas work for non-monotone P even if it may not have any linear terms but we only manage to prove polynomial sample complexity for finding w instead of polynomial time complexity.

In this section we will not correlate exactly with δ (z T x -t) but instead we will use this high level idea to estimate how fast the correlation with δ(z T x -t ) changes between two specific values as one changes t , to get an estimate for cot(α). Secondly since we can't to a smooth optimization over z, we will do a local search by using a random perturbation and iteratively check if the correlation has increased. We can assume that the polynomial P doesn't have a constant term c 0 as otherwise it can easily be determined and cancelled out6 .

We will refine the weights one by one. WLOG, let us assume that w * 1 = e 1 and we have z such that

z T w * 1 = z 1 = cos -1 (α 1 ). Let l(z, t, ) denote {x : z T x ∈ [t -, t]} for z ∈ S n-1 .

Algorithm 1 RefineEstimate 1: Run EstimateT anAlpha on z to get s = tan(α) where α is the angle between z and w * 1 . 2: Perturb current estimate z by a vector along the d -1 dimensional hyperplane normal to z with the distribution n(0, Θ(α/d)) d-1 to get z . 3: Run EstimateT anAlpha on z to get s = tan(α ) where α is the angle between z and w * 1 .

4: if α ≤ O(α/d) then 5: z ← z 6: Repeat till α ≤ . Algorithm 2 EstimateTanAlpha 1: Find t 1 and t 2 such that P r[sgn(f (x))|x ∈ l(z, t , )] at t 1 is 0.4 and at t 2 is 0.6. 2: Return t2-t1 Φ -1 (0.6)-Φ -1 (0.4) .

The algorithm (Algorithm 1) estimates the angle of the current estimate with the true vector and then subsequently perturbs the vector to get closer after each successful iteration. Theorem 5. Given a vector z ∈ S d-1 such that it is 1/poly(d)-close to the underlying true vector

w * 1 , that is cos -1 (z T w * 1 ) ≤ 1/poly(d), running Ref ineEstimate for O(T ) iterations outputs a vector z * ∈ S d-1 such that cos -1 ((z * ) T w * 1 ) ≤ 1 -c d T γ for some constant c > 0. Thus after O(d log(1/ )) iterations cos -1 ((z * ) T w * 1 ) ≤ .

We prove the correctness of the algorithm by first showing that EstimateT anAlpha gives a multiplicative approximation to tan(α). The following lemma captures this property. Lemma 3. EstimateT anAlpha(z) outputs y such that y = (1 ± O(η)) tan(α) where α is the angle between z and w * 1 .

Proof. We first show that the given probability when computed with sgn(x T w * 1 -t) is a well defined function of the angle between the current estimate and the true parameter up to multiplicative error. Subsequently we show that the computed probability is close to the one we can estimate using f (x) since the current estimate is close to one direction. The following two lemmas capture these properties.

Lemma 4. For t, t and ≤ 1/t , we have

P r[x T w * 1 ≥ t and x ∈ l(z, t , )|x ∈ l(z, t, )] = Φ c t -t * cos(α 1 ) | sin(α 1 )| ± O( )t Lemma 5. For t ∈ [0, t/ cos(α 1 )], we have P r[sgn(f (x))|x ∈ l(z, t , )] = P r[sgn((w * 1 ) T x -t)|x ∈ l(z, t, )] + de -Ω(t 2 ) .

Using the above, we can show that,

t 2 -t 1 = Φ -1 (0.6 -η 1 ± O( )t 1 ) -Φ -1 (0.4 -η 2 ± O( )t 2 ) tan(α) = Φ -1 (0.6) -Φ -1 (0.4) -(η 1 ± O( )t 1 )(Φ -1 ) (p 1 ) + (η 2 ± O( )t 2 )(Φ -1 ) (p 2 ) tan(α)

where η 1 , η 2 > 0 are the noise due to estimating using f and

p 1 ∈ [0.6 -η 1 ± O( )t 1 , 0.6] and p 2 ∈ [0.4 -η 2 ± O( )t 2 , 0.4] as long as t 1 , t 2 ∈ [0, t/ cos(α 1 )].

The following lemma bounds the range of t 1 and t 2 .

Lemma 6. We have

0 ≤ t 1 ≤ t 2 ≤ t cos(α1) .

Thus, we have,

t 2 -t 1 Φ -1 (0.6) -Φ -1 (0.4) = (1 ± O (η 1 + η 2 + t 2 )) tan(α)

as long as η 2 +O( )t 2 ≤ c for some constant c > 0. Thus, we can get a multiplicative approximation to tan(α) up to error η ( can be chosen to make its contribution smaller than η).

Finally we show (proof in Appendix ??) that with constant probability, a random perturbation reduces the angle by a factor of (1 -1/d) of the current estimate hence the algorithm will halt after O(d log(1/ν)) iterations.

Lemma 7. By applying a random Gaussian perturbation along the d -1 dimensional hyperplane normal to z with the distribution n(0, Θ(α/d)) d-1 and scaling back to the unit sphere, with constant probability, the angle α (< π/2) with the fixed vector decreases by at least Ω(α/d).

3. SAMPLE COMPLEXITY:

We extend the methods of the previous section to a broader class of polynomials but only to obtain results in terms of sample complexity. The main idea as in the previous section is to correlate with δ (z T x-t) (the derivative of the dirac delta function) and find

arg max ||z||2=1 E[f (x)δ (z T x-t)].

We will show that the correlation goes to infinity when z is one of w * i and bounded if it is far from all of them. From a practical standpoint we calculate δ (z T x -s) by measuring correlation with 1 2 (δ(z T x -s + ) -δ(z T x -s -). In the limit as → 0 this becomes δ (z T x -s). δ(z T x -s) in turn is estimated using 1 (sgn(z T x -s + )sgn(z T x -s)), as in the previous section, for an even smaller ; however, for ease of exposition, in this section, we will assume that correlations with δ(z T x -s) can be measured exactly.

Let us recall that:

f (x) = P (u((w * 1 ) T x), u((w * 2 ) T x), . . . , u((w * d ) T x)). Let C 1 (f, z, s) denote E[f (x)δ(z T x -s)] and let C 2 (f, z, s) denote E[f (x)(δ(z T x -s -) -δ(z T x -s + )].

If u = sgn then P has degree at most 1 in each X i . Let ∂P ∂Xi denote the symbolic partial derivative of P with respect to X i ; so, it drops monomials without X i and factors off X i from the remaining ones. Let us separate dependence on X i in P as follows:

P (X 1 , , .., X d ) = X i Q i (X 1 , ..X i-1 , X i+1 , .., X d ) + R 1 (X 1 , .X i-1 , X i+1 , .., X d ) then ∂P ∂Xi = Q i .

We will overload the polynomial P such that P [x] to denote the polynomial computed by substituting X i = u((w * 1 ) T x) and similarly for Q and R. Under this notation f (x) = P [x]. We will also assume that |P (X)| ≤ ||X|| O(1) = ||X|| c1 (say). By using simple correlations we will show: Theorem 6. If u is the sgn function, P (X) ≤ ||X|| c1 and for all i, E[

Q i [x]|(w * i ) T x = t] ≥ 3 then using poly( d 3 2

) samples one can determine the w * i 's within error 2 .7 

Note that if all the w * i 's are orthogonal then X i are independent and

E Q i [x] (w * i ) T x = t is just value of Q i evaluated

by setting X i = 1 and setting all the the remaining X j = µ where µ = E[X j ]. This is same as 1/µ times the coefficient of X i in P (µ(X 1 + 1), . . . , µ(X d + 1)).

Under review as a conference paper at ICLR 2019 Corollary 1. If u is the sgn function and w * i s are orthogonal then in sample complexity poly( d 3 2 ) one can determine W * within error 2 in each entry, if the coefficient of the linear terms in P (µ(X 1 + 1), µ(X 2 + 1), µ(X 3 + 1), ..) is larger than 3 µ, where µ = E[X i ].

The main point behind the proof of Theorem 6 is that the correlation is high when z is along one of w * i and negligible if it is not close to any of them. Lemma 8. Assuming P (X) 1) . Otherwise if all angles α i between z and w * i are at least 2 it is at most d O(1) / 2 .

< ||X|| c1 . If z = w * i then C 2 (f, z, t) = φ(t)E ∂P ∂Xi z T x = t + d O(

We will use the notation g(x) x=s to denote g(x) evaluated at x = s. Thus Cauchy's mean value theorem can be stated as

g(x + ) -g(x) = [g (s)](s = s ∈ [x, x + ]

). We will over load the notation a bit: φ(z T x = s) will denote the probability density that vz T x = s; so if z is a unit vector this is just φ(s); φ(z T 1 x = s 1 , z T 2 x = s 2 ) denotes the probability density that both z T 1 x = s 1 , z T 2 x = s 2 ; so again if z 1 , z 2 are orthonormal then this is just φ(s 1 )φ(s 2 ). The following claim interprets correlation with δ(z T x -s) as the expected value along the corresponding plane

z T x = s. Claim 3. E[f (x)δ(z T x -s)] = E[f (x)|z T x = s]φ(z T x = s).

The following claim computes the correlation of P with δ

(z T x -s). Claim 4. E[P [x]δ (z T x = s)] is equal to i | cot(α i )|φ(z T x = s, (w * i ) T x = t) E ∂P ∂Xi [x]|z T x = s, (w * i ) T x = t + φ (s)E[P [x]|z T x = s].

We use this to show that the correlation is bounded if all the angles are lower bounded.

Claim 5. If P (X) ≤ ||X|| c1 and if z has an angle of at least 2 with all the

w * i 's then C 2 (f, z, s) ≤ d O(1) / 2 .

Above claims can be used to prove main Lemma 8. Refer to the Appendix C for proofs.

Proof of Theorem 6. If we wish to determine w * i within an angle of accuracy 2 let us set to be O( 3 2 φ(t)d -c ). From Lemma 8, for some large enough c, this will ensure that if all α i > 2 the correlation is o(φ(t) 3 ). Otherwise it is φ(t) 3 (1±o(1)). Since φ(t) = poly(1/d), given poly( d 2 3 ) samples, we can test if a given direction is within accuracy 2 of a w * i or not.

4. STRONGER RESULTS UNDER STRUCTURAL ASSUMPTIONS:

Under additional structural assumptions on W * such as the weights being binary, that is, in {0, 1}, sparsity or certain restrictions on activation functions, we can give stronger recovery guarantees.

Proofs have been deferred to Appendix D.

Theorem 7. For activation u t (a) = e ρ(a-t) . Let the weight vectors w * i be 0, 1 vectors that select the coordinates of x. For each i, there are exactly d indices j such that w ij = 1 and the coefficient of the linear terms in P (µ(X 1 + 1), µ(X 2 + 1), µ(X 3 + 1), ..) for µ = e -ρt is larger than the coefficient of all the product terms (constant factor gap) then we can learn the W * .

In order to prove the above, we will construct a correlation graph over x 1 , . . . , x n and subsequently identify cliques in the graph to recover w * i 's. With no threshold, recovery is still possible for disjoint, low l 1 -norm vector. The proof uses simple correlations and shows that the optimization landscape for maximizing these correlations has local maximas being w * i 's. Theorem 8. For activation u(a) = e a . If all w * i ∈ {0, 1} n are disjoint, then we can learn w * i as long as P has all positive coefficients and product terms have degree at most 1 in each variable.

For even activations, it is possible to recover the weight vectors even when the threshold is 0. The technique used is the PCA like optimization using hermite polynomials as in Section 2. Denote C(S, µ)

= S⊆S ⊆[n] c S µ |S | .

Theorem 9. If the activation is even and for every i, j: C({i}, û0 ) + C({j}, û0 ) > 6û 2 2 û0 û4 C({i, j}, û0 ) then there exists an algorithm that can recover the underlying weight vectors.

5. CONCLUSION:

In this work we show how activations in a deep network that have a high threshold make it easier to learn the lowest layer of the network. We show that for a large class of functions that represent the upper layers, the lowest layer can be learned with high precision. Even if the threshold is low we show that the sample complexity is polynomially bounded. An interesting open direction is to apply these methods to learn all layers recursively. It would also be interesting to obtain stronger results if the high thresholds are only present at a higher layer based on the intuition we discussed.

A PREREQUISITES:

A.1 HERMITE POLYNOMIALS Hermite polynomials form a complete orthogonal basis for the gaussian distribution with unit variance. For more details refer to Wikipedia contributors (2018). Let h i be the normalized hermite polynomials. They satisfy the following,

Fact 0. E[h n (x)] = 0 for n > 0 and E[h 0 (x)] = 1. Fact 1. E a∼N (0,1) [h i (a)h j (a)] = δ ij where δ ij = 1 iff i = j.

This can be extended to the following: Fact 2. For a, b with marginal distribution N (0, 1) and correlation ρ,

E[h i (a)h j (b)] = δ ij ρ j .

Consider the following expansion of u into the hermite basis (h i ),

u(a) = ∞ i=0 ûi h i (a). Lemma 9. For unit norm vectors u, v, E[u(v T x)h j (w T x)] = ûj (v T w) j .

Proof. Observe that v T x and w T x have marginal distribution N (0, 1) and correlation v T w. Thus using Fact 2,

E[u(v T x)h j (w T x)] = ∞ i=1 ûi E[h i (v T x)h j (w T x)] = ∞ i=1 ûi δ ij (v T w) j = ûj (v T w) j .

For gaussians with mean 0 and variance σ 2 define weighted hermite polynomials H σ l (a) = |σ| l h l (a/σ). Given input v T x for x ∼ N (0, I), we suppress the superscript σ = ||v||. Corollary 2. For a non-zero vector v (not necessarily unit norm) and a unit norm vector w,

E[H i (v T x)h j (w T x)] = δ ij (v T w) j .

Proof. It follows as the proof of the previous lemma,

E[u(v T x)h j (w T x)] = ∞ i=1 ûi E[h i (v T x)h j (w T x)] = ∞ i=1 ûi δ ij (v T w) j = ûj (v T w) j . Fact 3. h n (x + y) = 2 -n 2 n k=0 n k h n-k (x √ 2)h k (y √ 2). Fact 4. h n (γx) = n 2 k=0 γ n-2k (γ 2 -1) k n 2k (2k)! k! 2 -k h n-2k (x). Fact 5. α(n, m, γ) = E[h m (x)h n (γx)] = γ n-2k (γ 2 -1) k n 2k (2k)! k! 2 -k for k = n-m 2 if k ∈ Z + else 0. A.2 PROPERTIES OF MATRICES Consider matrix A ∈ R m×m . Let σ i (A) to be the ith singular value of A such that σ 1 (A) ≥ σ 2 (A) ≥ . . . ≥ σ m (A) and set κ(A) = σ 1 (A)/σ m (A). Fact 6. |det(A)| = m i=1 σ i (A).

Fact 7. Let B be a (mk) × (mk) principal submatrix of A, then κ(B) ≤ κ(A).

A.3 ACTIVATION FUNCTIONS:

Lemma 10. For u being a high threshold ReLU, that is, u t (a) = max(0, a -t) we have for t ≥ C for large enough constant C > 0, E g∼N (0,σ 2 ) [u t (g)] ≤ e -t 2 2σ 2 . Also, û4 , û2 = t Θ(1) e -t 2 2 .

Proof. We have

E g∼N (0,σ 2 ) [u t (g)] = 1 √ 2πσ ∞ -∞ max(0, g -t)e -g 2 2σ 2 dg = 1 √ 2πσ ∞ t (g -t)e -g 2 2σ 2 dg ≤ 1 √ 2πσ ∞ t ge -g 2 2σ 2 dg = σ √ 2π ∞ t 2 2σ 2 e -h dh = σ √ 2π e -t 2 2σ 2 .

Also,

û4 = E g∼N (0,1) [u t (g)h 4 (g)] = 1 √ 2π ∞ -∞ max(0, g -t)(g 4 -6g 2 + 3)e -g 2 2 dg = 1 √ 2π ∞ t (g -t)(g 4 -6g 2 + 3)e -g 2 2 dg ≥ 1 √ 2π (t 4 -6t 2 ) 1 t e -t 2 2 -1-1 2t 2 ≥ Ω t 3 e -t 2 2 .

To upper bound,

û4 = 1 √ 2π ∞ -∞ max(0, g -t)(g 4 -6g 2 + 3)e -g 2 2 dg = 1 √ 2π ∞ t (g -t)(g 4 -6g 2 + 3)e -g 2 2 dg ≤ 1 √ 2π ∞ t 2g 5 e -g 2 2 dg = 1 √ 2π ∞ t 2 2 h 2 e -h dh = O t 4 e -t 2 2 .

Similar analysis holds for û2 .

Observe that sgn can be bounded very similarly replacing g -t by 1 which can affect the bounds up to only a polynomial in t factor. Lemma 11. For u being a high threshold sgn, that is, u t (a) = sgn(a -t) we have for t ≥ C for

large enough constant C > 0, E g∼N (0,σ 2 ) [u t (g)] ≤ e -t 2 2σ 2 . Also, û4 , û2 = t Θ(1) e -t 2 2 .

For sigmoid, the dependence varies as follows: Lemma 12. For u being a high threshold sigmoid, that is, u t (a) = 1 1+e -(a-t) we have for t ≥ C for large enough constant C > 0, E g∼N (0,σ 2 ) [u t (g)] ≤ e -t+ σ 2 2 . Also, û4 , û2 = Θ(e -t ).

Proof. We have

E g∼N (0,σ 2 ) [u t (g)] = 1 √ 2πσ ∞ -∞ 1 1 + e -(g-t) e -g 2 2σ 2 dg = e -t √ 2πσ ∞ -∞ 1 e -t + e -g e -g 2 2σ 2 dg ≤ e -t √ 2πσ ∞ -∞ e g e -g 2 2σ 2 dg = e -t e σ 2 2 √ 2πσ ∞ -∞ e -(g-σ 2 ) 2 2σ 2 dg = e -t e σ 2 2

Also,

û4 = E g∼N (0,1) [u t (g)h 4 (g)] = 1 √ 2π ∞ -∞ 1 1 + e -(g-t) e -g 2 2 dg = e -t √ 2π ∞ -∞ 1 e -t + e -g (g 4 -6g 2 + 3)e -g 2 2 dg ≥ e -t √ 2π ∞ 0 1 e -t + e -g (g 4 -6g 2 + 3)e -g 2 2 dg ≥ e -t √ 2π ∞ 0 1 2 (g 4 -6g 2 + 3)e -g 2 2 dg

= Ω(e -t ).

We can upper bound similarly and bound û2 .

B APPROXIMATE RECOVERY WITH LINEAR TERMS B.1 CONSTRAINED OPTIMIZATION VIEW OF LANDSCAPE DESIGN:

Let us consider the linear case with w * i 's are orthonormal. Consider the following maximization problem for even l ≥ 4, max

z∈S n-1 sgn(û l ) • E f (x) • H l z T x

where h l is the lth hermite polynomial. Then we have,

sgn(û l ) • E f (x) • h l z T x = sgn(û l ) • E k i=1 c i u t ((w * i ) T x) • h l z T x = sgn(û l ) • k i=1 c i E u t ((w * i ) T x) • h l z T x = |û l | k i=1 c i ((w * i ) T z) l .

It is easy to see that for z ∈ S n-1 , the above is maximized at exactly one of the w i 's (up to sign flip for even l) for l ≥ 3 as long as u l = 0. Thus, each w i is a local minima of the above problem.

Let

L(z) = - k i=1 c i z l i .

For constraint ||z|| 2 = 1, we have the following optimality conditions (see Nocedal & Wright (2006) for more details). This applied to our function gives us that for λ

= -i ciz l i ||z|| 2 (λ < 0), -lc i z l-1 i -2λz i = 0

The above implies that either z i = 0 or z l-2 i = -λ lci with ||z|| 2 = 1. For this to hold z is such that for some set S ⊆ [n], |S| > 1, only i ∈ S have z i = 0 and i∈S z 2 i = 1. This implies that for all i ∈ S, z l-2 i = -2λ lci . Second order: For all w = 0 such that w T z = 0, w T (∇ 2 L(z) -2λI)w ≥ 0.

For our function, we have:

∇ 2 L(z) = -l(l -1)diag(c • z) l-2 =⇒ (∇ 2 L(z)) ij = 2(l -1)λ if i = j and i ∈ S 0 otherwise.

The last follows from using the first order condition. For the second order condition to be satisfied we will show that |S| = 1. Suppose |S| > 2, then choosing w such that w i = 0 for i ∈ S and such that w T z = 0 (it is possible to choose such a value since |S| > 2), we get w T (∇ 2 L(z) -2λI)w = 2(l -2)λ||w|| 2 which is negative since λ < 0, thus these cannot be global minima. However, for |S| = 1, we cannot have such a w, since to satisfy w T z = 0, we need w i = 0 for all i ∈ S, this gives us w T (∇ 2 L(z) -2λI)w = -2λ||w|| 2 which is always positive. Thus z = ±e i are the only local minimas of this problem.

B.2 IMPORTANT RESULTS FROM GE ET AL. (2017):

Lemma 13 (Ge et al. (2017) ). If z is an ( , τ )-local minima of F (z) =i α i z 4 i + λ( i z 2 i -1) 2 for ≤ τ 3 /α min where α min = min i α i , then

• (Lemma 5.2) |z| 2nd ≤ τ αmin where |z| 2nd denotes the magnitude of the second largest entry in terms of magnitude of z.

• (Derived from Proposition 5.7) z max = ±1 ± O(dτ /α min ) ± O( /λ) where |z| max is the value of the largest entry in terms of magnitude of z.

B.3 OMITTED PROOFS FOR ONE-BY-ONE RECOVERY:

Proof of Lemma 1. Let O ∈ R d×d be the orthonormal basis (row-wise) of the subspace spanned by w * i for all i ∈ [d] generated using Gram-schmidt (with the procedure done in order with elements of |S| first). Now let O S ∈ R |S|×d be the matrix corresponding to the first S rows and let O ⊥ S ∈ R (d-|S|)×n be that corresponding to the remaining rows. Note that OW * (W * also has the same ordering) is an upper triangular matrix under this construction.

E   j∈S u t ((w * j ) T x)   = 1 (2π) n/2 x i∈S u t (x T w * i )e -||x|| 2 2 dx = 1 (2π) n/2 x i∈S u t ((O S w * i ) T O S x)e -||O S x|| 2 +||O ⊥ S x|| 2 2 dx = 1 (2π) |S| 2 x ∈R |S| i∈S u t ((O S w * i ) T x )e -||x || 2 2 dx 1 (2π) d-|S| 2 x ∈R d-|S| e -||x || 2 2 dx = 1 (2π) |S| 2 x ∈R |S| i∈S u t ((O S w * i ) T x )e -||x || 2 2 dx = |det(O S W * S )| -1 (2π) |S| 2 b∈R |S| i∈S u t (b i )e -||(O S W * S ) -T b|| 2 2 db

Now observe that O S W * S is also an upper triangular matrix since it is a principal sub-matrix of OW * . Thus using Fact 6 and 7, we get the last equality. Also, the single non-zero entry row has non-zero entry being 1 (||w * i || = 1 for all i). This gives us that the inverse will also have the single non-zero entry row has non-zero entry being 1. WLOG assume index 1 corresponds to this row. Thus we can split this as following

E   j∈S u t ((w * j ) T x)   ≤ |det(O S W * S )| -1 1 √ 2π b1 u t (b 1 )e -b 2 1 2 db 1   i∈S\{1} 1 √ 2π bi u t (b i )e - b 2 i 2||O S W * S || 2 db i   ≤ |det(O S W * S )| -1 1 √ 2π b1 u t (b 1 )e -b 2 1 2 db 1   i∈S\{1} 1 √ 2π bi u t (b i )e -b 2 i ||W * ||2 db i   ≤ ρ(t, 1) (κ(W * )ρ(t, ||W * ||))

|S|-1

Proof of Claim 1. Consider the SVD of matrix M = UDU T . Let W = UD -1/2 and y i = √ c i W T w * i for all i. It is easy to see that y i are orthogonal. Let F (z) = G(Wz):

F (z) = |û 4 | i c i (z T W T w * i ) 4 -λû 2 2 i c i (z T W T w * i ) 2 -1 2 = |û 4 | i 1 c i (z T y i ) 4 -λû 2 2 i (z T y i ) 2 -1 2 .

Since y i are orthogonal, for means of analysis, we can assume that y i = e i , thus the formulation reduces to

max z |û 4 | i 1 ci (z i ) 4 -λ ||z|| 2 -1 2 up to scaling of λ = λû 2 2 .

Note that this is of the form in Lemma 13 hence using that we can show that the approximate local minimas of F (z) are close to y i and thus the local maximas of G(z) are close to

Wy i = √ c i WW T w * i = √ c i M -1 w * i

due to the linear transformation. This can alternately be viewed as the columns of (TW * ) -1 since TW * M -1 (TW * ) T = I.

Proof of Theorem 4. Let Z be an ( , τ )-local minimum of A, then we have ||∇A(Z)|| ≤ and

λ min (∇ 2 A(Z)) ≥ -τ . Observe that ||∇B(Z)|| = ||∇(A + (B -A)(Z)|| ≤ ||∇A(Z)|| + ||∇(B -A)(Z)|| ≤ + ρ.

Also observe that

λ min (∇ 2 B(Z)) = λ min (∇ 2 (A + (B -A))(Z)) ≥ λ min (∇ 2 A(Z)) + λ min (∇ 2 (B -A)(Z)) ≥ -τ -||∇ 2 (B -A)(Z)|| ≥ -τ -γ

Here we use |λ min (M)| ≤ ||M|| for any symmetric matrix. To prove this, we have ||M|| = max x∈S n-1 ||Mx||. We have x = i x i v i where v i are the eigenvectors. Thus we have Mx =

i x i λ i (M)v i and x 2 i = 1. Which gives us that ||M|| = i x 2 i λ 2 i (M) ≥ |λ min (M)|.

Proof of Lemma 2. Expanding f , we have

E[|∆(x)|] = E   S⊆[d]:|S|>1 c S j∈S u t ((w * j ) T x)   ≤ S⊆[d]:|S|>1 |c S |E   j∈S u t ((w * j ) T x)   using Lemma 1 ≤ C S⊆[d]:|S|>1 ρ(t, 1) 1 σ min (W * ) ρ(t, ||W * ||) |S|-1 = C d i=1 d i ρ(t, 1) 1 σ min (W * ) ρ(t, ||W * ||) i-1 using d i ≤ d i ≤ C d i=1 dρ(t, 1) d σ min (W * ) ρ(t, ||W * ||) i-1

using assumption on t ≤ Cd 2 ρ(t, 1)

d σ min (W * ) ρ(t, ||W * ||) Lemma 14. For any function L such that ||L(z, x)|| ≤ C(z)||x|| O(1)

where C is a function that is not dependent on x, we have

||E[∆(x)L(x)]|| ≤ C(z)d -(1+p)η+3 O(log d).

Proof. We have 

||E[∆(x)L(x)]|| ≤ E[|∆(x)||L(x)||] ≤ E[|∆(x)C(z)||x|| O(1) ] = C(z) E[|∆(x)| ||x|| O(1) | ||x|| ≥ c]P r[||x|| ≥ c] + E[|∆(x)| ||x|| O(1) | ||x|| < c]P r[||x|| < c] ≤ C(z)(E[||x|| O(1) |||x|| ≥ c]P r[||x|| ≥ c] + cE[|∆(x)|]) = C(z)(c O(1) e -c 2 2 + c O(1) E[|∆(x)|]).

= Θ(|û 4 |/û 2 2 ) ≈ d η , ||∇G(z)|| ≥ Ω(1)d -η .

Proof. Let K = κ(W * ) which by assumption is θ(1). We will argue that local minima of G cannot have z with large norm. First lets argue this for G lin (z). We know that

G lin (z) = -α (z T w * i ) 4 + λβ 2 (( (z T w * i ) 2 ) -1) 2 where α = |û 4 | and β = û2 . We will argue that z T ∇G lin (z) is large if z is large. z T ∇G lin (z) = -4α (z T w * i ) 3 (z T w * i ) + 2λβ 2 (z T w * i ) 2 -1 2(z T w * i )(z T w * i ) = -4α (z T w * i ) 4 + 4λβ 2 (z T w * i ) 2 -1 (z T w * i ) 2 Let y = W * z then K||z|| ≥ ||y|| ≥ ||z||/K since K is the condition number of W * . Then this implies z T ∇G lin (z) = -4α y 4 i + 4λβ 2 (||y|| 2 -1)||y|| 2 = 4||y|| 2 ((-α + λβ 2 )||y|| 2 + λβ 2 ) ≥ ||y|| 4 (-α + λβ 2 ) ≥ Ω(1)d -η ||y|| 4

Since ||y|| ≥ ||z||/K = Ω(1) by assumptions on λ, z we have z T ∇G lin (z) ≥ Ω(λβ 2 ||y|| 4 ) = Ω(1)d -η ||z|| 4 . This implies ||∇G lin (z)|| = Ω(1)d -η ||z|| 3 .

Now we need to argue for:

G. G(z) -G lin (z) = -sgn(û 4 )E[(f lin (x) + ∆(x))H 4 (z T x)] + λ(E[(f lin (x) + ∆(x))H 2 (z T x)] -β) 2 + sgn(û 4 )E[(f lin (x))H 4 (z T x)] -λE[(f lin (x))H 2 (z T x)] -β] 2 = -sgn(û 4 )E[∆(x)H 4 (z T x)] + λE[∆(x)H 2 (z T x)] 2 + 2λE[∆(x)H 2 (z T x)]E[f lin (x)H 2 (z T x) -β] = -sgn(û 4 )||z|| 4 E[∆(x)h 4 (z T x/||z||)] + λ||z|| 4 E[∆(x)h 2 (z T x/||z||)] 2 + 2λ||z|| 4 E[∆(x)h 2 (z T x/||z||)]E[f lin (x)h 2 (z T x/||z||)] -2λβ||z|| 2 E[∆(x)h 2 (z T x/||z||)] Now h 4 (z T x/||z||) doesn't have a gradient in the direction of z so z T ∇h 4 (z T x/||z||) = 0. Simi- larly z T ∇h 2 (z T x/||z||) = 0. So z T ∇(G(z) -G lin (z)) = -4sgn(û 4 )||z|| 4 E[∆(x)h 4 (z T x/||z||)] + 4λ||z|| 4 (E[∆(x)h 2 (z T x/||z||)]) 2 + 8λ||z|| 4 E[∆(x)h 2 (z T x/||z||)]E[f lin (x)h 2 (z T x/||z||)] -4λβ||z|| 2 E[∆(x)h 2 (z T x/||z||)]

We know that E[f lin (x)h 2 (z T x/||z||)] has a factor of β giving us using Lemma 14:

|z T ∇(G(z) -G lin (z))| ≤ O(log d)d -(1+p)η+3 ||z|| 4 . So z T ∇G(z) is also Ω(||z|| 4 ). so ||∇G(z)|| ≥ Ω(1)d -η Proof of Claim 2. We have G -G lin as follows, G(z) -G lin (z) = -sgn(û 4 )E[(f lin (x) + ∆(x))H 4 (z T x)] + λ(E[(f lin (x) + ∆(x))H 2 (z T x)] -û2 ) 2 + sgn(û 4 )E[(f lin (x))H 4 (z T x)] -λ(E[(f lin (x))H 2 (z T x)] -û2 ) 2 = -sgn(û 4 )E[∆(x)H 4 (z T x)] + λ(E[∆(x)H 2 (z T x)]) 2 + 2λE[∆(x)H 2 (z T x)]E[f lin (x)H 2 (z T x) -û2 ] Thus we have, ∇(G(z) -G lin (z)) = -sgn(û 4 )E[∆(x)∇H 4 (z T x)] + 2λE[∆(x)H 2 (z T x)]E[∆(x)∇H 2 (z T x)] + 2λE[f lin (x)H 2 (z T x) -û2 ]E[∆(x)∇H 2 (z T x)] + 2λE[∆(x)H 2 (z T x)]E[f lin (x)∇H 2 (z T x)]

Observe that H 2 and H 4 are degree 2 and 4 (respectively) polynomials thus norm of gradient and hessian of the same can be bounded by at most O(||z||||x|| 4 ). Using Lemma 14 we can bound each term by roughly O(log d)d -(1+p)η+3 ||z|| 4 . Note that λ being large does not hurt as it is scaled appropriately in each term. Subsequently, using Lemma 15, we can show that ||z|| is bounded by a constant since ||G(z)|| ≤ d -2η . Similar analysis holds for the hessian too.

Now applying Theorem 4 gives us that:

z is an (O(log d)d -(1+p)η+3 , O(log d)d -(1+p)η+3 )- approximate local minima of G lin . This implies that it is also an ( := C log(d)d -(1+2p)η+3 , τ := C log(d)d -(1+2p/3)η+3 )-approximate local minima of G lin for large enough C > 0 by increasing τ .

Observe that 

τ 3 /|û 4 | = C 3/2 log 3/2 (d)d -(3/2+p)η+9/2 /d -η/2 = C 3/2 log 3/2 (d)d -(1+p)η+9/2 ≥ . Now

G lin (W) = E   f lin (x) j,k∈[d],j =k ψ(w j , w k , x)   -γE   f lin (x) j∈[d] H 4 (w T j x)   (1) + λ i E f lin (x)H 2 (w T i x) -û2 2 (2) where ψ(v, w, x) = H 2 (v T x)H 2 (w T x) + 2(v T w) 2 + 4(v T x)(w T x)v T w.

They gave the following result. We show that this minimization is robust. Let us consider the corresponding function G to G lin with the additional non-linear terms as follows:

G(W) = E   f (x) j,k∈[d],j =k ψ(w j , w d , x)   -γE   f (x) j∈[d] H 4 (w j , x)   + λ i (E [f (x)H 2 (w i , x)] -û2 ) 2

Now we can show that G and G lin are close as in the one-by-one case.

R(W):

:= G(W) -G lin (W) = E [∆(x)A(W, x)] -γE [∆(x)B(W, x)] + λ E [f (x)C(W, x)] 2 -E [f lin (x)C(W, x)] 2 = E [∆(x)A(W, x)] -γE [∆(x)B(W, x)] + λE [(∆(x)C(W, x)(f (x ) + f lin (x ))C(W, x )] = E [∆(x)A(W, x)] -γE [∆(x)B(W, x)] + λE [(∆(x)D(W, x)] = E [∆(x)(A(W, x) -γB(W, x) + λD(W, x))] = E [∆(x)L(W, x)]

where

A(W, x) = j,k∈[d],j =k ψ(w j , w d , x), B(W, x) = j∈[d] H 4 (w j , x), C(W, x) = i H 2 (w i , x), D(W, x) = C(W, x)E[(f (x ) + f lin (x ))C(W, x )] and L(W, x) = A(W, x) - γB(W, x) + λD(W, x).

Using similar analysis as the one-by-one case, we can show the required closeness. It is easy to see that ||∇L|| and ||∇ 2 L|| will be bounded above by a constant degree polynomial in O(log d)d -(1+p)η+3 max ||w i || 4 . No row can have large weight as if any row is large, then looking at the gradient for that row, it reduces to the one-by-one case, and there it can not be larger than a constant. Thus we have the same closeness as in the one-by-one case. Combining this with Theorem 10 and 4, we have the following theorem: Theorem 11. Let c be a sufficiently small universal constant (e.g. c = 0.01 suffices), and under Assumptions 1, 2 and 3. Assume γ ≤ c, λ = Θ(d η ), and W * be the true weight matrix. The function G satisfies the following 1. Any saddle point W has a strictly negative curvature in the sense that λ min (∇ 2 G lin (W)) ≥ -τ where τ 0 = O(log d)d -Ω(1) .

2. Suppose W is a (d -Ω(1) , d -Ω(1) )-approximate local minimum, then W can be written as

W -T = PDW * + E where D is a diagonal matrix with D ii ∈ {±1 ± O(γ) ± d -Ω(1)

)}, P is a permutation matrix, and the error term ||E|| ≤ O(log d)d -Ω(1) .

Using standard optimization techniques we can find a local minima.

B.5 APPROXIMATE TO ARBITRARY CLOSE:

Lemma 16. If u is the sign function then E[u(w T x)δ (z T x)] = c| cot(α)| where w, z are unit vectors and α is the angle between them and c is some constant.

Proof. WLOG we can work the in the plane spanned by z and w and assume that z is the vector i along and w = i cos α + j sin α. Thus we can replace the vector x by ix + jy where x, y are normally distributed scalars. Also note that u = δ (Dirac delta function). Proof of Lemma 4. Let us compute the probability of lying in the -band for any t:

E[u(w T x)δ (z T x)] = E[u(x cos α + y sin α)δ (x)]

P r[x ∈ l(z, t, )] = Pr[t -≤ z T x ≤ t] = Pr g∈N (0,||z|| 2 ) [t -≤ g ≤ t] = 1 √ 2π||z|| t g=t- e -g 2 2||z|| 2 dg = √ 2π||z|| e - t2 2||z|| 2

where the last equality follows from the mean-value theorem for some t ∈ [t -, t].

Next we compute the following:

P r[x T w * 1 ≥ t and x ∈ l(z, t , )] = 1 (2π) n 2 x sgn(x 1 -t)1[x ∈ l(z, t , )]e -||x|| 2 2 dx = 1 (2π) 1 2 ∞ x1=t e -x 2 1 2 1 (2π) n-1 2 x-1 1[x -1 ∈ l(z -1 , t -z 1 x 1 , )]e -||x -1 || 2 2 dx -1 dx 1 = 1 (2π) 1 2 ∞ x1=t e -x 2 1 2 P r[x -1 ∈ l(z -1 , t -z 1 x 1 , )]dx -1 = 2π||z -1 || t g=t - ∞ x1=t e -x 2 1 2 e - (g-z 1 x 1 ) 2 2||z -1 || 2 dx 1 dg = 1 2π||z -1 || t g=t - e -g 2 2||z|| 2 ∞ x1=t e - x 1 - gz 1 ||z|| 2 2 2 ||z -1 || 2 ||z|| 2 dx 1 dg = 1 √ 2π||z|| t g=t - e -g 2 2||z|| 2 Φ c t||z|| 2 -gz 1 ) ||z -1 | |||z|| dg = √ 2π e -t * 2 2 Φ c t -t * cos(α 1 ) | sin(α 1 )|

where the last equality follows from the mean-value theorem for some t * ∈ [t -, t ]. Combining, we get:

P r[x T w 1 * ≥ t and x ∈ l(z, t , )|x ∈ l(z, t, )] = e -t * 2 -t 2 2 Φ c t -t * cos(α 1 ) | sin(α 1 )| = Φ c t -t * cos(α 1 ) | sin(α 1 )| ± O( )t

for ≤ 1/t .

Proof of Lemma 5. Recall that P is monotone with positive linear term, thus for high threshold u (0 unless input exceeds t and positive after) we have sgn(f (x)) = ∨sgn(x T w * i -t). This is because, for any i, P applied to X i > 0 and ∀j = i, X j = 0 gives us c i which is positive. Also, P (0) = 0. Thus, sgn(P ) is 1 if any of the inputs are positive. Using this, we have,

P r[sgn(f (x))|x ∈ l(z, t , )] ≥ P r[sgn((w * 1 ) T x -t)|x ∈ l(z, t , )] Also, P r[sgn(f (x))|x ∈ l(z, t , )] ≤ P r[sgn(x T w * i -t)|x ∈ l(z, t , )] = P r[sgn((w * 1 ) T x -t)|x ∈ l(z, t , )] + i =1 P r[sgn(x T w * i -t)|x ∈ l(z, t , )] ≤ P r[sgn((w * 1 ) T x -t)|x ∈ l(z, t, )] + η where i =1 P r[sgn(x T w * i -t)|x ∈ l(z, t , )] ≤ η.

We will show that η is not large since a z is close to one of the vectors, it can not be close to the others thus α i will be large for all i = j. Let us bound η,

i =1 P r[sgn(x T w * i -t)|x ∈ l(z, t , )] ≤ i =1 Φ c t -t * i cos(α i ) | sin(α i )| + O( )t i ≤ i =1 Φ c t -t * i cos(α i ) | sin(α i )| + O( )t ≤ i =1 Φ c t -t cos(α i ) | sin(α i )| + O( )t ≤ i =1 1 √ 2πγ i e -γ 2 i 2 + O( )kt

where γ i = t-t cos(αi) | sin(αi)| . The above follows since γ i ≥ 0 by assumption on t . Under the assumption, let β = max i =1 cos(α i ) we have

γ i ≥ t 1 -β cos(α1) 1 -β 2 = Ω(t)

under our setting. Thus we have,

i =1 P r[sgn(x T w * i -t))|x ∈ l(z, t , )] ≤ de -Ω(t 2 ) + O( )dt = de -Ω(t 2 )

for small enough .

Proof of Lemma 6. Let us assume that < c/t for sufficiently small constant c, then we have that

0.6 = P r[sgn(f (x)) and x ∈ l(z, t 2 , )|x ∈ l(z, t 2 , )] ≥ P r[x T w * 1 ≥ t and x ∈ l(z, t 2 , )|x ∈ l(z, t 2 , )] ≥ Φ c t -t * cos(α 1 ) | sin(α 1 )| -0.1 =⇒ 0.7 ≥ Φ c t -t * cos(α 1 ) | sin(α 1 )| =⇒ (Φ c ) -1 (0.7) ≤ t -t * cos(α 1 ) | sin(α 1 )| =⇒ t 2 ≤ t -(Φ c ) -1 (0.7) sin(α 1 ) cos(α 1 ) + O(1) ≤ t cos(α) + O(1)

Similarly for t 1 . Now we need to argue that t 1 , t 2 ≥ 0. Observe that P r[sgn(f (x)) and x ∈ l(z, 0, )|x ∈ l(z, 0, )]

≤ P r[x T w * i ≥ t and x ∈ l(z, 0, )|x ∈ l(z, 0, )] = Φ c t -cos(α 1 ) | sin(α 1 )| + O( 2 )d ≤ de -Ω(t 2 ) < 0.4

Thus for sufficiently large t = Ω( √ log d), this will be less than 0.4. Hence there will be some t 1 , t 2 ≥ 0 with probability evaluating to 0.4 since the probability is an almost increasing function of t up to small noise in the given range (see proof of Lemma 5).

Proof of Lemma 7. Let V be the plane spanned by w * 1 and z and let v 1 = w * 1 and v 2 be the basis of this space. Thus, we can write z = cos(α)v 1 + sin(α)v 2 .

Let us apply a Gaussian perturbation ρ along the tangential hyperplane normal to z. Say it has distribution N (0, 1) along any direction tangential to the vector z. Let 1 be the component of ρ on to V and let 2 be the component perpendicular to it. We can write the perturbation as ρ = 1 (sin(α)v 1 -cos(α)v 2 ) + 2 v 3 where v 3 is orthogonal to both v 1 and v 2 . So the new angle α of z after the perturbation is given by

cos(α ) = v T 1 (z + ρ) ||z + ρ|| = cos(α) + 1 sin(α) 1 + ||ρ|| 2

Note that with constant probability 1 ≥ as ρ is a Gaussian variable with standard deviation . And with high probability ||ρ|| < O( √ d -1). We will set = Θ(sin(α)/d) = Θ(α/d). Thus with constant probability:

cos(α ) ≥ cos(α) + sin(α) 1 + O( 2 d) ≥ (cos(α) + sin(α))(1 -O( 2 d)) ≥ cos(α) + Ω( sin(α)) -O( 2 d) ≥ cos(α) + Ω( sin(α)).

Thus change in cos(α) is given by ∆ cos(α) ≥ Ω( sin(α)). Now change in the angle α satisfies by the Mean Value Theorem:

∆ cos(α) = ∆α d dx cos(x) x∈[α,α ] =⇒ -∆α = ∆ cos(α) 1 sin(x) x∈[α,α ] ≥ Ω( sin(α)) sin(α) = Ω( ) = Ω(α/d).

B.6 LEARNING UNION OF HALFSPACES FAR FROM THE ORIGIN Theorem 12. Given non-noisy labels from a union of halfspaces that are at a distance Ω( √ log d) and are each a constant angle apart, there is an algorithm to recover the underlying weights to closeness in polynomial time.

Proof. Observe that X i is equivalent to P (X 1 , •, X d ) = 1 -(1 -X i ). Thus f (x) = sgn(x T w * i -t) = 1 - (1 -sgn(x T w * i -t)).

Since P and sgn here satisfies our assumptions 1, 2, for t = Ω( √ log d) (see Lemma 11) we can apply Theorem 11 to recover the vectors w * i approximately. Subsequently, refining to arbitrarily close using Theorem 5 is possible due to the monotonicity. Thus we can recover the vectors to arbitrary closeness in polynomial time.

B.7 SIGMOID ACTIVATIONS:

Observe that for sigmoid activation, Assumption 2 is satisfied for ρ(t, σ) = e -t+σ 2 /2 . Thus to satisfy Assumption 3, we need t = Ω(η log d).

Note that for such value of t, the probability of the threshold being crossed is small. To avoid this we further assume that f is non-negative and we have access to an oracle that biases the samples towards larger values of f ; that after x is drawn from the Gaussian distribution, it retains the sample (x, f (x)) with probability proportional to f (x) -so P r[x] in the new distribution. This enables us to compute correlations even if E x Ñ (0,I [f (x)] is small. In particular by computing E[h(x)] from this distribution, we are obtaining E[f (x)h(x)]/E[f (x)] in the original distribution. Thus we can compute correlations that are scaled. We get our approximate theorem: Theorem 13. For t = Ω(log d), columns of (TW * ) -1 can be recovered within error 1/poly(d) using the algorithm in polynomial time.

B.8 POLYNOMIALS P WITH HIGHER DEGREE IN ONE VARIABLE:

In the main section we assumed that the polynomial has degree at most 1 in each variable. Let us give a high level overview of how to extend this to the case where each variable is allowed a large degree. P now has the following structure,

P (X 1 , . . . , X d ) = r∈Z d + c r d i=1

X ri i If P has a higher degree in X i then Assumption 2 changes to a more complex (stronger) condition. Let q i (x) = r∈Z d + |∀j =i,rj =0 c r x ri , that is q i is obtained by setting all X j for j = i to 0.

Assumption 4. E g∼N (0,σ 2 ) [|u t (g)| r ] ≤ ρ(t, σ) for all r ∈ Z + 8 . E[q i (u t (g))h k (g))] = t Θ(1) ρ(t, 1) for k = 2, 4. Lastly, for all d ≥ i > 1, r∈Z d + ||r||0=i |c r | ≤ d O(i) .

The last assumption holds for the case when the degree is a constant and each coefficient is upper bounded by a constant. It can hold for decaying coefficients.

Let us collect the univariate terms P uni (X) = d i=1 q i (X i ). Corresponding to the same we get f uni . This will correspond to the f lin we had before. Note that the difference now is that instead of being the same activation for each weight vector, now we have different ones q i for each. Using H 4 correlation as before, now we get that:

E[f uni (x)H 4 (z T x)] = d i=1 q i • u t4 (z T w * i ) 4 and E[f uni (x)H 2 (z T x)] = d i=1 q i • u t2 (z T w * i ) 2

where q i • u t are hermite coefficients for q i • u t . Now the assumption guarantees that these are positive which is what we had in the degree 1 case.

Second we need to show that even with higher degree, E[|f (x) -f uni (x)|] is small. Observe that Lemma 17. For r such that ||r|| 0 > 1, under Assumption 4 we have,

E d i=1 u t ((w * j ) T x) ri ≤ ρ(t, 1)O (ρ(t, ||W * ||)) ||r||0-1 .

The proof essentially uses the same idea, except that now the dependence is not Proof. We have,

E[|∆(x)|] = E      r∈Z d + ri≤D,||r||0>1 c r d i=1 u t ((w * j ) T x) ri      ≤ r∈Z d + ri≤D,||r||0>1 |c r |E d i=1 u t ((w * j ) T x) ri = r∈Z d + ri≤D,||r||0>1 |c r |ρ(t, 1) (ρ(t, ||W * ||)) ||r||0-1 ≤ C d i=1 d i D i ρ(t, 1) (ρ(t, ||W * ||)) i-1 ≤ d C d i=1 ρ(t, 1) d C ρ(t, ||W * ||) i-1 using Assumption 4 ≤ d 2C+1 ρ(t, 1)ρ(t, ||W * ||) since ρ(t, ||W * ||) ≤ d -C .

Thus as before, if we choose t appropriately, we get the required results. Similar ideas can be used to extend to non-constant degree under stronger conditions on the coefficients.

C SAMPLE COMPLEXITY:

Proof of Lemma 3.

C 1 (f, z, s) = E[f (x)δ(z T x -s)] = x f (x)δ(z T x -s)φ(x)

dx Let x 0 be the component of x along z and y be the component along z ⊥ . So x = x 0 ẑ + yz ⊥ . Interpreting x as a function of x 0 and y:

C 1 (f, z, s) = y x0 f (x)δ(x 0 -s)φ(x 0 )φ(y)dx 0 dy = y [f (x)] x0=s φ(y)dy = φ(s)E[f (x)|x 0 = s] = φ(z T x = s)E[f (x)|x 0 = s]

where the second equality follows from x δ(x -a)f (x) = [f (x)] x=a .

Proof of Claim 4. Let x 0 be the component of x along z and y be the component of x in the space orthogonal to z. Let ẑ denote a unit vector along z. We have x = x 0 ẑ + y and ∂x ∂x0 = ẑ. So, correlation can be computed as follows:

E[P [x].δ (z T x -s)] = y φ(y) x0 δ (x 0 -s)P [x]φ(x 0 )dx 0 dy Since x δ (x -a)f (x)dx = [ df dx ](x = a)

this implies:

E[P [x]δ (z T x -s)] = y ∂ ∂x 0 P [x]φ(x) x0=s dy = y φ(x 0 ) i ∂P ∂X i . ∂X i ∂x 0 + P [x]φ (x 0 )) x0=s φ(y)dy = i y φ(x 0 ) ∂P ∂X i . ∂X i ∂x 0 x0=s φ(y)dy + φ (s) y P [x]φ(y)dy Note that ∂Xi ∂x0 = ∂ ∂x0 u(x T w * i -t) = u (x T w * i -t) ẑT w * i .

If u is the sign function then u (x) = δ(x). So focusing on one summand in the sum we get

y φ(x 0 ) ∂P ∂X i . ∂X i ∂x 0 x0=s φ(y)dy = y [φ(x 0 )u (x T w * i -t)( ẑT w * i ) ∂P ∂X i ] x0=s φ(y)dy = y ( ẑT w * i )[φ(x 0 )δ( ẑT w * i x 0 + (w * i ) T y -t) ∂P ∂X i ] x0=s φ(y)dy = ( ẑT w * i ) y φ(s)δ(s(w * 1 ) T ẑ + (w * 1 ) T y -t) ∂P ∂X i φ(y)dy

Again let y = y 0 (w * i ) + z where z is perpendicular to w * i and z. And (w * i ) is perpendicular component of w * i to z. Interpreting x = t ẑ + y 0 (w * 1 ) + z as a function of y 0 , z we get:

= ( ẑT w * i ) z y0 φ(s)δ(s(w * 1 ) T ẑ + ((w * 1 ) T (w * 1 ) )y 0 -t)φ(y 0 )φ(z) ∂P ∂X i dy 0 dz Note that by substituting v = ax we get ∞ x=-∞ f (x)δ(ax -b)dx = ∞/a x=-∞/a f (x)δ(ax -b)dx = sgn(a) 1 a f ( b a ) = 1 |a| f ( b a )

. So this becomes:

= ẑT w * i |(w * i ) T (w * i ) | z φ(s)[φ(y 0 ) ∂P ∂X i ] y0= t-s ẑT w * i (w * i ) T (w * i ) φ(z)dz = ẑT w * i |(w * i ) T (w * i ) | z [φ(y 0 )φ(x 0 )φ(z) ∂P ∂X i ] x0=s,y0= t-s ẑT w * i (w * i ) T (w * i ) dz = ẑT w * i |(w * i ) T (w * i ) | φ y 0 = t -s ẑT w * i (w * i ) T (w * i ) φ(x 0 = t) z [φ(z) ∂P ∂X i ] x0=t,y0= t-s ẑT w * i (w * i ) T (w * i ) dz = ẑT w * i |(w * i ) T (w * i ) | z [φ(x) ∂P ∂X i ] x0=t,y0= t-s ẑT w * i (w * i ) T (w * i ) dz = ẑT w * i |(w * i ) T (w * i ) | x:z T x=s,x T w * i =t φ(x) ∂P ∂X i dx = ẑT w * i |(w * i ) T (w * i ) | φ(z T x = s, x T w * i = t)E[ ∂P ∂X i |x T z = s, x T w * i = t].

Let α i be the angle between z and w * i . Then this is

= | cot(α i )|φ(z T x = t, (w * i ) T x = t)E[ ∂P ∂X i ]|x T z = s, x T w * i = t]

Thus, overall correlation

= i∈S | cot(α i )|φ(z T x = s, x T w * i = t)E[ ∂P ∂X i |x T z = s, x T w * i = t] + φ (s)E[P [x]|x T z = s]

Proof of Claim 5. Note that for small α, | cot α| = O(1/α) ≤ O(1/ 2 ). Since P (X) ≤ ||X|| c1 , we have f (x) ≤ d c1 (as with sgn function each sgn((w * i ) T x) ≤ 1) and all

Q i [x], R i [x] are at most 2d c1 . By Cauchy's mean value theorem C 2 (f, z, t) = 2 [E[f (x)δ (z T x = s)]] s∈t± . Note that cot(α i )φ(z T x = t, (w * i ) T x = t) = cot(α i )φ(t tan(α i /2))

which is a decreasing function of α i in the range [0, π] So if all α i are upper bounded by 2 then by above corollary,

C 2 (f, z, s) ≤ 2 n cot( 2 )φ(z T x = t, (w * 1 ) T x = t)(2d c1 ) + (2d c1 ) = 2 n cot( 2 )φ(t tan( 2 /2))(2d c1 ) + (2d c1 ) ≤ 2 d O(1) .

Observe that the above proof does not really depend on P and holds for for any polynomial of u((w * i ) T x) as long as the polynomial is bounded and the w * i are far off from z.

Proof Of Lemma 8.

If z = w * i , then C 2 (f, z, t) = E[f (x)(δ(z T x -t -) -δ(z T x -t + ))] = E[u((w * i ) T x)Q i [x](δ(z T x -t -) -δ(z T x -t + ))] + E[R i [x](δ(z T x -t -) -δ(z T x -t + ))]

Since u((w * i ) T x) = 0 for z T x = tand 1 for z T x = t + , and using the Cauchy mean value theorem for the second term this is

= E[Q i [x]δ(z T x -t -)] + 2 [E[R i [x]δ (z T x -s 1 )] s1∈t± = E[Q i [x]δ(z T x -t)] + E[R i [x](δ(z T x -t -) -δ(z T x -t))] = φ(t)E[Q i [x]|z T x = t + ] + [C 2 (Q i , z, s 2 )] s2∈[t,t+ ] + 2 [C 2 (R i , z, s)] s∈t± = φ(t)E[Q i [x]|z T x = t] + d O(1)

The last step follows from Claim 5 applied on Q i and R i as all the directions of w * j are well separated from z = w * i and w * i is absent from both Q i and R i . Also the corresponding Q i and R i are bounded.

C.1 RELU ACTIVATION:

If u is the RELU activation, the high level idea is to use correlation with the second derivative δ of the Dirac delta function instead of δ . More precisely we will compute C

3 (f, z, s) = E[f.(δ (z T x- s -) -δ (z T x -s + )].

Although we show the analysis only for the RELU activation, the same idea works for any activation that has non-zero derivative at 0.

Note that now u = sgn and u = δ.

For ReLU activation, Lemma 8 gets replaced by the following Lemma. The rest of the argument is as for the sgn activation. We will need to assume that P has constant degree and sum of absolute value of all coefficients is poly(d) Lemma 19. Assuming polynomial P has constant degree, and sum of the magnitude of all coefficients is at most poly 1) . Otherwise if all angles α i between z and w i are at least 2 it is at most d O(1) / 2 .

(d), if z = w * i then C 3 (f, z, t) = E[φ(t) ∂ ∂Xi P + φ(t) j =i cos(α j )sgn(x T w * i -t) ∂ ∂Xj P + φ (t)P |x T w * i = t] + d O(

We will prove the above lemma in the rest of this section. First we will show that z is far from any of the w * i 's then E[P.δ (z T x -s)] is bounded. Lemma 20. If the sum of the absolute value of the coefficients of P is bounded by poly(d), its degree is at most constant,

α i > 2 then E[P.δ (z T x -s)] is d O(1) / 2 .

Proof. Let x 0 be the component of x along z and y be the component of x in the space orthogonal to z as before. We have x = x 0 ẑ + y and ∂x ∂x0 = ẑ. We will look at monomials M l in P = l M l . As before since

x δ (x -a)f (x)dx = d 2 f dx 2 x=a we get E[f (x).δ (z T x -s)] = y f (x)φ(x)δ (z T x -s)dy = y ∂ 2 ∂x 2 0 (P [x]φ(x)) x0=s dy = l y ∂ 2 ∂x 2 0 (M l [x]φ(x 0 )) x0=s φ(y)dy Now consider a monomial M = X i1 1 ..X i k k . Take the symbolic second derivative ∂ 2 ∂x 2 0 of M [x]φ(x 0 ) w.r.t x 0 . This will produce a polynomial involving X i 's, ∂Xi ∂x0 , ∂ 2 Xi ∂x 2 0 , φ, φ , φ . Let us examine each of these terms. ∂ ∂x 0 X i [x] = ∂ ∂x 0 u(x T w * i -t) = sgn(x T w * i -t)(w * i ) T ∂ ∂x 0 x = sgn(x T w * i -t)((w * i ) T z) = cos(α i )sgn(x T w * i -t) Thus ∂ ∂x0 X i [x] is a bounded function of x. We have ∂ 2 ∂x 2 0 X i (x) = ∂ ∂x 0 sgn(x T w * i -t)((w * i ) T z) = cos 2 (α i )δ(x T w * i -t).

Again as before y

[δ(x T w * i -t)g(x)φ(x 0 )φ(y)] x0=s dy = (1/| sin(α i )|)E[g(x)|x 0 = s, x T w * i = t]φ(x 0 = s, x T w * i = t)

Note that if the degree is bounded, since | sin(α i )| is at least 2 expected value of each monomial obtained is bounded. So the total correlation is poly(d)/ 2 .

Proof of Lemma 20. As in the case of sgn activation, if z = w * i ,

E[P [x]δ (x T w * i -s)] = x P [x]φ(x)δ (z T x -s)dx = y ∂ ∂x 0 (P [x]φ(x)) x0=s dy = y   j ∂P ∂X j [x] ∂X j [x] ∂x 0 φ(x) + P [x] ∂φ(x) ∂x 0   x0=s dy = y φ(y)   sgn(x 0 -t)φ(s) ∂P ∂X i [x] + j =i sgn(x T w * j -t) cos(α j )φ(s) ∂P ∂X j [x] + P φ (x 0 )   x0=s dy = y φ(y)   sgn(s -t)φ(s) ∂P ∂X i [x] + j =i sgn(x T w * j -t) cos(α j )φ(s) ∂P ∂X j [x] + P φ (x 0 )   x0=s dy For s = t + this is = y φ(y)   φ(s) ∂P ∂X i [x] + j =i cos(α j )sgn(x T w * j -t)φ(s) ∂P ∂X j [x] + P φ (s)   s=t+ dy =   y φ(y)   φ(t) ∂P ∂X i [x] + j =i cos(α j )sgn(x T w * j -t)φ(t) ∂P ∂X j [x] + P φ (t)   dy   + d O(1) =E   φ(t) ∂P ∂X i [x] + φ(t) j =i cos(α j )sgn(x T w * i -t) ∂P ∂X j [x] + φ (t)P x T w * i = t   + d O(1)

If z is away from every w * i by at least 2 then again E

[P C 3 (f, z, t)] = E[P δ (x T w * i -s)](s ∈ [t -, t + ]) = d O(1) / 2 .

D STRUCTURAL RESTRICTIONS HELPS LEARNING D.1 PROOF OF THEOREM 7:

To construct this correlation graph, we will run the following Algorithm 3 Denote T

i := {j : w ij = 1}. Let us compute E[f (x)x i x j ]: E[f (x)x i x j ] = S⊆[d] c S E     p∈S u(x T w * p )   x i x j   = S⊆[d]

c S e -ρt|S| E e ρ p∈S x T w * p x i x j Algorithm 3 ConstructCorrelationGraph 1: Let G be an undirected graph on n vertices each corresponding to the x i 's. 2: for every pair i, j do 3:

Compute

α ij = E[f (x)x i x j ]. 4: if α ij ≥ ρ then 5:

Add edge (i, j) to the graph.

= S⊆[d] c S e -ρt|S| E e ρ q∈∪ p∈S Tp xq x i x j = S⊆[d] c S e -ρt|S| 1[i, j ∈ ∪ p∈S T p ]E [x i e ρxi ] E [x j e ρxj ]

q∈∪ p∈S Tp\{i,j}

E [e ρxq ] = S⊆[d] c S e -ρt|S| 1[i, j ∈ ∪ p∈S T p ]ρ 2 e ρ 2 E [x j e ρxj ] q∈∪ p∈S Tp\{i,j} e ρ 2 /2 = S⊆[d] c S e -ρt|S| 1[i, j ∈ ∪ p∈S T p ]ρ 2 e ρ 2 | ∪ p∈S Tp | 2

By assumption, for all p, T p are disjoint. Now, if i, j ∈ T r for some r, we have

EQUATION

Similarly, if i ∈ T r1 and j ∈ T r2 with r 1 = r 2 , we have

EQUATION

It is easy to see that these correspond to coefficients of X r and X r1 X r2 (respectively) in the following polynomial:

Q(X 1 , . . . , X n ) = ρ 2 P (µ(X 1 + 1), . . . , µ(X n + 1))

for µ = e -ρ(t-ρd) 2

. For completeness we show that this is true. We have,

Q(X 1 , . . . , X n ) = ρ 2 S⊆[d] c S j∈S µ(X j + 1) = ρ 2 S⊆[d] c S µ |S| j∈S (X j + 1)

The coefficient of X r in the above form is clearly ρ 2 S⊆[d]:r∈S c S µ |S| (corresponds to picking the 1 in each product term). Similarly coefficient of X r1 X r2 is ρ 2 S⊆[d]:r1,r2∈S c S µ |S| . Now as in the assumptions, if we have a gap between these coefficients, then we can separate the large values from the small ones and form the graph of cliques. Each clique will correspond to the corresponding weight vector. Consider the following optimization problem:

max z g(z) -λ||z|| 1 -γ||z|| 2 2 h(z)

for λ, γ > 0 to be fixed later.

We can assume that z i ≥ 0 for all i at a local maxima else we can move in the direction of e i and this will not decrease g(z) since c S ≥ 0 for all S and will decrease ||z|| 1 and ||z|| 2 making h(z) larger. From now on, we assume this for any z at local maxima.

We will show that the local maximas of the above problem will have z such that most of the mass is equally divided among j ∈ S i for some i and close to 0 everywhere else. Lemma 21. There exists at most one i such that there exists j ∈ S i with |z j | ≥ β for γ satisfying 4γ < min i =j e β S:i∈S,j ∈S∨i ∈S,j∈S c S e

|∪ i∈S S i | 2

.

Proof. Let us prove by contradiction. Assume that there is a local maxima such that there are at least 2 indices say 1, 2 such that ∃j ∈ S 1 , k ∈ S 2 , |z j |, |z k | ≥ β. Now we will show that there exists a perturbation such that g(z) can be improved. Now consider the following perturbation, z + s e j -s e k for s ∈ {±1}. Observe that ||z|| 1 remains unchanged for < β also ||z|| 2 2 changes by 2s 2 2 + 2(z j -z k )s . We have For chosen value of γ, there will always be an improvement, hence can not be a local maxima.

Lemma 22. At the local maxima, for all i ∈ [n], z is such that for all j, k ∈ S i , z j = z k at local maxima.

Proof. We prove by contradiction. Suppose there exists j, k such that z j < z k . Consider the following perturbation: z + (z k -z j )(e j -e k ) for 1 ≤ > 0. Observe that g(z) depends on only r∈Si z r and since that remains constant by this update g(z) does not change. Also note that ||z|| 1 does not change. However ||z|| 2 2 decreases by 2 (1 -)(z k -z j ) 2 implying that overall h(z) increases. Thus there is a direction of improvement and thus it can not be a local maxima. Proof. We prove by contradiction. Suppose ||z|| 1 < α, consider the following perturbation, z + 1. For given λ there is a direction of improvement giving a contradiction that this is the local maxima.

Combining the above, we have that we can choose λ, γ = poly(n, 1/ , s) where s is a paramater that depends on structure of f such that at any local maxima there exists i such that for all j ∈ S i , z j ≥ 1 and for all k ∈ ∪ j∈Si , z k ≤ . Let us consider correlation with h 4 (z T x). This above can be further simplified by observing that when we correlate with h 4 , i∈S u (x i )h 4 (z T x) = 0 for |S | ≥ 2. Observe that h 4 (z T x) = d1,...,dn∈[4]: di≤4 c(d 1 , . . . , d n ) h di (x i ) for some coefficients c which are functions of z. Thus when we correlate i∈S u (x i )h 4 (z T x) for |S | ≥ 3 then we can only get a non-zero term if we have at least h 2k (x i ) with k ≥ 1 for all i ∈ S . This is not possible for |S | ≥ 3, hence, these terms are 0. Thus,

E i∈S u(x i )h 4 (z T x) = 2 k=0 û|S|-k 0 S ⊆S:|S |=k E i∈S u (x i )h 4 (z T x) .

Lets compute these correlations. 

E u (x i )h 4 (z T x) = E p>0 û2p h 2p (x i )h 4 (z i x i + z T -i x -i ) = 1 4 p>0 û2p E h 2p (x i ) 4 k=0 4 k h 4-k (z i x i √ 2)h k (z T -ij x -ij √ 2) = 1 4 p>0 û2p E h 2p (x i ) h 4 (z i x i √ 2) + 6h 2 (z i x i √ 2)(

2 k=0 2 k E [h 2p (x i )h 2q (x j )h 2-k (2z i x i )h k (2z j x j )] = 3 4 (2||z -ij || 2 -1) 2 k=0 2 k û2 û2 α(2 -k, 2, 2z i )α(k, 2, 2z j ) = 0.

Combining, we get E u (x i )u (x j )h 4 (z T x) = 6û 2 2 z 2 i z 2 j . Further, taking correlation with f , we get: 

E f (x)h 4 (z T x) = S⊆[n]

:

We suppress W when it is clear from context.

Here H4 and H2 are the fourth and second order hermite polynomials respectively.

Theorem 1 holds for sigmoid with t ≥ c log d.

We can handle ∈ [d -C , d C ] for some constant C by changing the scaling on t.

For similar bounds for sigmoid and sign refer to Appendix B.7.

for example with RELU activation, f will be c0 most of the time as other terms in P will never activate. So c0 can be set to say the median value of f .

The theorem can be extended to ReLU by correlating with the second derivative δ (see Appendix C.1).

For example, his would hold for any u bounded in [-1, 1] such as sigmoid or sign.

Figures and Tables:

z T ∇L(z) ||z|| 2 z = 0 and ||z|| 2 = 1.

Now using Lemma 2 to bound E[|∆(x)|], for c = Θ( √ η log d we get the required result. Lemma 15. For ||z|| = Ω(1) and λ

using Claim 1, we get the required result. B.4 SIMULTANEOUS RECOVERY Ge et al. (2017) also showed simultaneous recovery by minimizing the following loss function G lin defined below has a well-behaved landscape.

Ge et al. (2017)). Let c be a sufficiently small universal constant (e.g. c = 0.01 suffices), and suppose the activation function u satisfies û4 = 0. Assume γ ≤ c, λ ≥ Ω(|û 4 |/û 2 2 ), and W * be the true weight matrix. The function G lin satisfies the following:1. Any saddle point W has a strictly negative curvature in the sense that λ min (∇ 2 G lin (W)) ≥ -τ 0 where τ 0 = c min{γ|û 4 |/d, λû 2 2 }.2. Suppose W is an ( , τ 0 )-approximate local minimum, then W can be written as W -T = PDW * + E where D is a diagonal matrix with D ii ∈ {±1 ± O(γ|û 4 |/λû 2 2 ) ± O( /λ)}, P is a permutation matrix, and the error term ||E|| ≤ O( d/û 4 ).

cos α + y sin α)δ (x)φ(x)φ(y)dxdy Using the fact that x δ (x)h(x)dx = h (0) this becomes = y φ(y)[(∂/∂x)u(x cos α + y sin α)φ(x)] x=0 dy = y φ(y)[n(x)u (x cos α + y sin α) cos α + φ (x)u(x cos α + y sin α)] x=0 dy = ∞ y=-∞ φ(y)φ(0)δ(y sin α) cos αdy Substituting s = y sin α this becomes = ∞/ sin α s=-∞/ sin α φ(s/ sin α)φ(0)δ(s) cos α(1/ sin α)ds =sgn(sin α) cot(α)φ(0) s φ(s/ sin α)δ(s)ds =| cot(α)|φ(0)φ(0)

PROOF OF THEOREM 8 Consider f (x) = c S i∈S e x T w * i where w * i = j∈Si e j for S i ⊆ [n] such that for all i = j, S i ∩ S j = n and c S ≥ 0. Let us compute g(z) = e -||z|| 2 2 E f (x)e z T x where z = z i e i for some α i . E f (x)e z T x = E c S i∈S e j∈S i xj e z T x = c S E i∈S e p∈S i xp e q∈[n]

s [h(z + s e j -s e k ) -h(z)] E s e -s -1 -γE s 2 2 + 2(z j -z k )sThe inequality follows since E [e s ] = e +e - ,2 ∈S c S e |∪ i∈S S i | 2 e zj ≥ S:1∈S,2 ∈S c S e |∪ i∈S S i | 2 e β .

At the local maxima, ||z|| 1 ≥ α for λ < S c S |∪ i∈S Si| n e |∪ i∈S S i | 2-γ(2α + 1).

p∈∪ i∈S S i zp (e |∪ i∈S Si| -1) -nλ -nγ (2||z|| 1 + ) > c S e |∪ i∈S S i | 2 | ∪ i∈S S i | -nλ -nγ (2α + 1)

PROOF OF THEOREM 9 Let function f = S⊆[n] c S i∈S u((w * i ) T x) for orthonormal w * i . WLOG, assume w * i = e i . Let u(x) = ∞ i=0 û2i h 2i (x)where h i are hermite polynomials and u (x) = u(x) -û0 = ∞ i=1 û2i h 2i (x). This implies E[u (x)] = 0. Observe that,

i )h 4 (z T x)

on ||r|| 1 but only the number of non-zero entries (number of different weight vectors). With this bound, we can now bound the deviation in expectation. Lemma 18. Let ∆(x) = f (x)-f uni (x). Under Assumptions 4, if t is such that ρ(t, ||W * ||) ≤ d -C for large enough constant C > 0, we have, E[|∆(x)|] ≤ d O(1) ρ(t, 1)ρ(t, ||W * ||).

Table content: <table/>

2||z -i || 2 -1) + 3(2||z -i || 2 -1) 2 -1) + 6û 2 (2z 2 i )(2||z -i || 2 -1) = û4 z 4 i + 3û 2 z 2 i (2||z|| 2 -1) E u (x i )u (x j )h 4 (z T x) 2p (x i )h 2q (x j )h 4 (z i x i + z j x j + z T -ij x -ij ) 2p (x i )h 2q (x j ) h 4 ((z i x i + z j x j ) √ 2) + 6h 2 ((z i x i + z j x j ) √ 2)(2||z -ij || 2 -1) +3(2||z -ij || 2 -1) 2 2p (x i )h 2q (x j )h 4-k (2z i x i )h k (2z j x j )] 2p (x i )h 2q (x j )h 2-k (2z i x i )h k (2z j x j )] [h 2p (x i )h 4-k (2z i x i )] E [h 2q (x j )h k (2z j x j )]

Table content: <table><tr><td colspan="2">=</td><td colspan="2">1 4</td><td colspan="5">û4 α(4, 4, z i</td><td>√</td><td>2) + û2 α(4, 2, z i</td><td>√</td><td>2) + 6û 2 α(2, 2, z i</td><td>√</td><td>2)(2||z -i || 2 -1)</td></tr><tr><td colspan="10">= i = E 1 4 4û 4 z 4 i + 12û 2 z 2 i (2z 2 p,q&gt;0 û2p û2q h = 1 4 p,q&gt;0 û2p û2q E h 2p (x i )h 2q (x j )</td><td>4 k=0</td><td>4 k</td><td>h 4-k ((z i x i + z j x j ) √</td><td>2)h k (z T -ij x -ij</td><td>√</td><td>2)</td></tr><tr><td colspan="10">= û2p û2q E h = 1 4 p,q&gt;0 1 16 p,q ûp ûq 4 k=0 4 k</td><td>E [h 1</td></tr><tr><td>+</td><td colspan="2">3 4</td><td colspan="7">(2||z -ij || 2 -1)</td><td>p,q&lt;0</td><td>û2p û2q</td><td>2 k=0</td><td>2 k</td><td>E [h 2</td></tr><tr><td colspan="10">We will compute 1 and 2 :</td></tr><tr><td/><td/><td/><td/><td/><td colspan="5">1 = û2p û2q E = 1 16 4 k=0 4 k p,q&gt;0 1 16 4 k=0 4 p,q&gt;0 k û2p û2q α(4 -k, 2p, 2z i )α(k, 2q, 2z j )</td></tr><tr><td/><td/><td/><td/><td/><td/><td>=</td><td>3 8</td><td colspan="2">û2 2 α(2, 2, 2z i )α(2, 2, 2z j )</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="4">= 6û 2 2 z 2 i z 2 j</td></tr><tr><td colspan="4">Similarly,</td><td/><td/><td/><td/><td/></tr><tr><td/><td colspan="4">2 =</td><td>3 4</td><td colspan="4">(2||z -ij || 2 -1)</td><td>p,q&gt;0</td><td>û2p û2q</td></tr></table>

Peer review (iclr2019):
 This paper gives provable recovery guarantees for a class of neural networks which have high-threshold activation in the first layer, followed by a "well-behaved" polynomial, under Gaussian input. The algorithm is based on the approach by Ge et al. (2017), as well as an iterative refinement method. While this could be an interesting result, I have several concerns regarding the assumptions, correctness, and writing. 1) It is required that the threshold is at least sqrt{log d} (Thm. 1), where d is the number of hidden neurons in the first layer. It seems that this essentially zeros out almost all the neurons, since the maximum among d Gaussian random variables is roughly sqrt{log d}. The authors should explain what exactly this model is doing, i.e., what kind of functions it can compute, in order to justify why this is an interesting model. Furthermore, the authors claim that the studied model is a "deep" neural network, but I disagree. As I understand, the difference between this model and two-layer networks is that the second layer here is a polynomial instead of a linear function. This doesn*t make it a deep network since the (polynomial) part above the first layer is not modeled in a layer-wise fashion, not to mention that under the setting considered in the paper the polynomial behaves similar to a linear function. 2) It is stated at the end of Section 2 that the angle can be reduced by a factor of 1-1/d ***with constant probability***. How does this ensure you can succeed after O(d log(1/nu)) iterations? As far as I see you need the success probability in one iteration to be at least something like 1-1/Omega(d) so that you can apply a union bound. 3) Even if the issues of motivation and correctness are clarified, I find it very difficult to understand the overall intuition and main technical contributions in this paper. The writing needs to be significantly improved to reach the level of a top conference.