Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
11institutetext: School of Mathematics, Sun Yat-Sen University, China 22institutetext: School of Mathematics, Jiaying University, China
22email: {yanghr26@mail2, rchuanx@mail, luoyw28@mail2}.sysu.edu.cn

COD: Learning Conditional Invariant Representation for Domain Adaptation Regression

Hao-Ran Yang 11    Chuan-Xian Ren\orcidlink0000-0002-1861-3599 Corresponding author.11    You-Wei Luo\orcidlink0000-0002-3027-6679 1122
Abstract

Aiming to generalize the label knowledge from a source domain with continuous outputs to an unlabeled target domain, Domain Adaptation Regression (DAR) is developed for complex practical learning problems. However, due to the continuity problem in regression, existing conditional distribution alignment theory and methods with discrete prior, which are proven to be effective in classification settings, are no longer applicable. In this work, focusing on the feasibility problems in DAR, we establish the sufficiency theory for the regression model, which shows the generalization error can be sufficiently dominated by the cross-domain conditional discrepancy. Further, to characterize conditional discrepancy with continuous conditioning variable, a novel Conditional Operator Discrepancy (COD) is proposed, which admits the metric property on conditional distributions via the kernel embedding theory. Finally, to minimize the discrepancy, a COD-based conditional invariant representation learning model is proposed, and the reformulation is derived to show that reasonable modifications on moment statistics can further improve the discriminability of the adaptation model. Extensive experiments on standard DAR datasets verify the validity of theoretical results and the superiority over SOTA DAR methods.

Keywords:
Domain Adaptation Regression, Optimal Transport, Invariant Representation Learning, Kernel Embedding Theory.

1 Introduction

Deep learning has shown a promising ability for fitting and achieved great success in various application tasks. However, large-scale data with sufficient annotations plays a vital role in training high-quality deep models while the label annotations for unlabeled data are usually expensive. Therefore, aiming to reduce training costs by transferring knowledge, Unsupervised Domain Adaptation (UDA) [14, 28, 17, 2] has received increasing attention in recent years. Basically, UDA aims to transfer task-related information from a labeled source domain to a target domain without annotations, where the domain gap is characterized as the distribution shift across domains [43]. Following this line, considerable methods have been proposed to alleviate the UDA problem by shift correction and distribution alignment, such as adversarial training for cross-domain samples [14, 52] and explicitly distribution matching by statistical discrepancy minimization [28, 49]. From a theoretical perspective, recent literature [60, 6] show that the marginal distribution PXsubscript𝑃𝑋P_{X}italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT with covariate shift assumption is insufficient to characterize the cross-domain generalization error, while the conditional distribution alignment indeed serves as a sufficient condition to ensure small error and successful transfer. Mathematically, when the conditional distributions PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT for every y𝑦yitalic_y are aligned, a well-trained source predictor can perform well on the target domain, which is described in Fig. 1. Inspired by theoretical advances, recent works [41, 45, 23] consider the statistical methods for conditional alignment, and show promising performance in mitigating negative transfer and improving recognition accuracy. However, most of the existing results are designed for classification settings with discrete prior on label variables, which induces the infeasibility for regression problems. Consequently, it restricts the application of UDA methods for wider real-world problems, e.g., human pose estimation [36], product condition monitoring[27], depth estimation [42].

Considering the dilemma of UDA, pioneer works have explored the Domain Adaptation Regression (DAR) framework for solving the distribution shift in continuous label scenarios, where theoretical results [35, 7] and correction methods [15, 51, 47] are provided. To further relax the requirement of labeled training target data in pioneer DAR works, advanced works[4, 38] propose unsupervised DAR methods to learn transferable regression models without the help of labeled target data. Nonetheless, existing DAR theory and methods are directly facilitated by the covariate shift framework, where the marginal distribution alignment is considered. Therefore, it inevitably induces the negative transfer and misaligned local structure, which are shown in the results in classification settings [60, 6]. In fact, there is a significant obstacle for modeling conditional distribution alignment in regression scenarios, i.e., the conditional shift correction induces infinite matching problems for PX|Y=ysubscript𝑃conditional𝑋𝑌𝑦P_{X|Y=y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y = italic_y end_POSTSUBSCRIPT on all points y𝑦yitalic_y. The continuous label variable with infinite possible outcomes makes it infeasible to directly apply the existing conditional alignment methods for DAR. Thus, it is vital to understand the continuity problem in regression, and develop the generalization error theory and statistical methodology for DAR scenarios.

Refer to caption
Figure 1: Illustration of the conditional shift in regression setting, where label value y[0,1]𝑦01y\in[0,1]italic_y ∈ [ 0 , 1 ] and color gradients imply the continuity. (a) Before adaptation, distribution shift exists between source and target representations so the predictor trained on source domain cannot generalize to target domain. (b) After marginal distribution alignment, distributions of source and target domain are globally aligned. However, representations with different labels may be falsely aligned across domains, leading to significant inconsistency between cross-domain labeling rules. Thus, level sets provided by the source predictor are not suitable for target representations. (c) After conditional distribution alignment, label-wise matching is achieved, where the level sets of cross-domain labeling rules are consistent and the source predictor will be reliable.

In this work, we focus on the theory guarantee and conditional alignment methodology for successful DAR. Theoretically, we explore the generalization upper bound for continuous output scenarios and show the cross-domain joint error can be sufficiently upper-bounded by the conditional discrepancy. It ensures the sufficiency of conditional alignment still holds for regression scenarios and provides reliable guarantees for theory-driven modeling. Methodologically, to deal with the continuous problem in conditional alignment, we treat the conditional distribution as finite statistical moments in Hilbert space. By exploring the distribution embedding theory of Reproducing Kernel Hilbert Space (RKHS), a novel conditional metric called Conditional Operator Discrepancy (COD) is proposed, which can be interpreted as the conditional Wasserstein distance in RKHS. With the sufficiency theory and metric-based learning principle for regression, a COD-based conditional invariant representation learning method is proposed, which aims to correct the conditional shift on the pushforward measures in latent representation space. By reformulating the empirical COD estimation, a modified principle is developed to further improve the discriminability during knowledge transfer. Generally, our contributions can be summarized as follows:

  • The sufficiency theory for the conditional shift framework is established for DAR settings, which shows the cross-domain conditional discrepancy on continuous label variable serves as the major term in generalization error.

  • A novel COD metric for characterizing conditional distribution discrepancy with continuous and infinite conditioning variable is explored, providing a feasible solution for conditional discrepancy optimization in DAR.

  • A COD-based representation learning model with reformulation on moment statistics is proposed, which admits discriminability on task-related knowledge and shows superior performance in extensive experiments.

2 Related Work

Domain Adaptation Classification. Based on the distribution shift assumption, most of the UDA methods can be roughly categorized into three types: methods that align marginal distribution, conditional distribution and joint distribution. Most UDA methods consider only the marginal distribution alignment. Discrepancy-based methods attempt to directly minimize the distribution discrepancy and different types of discrepancy measurement greatly enrich this type of methodology. For example, these methods adopt moment matching [28, 49] or extension under manifold framework [17, 33]. In particular, Optimal transport (OT) and Wasserstein distance are introduced to UDA as a sound measurement for domain shift, based on which a large number of methods have been proposed [9, 8, 10, 46, 24]. Zhang et al. [59] further develop both theory and methodology to RKHS and Luo et al. [31] extends Kernel Bures metric to conditional distribution alignment. Adversarial-based methods [14, 20, 21] learn domain-invariant features via two-player games. Specifically, they alternatively optimize the feature generator and the domain discriminator to achieve domain confusion. In order to obtain better performance, advanced methods have explored how to utilize the label information to achieve conditional distribution alignment [41, 45, 61, 50, 53, 3, 44, 32] or joint distribution alignment [30, 29]. However, UDA methods for shift correction are proposed for discrete and finite conditioning variables, which are inapplicable for regression.

Domain Adaptation Regression. DA for regression has received relatively little attention compared with classification. Mansour et al. [35] and Cortes et al. [7] introduce theoretical analysis for DA in regression scenarios. Zhao et al. [62] discuss the tradeoffs between accuracy and invariance in learning invariant representations in both classification and regression settings. A range of methods in the shallow regime are proposed to tackle the DAR problem. These methods can be roughly categorized into two types: importance weighting [15, 37] and learning domain-invariant representation [40, 8, 39]. Besides, some attempts have been made in deep regime to tackle regression problem [51, 47, 25, 56]. Unfortunately, most of them belong to semi-supervised DA, which means they need to use a small number of target domain labels. Recently, deep UDA methods for regression that get rid of target label are proposed [4, 38, 57, 11]. RSD [4] states that deep neural networks are less robust to feature scaling in regression problems and aligning the deep features directly will impede task performance for altering feature scale. Thus, a strategy of aligning orthogonal bases of both domains is adopted in RSD. Based on RSD, DARE-GRAM [38] aligns the inverse Gram Matrix of deep features of both domains, and ensures promising improvements in accuracy and stability. Meanwhile, several DAR methods for specific tasks have been explored [22, 55, 1, 18, 26, 27]. However, the theory and method for conditional shift correction are still unexplored for DAR scenarios, which hinders the in-depth reasoning and development for practical modeling.

3 COD: Theory and Methodology

3.1 Preliminary

Setup. In this paper, we discuss UDA method dealing with regression problems. Let X and Y be covariate and label variable taken values on 𝒳𝒳\mathcal{X}caligraphic_X and 𝒴𝒴\mathcal{Y}caligraphic_Y. \mathbb{P}blackboard_P denotes the set of distributions and subscripts represent the corresponding random variable while the superscript represents the domain, e.g. PXYs𝒳𝒴superscriptsubscript𝑃𝑋𝑌𝑠subscript𝒳𝒴P_{XY}^{s}\in\mathbb{P}_{\mathcal{XY}}italic_P start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ∈ blackboard_P start_POSTSUBSCRIPT caligraphic_X caligraphic_Y end_POSTSUBSCRIPT implies the joint distribution on source domain. In UDA settings, we are provided with labeled source samples (xs,ys)PXYssimilar-tosuperscript𝑥𝑠superscript𝑦𝑠superscriptsubscript𝑃𝑋𝑌𝑠(x^{s},y^{s})\sim P_{XY}^{s}( italic_x start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_y start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) ∼ italic_P start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT and unlabeled target samples (xt,yt)PXYtsimilar-tosuperscript𝑥𝑡superscript𝑦𝑡superscriptsubscript𝑃𝑋𝑌𝑡(x^{t},y^{t})\sim P_{XY}^{t}( italic_x start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT , italic_y start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) ∼ italic_P start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT, which means the target labels are not available during training. Samples of source and target domain follow different distributions, i.e. PXYsPXYtsuperscriptsubscript𝑃𝑋𝑌𝑠superscriptsubscript𝑃𝑋𝑌𝑡P_{XY}^{s}\neq P_{XY}^{t}italic_P start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ≠ italic_P start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT. Denote g:XZ:𝑔𝑋𝑍g:X\to Zitalic_g : italic_X → italic_Z as the representation extractor, h:ZY:𝑍𝑌h:Z\to Yitalic_h : italic_Z → italic_Y as the predictor, εssubscript𝜀𝑠\varepsilon_{s}italic_ε start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and εtsubscript𝜀𝑡\varepsilon_{t}italic_ε start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT as the source and target prediction error, respectively. UDA methods aim to generalize the learner hg𝑔h{\circ}gitalic_h ∘ italic_g across domains and minimize target risk εtsubscript𝜀𝑡\varepsilon_{t}italic_ε start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT by mitigating the distribution shift problem. In contrast to the discrete label space in classification settings, unsupervised DAR discusses UDA problems with continuous labels, which means the number of classes is uncountable, leading to great difficulty in aligning conditional distribution PX|Ysubscript𝑃conditional𝑋𝑌P_{X|Y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT. The proofs are provided in supplementary materials.

OT in RKHS. Let (𝒳,)𝒳(\mathcal{X},\mathcal{B})( caligraphic_X , caligraphic_B ) be a measure space with Borel σlimit-from𝜎\sigma-italic_σ -field \mathcal{B}caligraphic_B. Denote (𝒳,k𝒳)subscript𝒳subscript𝑘𝒳(\mathcal{H}_{\mathcal{X}},k_{\mathcal{X}})( caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ) as the RKHS of 𝒳𝒳\mathcal{X}caligraphic_X, which is generated by the positive definite kernel k𝒳subscript𝑘𝒳k_{\mathcal{X}}italic_k start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT. Elements from 𝒳𝒳\mathcal{X}caligraphic_X are mapped to 𝒳subscript𝒳\mathcal{H}_{\mathcal{X}}caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT via feature map ϕ(x)=k𝒳(x,)italic-ϕ𝑥subscript𝑘𝒳𝑥\phi(x)=k_{\mathcal{X}}(x,\cdot)italic_ϕ ( italic_x ) = italic_k start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ( italic_x , ⋅ ), which is assumed to satisfy the reproducing properties ϕ(x),ϕ(x)𝒳=k𝒳(x,x)subscriptitalic-ϕ𝑥italic-ϕsuperscript𝑥subscript𝒳subscript𝑘𝒳𝑥superscript𝑥\left<\phi(x),\phi(x^{\prime})\right>_{\mathcal{H}_{\mathcal{X}}}=k_{\mathcal{% X}}(x,x^{\prime})⟨ italic_ϕ ( italic_x ) , italic_ϕ ( italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ start_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_k start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) and ϕ(x),f𝒳=f(x),f𝒳formulae-sequencesubscriptitalic-ϕ𝑥𝑓subscript𝒳𝑓𝑥for-all𝑓subscript𝒳\left<\phi(x),f\right>_{\mathcal{H}_{\mathcal{X}}}=f(x),\forall f\in\mathcal{H% }_{\mathcal{X}}⟨ italic_ϕ ( italic_x ) , italic_f ⟩ start_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_f ( italic_x ) , ∀ italic_f ∈ caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT. The mean element μXsubscript𝜇𝑋\mu_{X}italic_μ start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT in 𝒳subscript𝒳\mathcal{H}_{\mathcal{X}}caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT with law PXsubscript𝑃𝑋P_{X}italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT is given by μX=𝔼X[ϕ(X)]subscript𝜇𝑋subscript𝔼𝑋delimited-[]italic-ϕ𝑋\mu_{X}=\mathbb{E}_{X}\left[\phi(X)\right]italic_μ start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT = blackboard_E start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT [ italic_ϕ ( italic_X ) ], and the covariance operator 𝒞XXsubscript𝒞𝑋𝑋\mathcal{C}_{XX}caligraphic_C start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT is defined as 𝒞XX=EX(ϕ(X)μX)(ϕ(X)μX)subscript𝒞𝑋𝑋tensor-productsubscript𝐸𝑋italic-ϕ𝑋subscript𝜇𝑋italic-ϕ𝑋subscript𝜇𝑋\mathcal{C}_{XX}=E_{X}(\phi(X)-\mu_{X})\otimes(\phi(X)-\mu_{X})caligraphic_C start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT = italic_E start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( italic_ϕ ( italic_X ) - italic_μ start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ) ⊗ ( italic_ϕ ( italic_X ) - italic_μ start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ).

For any two distributions PXs,PXt𝒳superscriptsubscript𝑃𝑋𝑠superscriptsubscript𝑃𝑋𝑡subscript𝒳P_{X}^{s},P_{X}^{t}\in\mathbb{P}_{\mathcal{X}}italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ∈ blackboard_P start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT, the Kantorovitch optimal transport in RKHS 𝒦subscript𝒦\mathcal{H}_{\mathcal{K}}caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT is considered on pushforward measures ϕ#PXssubscriptitalic-ϕ#superscriptsubscript𝑃𝑋𝑠\phi_{\#}P_{X}^{s}italic_ϕ start_POSTSUBSCRIPT # end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT and ϕ#PXtsubscriptitalic-ϕ#superscriptsubscript𝑃𝑋𝑡\phi_{\#}P_{X}^{t}italic_ϕ start_POSTSUBSCRIPT # end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT:

dKW2(Xs,Xt)=infπ𝒦Π(ϕ#PXs,ϕ#PXt)𝒦×𝒦XsXt𝒦2𝑑π𝒦(Xs,Xt)subscriptsuperscript𝑑2KWsuperscript𝑋𝑠superscript𝑋𝑡subscriptinfimumsubscript𝜋𝒦Πsubscriptitalic-ϕ#superscriptsubscript𝑃𝑋𝑠subscriptitalic-ϕ#superscriptsubscript𝑃𝑋𝑡subscriptsubscript𝒦subscript𝒦superscriptsubscriptnormsuperscript𝑋𝑠superscript𝑋𝑡subscript𝒦2differential-dsubscript𝜋𝒦superscript𝑋𝑠superscript𝑋𝑡\displaystyle d^{2}_{\mathrm{KW}}(X^{s},X^{t})=\inf_{\pi_{\mathcal{K}}\in\Pi(% \phi_{\#}P_{X}^{s},\phi_{\#}P_{X}^{t})}\int_{\mathcal{H}_{\mathcal{K}}\times% \mathcal{H}_{\mathcal{K}}}\left\|X^{s}-X^{t}\right\|_{\mathcal{H}_{\mathcal{K}% }}^{2}d\pi_{\mathcal{K}}(X^{s},X^{t})italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_KW end_POSTSUBSCRIPT ( italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) = roman_inf start_POSTSUBSCRIPT italic_π start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT ∈ roman_Π ( italic_ϕ start_POSTSUBSCRIPT # end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_ϕ start_POSTSUBSCRIPT # end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT × caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT end_POSTSUBSCRIPT ∥ italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT - italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_π start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT ( italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) (1)

Zhang et al. [59] show that if the pushforward measures are Gaussian, the Kernel Gaussian Wasserstein (KGW) distance between them can be written as

dKGW2(Xs,Xt)superscriptsubscript𝑑KGW2superscript𝑋𝑠superscript𝑋𝑡\displaystyle d_{\mathrm{KGW}}^{2}(X^{s},X^{t})italic_d start_POSTSUBSCRIPT roman_KGW end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) =uXsuXt𝒦2+tr(𝒞XXss+𝒞XXtt2𝒞XXst)absentsubscriptsuperscriptnormsubscriptsuperscript𝑢𝑠𝑋subscriptsuperscript𝑢𝑡𝑋2subscript𝒦𝑡𝑟subscriptsuperscript𝒞𝑠𝑠𝑋𝑋subscriptsuperscript𝒞𝑡𝑡𝑋𝑋2subscriptsuperscript𝒞𝑠𝑡𝑋𝑋\displaystyle=\|u^{s}_{X}-u^{t}_{X}\|^{2}_{\mathcal{H}_{\mathcal{K}}}+tr(% \mathcal{C}^{ss}_{XX}+\mathcal{C}^{tt}_{XX}-2\mathcal{C}^{st}_{XX})= ∥ italic_u start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT - italic_u start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ∥ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_t italic_r ( caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT + caligraphic_C start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT - 2 caligraphic_C start_POSTSUPERSCRIPT italic_s italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT ) (2)

where 𝒞XXst=𝒞XXss𝒞XXtt𝒞XXsssubscriptsuperscript𝒞𝑠𝑡𝑋𝑋subscriptsuperscript𝒞𝑠𝑠𝑋𝑋subscriptsuperscript𝒞𝑡𝑡𝑋𝑋subscriptsuperscript𝒞𝑠𝑠𝑋𝑋\mathcal{C}^{st}_{XX}=\sqrt{\sqrt{\mathcal{C}^{ss}_{XX}}\mathcal{C}^{tt}_{XX}% \sqrt{\mathcal{C}^{ss}_{XX}}}caligraphic_C start_POSTSUPERSCRIPT italic_s italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT = square-root start_ARG square-root start_ARG caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT end_ARG caligraphic_C start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT square-root start_ARG caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT end_ARG end_ARG and 𝒞XXsssubscriptsuperscript𝒞𝑠𝑠𝑋𝑋\mathcal{C}^{ss}_{XX}caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT, 𝒞XXttsubscriptsuperscript𝒞𝑡𝑡𝑋𝑋\mathcal{C}^{tt}_{XX}caligraphic_C start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT are the covariance operators of PXssuperscriptsubscript𝑃𝑋𝑠P_{X}^{s}italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT and PXtsuperscriptsubscript𝑃𝑋𝑡P_{X}^{t}italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT on 𝒦subscript𝒦\mathcal{H}_{\mathcal{K}}caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT, respectively.

Remark 1

Note that KGW consists of two terms, i.e., the Maximum Mean Discrepancy (MMD) on first-order statistics and kernel Bures metric on second-order statistics. It provides a basic innovation for exploring conditional discrepancy metric by considering the conditional statistics in RKHS, i.e., first-order statistics in literature [48] and second-order statistics in literature [13].

3.2 Generalization Error Analysis for DAR

Most of DAR methods focus on aligning the marginal distribution PXsubscript𝑃𝑋P_{X}italic_P start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT. However, Zhao et al. [60] have shown that in the presence of label shift, only aligning the marginal distribution does not guarantee satisfying generalization performance on target domain. In fact, [60, 6] have proven the following theorem (3.1) for classification problems:

Theorem 3.1 ([60])

For representation variable Z=g(X)𝑍𝑔𝑋Z=g(X)italic_Z = italic_g ( italic_X ), suppose inequality dJS(PYs,PYt)dJS(PZs,PZt)subscript𝑑𝐽𝑆superscriptsubscript𝑃𝑌𝑠superscriptsubscript𝑃𝑌𝑡subscript𝑑𝐽𝑆superscriptsubscript𝑃𝑍𝑠superscriptsubscript𝑃𝑍𝑡d_{JS}(P_{Y}^{s},P_{Y}^{t})\geq d_{JS}(P_{Z}^{s},P_{Z}^{t})italic_d start_POSTSUBSCRIPT italic_J italic_S end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) ≥ italic_d start_POSTSUBSCRIPT italic_J italic_S end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) holds, then for predictor h:ZY:𝑍𝑌h:Z\to Yitalic_h : italic_Z → italic_Y, we have

εs(hg)+εt(hg)12(dJS(PYs,PYt)dJS(PZs,PZt))2,subscript𝜀𝑠𝑔subscript𝜀𝑡𝑔12superscriptsubscript𝑑𝐽𝑆superscriptsubscript𝑃𝑌𝑠superscriptsubscript𝑃𝑌𝑡subscript𝑑𝐽𝑆superscriptsubscript𝑃𝑍𝑠superscriptsubscript𝑃𝑍𝑡2\displaystyle\varepsilon_{s}(h{\circ}g){+}\varepsilon_{t}(h{\circ}g){\geq}% \frac{1}{2}\left(\sqrt{d_{JS}(P_{Y}^{s},P_{Y}^{t})}-\sqrt{d_{JS}(P_{Z}^{s},P_{% Z}^{t})}\right)^{2},italic_ε start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_h ∘ italic_g ) + italic_ε start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_h ∘ italic_g ) ≥ divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( square-root start_ARG italic_d start_POSTSUBSCRIPT italic_J italic_S end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) end_ARG - square-root start_ARG italic_d start_POSTSUBSCRIPT italic_J italic_S end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (3)

where dJSsubscript𝑑𝐽𝑆d_{JS}italic_d start_POSTSUBSCRIPT italic_J italic_S end_POSTSUBSCRIPT denotes the Jensen-Shannon divergence.

Assuming the label distributions differ between source and target domain, i.e., dJS(PYs,PYt)>0subscript𝑑𝐽𝑆superscriptsubscript𝑃𝑌𝑠superscriptsubscript𝑃𝑌𝑡0d_{JS}(P_{Y}^{s},P_{Y}^{t})>0italic_d start_POSTSUBSCRIPT italic_J italic_S end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) > 0, the lower bound shows that good alignment of representation distributions leads to bad joint error. In other words, when label shift exists, which is rather common in realistic scenarios, only aligning the marginal distribution is not enough for low target error.

To deal with the insufficiency problem shown in the lower bound analysis, it is intuitive to introduce label information into generalization analysis, e.g., the balanced error rate for classification scenarios [6]. Analogously, for the continuous regression scenarios, we first define the continuous balanced error rate for quantifying maximum conditional risk in DAR.

Definition 1

The continuous balanced error rate of predictor hhitalic_h on distribution P𝑃Pitalic_P is defined as:

cBERP(h,Y):=sup(x,y)𝒳𝒴P(h(x)YY=y).assignsubscriptcBERP𝑌subscriptsupremum𝑥𝑦direct-sum𝒳𝒴𝑃𝑥conditional𝑌𝑌𝑦\displaystyle\mathrm{cBER}_{\mathrm{P}}(h,Y):=\sup_{(x,y)\in\mathcal{X}\oplus% \mathcal{Y}}P(h(x)\neq Y\mid Y=y).roman_cBER start_POSTSUBSCRIPT roman_P end_POSTSUBSCRIPT ( italic_h , italic_Y ) := roman_sup start_POSTSUBSCRIPT ( italic_x , italic_y ) ∈ caligraphic_X ⊕ caligraphic_Y end_POSTSUBSCRIPT italic_P ( italic_h ( italic_x ) ≠ italic_Y ∣ italic_Y = italic_y ) . (4)

Intuitively, the continuous balanced error rate can be regarded as the maximum potential risk for predicting sample x𝑥xitalic_x with label y𝑦yitalic_y. Based on cBERcBER\mathrm{cBER}roman_cBER and conditional shift correction with representation transformation Z=g(X)𝑍𝑔𝑋Z=g(X)italic_Z = italic_g ( italic_X ), a sufficient condition for DAR can be induced as follows.

Theorem 3.2

For representation variable Zs=g(Xs)superscript𝑍𝑠𝑔superscript𝑋𝑠Z^{s}=g(X^{s})italic_Z start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT = italic_g ( italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) and Zt=g(Xt)superscript𝑍𝑡𝑔superscript𝑋𝑡Z^{t}=g(X^{t})italic_Z start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT = italic_g ( italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ), if g()𝑔g(\cdot)italic_g ( ⋅ ) satisfies conditional invariant property, i.e.,

PZ|Y=ys=PZ|Y=yt,y𝒴,formulae-sequencesuperscriptsubscript𝑃conditional𝑍𝑌𝑦𝑠superscriptsubscript𝑃conditional𝑍𝑌𝑦𝑡for-all𝑦𝒴\displaystyle P_{Z|Y=y}^{s}=P_{Z|Y=y}^{t},~{}~{}\forall y\in\mathcal{Y},italic_P start_POSTSUBSCRIPT italic_Z | italic_Y = italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT = italic_P start_POSTSUBSCRIPT italic_Z | italic_Y = italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT , ∀ italic_y ∈ caligraphic_Y , (5)

then for any predictor h:ZY:𝑍𝑌h:Z\to Yitalic_h : italic_Z → italic_Y, we have

εs(h)+εt(h)2cBERPs(h,Y).subscript𝜀𝑠subscript𝜀𝑡2cBEsubscriptRsuperscriptPs𝑌\displaystyle\varepsilon_{s}(h)+\varepsilon_{t}(h)\leq 2\mathrm{cBER}_{\mathrm% {P^{s}}}(h,Y).italic_ε start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_h ) + italic_ε start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_h ) ≤ 2 roman_c roman_B roman_E roman_R start_POSTSUBSCRIPT roman_P start_POSTSUPERSCRIPT roman_s end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_h , italic_Y ) . (6)

Therefore, the key for DA, no matter for classification or regression problems, is to achieve conditional distribution alignment. In classification settings, conditional alignment can be achieved by conducting class-wise marginal distribution alignment with the help of target pseudo label. However, in regression scenarios, the label variable is continuous, which means previous methods correcting conditional shift in classification settings are no longer applicable. Obviously, it is impossible to align the PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT for every single y𝑦yitalic_y when y𝑦yitalic_y is a continuous variable and the possible number of y is infinite. Thus, it is important and meaningful to explore a new framework that is capable of conditional distribution alignment for continuous label variable.

3.3 COD metric for DAR

Theory for COD. Due to the continuity of label variable, it is impractical to correct conditional shift by aligning marginal distribution of PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT for every single y𝑦yitalic_y. For this inevitable obstacle, an intuitive idea is can we align all PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT as a whole? i.e., PX|Ysubscript𝑃conditional𝑋𝑌P_{X|Y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT, rather than align PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT one by one? It is more than attractive that if we can characterize PX|Ysubscript𝑃conditional𝑋𝑌P_{X|Y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT with a well-defined discrepancy, then all PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT can be aligned by aligning the whole PX|Ysubscript𝑃conditional𝑋𝑌P_{X|Y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT.

In kernel embedding theory, a distribution can be characterized by a point in RKHS. Then PX|Ysubscript𝑃conditional𝑋𝑌P_{X|Y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT can be seen as a point set in RKHS in which every element corresponds to a distribution PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT for a fixed y𝑦yitalic_y. From another perspective, we can consider PX|Ysubscript𝑃conditional𝑋𝑌P_{X|Y}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT as the distribution in RKHS and measure the cross-domain conditional discrepancy by defining the COD distance on their statistical moments in RKHS. Following the definitions of conditional statistics in pioneer works, i.e., the conditional mean operator as 𝒰X|Y:=𝒞XY𝒞YY1assignsubscript𝒰conditional𝑋𝑌subscript𝒞𝑋𝑌superscriptsubscript𝒞𝑌𝑌1\mathcal{U}_{X|Y}:=\mathcal{C}_{XY}\mathcal{C}_{YY}^{-1}caligraphic_U start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT := caligraphic_C start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT caligraphic_C start_POSTSUBSCRIPT italic_Y italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT [48] and conditional covariance operator as 𝒞XX|Y=𝒞XX𝒞XY𝒞YY1𝒞YXsubscript𝒞conditional𝑋𝑋𝑌subscript𝒞𝑋𝑋subscript𝒞𝑋𝑌superscriptsubscript𝒞𝑌𝑌1subscript𝒞𝑌𝑋\mathcal{C}_{XX|Y}=\mathcal{C}_{XX}-\mathcal{C}_{XY}\mathcal{C}_{YY}^{-1}% \mathcal{C}_{YX}caligraphic_C start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT = caligraphic_C start_POSTSUBSCRIPT italic_X italic_X end_POSTSUBSCRIPT - caligraphic_C start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT caligraphic_C start_POSTSUBSCRIPT italic_Y italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT caligraphic_C start_POSTSUBSCRIPT italic_Y italic_X end_POSTSUBSCRIPT [13], we present the definition of COD as follows.

Definition 2 (COD)

Given conditional distributions PX|Yssuperscriptsubscript𝑃conditional𝑋𝑌𝑠P_{X|Y}^{s}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT and PX|Ytsuperscriptsubscript𝑃conditional𝑋𝑌𝑡P_{X|Y}^{t}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT, denote 𝒞XX|Ysssubscriptsuperscript𝒞𝑠𝑠conditional𝑋𝑋𝑌\mathcal{C}^{ss}_{XX|Y}caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT and 𝒞XX|Yttsubscriptsuperscript𝒞𝑡𝑡conditional𝑋𝑋𝑌\mathcal{C}^{tt}_{XX|Y}caligraphic_C start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT as their corresponding covariance operators. The COD is defined as

dCOD2(PX|Ys,PX|Yt)=superscriptsubscript𝑑COD2superscriptsubscript𝑃conditional𝑋𝑌𝑠superscriptsubscript𝑃conditional𝑋𝑌𝑡absent\displaystyle d_{\mathrm{COD}}^{2}({P}_{X|Y}^{s},{P}_{X|Y}^{t})=italic_d start_POSTSUBSCRIPT roman_COD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) = 𝒰X|Ys𝒰X|Yt𝒦2+tr(𝒞XX|Yss+𝒞XX|Ytt2𝒞XX|Yst),subscriptsuperscriptnormsubscriptsuperscript𝒰𝑠conditional𝑋𝑌subscriptsuperscript𝒰𝑡conditional𝑋𝑌2subscript𝒦trsubscriptsuperscript𝒞𝑠𝑠conditional𝑋𝑋𝑌subscriptsuperscript𝒞𝑡𝑡conditional𝑋𝑋𝑌2subscriptsuperscript𝒞𝑠𝑡conditional𝑋𝑋𝑌\displaystyle\|\mathcal{U}^{s}_{X|Y}-\mathcal{U}^{t}_{X|Y}\|^{2}_{\mathcal{H}_% {\mathcal{K}}}+\mathrm{tr}(\mathcal{C}^{ss}_{XX|Y}+\mathcal{C}^{tt}_{XX|Y}-2% \mathcal{C}^{st}_{XX|Y}),∥ caligraphic_U start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT - caligraphic_U start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT ∥ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT caligraphic_K end_POSTSUBSCRIPT end_POSTSUBSCRIPT + roman_tr ( caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT + caligraphic_C start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT - 2 caligraphic_C start_POSTSUPERSCRIPT italic_s italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT ) , (7)

where 𝒞XX|Yst=𝒞XX|Yss𝒞XX|Ytt𝒞XX|Ysssubscriptsuperscript𝒞𝑠𝑡conditional𝑋𝑋𝑌subscriptsuperscript𝒞𝑠𝑠conditional𝑋𝑋𝑌subscriptsuperscript𝒞𝑡𝑡conditional𝑋𝑋𝑌subscriptsuperscript𝒞𝑠𝑠conditional𝑋𝑋𝑌\mathcal{C}^{st}_{XX|Y}=\sqrt{\sqrt{\mathcal{C}^{ss}_{XX|Y}}\mathcal{C}^{tt}_{% XX|Y}\sqrt{\mathcal{C}^{ss}_{XX|Y}}}caligraphic_C start_POSTSUPERSCRIPT italic_s italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT = square-root start_ARG square-root start_ARG caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT end_ARG caligraphic_C start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT square-root start_ARG caligraphic_C start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT end_ARG end_ARG.

Refer to caption
Figure 2: Illustration of COD metric. (a) Classification DA methods usually correct conditional shift by class-wise computations on discrete clusters PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT, which is infeasible for regression due to the infinite slices PX|ysubscript𝑃conditional𝑋𝑦P_{X|y}italic_P start_POSTSUBSCRIPT italic_X | italic_y end_POSTSUBSCRIPT. (b) In COD, the continuous conditional distributions are embedded into the RKHS and characterized by finite statistical moments in kernel spaces, which do not rely on specific conditions y𝑦yitalic_y. Under the guarantees of distribution embedding property, the conditional alignment is equivalent to the matching on first-order and second-order statistics, i.e., mean embedding operator 𝒰X|Ysubscript𝒰conditional𝑋𝑌\mathcal{U}_{X|Y}caligraphic_U start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT and covariance operator 𝒞XX|Ysubscript𝒞conditional𝑋𝑋𝑌\mathcal{C}_{XX|Y}caligraphic_C start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT. When COD is minimized, conditional alignment is achieved, and the level sets of cross-domain labeling rules are then aligned.

Theoretically, it can be proved that COD indeed admits the metric property on conditional distributions, i.e., two conditional distributions are the same if and only if the COD metric between them is zero. An intuitive illustration of COD is shown in Fig. 2. Formally, COD admits two key properties.

Proposition 1

(a) COD defines a metric between Gaussian measures.
(b) Let (𝒳,𝒳)𝒳subscript𝒳(\mathcal{X},\mathcal{B}_{\mathcal{X}})( caligraphic_X , caligraphic_B start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ) be the locally compact and Hausdorff measurable space and k be c0subscript𝑐0c_{0}italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT-universal kernel. Assuming that (ϕ(X),ψ(Y))italic-ϕ𝑋𝜓𝑌(\phi(X),\psi(Y))( italic_ϕ ( italic_X ) , italic_ψ ( italic_Y ) ) is a Gaussian random variable in 𝒳𝒴direct-sumsubscript𝒳subscript𝒴\mathcal{H}_{\mathcal{X}}\oplus\mathcal{H}_{\mathcal{Y}}caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ⊕ caligraphic_H start_POSTSUBSCRIPT caligraphic_Y end_POSTSUBSCRIPT. For any PX|Ys,PX|Yt𝒳|𝒴superscriptsubscript𝑃conditional𝑋𝑌𝑠superscriptsubscript𝑃conditional𝑋𝑌𝑡subscriptconditional𝒳𝒴P_{X|Y}^{s},P_{X|Y}^{t}\in\mathbb{P}_{\mathcal{X}|\mathcal{Y}}italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ∈ blackboard_P start_POSTSUBSCRIPT caligraphic_X | caligraphic_Y end_POSTSUBSCRIPT, we have

dCOD2(PX|Ys,PX|Yt)=0PX|Ys=PX|Yt.superscriptsubscript𝑑COD2superscriptsubscript𝑃conditional𝑋𝑌𝑠superscriptsubscript𝑃conditional𝑋𝑌𝑡0superscriptsubscript𝑃conditional𝑋𝑌𝑠superscriptsubscript𝑃conditional𝑋𝑌𝑡\displaystyle d_{\mathrm{COD}}^{2}(P_{X|Y}^{s},P_{X|Y}^{t})=0\implies P_{X|Y}^% {s}=P_{X|Y}^{t}.italic_d start_POSTSUBSCRIPT roman_COD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) = 0 ⟹ italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT = italic_P start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT . (8)

Empirical Estimation of COD. Based on the above properties, we can perform conditional distribution alignment by minimizing the COD metric between source and target domain. To do this, we shall derive the empirical estimation of COD, which can be obtained with the help of the kernel trick.

Let (𝐱is,𝐲is)i=1nssuperscriptsubscriptsuperscriptsubscript𝐱𝑖𝑠superscriptsubscript𝐲𝑖𝑠𝑖1subscript𝑛𝑠{(\mathbf{x}_{i}^{s},\mathbf{y}_{i}^{s})}_{i=1}^{n_{s}}( bold_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , bold_y start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUPERSCRIPT and (𝐱it,𝐲it)i=1ntsuperscriptsubscriptsuperscriptsubscript𝐱𝑖𝑡superscriptsubscript𝐲𝑖𝑡𝑖1subscript𝑛𝑡{(\mathbf{x}_{i}^{t},\mathbf{y}_{i}^{t})}_{i=1}^{n_{t}}( bold_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT , bold_y start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUPERSCRIPT be two set of samples drawn i.i.d. from source and target domain. For simplicity, nssubscript𝑛𝑠n_{s}italic_n start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and ntsubscript𝑛𝑡n_{t}italic_n start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT are both set to n𝑛nitalic_n. In kernel method, samples are mapped to RKHS 𝒳𝒴direct-sumsubscript𝒳subscript𝒴\mathcal{H}_{\mathcal{X}}\oplus\mathcal{H}_{\mathcal{Y}}caligraphic_H start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ⊕ caligraphic_H start_POSTSUBSCRIPT caligraphic_Y end_POSTSUBSCRIPT by the implicit feature map (ϕ,ψ)italic-ϕ𝜓(\phi,\psi)( italic_ϕ , italic_ψ ). Denote the feature map matrices by (𝚽s,𝚿s)superscript𝚽𝑠superscript𝚿𝑠(\mathbf{\Phi}^{s},\mathbf{\Psi}^{s})( bold_Φ start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , bold_Ψ start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) and (𝚽t,𝚿t)superscript𝚽𝑡superscript𝚿𝑡(\mathbf{\Phi}^{t},\mathbf{\Psi}^{t})( bold_Φ start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT , bold_Ψ start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ), the explicit kernel matrices 𝐊Xss=𝚽s𝚽sTsubscriptsuperscript𝐊𝑠𝑠𝑋superscript𝚽𝑠superscriptsuperscript𝚽𝑠𝑇\mathbf{K}^{ss}_{X}=\mathbf{\Phi}^{s}{\mathbf{\Phi}^{s}}^{T}bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT = bold_Φ start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT bold_Φ start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT, 𝐊Yss=𝚿s𝚿sTsubscriptsuperscript𝐊𝑠𝑠𝑌superscript𝚿𝑠superscriptsuperscript𝚿𝑠𝑇\mathbf{K}^{ss}_{Y}=\mathbf{\Psi}^{s}{\mathbf{\Psi}^{s}}^{T}bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT = bold_Ψ start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT bold_Ψ start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT are computed as (𝐊Xss)ij=k𝒳(𝐱is,𝐱js)subscriptsubscriptsuperscript𝐊𝑠𝑠𝑋𝑖𝑗subscript𝑘𝒳superscriptsubscript𝐱𝑖𝑠superscriptsubscript𝐱𝑗𝑠(\mathbf{K}^{ss}_{X})_{ij}=k_{\mathcal{X}}(\mathbf{x}_{i}^{s},\mathbf{x}_{j}^{% s})( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = italic_k start_POSTSUBSCRIPT caligraphic_X end_POSTSUBSCRIPT ( bold_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , bold_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ), (𝐊Yss)ij=k𝒴(𝐲is,𝐲js)subscriptsubscriptsuperscript𝐊𝑠𝑠𝑌𝑖𝑗subscript𝑘𝒴superscriptsubscript𝐲𝑖𝑠superscriptsubscript𝐲𝑗𝑠(\mathbf{K}^{ss}_{Y})_{ij}=k_{\mathcal{Y}}(\mathbf{y}_{i}^{s},\mathbf{y}_{j}^{% s})( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = italic_k start_POSTSUBSCRIPT caligraphic_Y end_POSTSUBSCRIPT ( bold_y start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , bold_y start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ), respectively. And so as 𝐊Xttsubscriptsuperscript𝐊𝑡𝑡𝑋\mathbf{K}^{tt}_{X}bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT, 𝐊Xstsubscriptsuperscript𝐊𝑠𝑡𝑋\mathbf{K}^{st}_{X}bold_K start_POSTSUPERSCRIPT italic_s italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT and 𝐊Yttsubscriptsuperscript𝐊𝑡𝑡𝑌\mathbf{K}^{tt}_{Y}bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT. Let 𝐈nsubscript𝐈𝑛\mathbf{I}_{n}bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT be the n-dimensional identity matrix, 𝟏nsubscript1𝑛\mathbf{1}_{n}bold_1 start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT be the n-dimensional vector with all elements equal to 1. The n𝑛nitalic_n-dimensional centering matrix can be defined as 𝐇n=𝐈n1n𝟏n𝟏nTsubscript𝐇𝑛subscript𝐈𝑛1𝑛subscript1𝑛superscriptsubscript1𝑛𝑇\mathbf{H}_{n}=\mathbf{I}_{n}-\frac{1}{n}\mathbf{1}_{n}\mathbf{1}_{n}^{T}bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG italic_n end_ARG bold_1 start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_1 start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT. Then the centralized kernel matrix is defined as 𝐆=𝐇n𝐊𝐇n𝐆subscript𝐇𝑛subscript𝐊𝐇𝑛\mathbf{G}=\mathbf{H}_{n}\mathbf{K}\mathbf{H}_{n}bold_G = bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_KH start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT.

With above notations, the covariance operator 𝒞XYsubscript𝒞𝑋𝑌\mathcal{C}_{XY}caligraphic_C start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT can be empirically estimated as 𝒞^XY=1n(𝚽μ^X𝟏T)(𝚿μ^Y𝟏T)T=1n𝚽𝐇𝚿Tsubscript^𝒞𝑋𝑌1𝑛𝚽subscript^𝜇𝑋superscript1𝑇superscript𝚿subscript^𝜇𝑌superscript1𝑇𝑇1𝑛𝚽𝐇superscript𝚿𝑇\hat{\mathcal{C}}_{XY}=\frac{1}{n}(\mathbf{\Phi}-\hat{\mu}_{X}\mathbf{1}^{T})(% \mathbf{\Psi}-\hat{\mu}_{Y}\mathbf{1}^{T})^{T}=\frac{1}{n}\mathbf{\Phi}\mathbf% {H}\mathbf{\Psi}^{T}over^ start_ARG caligraphic_C end_ARG start_POSTSUBSCRIPT italic_X italic_Y end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_n end_ARG ( bold_Φ - over^ start_ARG italic_μ end_ARG start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT bold_1 start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) ( bold_Ψ - over^ start_ARG italic_μ end_ARG start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT bold_1 start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_n end_ARG bold_Φ bold_H bold_Ψ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT, then the conditional moment statistics can be estimated as

𝒰^X|Y=(1n𝚽𝐇𝚿T)(1n𝚿𝐇𝚿T)1,𝒞^XX|Y=1n𝚽𝐇n𝐀(𝚽𝐇n𝐀)T,formulae-sequencesubscript^𝒰conditional𝑋𝑌1𝑛𝚽𝐇superscript𝚿𝑇superscript1𝑛𝚿𝐇superscript𝚿𝑇1subscript^𝒞conditional𝑋𝑋𝑌1𝑛𝚽subscript𝐇𝑛𝐀superscript𝚽subscript𝐇𝑛𝐀𝑇\displaystyle\hat{\mathcal{U}}_{X|Y}=(\frac{1}{n}\mathbf{\Phi}\mathbf{H}% \mathbf{\Psi}^{T})(\frac{1}{n}\mathbf{\Psi}\mathbf{H}\mathbf{\Psi}^{T})^{-1},~% {}~{}~{}\hat{\mathcal{C}}_{XX|Y}=\frac{1}{n}{\mathbf{\Phi}}\mathbf{H}_{n}{% \mathbf{A}}\left({\mathbf{\Phi}}\mathbf{H}_{n}{\mathbf{A}}\right)^{T},over^ start_ARG caligraphic_U end_ARG start_POSTSUBSCRIPT italic_X | italic_Y end_POSTSUBSCRIPT = ( divide start_ARG 1 end_ARG start_ARG italic_n end_ARG bold_Φ bold_H bold_Ψ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) ( divide start_ARG 1 end_ARG start_ARG italic_n end_ARG bold_Ψ bold_H bold_Ψ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , over^ start_ARG caligraphic_C end_ARG start_POSTSUBSCRIPT italic_X italic_X | italic_Y end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_n end_ARG bold_Φ bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_A ( bold_Φ bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_A ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT , (9)

where 𝐀𝐀\mathbf{A}bold_A satisfies that 𝐁=𝐀𝐀T𝐁superscript𝐀𝐀𝑇\mathbf{B}={\mathbf{A}}{{\mathbf{A}}}^{T}bold_B = bold_AA start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT and B𝐈n1nε[𝐆Y𝐆Y(𝐆Y+εn𝐈n)1𝐆Y]𝐵subscript𝐈𝑛1𝑛𝜀delimited-[]subscript𝐆𝑌subscript𝐆𝑌superscriptsubscript𝐆𝑌𝜀𝑛subscript𝐈𝑛1subscript𝐆𝑌B\triangleq\mathbf{I}_{n}-\frac{1}{n\varepsilon}\left[{\mathbf{G}}_{Y}-{% \mathbf{G}}_{Y}\left({\mathbf{G}}_{Y}+\varepsilon n{\mathbf{I}}_{n}\right)^{-1% }{\mathbf{G}}_{Y}\right]italic_B ≜ bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG italic_n italic_ε end_ARG [ bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT - bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT + italic_ε italic_n bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ]. Based on these estimations, the closed-form empirical estimation of the COD metric can be derived as follows.

Proposition 2

The empirical Conditional Operator Discrepancy distance is

d^COD2=superscriptsubscript^𝑑COD2absent\displaystyle\widehat{d}_{\mathrm{COD}}^{2}=over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT roman_COD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = tr(𝐊Yss(𝐊Yss)1𝐊Xss(𝐊Yss)1)+tr(𝐊Ytt(𝐊Ytt)1𝐊Xtt(𝐊Ytt)1)trsubscriptsuperscript𝐊𝑠𝑠𝑌superscriptsubscriptsuperscript𝐊𝑠𝑠𝑌1subscriptsuperscript𝐊𝑠𝑠𝑋superscriptsubscriptsuperscript𝐊𝑠𝑠𝑌1trsubscriptsuperscript𝐊𝑡𝑡𝑌superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1subscriptsuperscript𝐊𝑡𝑡𝑋superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1\displaystyle\ \mathrm{tr}(\mathbf{K}^{ss}_{Y}(\mathbf{K}^{ss}_{Y})^{-1}% \mathbf{K}^{ss}_{X}(\mathbf{K}^{ss}_{Y})^{-1})+\mathrm{tr}(\mathbf{K}^{tt}_{Y}% (\mathbf{K}^{tt}_{Y})^{-1}\mathbf{K}^{tt}_{X}(\mathbf{K}^{tt}_{Y})^{-1})roman_tr ( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) + roman_tr ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) (10)
2tr(𝐊Yts(𝐊Yss)1𝐊XsXt(𝐊Ytt)1)2trsubscriptsuperscript𝐊𝑡𝑠𝑌superscriptsubscriptsuperscript𝐊𝑠𝑠𝑌1subscript𝐊superscript𝑋𝑠superscript𝑋𝑡superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1\displaystyle-2\mathrm{tr}(\mathbf{K}^{ts}_{Y}(\mathbf{K}^{ss}_{Y})^{-1}% \mathbf{K}_{X^{s}X^{t}}(\mathbf{K}^{tt}_{Y})^{-1})- 2 roman_t roman_r ( bold_K start_POSTSUPERSCRIPT italic_t italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUBSCRIPT italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT )
+εtr[𝐆Xs(εn𝐈n+𝐆Ys)1]+εtr[𝐆Xt(εn𝐈n+𝐆Yt)1]𝜀trdelimited-[]superscriptsubscript𝐆𝑋𝑠superscript𝜀𝑛subscript𝐈𝑛superscriptsubscript𝐆𝑌𝑠1𝜀trdelimited-[]superscriptsubscript𝐆𝑋𝑡superscript𝜀𝑛subscript𝐈𝑛superscriptsubscript𝐆𝑌𝑡1\displaystyle+\varepsilon\mathrm{tr}\left[{\mathbf{G}}_{X}^{s}\left(% \varepsilon n{\mathbf{I}}_{n}+{\mathbf{G}}_{Y}^{s}\right)^{-1}\right]+% \varepsilon\mathrm{tr}\left[\mathbf{G}_{X}^{t}\left(\varepsilon n{\mathbf{I}}_% {n}+\mathbf{G}_{Y}^{t}\right)^{-1}\right]+ italic_ε roman_tr [ bold_G start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ( italic_ε italic_n bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] + italic_ε roman_tr [ bold_G start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ( italic_ε italic_n bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ]
2n(𝐇n𝐀t)T𝐊Xts(𝐇n𝐀s),2𝑛subscriptnormsuperscriptsubscript𝐇𝑛superscript𝐀𝑡𝑇subscriptsuperscript𝐊𝑡𝑠𝑋subscript𝐇𝑛superscript𝐀𝑠\displaystyle-\frac{2}{n}\left\|\left(\mathbf{H}_{n}{\mathbf{A}}^{t}\right)^{T% }{\mathbf{K}}^{ts}_{X}\left(\mathbf{H}_{n}{\mathbf{A}}^{s}\right)\right\|_{*},- divide start_ARG 2 end_ARG start_ARG italic_n end_ARG ∥ ( bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_A start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_A start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) ∥ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ,

where \left\|\cdot\right\|_{*}∥ ⋅ ∥ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is the nuclear norm. Empirically, 𝐊Y1superscriptsubscript𝐊𝑌1\mathbf{K}_{Y}^{-1}bold_K start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT is approximated by the regularized formulation (𝐊Y+λ𝐈)1superscriptsubscript𝐊𝑌𝜆𝐈1(\mathbf{K}_{Y}+\lambda\mathbf{I})^{-1}( bold_K start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT + italic_λ bold_I ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT.

Discriminability of empirical COD. As discussed above, the first term of COD measures the distance of the first-order statistic between conditional distributions in RKHS. However, recent study [5, 54] shows that the mean matching may reduce the discriminability when bringing the two distributions close, which undoubtedly degrades the performance of the algorithm. Therefore, it is worth making an in-depth analysis of the discriminability of COD as well.

Without loss of generality, we first consider Kronechker Delta kernel on conditioning variable Y𝑌Yitalic_Y, i.e., k(𝐲1,𝐲2)=δ(𝐲1,𝐲2)𝑘subscript𝐲1subscript𝐲2𝛿subscript𝐲1subscript𝐲2k(\mathbf{y}_{1},\mathbf{y}_{2})=\delta(\mathbf{y}_{1},\mathbf{y}_{2})italic_k ( bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = italic_δ ( bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ). More explicitly, k(𝐲1,𝐲2)𝑘subscript𝐲1subscript𝐲2k(\mathbf{y}_{1},\mathbf{y}_{2})italic_k ( bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) equals to 1111 if 𝐲1=𝐲2subscript𝐲1subscript𝐲2\mathbf{y}_{1}=\mathbf{y}_{2}bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT or 00 otherwise. For empirical scenarios with finite samples, the numbers of samples and observed labels are finite, which implies the number of distinct label c𝑐citalic_c will be obviously no more than sample size n𝑛nitalic_n. For convenience, we call two samples are in the same class when their label values are the same. Consequently, the first term in COD can be reformulated as follows.

Proposition 3

Under the conditions mentioned above, the empirical estimation of dCMMDsubscript𝑑𝐶𝑀𝑀𝐷d_{CMMD}italic_d start_POSTSUBSCRIPT italic_C italic_M italic_M italic_D end_POSTSUBSCRIPT can be written as

d^CMMD2=superscriptsubscript^𝑑CMMD2absent\displaystyle\widehat{d}_{\mathrm{CMMD}}^{2}=over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT roman_CMMD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = p=1c1(λ+np)2(i,jk(𝐱(p),is,𝐱(p),js))superscriptsubscript𝑝1𝑐1superscript𝜆subscript𝑛𝑝2subscript𝑖𝑗𝑘superscriptsubscript𝐱𝑝𝑖𝑠superscriptsubscript𝐱𝑝𝑗𝑠\displaystyle\sum_{p=1}^{c}\frac{1}{(\lambda+n_{p})^{2}}\left(\sum_{i,j}k(% \mathbf{x}_{(p),i}^{s},\mathbf{x}_{(p),j}^{s})\right)∑ start_POSTSUBSCRIPT italic_p = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG ( italic_λ + italic_n start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( ∑ start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT italic_k ( bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) ) (11)
+p=1c1(λ+np)2(i,jk(𝐱(p),it,𝐱(p),jt))superscriptsubscript𝑝1𝑐1superscript𝜆subscript𝑛𝑝2subscript𝑖𝑗𝑘superscriptsubscript𝐱𝑝𝑖𝑡superscriptsubscript𝐱𝑝𝑗𝑡\displaystyle+\sum_{p=1}^{c}\frac{1}{(\lambda+n_{p})^{2}}\left(\sum_{i,j}k(% \mathbf{x}_{(p),i}^{t},\mathbf{x}_{(p),j}^{t})\right)+ ∑ start_POSTSUBSCRIPT italic_p = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG ( italic_λ + italic_n start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( ∑ start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT italic_k ( bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT , bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) )
2p=1c1(λ+np)2(i,jk(𝐱(p),is,𝐱(p),jt)),2superscriptsubscript𝑝1𝑐1superscript𝜆subscript𝑛𝑝2subscript𝑖𝑗𝑘superscriptsubscript𝐱𝑝𝑖𝑠superscriptsubscript𝐱𝑝𝑗𝑡\displaystyle-2\sum_{p=1}^{c}\frac{1}{(\lambda+n_{p})^{2}}\left(\sum_{i,j}k(% \mathbf{x}_{(p),i}^{s},\mathbf{x}_{(p),j}^{t})\right),- 2 ∑ start_POSTSUBSCRIPT italic_p = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG ( italic_λ + italic_n start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( ∑ start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT italic_k ( bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT , bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) ) ,

where nisubscript𝑛𝑖n_{i}italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT denotes the sample-size of i𝑖iitalic_i-th label such that n=i=1cni𝑛superscriptsubscript𝑖1𝑐subscript𝑛𝑖n=\sum_{i=1}^{c}n_{i}italic_n = ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, 𝐱(p),isubscript𝐱𝑝𝑖\mathbf{x}_{(p),i}bold_x start_POSTSUBSCRIPT ( italic_p ) , italic_i end_POSTSUBSCRIPT denotes the i𝑖iitalic_i-th sample of p𝑝pitalic_p-th label.

Eq. (11) contains three terms. The third term characterizes the cross-domain intra-class similarity. When Eq. (11) is minimized, the third term is maximized so that the cross-domain sample with the same label will be aligned, leading to conditional alignment as expected. However, the first two terms of the equation represent the intra-class similarity of both domains. Minimizing them will decrease the intra-class similarity, which means intra-class samples will be pushed away from each other and the discriminability of representations is degraded.

Remark 2

Note the Kronechker Delta kernel is defined as

k(𝐲1,𝐲2)={1 if 𝐲1=𝐲20 otherwise.𝑘subscript𝐲1subscript𝐲2cases1 if subscript𝐲1subscript𝐲2otherwise0 otherwiseotherwise\displaystyle k(\mathbf{y}_{1},\mathbf{y}_{2})=\begin{cases}1\text{ if }% \mathbf{y}_{1}=\mathbf{y}_{2}\\ 0\text{ otherwise}\end{cases}.italic_k ( bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = { start_ROW start_CELL 1 if bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL 0 otherwise end_CELL start_CELL end_CELL end_ROW . (12)

When it comes to the Gaussian kernel for empirical modeling, the case of 𝐲1=𝐲2subscript𝐲1subscript𝐲2\mathbf{y}_{1}=\mathbf{y}_{2}bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is the same; for other cases, the kernel values can be arbitrarily close to zero, i.e., for any ϵ>0italic-ϵ0\epsilon>0italic_ϵ > 0, there always width σ0subscript𝜎0\sigma_{0}italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, such that

k(𝐲1,𝐲2)=exp(𝐲1𝐲222/σ02)<ϵ.𝑘subscript𝐲1subscript𝐲2superscriptsubscriptnormsubscript𝐲1subscript𝐲222superscriptsubscript𝜎02italic-ϵ\displaystyle k(\mathbf{y}_{1},\mathbf{y}_{2})=\exp\left(-\|{\mathbf{y}_{1}}-{% \mathbf{y}_{2}}\|_{2}^{2}/\sigma_{0}^{2}\right)<\epsilon.italic_k ( bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = roman_exp ( - ∥ bold_y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ∥ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) < italic_ϵ . (13)

Thus, it is generally reasonable to use Gaussian kernel, which admits a similar property to Kronechker Delta kernel, for empirical modeling.

Remark 3

Note that although we treat the number of distinct labels as finite numbers in mini-batch training, it does not mean class-wise matching can be conducted. In fact, even though the pseudo target label has finite range in mini-batch, it still admits the continuous property, which means we can not cluster samples as usually done in classification methods. When two samples have predicted labels with little numerical difference, they can be regarded as samples of the same class in classification settings. Nonetheless, this kind of clustering is infeasible in regression settings since the labels are continuous. Therefore, class-wise alignment is not practical even in mini-batch training.

With these observation, it is reasonable to reformulate Eq. (11) to promote the discriminability of representations:

d^CMMDmod=subscript^𝑑subscriptCMMDmodabsent\displaystyle\widehat{d}_{\mathrm{CMMD_{mod}}}=over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT roman_CMMD start_POSTSUBSCRIPT roman_mod end_POSTSUBSCRIPT end_POSTSUBSCRIPT = tr(𝐊Ytt(𝐊Ytt)1𝐊Xtt(𝐊Ytt)1)tr(𝐊Ytt(𝐊Ytt)1𝐊Xtt(𝐊Ytt)1)trsubscriptsuperscript𝐊𝑡𝑡𝑌superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1subscriptsuperscript𝐊𝑡𝑡𝑋superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1trsubscriptsuperscript𝐊𝑡𝑡𝑌superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1subscriptsuperscript𝐊𝑡𝑡𝑋superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1\displaystyle-\mathrm{tr}(\mathbf{K}^{tt}_{Y}(\mathbf{K}^{tt}_{Y})^{-1}\mathbf% {K}^{tt}_{X}(\mathbf{K}^{tt}_{Y})^{-1})-\mathrm{tr}(\mathbf{K}^{tt}_{Y}(% \mathbf{K}^{tt}_{Y})^{-1}\mathbf{K}^{tt}_{X}(\mathbf{K}^{tt}_{Y})^{-1})- roman_tr ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) - roman_tr ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) (14)
2tr(𝐊YtYs(𝐊Yss)1𝐊XsXt(𝐊Ytt)1)2trsubscript𝐊superscript𝑌𝑡superscript𝑌𝑠superscriptsubscriptsuperscript𝐊𝑠𝑠𝑌1subscript𝐊superscript𝑋𝑠superscript𝑋𝑡superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1\displaystyle-2\mathrm{tr}(\mathbf{K}_{Y^{t}Y^{s}}(\mathbf{K}^{ss}_{Y})^{-1}% \mathbf{K}_{X^{s}X^{t}}(\mathbf{K}^{tt}_{Y})^{-1})- 2 roman_t roman_r ( bold_K start_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT italic_Y start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUBSCRIPT italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT )

Now minimizing Eq. (14) will maximize intra-class similarity regardless of domain, yielding conditional matching with great discriminability. Then the empirical COD is modified as

d^CODmod=subscript^𝑑subscriptCODmodabsent\displaystyle\widehat{d}_{\mathrm{COD_{mod}}}=over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT roman_COD start_POSTSUBSCRIPT roman_mod end_POSTSUBSCRIPT end_POSTSUBSCRIPT = tr(𝐊Ytt(𝐊Ytt)1𝐊Xtt(𝐊Ytt)1)tr(𝐊Ytt(𝐊Ytt)1𝐊Xtt(𝐊Ytt)1)trsubscriptsuperscript𝐊𝑡𝑡𝑌superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1subscriptsuperscript𝐊𝑡𝑡𝑋superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1trsubscriptsuperscript𝐊𝑡𝑡𝑌superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1subscriptsuperscript𝐊𝑡𝑡𝑋superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1\displaystyle-\mathrm{tr}(\mathbf{K}^{tt}_{Y}(\mathbf{K}^{tt}_{Y})^{-1}\mathbf% {K}^{tt}_{X}(\mathbf{K}^{tt}_{Y})^{-1})-\mathrm{tr}(\mathbf{K}^{tt}_{Y}(% \mathbf{K}^{tt}_{Y})^{-1}\mathbf{K}^{tt}_{X}(\mathbf{K}^{tt}_{Y})^{-1})- roman_tr ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) - roman_tr ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) (15)
2tr(𝐊YtYs(𝐊Yss)1𝐊XsXt(𝐊Ytt)1)2trsubscript𝐊superscript𝑌𝑡superscript𝑌𝑠superscriptsubscriptsuperscript𝐊𝑠𝑠𝑌1subscript𝐊superscript𝑋𝑠superscript𝑋𝑡superscriptsubscriptsuperscript𝐊𝑡𝑡𝑌1\displaystyle-2\mathrm{tr}(\mathbf{K}_{Y^{t}Y^{s}}(\mathbf{K}^{ss}_{Y})^{-1}% \mathbf{K}_{X^{s}X^{t}}(\mathbf{K}^{tt}_{Y})^{-1})- 2 roman_t roman_r ( bold_K start_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT italic_Y start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_K start_POSTSUBSCRIPT italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( bold_K start_POSTSUPERSCRIPT italic_t italic_t end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT )
+εtr[𝐆Xs(εn𝐈n+𝐆Ys)1]+εtr[𝐆Xt(εn𝐈n+𝐆Yt)1]𝜀trdelimited-[]superscriptsubscript𝐆𝑋𝑠superscript𝜀𝑛subscript𝐈𝑛superscriptsubscript𝐆𝑌𝑠1𝜀trdelimited-[]superscriptsubscript𝐆𝑋𝑡superscript𝜀𝑛subscript𝐈𝑛superscriptsubscript𝐆𝑌𝑡1\displaystyle+\varepsilon\mathrm{tr}\left[{\mathbf{G}}_{X}^{s}\left(% \varepsilon n{\mathbf{I}}_{n}+{\mathbf{G}}_{Y}^{s}\right)^{-1}\right]+% \varepsilon\mathrm{tr}\left[\mathbf{G}_{X}^{t}\left(\varepsilon n{\mathbf{I}}_% {n}+\mathbf{G}_{Y}^{t}\right)^{-1}\right]+ italic_ε roman_tr [ bold_G start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ( italic_ε italic_n bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] + italic_ε roman_tr [ bold_G start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ( italic_ε italic_n bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + bold_G start_POSTSUBSCRIPT italic_Y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ]
2n(𝐇n𝐀t)T𝐊Xts(missingHn𝐀s)2𝑛subscriptnormsuperscriptsubscript𝐇𝑛superscript𝐀𝑡𝑇subscriptsuperscript𝐊𝑡𝑠𝑋missingsubscript𝐻𝑛superscript𝐀𝑠\displaystyle-\frac{2}{n}\left\|\left(\mathbf{H}_{n}{\mathbf{A}}^{t}\right)^{T% }{\mathbf{K}}^{ts}_{X}\left(\mathbf{\mathbf{missing}}{H}_{n}{\mathbf{A}}^{s}% \right)\right\|_{*}- divide start_ARG 2 end_ARG start_ARG italic_n end_ARG ∥ ( bold_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_A start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT bold_K start_POSTSUPERSCRIPT italic_t italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_X end_POSTSUBSCRIPT ( roman_missing italic_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_A start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ) ∥ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT

3.4 COD-based DAR Modeling

According to Theorem 3.2 , when the conditional distributions PXsYsubscript𝑃conditionalsuperscript𝑋𝑠𝑌P_{X^{s}\mid Y}italic_P start_POSTSUBSCRIPT italic_X start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ∣ italic_Y end_POSTSUBSCRIPT and PXtYsubscript𝑃conditionalsuperscript𝑋𝑡𝑌P_{X^{t}\mid Y}italic_P start_POSTSUBSCRIPT italic_X start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ∣ italic_Y end_POSTSUBSCRIPT are aligned, a well-trained source predictor can gain promising performance on target domain. However, when the shift is significant, it is difficult for predictor to produce high-quality pseudo labels at the beginning. For this reason, the model is likely to miss the optimal optimization direction in the early training stage and thus fall into sub-optimal solutions. Therefore, the marginal alignment can improve the reliability of initial pseudo labels, and further benefit the conditional alignment. Overall, the learning objective of COD can be formulated as

ming,hsrc+λ1d^CODmod+λ2d^KGW,subscript𝑔subscriptsrcsubscript𝜆1subscript^𝑑subscriptCODmodsubscript𝜆2subscript^𝑑KGW\displaystyle\min_{g,h}~{}\mathcal{L}_{\mathrm{src}}+\lambda_{1}\widehat{d}_{% \mathrm{COD_{mod}}}+\lambda_{2}\widehat{d}_{\mathrm{KGW}},roman_min start_POSTSUBSCRIPT italic_g , italic_h end_POSTSUBSCRIPT caligraphic_L start_POSTSUBSCRIPT roman_src end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT roman_COD start_POSTSUBSCRIPT roman_mod end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT roman_KGW end_POSTSUBSCRIPT , (16)

where λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, λ2subscript𝜆2\lambda_{2}italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are the trade-off parameters. For regression task, the source risk srcsubscriptsrc\mathcal{L}_{\mathrm{src}}caligraphic_L start_POSTSUBSCRIPT roman_src end_POSTSUBSCRIPT is commonly set as Mean Square Error (MSE) over continuous outputs:

src=1nsi=1ns𝐲is^𝐲is22subscriptsrc1superscript𝑛𝑠superscriptsubscript𝑖1superscript𝑛𝑠superscriptsubscriptdelimited-∥∥^subscriptsuperscript𝐲𝑠𝑖superscriptsubscript𝐲𝑖𝑠22\displaystyle\mathcal{L}_{\mathrm{src}}=\frac{1}{n^{s}}\sum_{i=1}^{n^{s}}\left% \lVert\hat{\mathbf{y}^{s}_{i}}-\mathbf{y}_{i}^{s}\right\rVert_{2}^{2}caligraphic_L start_POSTSUBSCRIPT roman_src end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_n start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ∥ over^ start_ARG bold_y start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG - bold_y start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (17)

where 𝐲is^^subscriptsuperscript𝐲𝑠𝑖\hat{\mathbf{y}^{s}_{i}}over^ start_ARG bold_y start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG is the predicted value of the i𝑖iitalic_i-th source sample. Note that the marginal alignment term can be replaced by any other discrepancy metric, while we adopt KGW here to preserve the consistency with OT framework.

4 Experiments

We evaluate our method with several state-of-the-art domain adaptation methods on three regression benchmark datasets, implementation details and additional illustrations are provided in supplementary material.

dSprites [19] is a standard 2D synthetic dataset for deep representation learning, which contains three domains each with 737,280 images: Color (C), Noisy (N) and Scream (S). Following [4], we consider three factors for regression tasks: scale and two plane coordinates, while the orientation factor is excluded for its complexity. We evaluate all methods on six transfer tasks: CNCN\textbf{C}\to\textbf{N}C → N, CSCS\textbf{C}\to\textbf{S}C → S, NCNC\textbf{N}\to\textbf{C}N → C, NSNS\textbf{N}\to\textbf{S}N → S, SCSC\textbf{S}\to\textbf{C}S → C and SNSN\textbf{S}\to\textbf{N}S → N. Following previous works [4, 38], the sum of Mean Absolute Error (MAE) on all three sub-regression tasks is reported as the evaluation metric.

Biwi kinect [12] is a real-world dataset containing two domains according to gender: Female (F) with 5874 images and Male (M) with 9804 images. The three factors that can be employed for regression tasks are pitch, yaw and roll. We evaluate all methods on two transfer tasks: FMFM\textbf{F}\to\textbf{M}F → M and MFMF\textbf{M}\to\textbf{F}M → F. The sum of MAE on three sub-regression tasks is reported.

MPI3D [16] is a simulation-to-real dataset of 3D objects which consist of three different domain: Toy (T), RealistiC (RC) and ReaL (RL). Each domain contains 1,036,800 images. There are two factors in MPI3D that can be employed for regression tasks: a rotation about a vertical axis and a second rotation about a horizontal axis. We evaluate all methods on six transfer tasks: RLRCRLRC\textbf{RL}\to\textbf{RC}RL → RC, RLTRLT\textbf{RL}\to\textbf{T}RL → T, RCTRCT\textbf{RC}\to\textbf{T}RC → T, RCRLRCRL\textbf{RC}\to\textbf{RL}RC → RL, TRLTRL\textbf{T}\to\textbf{RL}T → RL and TRCTRC\textbf{T}\to\textbf{RC}T → RC. The sum of MAE on two sub-regression tasks is reported.

Comparison methods. SOTA methods for UDA and DAR are selected: TCA [40], MCD [45], DAN [28], DANN [14], JDOT [8], AFN [58], RSD [4], DARE-GRAM [38]. Following previous works [4, 38], Mean Absolute Error (MAE) is used as our evaluation metric across all regression tasks. In comparison experiment, the proposed DAR model in Eq. (16) is denoted as COD for convenience.

4.1 Results and Analysis

Comparison on dSprites. As shown in Tab. 1, COD achieves the best performance among all competing methods in nearly all tasks of dSprites. In CNCN\textbf{C}\to\textbf{N}C → N and CSCS\textbf{C}\to\textbf{S}C → S, in which the Resnet-18 and UDA methods exhibit terrible performance because of the significant domain gap. Although SOTA DAR methods have greatly boosted the performance compared with classification methods, COD still achieves a significant improvement, reducing the error to only 40%percent4040\%40 % of that produced by SOTA methods. These results demonstrate that the conditional invariant representations learned by COD indeed ensure better discriminability.

Comparison on Biwi Kinect. Since Biwi Kinect is the closest to real-world scenarios among the three benchmark datasets, the performance improvement of model is more challenging. In LABEL:{tab:result_BiwiKinect}, the performance gap between UDA methods and DAR methods is relatively smaller, which implies the SOTA DAR methods with marginal alignment indeed show limitations in the difficult transfer tasks. In such a scenario, the sufficiency of condition alignment is directly validated, where the MAE of COD is significantly lower than other SOTA methods. Thus, for real-world regression applications with complex shift, the COD model indeed shows superiority in both theoretical and methodological aspects.

Comparison on MPI3D. In MPI3D, the resolution remains the same with dSprites, while the main challenge is the large sample size. In Tab. 3, the MAE of COD is competitive with the SOTA method DARE-GRAM in the first four tasks, while higher in the last two tasks. However, note that COD still achieves the second-best result in averaged MAE, and obtains the lowest MAE in task RL\toT. In conclusion, the overall MAE of COD on three benchmarks is still lower, which demonstrates COD is generally superior to SOTA methods.

Table 1: Domain adaptation regression results (MAE) on dSprites (ResNet-18).
Method C\to N C\to S N\to C N\to S S\to C S\to N Avg.
Resnet-18 0.94 0.90 0.16 0.65 0.08 0.26 0.498
TCA 0.94 0.87 0.19 0.66 0.10 0.23 0.498
MCD 0.81 0.81 0.17 0.65 0.07 0.19 0.450
JDOT 0.86 0.79 0.19 0.64 0.10 0.23 0.468
AFN 1.00 0.96 0.16 0.62 0.08 0.32 0.523
DAN 0.70 0.77 0.12 0.50 0.06 0.11 0.377
DANN 0.47 0.46 0.16 0.65 0.05 0.10 0.315
RSD 0.31 0.31 0.12 0.53 0.07 0.08 0.237
DARE-GRAM 0.30 0.20 0.11 0.25 0.05 0.07 0.164
COD 0.12 0.16 0.10 0.23 0.08 0.12 0.134
Table 2: Domain adaptation regression results (MAE) on Biwi Kinect (ResNet-18).
Method M\to F F\to M Avg.
Resnet-18 0.29 0.38 0.335
TCA 0.31 0.39 0.350
MCD 0.31 0.37 0.340
JDOT 0.29 0.39 0.340
AFN 0.32 0.41 0.365
DAN 0.28 0.37 0.325
DANN 0.30 0.37 0.335
RSD 0.26 0.30 0.280
DARE-GRAM 0.23 0.29 0.260
COD 0.20 0.21 0.205
Table 3: Domain adaptation regression results (MAE) on MPI3D (ResNet-18).
Method RL\to RC RL\to T RC\to RL RC\to T T\to RL T\to RC Avg.
Resnet-18 0.17 0.44 0.19 0.45 0.51 0.50 0.377
TCA 0.17 0.42 0.19 0.42 0.50 0.50 0.373
MCD 0.13 0.40 0.15 0.45 0.52 0.50 0.358
JDOT 0.16 0.41 0.16 0.41 0.47 0.47 0.353
AFN 0.18 0.45 0.20 0.46 0.53 0.53 0.390
DAN 0.12 0.35 0.12 0.27 0.40 0.41 0.278
DANN 0.09 0.24 0.11 0.41 0.48 0.37 0.283
RSD 0.09 0.19 0.08 0.15 0.36 0.36 0.205
DARE-GRAM 0.09 0.15 0.10 0.14 0.24 0.24 0.160
COD 0.10 0.11 0.13 0.15 0.37 0.29 0.192
Table 4: Ablation study on different metrics for DAR.
Objectives dSprites Biwi Kinect MPI3D
MSE dKGWsubscript𝑑KGWd_{\mathrm{KGW}}italic_d start_POSTSUBSCRIPT roman_KGW end_POSTSUBSCRIPT dCODsubscript𝑑CODd_{\mathrm{COD}}italic_d start_POSTSUBSCRIPT roman_COD end_POSTSUBSCRIPT dCODmodsubscript𝑑subscriptCOD𝑚𝑜𝑑{d}_{\mathrm{COD}_{mod}}italic_d start_POSTSUBSCRIPT roman_COD start_POSTSUBSCRIPT italic_m italic_o italic_d end_POSTSUBSCRIPT end_POSTSUBSCRIPT C\to S N\to S T\to RL T\to RC M\to F F\to M
\checkmark 0.90 0.65 0.29 0.38 0.51 0.50
\checkmark \checkmark 0.18 0.26 0.21 0.25 0.41 0.34
\checkmark \checkmark 0.26 0.45 0.21 0.22 0.40 0.31
\checkmark \checkmark \checkmark 0.17 0.25 0.21 0.22 0.41 0.29
\checkmark \checkmark \checkmark 0.16 0.23 0.20 0.21 0.37 0.29

Ablation study. The ablation experiment was conducted on the most challenging tasks from all three datasets: CSCS\textbf{C}\to\textbf{S}C → S and NSNS\textbf{N}\to\textbf{S}N → S from dSprites, TRLTRL\textbf{T}\to\textbf{RL}T → RL and TRCTRC\textbf{T}\to\textbf{RC}T → RC from MPI3D, and both MFMF\textbf{M}\to\textbf{F}M → F and FMFM\textbf{F}\to\textbf{M}F → M from Biwi Kinect. The major components of proposed method, i.e., MSE objective srcsubscript𝑠𝑟𝑐\mathcal{L}_{src}caligraphic_L start_POSTSUBSCRIPT italic_s italic_r italic_c end_POSTSUBSCRIPT, KGW distance dKGWsubscript𝑑KGWd_{\mathrm{KGW}}italic_d start_POSTSUBSCRIPT roman_KGW end_POSTSUBSCRIPT, COD distance dCODsubscript𝑑CODd_{\mathrm{COD}}italic_d start_POSTSUBSCRIPT roman_COD end_POSTSUBSCRIPT and modified COD dCODmodsubscript𝑑subscriptCOD𝑚𝑜𝑑{d}_{\mathrm{COD}_{mod}}italic_d start_POSTSUBSCRIPT roman_COD start_POSTSUBSCRIPT italic_m italic_o italic_d end_POSTSUBSCRIPT end_POSTSUBSCRIPT, are evaluated, where the results are shown in Tab. 4.

From the 2nd and 3rd rows, it can be observed that both OT-based metrics show promising results on DAR datasets, where the conditional alignment with COD shows better performance than marginal alignment with KGW on MPI3D and Biwi Kinect datasets. On the other hand, since dSprites induce more difficult transfer problems, the baseline performance is generally lower, and the pseudo labels for conditional alignment are less reliable. Therefore, the combination of KGW and COD is necessary to ensure higher performance on challenging tasks, where the results in the 4th row validate the mutually beneficial relation between these two objectives. Finally, the 5th row shows the effectiveness of rethinking the discriminability of COD, which implies the modifications on the moment statistics indeed ensure lower MAE values on all datasets. Overall, these results demonstrate that conditional alignment is indeed crucial for successful DAR, and the theoretical results and proposed method are generally valid.

Refer to caption
Figure 3: t-SNE visualization of learned representations. Label values are denoted by color gradients. Eight label values are selected from the range of the variable for visualization. ’+++’: source samples; ’\circ’: target samples.

Feature Visualization. t-SNE [34] is employed to visualize the 2-D features of different alignments on dSprites C\toN task. (a) Before adaptation, distributions of source and target domain differ severely, as shown in Fig. 3(a), which implies predictor is infeasible for the target samples with a significant domain gap. (b) After marginal distribution alignment via the SOTA DARE-GRAM method, source and target representations are aligned globally, and the cross-domain distributions are overlapped on the whole. However, there are still misaligned local patterns where source and target samples of different colors are matched. Therefore, it is hard to fit clear level sets in Fig. 3 (b). (c) After conditional distribution alignment via COD, source and target distributions are aligned not only globally but also locally. As shown in Fig. 3 (c), representations are matched according to label values so the change of color shows a clear gradient and level sets suitable for both domains can be fitted.

Refer to caption
Figure 4: MAE values under different settings of hyper-parameters.

Hyper-parameter. We investigate the sensitivity of hyper-parameters on dSprites C\toN. MAE values by varying λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and λ2subscript𝜆2\lambda_{2}italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, i.e., parameters of modified COD and KGW, are shown in Fig. 4. The results show that the performance of the COD-based DAR model is generally robust to the different choice of parameters. In fact, the MAE values of the COD-based model are stable in random experiments, which implies the proposed method is both robust from the view of randomness and parameter setting. In conclusion, the stable performance validates that the proposed model is also empirically reliable.

5 Conclusion

In this work, we studied the limitations of existing DAR research, i.e., the misaligned intrinsic discriminant structure and negative transfer of marginal distribution matching framework. To deal with the challenges of conditional shift with continuous output space, we provide feasible solutions from both theoretical and methodological views, where the main results show that: 1) the cross-domain generalization regression error can be sufficiently bounded by the conditional discrepancy on continuous labels; 2) the proposed COD serves as a well-defined metric for conditional discrepancy measures and empirical modeling. With these guarantees, a COD-based DAR method is proposed with a discriminability enhancement mechanism. Numerical validations on standard DAR datasets validate the superiority of COD.

Acknowledgements

This work is supported in part by National Natural Science Foundation of China (Grant No. 62376291), in part by Guangdong Basic and Applied Basic Research Foundation (2023B1515020004), in part by Science and Technology Program of Guangzhou (2024A04J6413), in part by the Fundamental Research Funds for the Central Universities, Sun Yat-sen University (24xkjc013), in part by Guangdong Province Key Laboratory of Computational Science at Sun Yat-sen University (2020B1212060032), and in part by Key Laboratory of Machine Intelligence and Advanced Computing, Ministry of Education.

References

  • [1] Bao, Y., Liu, Y., Wang, H., Lu, F.: Generalizing gaze estimation with rotation consistency. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 4207–4216 (2022)
  • [2] Ben-David, S., Blitzer, J., Crammer, K., Kulesza, A., Pereira, F., Vaughan, J.W.: A theory of learning from different domains. Mach. Learn. 79, 151–175 (2010)
  • [3] Chen, L., Chen, H., Wei, Z., Jin, X., Tan, X., Jin, Y., Chen, E.: Reusing the task-specific classifier as a discriminator: Discriminator-free adversarial domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 7181–7190 (2022)
  • [4] Chen, X., Wang, S., Wang, J., Long, M.: Representation subspace distance for domain adaptation regression. In: Int. Conf. Mach. Learn. pp. 1749–1759 (2021)
  • [5] Chen, Y., Song, S., Li, S., Wu, C.: A graph embedding framework for maximum mean discrepancy-based domain adaptation algorithms. IEEE Trans. Image Process. 29, 199–213 (2019)
  • [6] Tachet des Combes, R., Zhao, H., Wang, Y.X., Gordon, G.J.: Domain adaptation with conditional distribution matching and generalized label shift. In: Adv. Neural Inform. Process. Syst. vol. 33, pp. 19276–19289 (2020)
  • [7] Cortes, C., Mohri, M.: Domain adaptation in regression. In: ICALT. pp. 308–323 (2011)
  • [8] Courty, N., Flamary, R., Habrard, A., Rakotomamonjy, A.: Joint distribution optimal transportation for domain adaptation. Adv. Neural Inform. Process. Syst. 30 (2017)
  • [9] Courty, N., Flamary, R., Tuia, D.: Domain adaptation with regularized optimal transport. In: ECML PKDD. pp. 274–289 (2014)
  • [10] Damodaran, B.B., Kellenberger, B., Flamary, R., Tuia, D., Courty, N.: Deepjdot: Deep joint distribution optimal transport for unsupervised domain adaptation. In: Eur. Conf. Comput. Vis. pp. 447–463 (2018)
  • [11] Dhaini, M., Berar, M., Honeine, P., Van Exem, A.: Unsupervised domain adaptation for regression using dictionary learning. Knowl.-Based Syst. 267, 110439 (2023)
  • [12] Fanelli, G., Dantone, M., Gall, J., Fossati, A., Van Gool, L.: Random forests for real time 3d face analysis. Int. J. Comput. Vis. 101, 437–458 (2013)
  • [13] Fukumizu, K., Bach, F.R., Jordan, M.I.: Kernel dimension reduction in regression. Ann. Stat. 37(4), 1871–1905 (2009)
  • [14] Ganin, Y., Ustinova, E., Ajakan, H., Germain, P., Larochelle, H., Laviolette, F., March, M., Lempitsky, V.: Domain-adversarial training of neural networks. J. Mach. Learn. Res. 17(59), 1–35 (2016)
  • [15] Geng, X., Zhou, Z.H., Smith-Miles, K.: Automatic age estimation based on facial aging patterns. IEEE Trans. Pattern Anal. Mach. Intell. 29(12), 2234–2240 (2007)
  • [16] Gondal, M.W., Wuthrich, M., Miladinovic, D., Locatello, F., Breidt, M., Volchkov, V., Akpo, J., Bachem, O., Schölkopf, B., Bauer, S.: On the transfer of inductive bias from simulation to the real world: a new disentanglement dataset. Adv. Neural Inform. Process. Syst. 32 (2019)
  • [17] Gong, B., Shi, Y., Sha, F., Grauman, K.: Geodesic flow kernel for unsupervised domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 2066–2073 (2012)
  • [18] Goswami, S., Kontolati, K., Shields, M.D., Karniadakis, G.E.: Deep transfer operator learning for partial differential equations under conditional shift. Nat. Mach. Intell. 4(12), 1155–1164 (2022)
  • [19] Higgins, I., Matthey, L., Pal, A., Burgess, C., Glorot, X., Botvinick, M., Mohamed, S., Lerchner, A.: beta-vae: Learning basic visual concepts with a constrained variational framework. In: Int. Conf. Learn. Represent. (2016)
  • [20] Hoffman, J., Tzeng, E., Park, T., Zhu, J.Y., Isola, P., Saenko, K., Efros, A., Darrell, T.: Cycada: Cycle-consistent adversarial domain adaptation. In: Int. Conf. Mach. Learn. pp. 1989–1998 (2018)
  • [21] Hu, J., Zhong, H., Yang, F., Gong, S., Wu, G., Yan, J.: Learning unbiased transferability for domain adaptation by uncertainty modeling. In: Eur. Conf. Comput. Vis. pp. 223–241 (2022)
  • [22] Jiang, J., Ji, Y., Wang, X., Liu, Y., Wang, J., Long, M.: Regressive domain adaptation for unsupervised keypoint detection. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 6780–6789 (2021)
  • [23] Kirchmeyer, M., Rakotomamonjy, A., de Bezenac, E., Gallinari, P.: Mapping conditional distributions for domain adaptation under generalized target shift. In: Int. Conf. Learn. Represent. (2022)
  • [24] Lee, C.Y., Batra, T., Baig, M.H., Ulbricht, D.: Sliced wasserstein discrepancy for unsupervised domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 10285–10295 (2019)
  • [25] Li, B., Wang, Y., Zhang, S., Li, D., Keutzer, K., Darrell, T., Zhao, H.: Learning invariant representations and risks for semi-supervised domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 1104–1113 (2021)
  • [26] Li, X., Li, J., Zuo, L., Zhu, L., Shen, H.T.: Domain adaptive remaining useful life prediction with transformer. IEEE Trans. Instrum. 71, 1–13 (2022)
  • [27] Li, Z., Zhang, K., Liu, Y., Zou, Y., Ding, G.: A novel remaining useful life transfer prediction method of rolling bearings based on working conditions common benchmark. IEEE Trans. Instrum. 71,  1–9 (2022)
  • [28] Long, M., Cao, Y., Wang, J., Jordan, M.: Learning transferable features with deep adaptation networks. In: Int. Conf. Mach. Learn. pp. 97–105 (2015)
  • [29] Long, M., Cao, Z., Wang, J., Jordan, M.I.: Conditional adversarial domain adaptation. Adv. Neural Inform. Process. Syst. 31 (2018)
  • [30] Long, M., Zhu, H., Wang, J., Jordan, M.I.: Deep transfer learning with joint adaptation networks. In: Int. Conf. Mach. Learn. pp. 2208–2217 (2017)
  • [31] Luo, Y.W., Ren, C.X.: Conditional bures metric for domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 13989–13998 (2021)
  • [32] Luo, Y.W., Ren, C.X.: Mot: Masked optimal transport for partial domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 3531–3540 (2023)
  • [33] Luo, Y.W., Ren, C.X., Dai, D.Q., Yan, H.: Unsupervised domain adaptation via discriminative manifold propagation. IEEE Trans. Pattern Anal. Mach. Intell. 44(3), 1653–1669 (2022)
  • [34] Van der Maaten, L., Hinton, G.: Visualizing data using t-sne. J. Mach. Learn. Res. 9(11) (2008)
  • [35] Mansour, Y., Mohri, M., Rostamizadeh, A.: Domain adaptation: Learning bounds and algorithms. arXiv preprint arXiv:0902.3430 (2009)
  • [36] Mao, W., Ge, Y., Shen, C., Tian, Z., Wang, X., Wang, Z., den Hengel, A.v.: Poseur: Direct human pose regression with transformers. In: Eur. Conf. Comput. Vis. pp. 72–88 (2022)
  • [37] de Mathelin, A., Richard, G., Deheeger, F., Mougeot, M., Vayatis, N.: Adversarial weighting for domain adaptation in regression. In: ICTAI. pp. 49–56 (2021)
  • [38] Nejjar, I., Wang, Q., Fink, O.: Dare-gram: Unsupervised domain adaptation regression by aligning inverse gram matrices. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 11744–11754 (2023)
  • [39] Nikzad-Langerodi, R., Zellinger, W., Saminger-Platz, S., Moser, B.A.: Domain adaptation for regression under beer–lambert’s law. Knowl.-Based Syst. 210, 106447 (2020)
  • [40] Pan, S.J., Tsang, I.W., Kwok, J.T., Yang, Q.: Domain adaptation via transfer component analysis. IEEE Trans. Neural Netw. 22(2), 199–210 (2010)
  • [41] Pei, Z., Cao, Z., Long, M., Wang, J.: Multi-adversarial domain adaptation. In: AAAI. vol. 32 (2018)
  • [42] Peng, R., Wang, R., Wang, Z., Lai, Y., Wang, R.: Rethinking depth estimation for multi-view stereo: A unified representation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 8645–8654 (2022)
  • [43] Quinonero-Candela, J., Sugiyama, M., Schwaighofer, A., Lawrence, N.D.: Dataset shift in machine learning. Mit Press (2008)
  • [44] Ren, C.X., Luo, Y.W., Dai, D.Q.: Buresnet: Conditional bures metric for transferable representation learning. IEEE Trans. Pattern Anal. Mach. Intell. 45(4), 4198–4213 (2023)
  • [45] Saito, K., Watanabe, K., Ushiku, Y., Harada, T.: Maximum classifier discrepancy for unsupervised domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 3723–3732 (2018)
  • [46] Shen, J., Qu, Y., Zhang, W., Yu, Y.: Wasserstein distance guided representation learning for domain adaptation. In: AAAI. vol. 32 (2018)
  • [47] Singh, A., Chakraborty, S.: Deep domain adaptation for regression. Development and Analysis of Deep Learning Architectures pp. 91–115 (2020)
  • [48] Song, L., Huang, J., Smola, A., Fukumizu, K.: Hilbert space embeddings of conditional distributions with applications to dynamical systems. In: Int. Conf. Mach. Learn. pp. 961–968 (2009)
  • [49] Sun, B., Feng, J., Saenko, K.: Return of frustratingly easy domain adaptation. In: AAAI. vol. 30 (2016)
  • [50] Tang, H., Jia, K.: Discriminative adversarial domain adaptation. In: AAAI. vol. 34, pp. 5940–5947 (2020)
  • [51] Teshima, T., Sato, I., Sugiyama, M.: Few-shot domain adaptation by causal mechanism transfer. In: Int. Conf. Mach. Learn. pp. 9458–9469 (2020)
  • [52] Tzeng, E., Hoffman, J., Saenko, K., Darrell, T.: Adversarial discriminative domain adaptation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 7167–7176 (2017)
  • [53] Wang, H., Shen, T., Zhang, W., Duan, L.Y., Mei, T.: Classes matter: A fine-grained adversarial approach to cross-domain semantic segmentation. In: Eur. Conf. Comput. Vis. pp. 642–659 (2020)
  • [54] Wang, W., Li, H., Ding, Z., Nie, F., Chen, J., Dong, X., Wang, Z.: Rethinking maximum mean discrepancy for visual domain adaptation. IEEE Trans. Neural Netw. Learn. Syst. (2021)
  • [55] Wang, Y., Jiang, Y., Li, J., Ni, B., Dai, W., Li, C., Xiong, H., Li, T.: Contrastive regression for domain adaptation on gaze estimation. In: IEEE Conf. Comput. Vis. Pattern Recog. pp. 19376–19385 (2022)
  • [56] Wei, P., Ke, Y., Ong, Y.S., Ma, Z.: Adaptive transfer kernel learning for transfer gaussian process regression. IEEE Trans. Pattern Anal. Mach. Intell. (2022)
  • [57] Wu, J., He, J., Wang, S., Guan, K., Ainsworth, E.: Distribution-informed neural networks for domain adaptation regression. Adv. Neural Inform. Process. Syst. 35, 10040–10054 (2022)
  • [58] Xu, R., Li, G., Yang, J., Lin, L.: Larger norm more transferable: An adaptive feature norm approach for unsupervised domain adaptation. In: Int. Conf. Comput. Vis. pp. 1426–1435 (2019)
  • [59] Zhang, Z., Wang, M., Nehorai, A.: Optimal transport in reproducing kernel hilbert spaces: Theory and applications. IEEE Trans. Pattern Anal. Mach. Intell. 42(7), 1741–1754 (2019)
  • [60] Zhao, H., Combes, R.T.D., Zhang, K., Gordon, G.: On learning invariant representations for domain adaptation. In: Int. Conf. Mach. Learn. pp. 7523–7532 (2019)
  • [61] Zhao, H., Coston, A., Adel, T., Gordon, G.J.: Conditional learning of fair representations. In: Int. Conf. Learn. Represent. (2020)
  • [62] Zhao, H., Dan, C., Aragam, B., Jaakkola, T.S., Gordon, G.J., Ravikumar, P.: Fundamental limits and tradeoffs in invariant representation learning. J. Mach. Learn. Res. 23(340), 1–49 (2022)