On Orthogonal Projections for Dimension Reduction and Applications in Augmented Target Loss Functions for Learning Problems

The use of orthogonal projections on high-dimensional input and target data in learning frameworks is studied. First, we investigate the relations between two standard objectives in dimension reduction, preservation of variance and of pairwise relative distances. Investigations of their asymptotic correlation as well as numerical experiments show that a projection does usually not satisfy both objectives at once. In a standard classification problem, we determine projections on the input data that balance the objectives and compare subsequent results. Next, we extend our application of orthogonal projections to deep learning tasks and introduce a general framework of augmented target loss functions. These loss functions integrate additional information via transformations and projections of the target data. In two supervised learning problems, clinical image segmentation and music information classification, the application of our proposed augmented target loss functions increases the accuracy.


Introduction
Linear dimension reduction is commonly used for preprocessing of high-dimensional data in complicated learning frameworks to compress and weight important data features. In contrast to nonlinear approaches, the use of orthogonal projections is computationally cheap, since it corresponds to a simple matrix multiplication. Conventional approaches apply specific projections that preserve essential information and complexity within a more compact representation. The projector is usually selected by optimizing distinct objectives, such as information preservation of the sample variance or of pairwise relative distances. Widely used orthogonal projections for dimension reduction are variants of the prin-cipal component analysis (PCA) that maximize the variance of the projected data [37]. Preservation of relative pairwise distances asks for a near-isometric embedding, and random projections guarantee this embeddings with high probability, cf. [5,15] and see also [1,6,12,27,30,35]. The use of random projections is especially favorable for large, highdimensional data [48], since the computational complexity is just O(dkm), e.g., using the construction in [1], with d, k ∈ N being the original and lower dimensions and m ∈ N the number of samples. In contrast, PCA needs O(d 2 m) + O(d 3 ) operations [24]. Moreover, tasks that do not have all data available at once, e.g., data streaming, ask for dimension reduction methods that are independent of the data.
In the present manuscript, we study orthogonal projections regarding the interplay between (O1) preservation of variance, (O2) preservation of pairwise relative distances, aiming for a sufficient lower-dimensional data representation. We shall consider the Euclidean distance exclusively since it is most widely used in applications, especially for error estimation. On manifolds, the geodesic distance is locally equivalent to the Euclidean distance. The two objectives (O1) and (O2) are directly addressed by PCA (O1) and random projections (O2). We achieve the following goals: First, we clarify mathematically and numerically that the two objectives are competing, i.e., PCA and random projections preserve different kinds of information. Depending on the objectives, we discuss beneficial choices of orthogonal projections and numerically find a balancing projector for a given data set. Finally, we define a general framework of augmented target (AT) loss functions for deep neural networks that integrate information about target characteristics via features and projections. We observe that our proposed methodology can increase the accuracy in two deep learning problems.
In contrast to conventional approaches, we study the joint behavior of the two objectives with respect to the entire set of orthogonal projectors. By analyzing the correlation between the variance and pairwise relative distances of projected data, we observe that (O1) and (O2) are competing and usually cannot be reached at the same time. In classification experiments with support vector machine and shallow neural networks, we investigate heuristic choices of projections applied to the input features.
In view of learning frameworks, we utilize features and projections on target data. The class of augmented target loss functions incorporates suitable transformations and projections that provide beneficial representations of the target space. It is applied in two supervised deep learning problems dealing with real-world data.
The first experiment is a clinical image segmentation problem in optical coherence tomography (OCT) data of the human retina. Related principles of dimension reduction for other clinical classification problems in OCT have already been successfully applied in [9]. In the second experiment, we aim to categorize musical instruments based on their spectrogram; see [19] for related results. Our utilized augmented target loss functions can increase the accuracy in both experiments.
The outline is as follows. In Sect. 2, we address the analysis of the competing objectives and Theorem 2.5 yields the asymptotic correlation between variance and pairwise relative distances of projected data. Section 3 prepares for the numerical investigations by recalling t-designs as considered in [10], enabling subsequent numerics. Heuristic investigations on projected input used in a straightforward classification task are presented in Sect. 4. Our framework of augmented target loss functions as modified standard loss functions for deep learning is introduced in Sect. 5. Finally, in Sects. 6 and 7 we present classification experiments on OCT images and musical instruments using aligned augmented target loss functions.

Dimension Reduction with Orthogonal Projections
To reduce the dimension of a high-dimensional data set denotes the Grassmannian, so that the lower-dimensional data representation is with range( p) = V . A suitable choice of p within G k,d depends on further objectives, i.e., which kind of information preservation shall be favored for subsequent analysis tasks. In the following, we consider two objectives associated with popular choices of orthogonal projectors for dimension reduction, in particular, random projectors and PCA. We will first observe that the two objectives are competing, especially in high dimensions, and then discuss consequences.

Objective (O1)
The total sample variance tvar(x) of x = {x i } m i=1 ⊂ R d is the sum of the corrected variances along each dimension: (2.2) PCA aims to construct p ∈ G k,d , such that the total sample variance of (2.1) is maximized among all projectors in G k,d . For other equivalent optimality criteria, we refer to [49]. The total sample variance of px = {px i } m i=1 ⊂ V coincides with the one of (2.1) and satisfies tvar( px) ≤ tvar(x) for all p ∈ G k,d . Thus, PCA achieves optimal variance preservation. The total variance (2.2) can also be expressed via pairwise absolute distances: Equally, it holds that which reveals that PCA maximizes the sample mean of the projected pairwise absolute distances.

Objective (O2)
In contrast to pairwise absolute distances, the Johnson-Lindenstrauss lemma targets the global property of preservation of pairwise relative distances: Lemma 2.1 (Johnson-Lindenstrauss, cf. [15,35]). For any holds for all i < j.
For small > 0, the projector p in Lemma 2.1 yields that all of the m(m−1) are close to 1, i.e., the projection p preserves all scaled pairwise relative distances well. A good choice of p in Lemma 2.1 is based on random projectors 1 P ∼ λ k,d , where λ k,d denotes the unique orthogonally invariant probability measure on G k,d . The following theorem is essentially proved by following the lines of the proof of Lemma 2.1 in [15] after replacing the constant 4 with (2 + τ )2 in the respective bound on k.
Theorem 2.2 For any 0 < < 1, any k ≤ d, m ∈ N and any 0 < τ with with probability at least 1 − 1 m τ + 1 m τ +1 . The theorem tells that the preservation property of pairwise relative distances is achieved with high probability when choosing a random projection according to k;d. Note that the random choice is completely independent from the actual data set.

Competing Objectives
A projector p satisfying the near-isometry property (2.5) implies so that the total variance of the projected data px is not preserved for k < d. In particular, with high probability a random projector P ∼ λ k,d does not suit the objective of maximizing the total variance, and we even observe E tvar(Px) = k d tvar(x); see (A.2) in the Appendix. PCA does not guarantee any local geometric property, and distances between pairs of points can be arbitrarily distorted [1]; see [39] for more robust PCA. The preservation of larger distances is favored since PCA maximizes (2.4) among all p ∈ G k,d and p(x i ) − p(x j ) ≤ x i − x j holds for all i < j. Close but distinct points could even be projected onto a single point, which violates the preservation of pairwise relative distances; see Fig. 1.
To more quantitatively understand the relation between the two competing objectives, we consider the sample mean (2.9) Recall that good preservation of the relative pairwise distances in (2.6) asks for M( p, x) being close to 1 and the variance V( p, x) being small. In the following, we analyze tvar( px), M( p, x) and V( p, x) and their expectations for random P ∼ λ k,d . In Fig. 2, we see a simple numerical experiment, where we first create an independent, normally distributed fixed data . . , m and m = 100, d = 50. We then compute PCA, for k = 10, 20, 30, 40, as well as n = 10,000 random projections p distributed according to λ k, 50 . In Fig. 2a-d, we can see that the more the k differs from d, the more the PCA and random projections differ concerning tvar( px) and M( p, x). Those differences may lead to diverse behavior in subsequent data analysis. Moreover, we compare M( p, x) and V( p, x) in Fig. 2e-h for the different k. We can see that again when k is much smaller than d, random projections and PCA differ more concerning the variance of pairwise distances V( p, x). For k = 10, the variance for PCA is higher in comparison with random projections (Fig. 2e); for k = 40 vice versa (Fig. 2h). Note that the theoretical bounds stated in Theorem 2.1 are much higher than the dimensions k used in the experiments, but the projections still preserve relative pairwise distances very well. In [7], similar observations were made on empiric experiments with image and text data.
The amount of variance kept in the principal components comparing real-world and random data has been experimentally studied, e.g., in [29] and [46]. Both studies determine that the difference occurs mainly in the first principal component.

Remark 2.3
In the numerical example, we compare random projections and PCA directly, serving as the corresponding projections to the objectives (O1) and (O2). We observe that even for not so high-dimensional (d = 50) data x and k d/2, PCA severely loses information in terms of total variance, i.e., more than 50% for k = 10, and more importantly loses much more information on pairwise relative distances than random projections. If both types of information are of interest, pairwise relative distances and high total variance, one should therefore favor random projections over PCA for k d/2 to balance the two objectives (O1) and (O2) and vice versa. Note that with a large amount of data one might still want to favor random projectors since their construction is computationally much cheaper and independent from the data. On the other hand, if objective (O2) is negligible, e.g., tasks with very noisy data, then PCA would be the favorable choice for all k.
Information of data can be quantified and expressed in different ways. One crucial part in dimension reduction is the decision of what kind of information shall be kept, which depends on several parameters including the quality of the data and the analysis task. Variants of PCA, focusing on the preservation of variance, have been widely used in real-world problems with big success, especially in denoising, when the preservation of all pairwise relative distances may be counterproductive, e.g., in dMRI imaging [51] and color filter array images [56]. Drawbacks are the necessity for all data being available from the start and the high computational costs. For very high-dimensional and large data sets, the computation of PCA is often not feasible. Besides the huge benefit of data independence and low computational cost when using random projections, the near-isometry property often allows to establish that the solution found in the low-dimensional space is a good approximation to the solution in the original space [1,34].
Algorithms in machine learning often need or benefit from sufficient estimates of pairwise distances, e.g., approximate nearest-neighbor problems, supervised classification [27] and subspace clustering [26]. In [32], algorithmic applications of near-isometry embeddings have been introduced. In [7], random projections have been successfully applied to noisy and noiseless text and image data. The experimental studies include the comparison of preservation of pairwise distances between random projections and PCA. The results coincide with our observations that for k > d/2 PCA is able to preserve the pairwise distances sufficiently, whereas for k < d/2 PCA distorts them. The smaller the k, the worse the distortion, whereas random projections preserve similarities still well for very small k, while being computationally much cheaper than PCA. One should point out again that favoring preservation of pairwise distances relies on the accuracy of the original distances.
PCA and random projections are orthogonal projections favoring two different aims. We want to study in the context of the whole set of orthogonal projections if the two objectives (O1) and (O2) could be reached at the same time. We will see that the objectives act competing, and therefore we suggest a balancing projector for tasks that benefit from both objectives.

Covariances and Correlation Between Competing Objectives
For further mathematical analysis, we first introduce a more general class of probability measures on G k,d that resemble λ k,d sufficiently well: where Pol t (R d 2 ) denotes the set of multivariate polynomials of total degree t in d 2 variables.
Existence of cubature measures is studied, for instance, in [17]. For random P, we now determine the expectation values for our three quantities of interest: tvar(Px), M(P, x) and V(P, x). If P ∼ λ and λ is a cubature measure of strength at least 2, the identities (A.2) and (A.3) in the Appendix and a short calculation yield This asymptotic bound relates to Theorem 2.2 and alludes to a near-isometry property of the type (2.7) for k sufficiently large.
The following theorem provides a lower bound for random P on the population correlation ) . (2.13) (2.14) are random points, whose entries are independent, identically distributed with finite 4-th moments that are uniformly bounded in d, then (2.14) converges towards 1 in probability for d → ∞.
The strong correlation for large dimensions d in the second part of Theorem 2.5 suggests that increasing tvar(Px) may also lead to increasing M(P, x); see Fig. 3 for illustration. Thus, large projected total variance tvar(Px) and the preservation of scaled pairwise distances, i.e., M(P, x) being close to 1, are competing properties. As discussed in Sect. 2.3, the choice of which kind of information is favorable to preserve depends on the data and the task, e.g., denoising (O1) and nearest-neighbor classification (O2). PCA and random projections are extreme in preserving either (O1) or (O2). We will heuristically study the behavior of orthogonal projections balancing both objectives in the next section and will state a numerical experiment where a balancing projector yields the highest classification accuracy.

Remark 2.6
The second part of Theorem 2.5 relates to the well-known fact that random vectors in high dimensions are almost orthogonal [4], and standard concentration of measure arguments may lead to more quantitative statements, cf. [52].

Optimal Covering Sequences
Let the covering radius of a set { p l } n l=1 ⊂ G k,d be denoted by where · F is the Frobenius norm. The smaller the covering radius, the better the set { p l } n l=1 represents the entire space G k,d , i.e., there are smaller holes and the points { p l } n l=1 are better distributed within G k,d . Following Lemma 2.1, we can connect finite sets of projections and their covering radius to the near-isometry property: Proof Given an arbitrary projector p ∈ G k,d , there is an index l 0 ∈ {1, . . . , n} such that From here, standard computations imply Lemma 3.1. We omit the details.
The accuracy of the near-isometry property in (3.2) depends on the covering radius. Therefore, a set { p l } n l=1 ∈ G k,d with a small covering radius is more likely to contain a projector with better preservation of pairwise relative distances. According to [11], it holds that 2 n − 1 k(d−k) and we shall see next how to achieve this lower bound.
A set of projectors { p l } n l=1 ⊂ G k,d is called a t-design if the associated normalized atomic measure 1 n n l=1 δ p l is a cubature measure of strength t (see Definition 2.4); see [44] for general existence results. Any sequence of t i -designs 3) 2 We use the symbols and to indicate that the corresponding inequalities hold up to a positive constant factor on the respective righthand side. The notation means that both relations and hold. and moreover, the bound n i t k(d−k) i holds, cf. [11,17]. To For their existence and numerical constructions, we refer to [21] and [10,11]. According to [11], see also (3.3) and (3.4), any low-cardinality design sequence { p i l } n i l=1 covers asymptotically optimal, i.e., Benefiting from the covering property, we will use lowcardinality design sequences as a representation of the overall space of orthogonal projectors G k,d .

Linear Least Squares Fit
With the linear least squares fit, we can directly gain information about the relation between M( p, x) and tvar( px) for a given data set where s and γ are determined by the sample variances and the sample covariance. If { p l } n l=1 is a 2-design, then the sample (co)variances coincide with the respective population (co)variances for P ∼ λ k,d ; see Appendix A.3 for further details. It follows that The quantities s and γ can be directly computed, where tvar(x) is given by (2.2) and the covariances are stated in Corollary A.1. Note that (3.6) and (3.7) are now independent of the particular choice of { p l } n l=1 . The correlation between the two samples (3.5) yields additional information about their relation. As before, if { p l } n l=1 is a 2-design, then the sample correlation coincides with the population correlation (2.13) for P ∼ λ k,d , cf. Appendix A.3. High correlation for a specific data set x suggests that random projections and PCA preserve competing properties, whose benefits need to be assessed for the specific subsequent task.

Numerical Experiments in Pattern Recognition
We investigate the impact on classification accuracy when applying specific orthogonal projections to input data. The chosen real-world data yields a straightforward classification task, serving as a toy example for comparing the accuracy of several projected input data in simple learning frameworks.

Choice of Orthogonal Projection
In the experiment, we use projections p ∈ G 2,4 reducing the original dimension from d = 4 to k = 2. As a finite representation of the overall space, we use a t-design of strength 14 from a low-cardinality sequence (see Sect. 3.1) consisting of 8475 orthogonal projectors. Note that the dimension reduction in practice takes place by applying q ∈ V k,d with q q = p ∈ G k,d , where V k,d := {q ∈ R k×d : qq = I k } denotes the Stiefel manifold. When taking norms, p and q are interchangeable, i.e., q(x) 2 = p(x) 2 , for all x ∈ R d . Therefore, we can use w.l.o.g. the theory developed for p.
The projections are chosen in a deterministic manner viewing the previously described competing properties. In

Results
In Fig. 4b, we see the linear least squares fitting line, computed directly and via the slope and intercept as stated in (3.6) and (3.7). The correlation coefficient (2.13) is 0.98, which suggests that preserving the two properties is highly competing and needs to be balanced.
In Table 1, the classification results of the iris data are presented. We can see that in this comparison the projector p ♦ , which corresponds to preserving M( p, x) ≈ 1 and maximizing tvar( px), yields the highest and most robust results. It even yields better results than working with the original input data. The projections that preserve M( p, x) ≈ 1 but do not take care of the magnitude of the total variance yield much worse results. On the other hand, the projections that just focus on high total variance still do not yield as high results as the projection p ♦ that balances both properties.

Augmented Target Loss Functions
In the previous section, projectors were applied to input features of shallow neural networks. In more complex architectures, such as deep neural networks, the adaption of weights can be viewed as optimization of input features, e.g., arising features can be used for transfer learning [54]. Whereas the input data are processed and optimized in each iteration, the target data stay usually unchanged during the whole learning process, serving as a measure of accuracy. The representation of the target data is one key property for successful approximation with neural networks. Here, we will introduce a general class of loss functions, i.e., augmented target (AT) loss functions, that use projections and features to yield beneficial representations of the target space, emphasizing important characteristics. In optimization problems, additional penalty terms are used for regularization or to enforce other constraints. In deep learning, weight decay (i.e., Tikhonov regularization) is a standard adaption of the loss function to that effect. Incorporating additional underlying information via features of the output/target data has been studied in diverse settings tailored to particular imaging applications. Perceptual loss functions have been used in [31] for image super-resolution, incorporating the comparison of high-level image features that arise from pretrained convolutional neural networks, i.e., the VGG network [45]. Deep perceptual similarity metrics have been proposed in [20] for generating images, comparing image features instead of the original images. In [28], a similar approach was successfully used for style transfer and super-resolution, adding a network that defines loss functions. Anatomically constrained neural networks (ACNN) have been introduced in [40] and applied to cardiac image enhancement and segmentation. Their loss functions incorporate structural information by using autoencoders to gain features about lower-dimensional parametrization of the segmentation. Brain segmentation was studied in [22], where information about the desired structure has been added in the loss function via an adjacency matrix. It was used for finetuning the supervised learned network with unlabeled data, reducing the number of abnormalities in the segmentation.
The information of certain target characteristics can be very powerful and even replace the need of annotations in some tasks. In [47], label-free learning is approached by using just structural information of the desired output in the loss function instead of annotated target values.
In the following, we will define a general framework of loss functions that add information of target characteristics via features and projections in supervised learning tasks.

General Framework
Let the training data be input vectors We consider training a neural network where θ ∈ R N corresponds to the vector of all free parameters of a fixed architecture. In each optimization step for θ , the network's output via an underlying loss function L. In contrast to ordinary learning problems with highly accurate target data, complicated learning tasks arising in many real-world problems do not yield sufficient results when optimizing neural networks with standard loss functions L, such as the widely used mean least squares error The training data may include important information that is obvious for humans, but poorly represented within the original target data and therefore lacks consideration in the learning process. To overcome this issue, we propose to add information tailored to the particular learning problem represented by additional features of the outputs and targets. First, we select transformations to enable error estimation in transformed output/target spaces. Note that the transformations T j are not required to be linear. However, they should be piecewise differentiable to enable subsequent optimization of the loss function with gradient methods. We shall allow for additional weighting of the transformations T 1 , . . . , T d to facilitate the selection of features for a specific learning problem. The previous sections suggest that orthogonal projections can provide favorable feature combinations, which essentially turns into a weighting procedure.
To enable suitable projections, we stack the d output/target features . . .
so that applying a projector p ∈ G k,d to each column of T (y i ) yields p(T (y i )) ∈ R d×t . We now define the augmented target loss function with projections by where α > 0 and L andL correspond to conventional loss functions. Apparently, L p depends on the choice of p ∈ G k,d . The projection p(T (y i )) weighs the previously chosen feature transformations T (y i ). Standard choices of L andL are L MSE , in which case L p becomes

Remark 5.1
For k = d, the projector p is the identity. In this case, the transformations can map onto different spaces, i.e., and we can now write the standard augmented target loss function by where T 1 corresponds to the identity function, L 1 , . . . , L d are common loss functions and α 1 , . . . , α d > 0 are weighting parameters.
It should be mentioned that α resembles a regularization parameter. The actual minimization of (5.1) among θ is usually performed through Tikhonov-type regularization in many standard deep neural network implementations. The formulation (5.2) adds one further variational step for beneficial output data representation.

Remark 5.2
Our proposed structure with target feature maps T 1 , . . . , T d as in (5.4) relates to multitask learning, which has been successfully used in deep neural networks [13]. It handles multiple learning problems with different outputs at the same time. In contrast to multitask learning, we aim to solve a single problem but also penalize the error in transformed spaces enhancing certain target characteristics.
For the projected feature transformations in the augmented target loss function, it is not possible to identify a balancing projection p heuristically (such as p ♦ in Sect. 4), because the output y changes in each iteration when the loss function is called. In the following clinical numerical experiment we overcome this issue by choosing random projections in each optimization step and compare it to prior deterministic choices of projections, including PCA.

Application to Clinical Image Data
The first experiment is a clinical problem in retinal image analysis of the human eye, where the disruptions of the socalled photoreceptor layers need to be quantified in optical coherence tomography (OCT) images. The photoreceptors have been identified as the most important retinal biomarker for prediction of vision from OCT in various clinical publications, see e.g., [23]. As OCT technology advances, clinicians are not able to look at each slice of OCT themselves. (In mean, they get 250 slices per patient and have 3-5 minutes/patients including their clinical examination.) Therefore, automated classification of, for example, photoreceptor status is necessary for clinical guidance.

Data and Objective
In this application, OCT images of different retinal diseases (diabetic macular edema and retinal vein occlusion) were provided by the Vienna Reading Center recorded with the Spectralis OCT device (Heidelberg Engineering, Heidelberg, Germany). Each patient's OCT volume consists of 49 cross sections/slices (496 × 512 pixels) recorded in an area of 6 × 6 mm in the center of the human retina, which is the part of the retina responsible for vision. Each of the slices was manually annotated by a trained grader of the reading center. This is a challenging and time-consuming procedure that is not feasible in clinical routine but only in a research setting. The binary pixelwise annotations serve as target values, enabling a supervised learning framework (Fig. 5).
The objective is to accurately detect the photoreceptor layers and their disruptions pixelwise in each OCT slice by training a deep convolutional neural network with a suitable loss function. The learning problem is complicated by potentially inaccurate target annotations, as studies have shown that inconsistencies between trained graders are common, cf. [50]. Moreover, the learning task is unbalanced in the sense that there are many more slices showing none or very little disruptions. We shall observe that optimization with respect to standard loss functions performs poorly in regards to detecting disruptions. The augmented target loss function proposed in the previous section can enhance the detection.

Convolutional Neural Network Learning
We implemented our experiments using Python 3.6 with Pytorch 1.0.0. A deep convolutional neural network f θ is trained by applying the U-Net architecture reported in [43] with a sigmoid activation function and Tikhonov regularization. A set of 20 OCT volumes (980 slices) from different patients with corresponding annotations are used for training, where four volumes were used for calibration (validation set). Another two independent test volumes were identified for evaluating the results, one without any disruptions in the photoreceptor layers, whereas the other one includes a high number of disruptions. Each OCT slice is represented by a vector x i ∈ R r with r = 496 · 512. The collection {x i } m i=1 corresponds to all slices from the training volumes, i.e., m = 20 · 49. Further matching the notation of the previous section, we have r = s and f θ : R r → R r with binary target vectors y i ∈ {0, 1} r . We observe that disruptions are not identified reliably when using the least squared loss function (5.1). To overcome this issues, we use the proposed augmented target loss function with least squared losses as stated in (5.3).
To enhance disruptions within the output/target space, we heuristically choose d = 4 local features of the original rep- Bold values correspond to the highest AUC value, which serves us as a measure of accuracy resentation. They are derived from convolutions with two edge filters, T 1 (Prewitt) and T 2 (Laplacian of Gaussian), and from two Gaussian high-pass filters, yielding T 3 and T 4 . Note that these feature transformations keep the same size, i.e., T j : R r → R r for j = 1, . . . , d. See Fig. 6 for example images.
We can derive several augmented target loss functions L p by choosing different p ∈ G k,d . In this experiment, we use the following projections: l=1 , all projections from a t-design of strength 2 ⊂ G 2,4 (see [10]), • p PCA ∈ G 2,4 , projection determined by PCA on the training data, • p λ 2,4 , random projection chosen according to λ 2,4 in each mini-batch.

Results
Since the detection problem is highly unbalanced, we use precision/recall curves [16] for evaluating the overall performance of each loss function model. The area under the curve (AUC) was used as a numerical indicator of the success rate [41]. The higher the AUC, the better the classification. The results of the different loss functions on the independent test set are stated in Table 2. Due to the imbalance within the data, the photoreceptor region is identified well, but disruptions are not identified reliably when using the least squared loss function (5.1). For α = 0.1, all proposed augmented target loss functions L p clearly increase the success rate of the disruption quantification. Note that all projections are independent from the actual data set, except PCA that was computed beforehand on the training data.
The features itself (i.e., p = I 4 ) improve the quantification, and weighting them by projections increases the results even more: using the fixed projection p 12 from the t-design sequence { p l } 15 l=1 on the output/target features yields the highest accuracy for photoreceptors and disruptions. This corresponds to the results of the previous sections, stating that depending on the particular data there are projections in the overall space acting beneficially. Since this projection generally cannot be found beforehand, using random projections in each loss function's evaluation step is easier, possible in practice, and independent from the data. The computation is efficient and randomization can have regularization effects that yield more robust results, cf. [34]. In the following, we will view a second classification problem based on spectrograms, where augmented target loss functions with random projections can improve the accuracy.

Application to Musical Data
Here, the learning task is a prototypical problem in Music Information Retrieval, namely multi-class classification of musical instruments. In analogy to the MNIST problem in image recognition, this classification problem is commonly used as a basis of comparison for innovative methods, since the ground truth is unambiguous and sufficiently many anno-tated data are available. The input to the neural network is spectrograms of audio signals, which is the standard choice in audio machine learning. Spectrograms are calculated from the time signal using a short-time Fourier transform and taking the absolute value squared of the resulting spectra, thus yielding a vector for each time step and a two-dimensional array, like an image, cf. [18].
Reproducible code and more detailed information of our computational experiments can be found in the online repository [25].

Data and Objective
The publicly available GoodSounds data set [42] contains recordings of single notes and scales played by several single instruments. To gain equally balanced input classes, we restrict the classification problem to six instruments: clarinet, flute, trumpet, violin, alto saxophone and cello. Note that the recordings are monophonic, so that each recording yields one spectrogram that we aim to correctly assign to one of the six instruments.
After removing the silence [3,38], segments from the raw audio files are transformed into log-mel spectrograms [36], so that we obtain images of time-frequency representations with size 100×100. One example spectrogram for each class of instruments is depicted in Fig. 7.

Convolutional Neural Network Learning
We implemented a fully convolutional neural network f θ : R r → [0, 1] s , cf. [33], where r = 100 × 100 and s = 6, in Python 3.6 using Keras 2.2.4 framework [14] and trained it on the Nvidia GTX 1080 Ti GPU. The data are split into 1, 40, 722 training, 36,000 validation and 36,000 independent test samples. We heuristically choose d = 16 output features arising directly from the particular output class. The transformations T 1 , . . . , T 16 , with T j : R 6 → R for j = 1, . . . , 16, are then given by the inner product of the output/target and the feature vectors. Among others, the features are chosen from the enhanced scheme of taxonomy [53] and from the table of frequencies, harmonics and under tones [57]. We use the proposed augmented target loss function L p (5.2), where L 1 corresponds to the categorical cross-entropy  The standard inbuilt Tikhonov regularization ( 2 -norm of θ) is weighted by β. For α > 0, the feature transformations {T j } 16 j=1 are used in the loss function, either directly or weighted by a random projection p λ6, 16 . The accuracy of the model is measured by the number of correctly classified samples divided by the number of all samples Bold values correspond to the highest classification accuracy loss [55] and L 2 to the mean squared error as in (5.3). We consider here two choices of p: the identity I 16 and random projectors p ∼ λ 6,16 in G 6,16 .
The deep learning model is sensitive to various hyperparameters, including α and p, in addition to conventional parameters, such as the number of convolutional kernels, learning rate and the parameter β for Tikhonov regularization. To find the best choices in a fair trial, we utilize a random hyper-parameter search approach, where we train 60 models and select the three best ones for a more precise search over different α in the augmented target loss function and β for Tikhonov regularization. This results in 212 models that are evaluated on the training and validation set. Finally, we select the best model based on the accuracy of the validation set and evaluate it on the independent test set. For comparison, we also evaluate this model with no Tikhonov regularization, i.e., β = 0; see Table 3. Table 3 shows that no regularization and no features provide the poorest results. It seems that adding features with random projections has a regularizing effect and improves the results significantly. As expected, it is important to include Tikhonov regularization on θ . Further enhancement happens by adding features via the modified augmented target loss function with or without additional weighting from projections. All results are very stable and are generalizing very well from training to the independent test set; see [25] for further details. For fixed parameters μ > 0, σ 2 > 0, that do not depend on d, let Y 1 ∈ R d be a random vector, whose squared entries are independent, identically distributed with mean EY 2 1,l = μ and variance Var(Y 2 1,l ) = σ 2 , for l = 1, . . . , d. We immediately observe

Results
For any c > 0, Chebyshev's inequality yields

M. Dörfler obtained her PhD in
Mathematics from the University of Vienna and is a research assistant at the Mathematics Department of the University of Vienna. She is heading her own research group and is interested in the interplay between aspects of harmonic analysis and the design of learning algorithms based on convolutional neural networks. She has worked on audio signal processing and, being a musician herself, has been particularly interested in developing representations adapted to music applications.