Title: AN INFORMATION RETRIEVAL APPROACH FOR FIND-ING DEPENDENT SUBSPACES OF MULTIPLE VIEWS

Abstract: Finding relationships between multiple views of data is essential both in exploratory analysis and as pre-processing for predictive tasks. A prominent approach is to apply variants of Canonical Correlation Analysis (CCA), a classical method seeking correlated components between views. The basic CCA is restricted to maximizing a simple dependency criterion, correlation, measured directly between data coordinates. We introduce a new method that finds dependent subspaces of views directly optimized for the data analysis task of neighbor retrieval between multiple views. We optimize mappings for each view such as linear transformations to maximize cross-view similarity between neighborhoods of data samples. The criterion arises directly from the well-defined retrieval task, detects nonlinear and local similarities, measures dependency of data relationships rather than only individual data coordinates, and is related to well understood measures of information retrieval quality. In experiments the proposed method outperforms alternatives in preserving cross-view neighborhood similarities, and yields insights into local dependencies between multiple views.

1. INTRODUCTION:

Finding dependent subspaces across views (subspaces where some property of data is statistically related or similar across views) is a common data analysis need, where Canonical Correlation Analysis (CCA) (Hotelling, 1936 ) is a standard unsupervised tool. Preprocessing to find dependent subspaces is useful both for prediction and for analysis: in predictive tasks, such subspaces help if non-dependent parts of each view may arise from noise and distortions. In some data analysis tasks, finding the dependent subspaces may itself be the main goal; for example in bioinformatics domains dependency seeking projections have been used to identify relationships between different views of cell activity (Tripathi et al., 2008; Klami et al., 2013) ; in signal processing a similar task could be identifying optimal filters for dependent signals of different nature, e.g., speech and the corresponding tongue movements of the speakers as in Westbury (1994) .

Methods like CCA maximize simple correlations between data point coordinate features across the projected subspaces. However, in many data domains the coordinates may not be of main interest but rather the data relationships that they reveal. It is then of great interest to develop dependency seeking methods that directly focus on the data relationships. For example, consider a database of scientists, defined in one view by their level of interest in various research topics, and in another view by their level of interest in various hobbies. In a database like this, finding relationships of people is the common interest, e.g. to find nearest colleagues for a scientist, having the most similar (neighboring) research interests; or to find hobby partners having the most similar (neighboring) hobby interests; the question is then, can we predict the research colleagues from hobby partners or vice versa? Research topics and hobbies are very dissimilar views, and not all of their variation will be related, but we can try to find subspaces of research and hobby interests, so that research neighbors and hobby neighbors are as highly related as possible in those subspaces.

In this paper we propose a method that solves this task: we present a novel method for seeking dependent subspaces across multiple views, preserving neighborhood relationships of data. Our method directly maximizes the between-view similarity of neighborhoods of data samples, a natural measure for similarity of data relationships among the views. The method detects nonlinear and local dependencies, has strong invariance properties, is related to an information retrieval task of the analyst, and performs well in experiments.

Relating data items is one of the main elementary tasks in Shneiderman's taxonomy of tasks in visual data analysis (Shneiderman, 1996) . Our method is optimized for finding related (neighboring) data items, formulated as optimizing an information retrieval task. Since our method directly serves the task of relating data items (across views) in Shneiderman's taxonomy, in this sense it arguably comes closer to needs of data analysts than maximizing some variant of coordinate correlation.

We find linear projections (linear subspaces) of views. Linear projections have advantages of simplicity and easy interpretability with respect to original data features. Even if projections are linear, the dependency criterion we optimize is flexible and detects nonlinear dependencies across views.

We present our method in Section 2, properties and extensions in Section 3, related work in Section 4, experiments in Section 5, and conclusions in Section 6.

2. METHOD: DEPENDENT NEIGHBORHOODS OF VIEWS:

Our method finds similar neighborhood relationships between views. We define the neighborhood relationships and then discuss how to measure their cross-view similarity. Instead of hard neighborhoods where two points are or are not neighbors, we use more realistic probabilistic neighborhoods. Assume input data items x i = (x i,1 , . . . , x i,N Views ) have paired features x i,V in each view V . We consider transformations of each view by a mapping f V which is typically a dimensionality reducing transformation to a subspace of interest; in this paper, for simplicity and interpretability we use linear mappings f

V (x i,V ) = W T V x i,V with W T V ∈ R dim low (V )×dim orig (V )

as the to-be-optimized parameters, where dim orig (V ) and dim low (V ) are the number of dimensions of V and its subspace respectively. The local neighborhood of a data item i in any transformation of view V can be represented by the conditional probability distribution p i,V = {p V ( j|i; f V )} where p V ( j|i; f V ) tells the probability that data item j = i is picked as a representative neighbor of i; that is, the probability that an analyst who inspected item i will next pick j for inspection. The p V ( j|i; f V ) can be defined in several ways, as a decreasing function of distance d V (i, j; f V ) between features of i and j in view V . Here we define it by a simple exponential falloff with respect to squared distance of i and j, as

p V ( j|i; f V ) = exp(-d 2 V (i, j; f V )/σ 2 i,V ) • ( ∑ k =i exp(-d 2 V (i, k; f V )/σ 2 i,V )) -1 (1)

where σ i,V sets the falloff rate around i in the view. We tried two simple ways to set the σ i,V : one is as a fraction of maximal pairwise distance so σ i,V = 0.05 • max j,k x j,Vx k,V , or alternatively, set

EQUATION

i.e., calculate the average distance between x j,V and its k-th nearest neighbor x l,V , then give the average a heuristic correction factor dim low (V )/dim orig (V ) since the average distance is obtained from the original space yet σ i,V is used in a subspace. We use the first simple σ i,V for artificial data experiments and the more data-driven second σ i,V from (2) with k = 5 for the other experiments. Both choices give good results. Other local choices to e.g. achieve a desired entropy are possible, see Venna et al. (2010) . With linear mappings the probabilities become

EQUATION

where the matrix W V defines the subspace of interest for the view and also the distance metric within the subspace. Our method learns the mapping parameters W V for each view.

2.1. COMPARISON OF NEIGHBORHOODS ACROSS VIEWS:

Neighborhoods represented as probability distributions can be compared by difference measures. We discuss two measures for different purposes, and their information retrieval interpretations.

Kullback-Leibler divergence. For two distributions p = {p( j)} and q = {q( j)}, the Kullback-Leibler (KL) divergence is an information-theoretic asymmetric difference measure defined as D KL (p, q) = ∑ j p( j)(log p( j)/q( j)) .

(4)

The KL divergence is nonnegative and zero if and only if p = q. Traditionally it is interpreted to measure the amount of extra coding length needed when coding examples with codes generated for distribution q when the samples actually come from distribution p. We treat views symmetrically and compute the symmetrized divergence (D KL (p, q) + D KL (q, p))/2.

KL divergence is related to an information retrieval criterion: D KL (p, q) is the cost of misses in information retrieval of neighbors, when neighbors are retrieved using retrieval distribution q but they actually follow distribution p. D KL (p, q) is also the cost of false neighbors when neighbors are retrieved using p but they actually follow q. The relationships were shown in Venna et al. (2010) and used to compare a reduced-dimensional neighborhood to an original one; we use it in a novel way to compare neighborhoods across (transformed) views of data. The symmetrized divergence is the total cost of both misses and false neighbors when neighbors following the distribution in one transformed view are retrieved from the other transformed view with its neighbor distribution.

The value of the KL divergence can depend highly on differences between individual probabilities p( j) and q( j). A single missed neighbor can yield a high divergence value: for any index j if p( j) > ε for some ε > 0, D KL (p, q) → ∞ as q( j) → 0. In real-life multi-view data differences between views may be unavoidable, so we prefer a less strict measure focusing more on overall similarity of neighborhoods than severity of individual misses. We discuss such a measure below.

Angle cosine. A simple similarity measure between discrete distributions is the angle cosine between the distributions as vectors, that is, Cos(p, q) = (∑ j p( j)q( j))/ (∑ j (p( j))2 )(∑ j (q( j)) 2 ), which can be seen as the Pearson correlation coefficient between elements of p and q; it is thus a neighborhood correlation-a neighborhood based analogue of the coordinate correlation cost function of CCA.1 The angle cosine is bounded above and below: it has highest value 1 if and only if p = q and lowest value 0 if supports of p and q are nonoverlapping.

Similarity of neighborhoods by itself is not enough. The KL divergence and angle cosine (neighborhood correlation) measures only compare similarity of neighborhoods but not potential usefulness of the found subspaces. In high-dimensional data it is often possible to find subspaces where neighborhoods are trivially similar. For example, in data with sparse features it is often possible to find two dimensions where all data is reduced to a single value; in such dimensions neighborhood distributions would become uniform across all data since, hence any two such dimensions appear similar. To avoid discovering trivial similarities we wish to complement the measures of similarity between neighborhoods with terms favoring nontrivial (sparse) neighborhoods. A simple way to prefer sparse neighborhoods is to omit the normalization from neighborhood correlation, yielding

EQUATION

which is the inner product between the vectors of neighborhood probabilities. Unlike Cos(p, q), Sim(p, q) favors sparse neighborhoods: it has highest value 1 if and only if p = q and p( j) = q( j) = 1 for only one element j, and lowest value 0 if the supports of p and q are nonoverlapping.

The information retrieval interpretation is: Sim(p, q) is a proportional count of true neighbors from p retrieved from q or vice versa. If p has K neighbors with near-uniform high probabilities p( j) ≈ 1/K and other neighbors have near-zero probabilities, and q has L neighbors with high probability q( j) ≈ 1/L, then Sim(p, q) ≈ M/KL where M is the number of neighbors for which both p and q are high (retrieved true neighbors). Thus Sim(p, q) rewards matching neighborhoods and favors sparse neighborhoods (small K and L). One advantage of this formulation is to avoid matching two neighborhoods that seem to match only because they are highly uninformative: for example if p and q are both uniform over all neighbors, they have the same probability values and would be "similar" in a naive comparison of probability values, but both are actually simply uninformative about the choice of neighbors. Sim(p, q) would prefer a more sparse, more informative match, as desired.

2.2. FINAL COST AND OPTIMIZATION TECHNIQUE:

We wish to evaluate similarity of neighborhoods between subspaces of each view, and optimize the subspaces to maximize the similarity, while favoring subspaces having sparse (informative) neighborhoods for data items. We then evaluate similarities as Sim(p i,V , p i,U ) where p i,V = {p V ( j|i; f V )} is the neighborhood distribution around data item i in the dependent subspace of view V and f V is the mapping (parameters) of the subspace, and p i,U = {p U ( j|i; f U )} is the corresponding neighborhood distribution in the dependent subspace of view U having the mapping f U . As the objective function for finding dependent projections, we sum the above over each pair of views (U,V ) and over the neighborhoods of each data item i, yielding

C( f 1 , . . . , f N views ) = N views ∑ V =1 N views ∑ U=1,U =V N data ∑ i=1 N data ∑ j=1, j =i p V ( j|i; f V )p U ( j|i; f U ) (6)

where, in the setting of linear mappings and neighborhoods with Gaussian falloffs, p V is defined by (3) and is parameterized by the projection matrix W V of the linear mapping.

Optimization. The function C( f 1 , . . . , f N views ) is a well-defined objective for dependent projections and can be maximized with respect to the projection matrices W V of each view. We use gradient techniques for optimization, specifically limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS). Even with L-BFGS, ( 6) can be hard to optimize due to several local optima. To find a good local optimum, we optimize over L-BFGS rounds with a shrinking added penalty driving the objective away from the worst local optima during the first rounds; we use the optimum in each round as initialization of the next. For the penalty we use KL divergence based dissimilarity between neighborhoods, summed over neighborhoods of all data items i and view pairs (U,V ), giving

C Penalty ( f 1 , . . . , f N views ) = N views ∑ V =1 N views ∑ U=1,U =V N data ∑ i=1 (D KL (p i,V , p i,U ) + D KL (p i,U , p i,V )) 2 (7)

which is a function of all mapping parameters and can be optimized by L-BFGS; (7) penalizes severe misses (pairs (i, j) with nonzero neighborhood probability in one view but near-zero in another) driving the objective away from bad local optima. In practice KL divergence is too strict about misses; we use two remedies below.

Bounding KL divergence by neighbor distribution smoothing. To bound the KL divergence, one way is to give the neighbor distributions (1) a positive lower bound. In the spirit of the well-known Laplace smoothing in statistics, we revise the neighbor probabilities (1) as

p V ( j|i; f V ) = (exp(-d 2 V (i, j; f V )/σ 2 i,V ) + ε) • ( ∑ k =i exp(-d 2 V (i, k; f V )/σ 2 i,V ) + (N data -1)ε) -1 (8)

where ε > 0 is a small positive number. To keep notations simple, we still denote this smoothed neighbor distribution as p V ( j|i; f V ). To avoid over-complicated formulation and for consistency, we also use this version of neighbor probabilities in our objective function ( 6), even though the value of the objective is bounded by itself. We simply set ε = 1e -6 which empirically works well.

Shrinking the penalty. Even with bounded KL divergence, optimization stages need different amounts of penalty. At end of optimization, nearly no penalty is preferred, as views may not fully agree even with the best mapping. We shrink the penalty during optimization; the objective becomes

EQUATION

where γ controls the penalty. We initially set γ so the two parts of the objective function are equal for the initial mappings, C( f 1 , . . . , f N views ) = γC Penalty ( f 1 , . . . , f N views ), and multiply γ by a small factor (0.9 in experiments) at the start of each L-BFGS round to yield exponential shrinkage.

Time complexity. We calculate the neighbor distributions for all views, and optimize the objective function involving each pairs of views, thus the naive implementation takes O(dN 2 data N 2 views ) time, with d the maximal dimensionality among views. Acceleration techniques (Yang et al., 2013; Van Der Maaten, 2014; Vladymyrov & Carreira-Perpinán, 2014) from neighbor embedding could be adopted to reduce time complexity of a single view from O(N 2 data ) to O(N data log N data ) or even O(N data ). But scalability is not our first concern in this paper, so we use the naive O(N 2 data ) implementation for calculating the neighbor distributions for each view involved.

3. PROPERTIES OF THE METHOD AND EXTENSIONS:

Information retrieval. Our objective measures success in a neighbor retrieval task of the analyst: we maximize count of retrieved true neighbors across views, and penalize by severity of misses. Invariances. For any subspace of any view, (1) and (3) depend only on pairwise distances and are thus invariant to global translation, rotation, and mirroring of data in that subspace. The cost is then invariant to differences of global translation, rotation, and mirroring between views and finds view dependencies despite such differences. If in any subspace the data has isolated subsets (where data pairs from different subsets have zero neighbor probability) invariance holds for local translation/rotation/mirroring of the subsets as long as they preserve the isolation. Dependency is measured between whole subspaces. Unlike CCA where each canonical component of one view has a particular correlated pair in the other view, we maximize dependency with respect to the entire subspaces (transformed representations) of each view, as neighborhoods of data depend on all coordinates within the dependent subspace. Our method thus takes into account within-view feature dependencies when measuring dependency. Moreover, dependent subspaces do not need to be same-dimensional, and in some views we can choose not to reduce dimensionality but to learn a metric (full-rank linear transformation). Finding dependent neighborhoods between feature-based views and views external neighborhoods. In some domains, some data views may directly provide neighborhood relationships or similarities between data items, e.g., friendships in a social network, followerships in Twitter, or citations between scientific papers. Such relationships or similarities can be used in place of the feature-based neighborhood probabilities p V ( j|i; f V ) above. This shows an interesting similarity to a method (Peltonen, 2009) used to find similarities of one view to an external neighborhood definition; our method contains this task as one special case. Other falloffs. Exponential falloff in ( 1) and ( 3) can be replaced with other forms like t-distributed neighborhoods (van der Maaten & Hinton, 2008) . Such replacement preserves the invariances. Other transformations. Our criterion is extensible to nonlinear transformations in future work; replace linear projections by another parametric form, e.g. neural networks, optimize (9) with respect to its parameters; the transformation can be chosen on a view-by-view basis. Optimization difficulty of transformations varies; the best form of nonlinear transformation is outside the paper scope.

4. RELATED WORK:

In general, multi-view learning (Xu et al., 2013) denotes learning models by leveraging multiple potentially dependent data views; such models could be built either for unsupervised tasks based on the features in the views or for supervised tasks that involve additional annotations like categories of samples. In this paper we concentrate on unsupervised multi-view learning, and our prediction tasks of interest are predicting neighbors across views.

The standard Canonical Correlation Analysis (CCA) (Hotelling, 1936) iteratively finds component pairs maximizing correlation between data points in the projected subspaces. Such correlation is a simple restricted measure of linear and global dependency. To measure dependency in a more flexible way and handle nonlinear local dependency, linear and nonlinear CCA variants have been proposed: Local CCA (LCCA, Wei & Xu 2012) seeks linear projections for local patches in both views that maximize correlation locally, and aligns local projections into a global nonlinear projection; its variant Linear Local CCA (LLCCA) finds a linear approximation for the global nonlinear projection; Locality Preserving CCA (LPCCA, Sun & Chen 2007) maximizes reweighted correlation between data coordinate differences in both views. In experiments we compare to the well known traditional CCA and LPCCA as an example of recent state of the art.

As a more general framework, Canonical Divergence Analysis (Nguyen & Vreeken, 2015) minimizes a general divergence between distributions of data coordinates in the projected subspace.

The methods mentioned above work on data coordinates in the original spaces. There are also nonlinear CCA variants (e.g., Lai & Fyfe 2000; Bach & Jordan 2003; Verbeek et al. 2003; Andrew et al. 2013; Wang et al. 2015; Hodosh et al. 2013) for detecting nonlinear dependency between multiple views. Although some variants above are locality-aware, they introduce locality from the original space before maximizing correlation or other similarity measures in the low-dimensional subspaces. Since locality in the original space may not reflect locality in the subspaces due to noise or distortions, such criteria may not be suited for finding local dependencies in subspaces.

The CODE method (Globerson et al., 2007) creates an embedding of co-occurrence data of pairs of original categorical variables, mapping both variables into a shared space. Our method is not restricted to categorical inputs -its main applicability is to high-dimensional vectorial data, with several other advantages. In contrast to CODE, we find dependent subspaces (mappings) from multiple high-dimensional real-valued data views. Instead of restricting to a single mapping space we find several mappings, one from each view, which do not need to go into the same space; our output spaces can even be different dimensional for each view. Unlike CODE our method is not restricted to maximizing coordinate similarity: we only need to make neighborhoods similar which is more invariant to various transformations.

The above methods and several in Xu et al. (2013) all maximize correlation or alternative dependency measures between data coordinates across views. As we pointed out, in many domains coordinates are not of main interest but rather the data relationships they reveal; we consider neighborhood relationships and our method directly finds subspaces having similar neighborhoods.

5. EXPERIMENTS:

We demonstrate our method on artificial data with multiple dependent groups between views, and three real data sets: a variant of MNIST digits (LeCun & Cortes, 2010) , video data, and stock prices. In this paper we restrict our method to find linear subspaces, important in many applications for interpretability, and compare with two prominent linear subspace methods for multiple views, CCA and LPCCA. To our knowledge, no other information retrieval based approaches for finding linear subspaces is known up to the time when we did the experiment. Future work could compare methods for nonlinear mappings (Xu et al., 2013) to variants of ours for the same mapping; we do not focus on the mapping choice, and focus on showing the benefit or our neighborhood based objective.

On the artificial data set, we measure performance by correspondence between found projections and the ground truth. On the real data we use mean precision-mean recall curves, a natural performance measure for information retrieval tasks, and a measure for dependency as argued in Section 2.

Experiment on artificial data sets. We generate an artificial data set with 2 views with multiple dependent groups in each pair of corresponding dimensions as follows. Denote view V (∈ {1, 2}) as X (V ) ∈ R 5×1000 , and its i-th dimension as X (V ) i . For each i, we divide the 1000 data points in that dimension into 20 groups {g i j } 20 j=1 with 50 data points each. For each g i j and view V , we let x(V)

i jk F i j m (V ) i j + ε i jk (1 ≤ k ≤ 50), with m (V )

i j ∼ N (0, 5), ε i jk ∼ U[-0.5, 0.5] and F i j ∈ {-1, 1} a random variable allowing positive or negative correlation inside the group. We assemble x(V) i jk into X(V) ∈ R 5×1000 , and randomly permute entries of X(1) i and X(2) i in the same way but differently for different i, to ensure cross-dimension independency. Lastly we perform a PCA between X(1) i and X(2) i for each i, to remove cross-dimension correlation. We assemble the resulting

X (V ) i into X (V ) .

We pursue 2 transformations mapping from the 5D original space to a 1D latent space for each of the two views. Ground truth projections for both views will then be W (i) = (δ i j ) 5 j=1 ∈ R 1×5 . Results are in Fig. 1 : compared with CCA, our method successfully finds one of the ground truth transformations (= the 5th one), despite mirroring and scale, recovering the between-view dependency.

We measure performance by correspondence between the found projections and the ground truth transformation: let W 1 ,W 2 ∈ R 1×5 be projections found by a method, define

EQUATION

as the correspondence score. High score means good match between the found projections and ground truth. We repeat the experiment calculating correspondence on 20 artificial data sets generated in the same way. The table in Figure 1 summarizes the statistics. Our method outperforms CCA and LPCCA (with k = 5), finding the dependency on all 20 data sets. Experiment on real data sets with two views. We show our method helps match neighbors between the subspaces of two views after transformation. We use the following 3 data sets.

MNIST handwritten digit database (MNIST). MNIST contains 70000 gray-scale hand-written digit images of size 28 × 28. We create a training set and a testing set with 2000 images each. In the training set, we randomly choose 200 images from each digit to balance the distribution, while the testing set is another 2000 random images without balancing. We apply nonlinear dimensionality algorithm on both the left half and the right half of the images to create the two views to simulate a scenario where views have nonlinear dependency; we use Neighbor Retrieval Visualizer (NeRV, Venna et al. 2010) embedding to 5 dimensions with λ = 0.1 and λ = 0.9 respectively. The experiment is repeated 17 times, covering 68000 images.

Image patches from video (Video). We take image patches from a subset of frames in a video (db.tt/rcAS5tII). Starting from frame 50, we take 5 consecutive frames as a short clip at every 50 frames until frame 5200, then create two views from image patches in two fixed rectangles in those frames, rect 1 = [800, 900] × [250, 350] and rect 2 = [1820, 1920] × [800, 900]. We measure 5-fold cross-validation performance after randomly permuting the clips.

Stock prices (Stock), from the Kaggle "Winton stock market challenge" (goo.gl/eqdhKK). It contains prices of a stock at different times. We split the time series in the given training set into two halves, and let view 1 be the amplitudes from the Fourier transform results of the first half, and view 2 be the phases from the Fourier transform results of the second half.

For each data set we seek a pair of transformations onto 2D subspaces for the views. We measure performance by mean precision-mean recall curves of neighbor retrieval between 1) the two subspaces from the transformations, and 2) one of the original views and the subspace from the transformation for the other view. The better the performance is, the more to the top (better mean precision) and right (better mean recall) the curve will be in the figure. We set the number of neighbors in the ground truth as 5 for MNIST and Stock, 4 for Video, and let the number of retrieved neighbors vary from 1 to 10 as we focus on the performance of the matching for the nearest neighbors. We compare the performance with CCA and LPCCA. Figure 2 (column 1-3) shows the results.

We now show our method can find dependent subspaces for multiple (more than two) views. In this experiment we use Cell-Cycle data with five views. The views are from different measurements of cell cycle-regulated gene expression for the same set of 5670 genes (Spellman et al., 1998) . We preprocess data as in Tripathi et al. (2008) with an extra feature normalization step. We seek five two-dimensional subspaces from the five views, comparing to the PCA baseline with 2 components. We again use mean precision-mean recalls curves as the performance measure, additionally average the curves across the 10 view pairs or view-transformed coordinate pairs, besides averaging over the five repetitions in 5-fold cross-validation. Figure 2 (column 4) shows we outperform the baseline.

Finding subspaces with different dimensions. We show our method can find dependent subspaces with different dimensions. We create three two-dimensional Lissajous curves L k (t) = (cos √ 2k -1t + 2π(k -1)/3, cos √ 2kt + 2π(k -1)/3), k = 1, 2, 3. We create the first view X (1) ∈ R 6×1000 as X

(1)

1,1:1000 = (0, • • • , 999) and X

(1) d≥2,1:1000 i.i.d

∼ N (0, 1), and the second view X (2) ∈ R 6×1000 as the concatenation of the coordinates in the Lissajous curves. We seek a one-dimensional subspace Figure 2 : Mean precision-mean recalls curves from different real data on the test sets. Top row: view 1 as the ground truth; bottom row: subspace from view 1 as the ground truth. We can see curves from our method are to the top and/or right of the curves from other methods in most parts of all sub-figures, meaning our method achieves better precision and recall on average. Column 1-3: our method outperforms CCA and LPCCA; column 4: our method outperforms PCA. from X (1) , and a two-dimensional subspace from X (2) ; the aim is to find the nonlinear dependency between one-dimensional timestamps, and a two-dimensional representation for the three trajectories summarizing the two-dimensional movements of the three points along Lissajous curves. Figure 3 shows the Lissajous curves, found subspaces, and optimized projection pair. Our method successfully finds the informative feature in X (1) , and keeps transformed coordinates of X (2) smooth, with roughly the same amount of "quick turns" as in original Lissajous curves. The magnitudes in the optimized projections also suggest they capture the correct dependency. 

:

To make the connection exact, typically correlation is computed after substracting the mean from coordinates; for neighbor distributions of n data items, the mean neighborhood probability is the data-independent value 1/(n -1

) 2 which can be substracted from each sum term if an exact analogue to correlation is desired.

CONCLUSIONSWe presented a novel method for seeking dependent subspaces across multiple views, preserving neighborhood relationships of data. It has strong invariance properties, detects nonlinear dependencies, is related to an information retrieval task of the analyst, and performs well in experiments.

Figures and Tables:

Figure 3: Lissajous curves (a) and found subspaces from our method. (b) -(d) show we find the correct dependency: (b): perfect linear correlation shows the time dimension was found. (c): the number of "quick turns" (14 in total) in the smooth curves roughly matches that in the original curves. (d): projection weights, darker color means smaller magnitude; high magnitude of P (1) 's first entry and the complementary pattern in P (2) suggest we capture the dependency correctly.

Result for artificial data with dependent groups. Figures (left to right): one of the ground truths; 1D subspace pair recovered by our method; 1D subspace pair recovered by CCA. Our method recovers the dependency between views in the 5th dimension despite mirroring and scale differences; CCA fails to do so. Table:means and standard deviations of the correspondence measure (10); our method outperforms CCA and LPCCA, recovering the dependency in all artificial data sets.

Table content: <table><tr><td/><td colspan="2">Transformed coordinates found</td><td/><td colspan="3">Transformed coordinates</td></tr><tr><td>transformed coordinate of view 2</td><td colspan="2">by the proposed method</td><td>-0.10 0.10 -0.05 0.00 0.05 transformed coordinate of view 2</td><td colspan="2">found by CCA</td><td>Our method CCA LPCCA</td><td>Mean 1.00 0.51 0.51</td><td>Std 0.00 0.043 0.043</td></tr><tr><td/><td>-5000</td><td>0</td><td>-0.10</td><td>-0.05</td><td>0.00</td><td>0.05</td></tr><tr><td/><td colspan="2">transformed coordinate of view 1</td><td colspan="4">transformed coordinate of view 1</td></tr><tr><td>Figure 1:</td><td/><td/><td/><td/><td/><td/></tr></table>

Peer review (iclr2017):
This paper proposes a multiview learning approach to finding dependent subspaces optimized for maximizing cross-view similarity between neighborhoods of data samples. The motivation comes from information retrieval tasks. Authors position their work as an alternative to CCA-based multiview learning; note, however, that CCA based techniques have very different purpose and are rather broadly applicable than the setting considered here. Main points: - I am not sure what authors mean by time complexity. It would appear that they simply report the computational cost of evaluating the objective in equation (7). Is there a sense of how many iterations of the L-BFGS method? Since that is going to be difficult given the nature of the optimization problem, one would appreciate some sense of how hard or easy it is in practice to optimize the objective in (7) and how that varies with various problem dimensions. Authors argue that scalability is not their first concern, which is understandable, but if they are going to make some remarks about the computational cost, it better be clarified that the reported cost is for some small part of their overall approach rather than “time complexity”. - Since authors position their approach as an alternative to CCA, they should remark about how CCA, even though a nonconvex optimization problem, can be solved exactly with computational cost that is linear in the data size and only quadratic with dimensionality even with a naive implementation. The method proposed in the paper does not seem to be tractable, at least not immediately. - The empirical results with synthetic data are a it confusing. First of all the data generation procedure is quite convoluted, I am not sure why we need to process each coordinate separately in different groups, and then permute and combine etc. A simple benchmark where we take different linear transformations of a shared representation and add independent noise would suffice to confirm that the proposed method does something reasonable. I am also baffled why CCA does not recover the true subspace - arguably it is the level of additive noise that would impact the recoverability - however the proposed method is nearly exact so the noise level is perhaps not so severe. It is also not clear if authors are using regularization with CCA - without regularization CCA can be have in a funny manner. This needs to be clarified.