Fewshot Domain Adaptation by Causal Mechanism Transfer
Abstract
We study fewshot supervised domain adaptation (DA) for regression problems, where only a few labeled target domain data and many labeled source domain data are available. Many of the current DA methods base their transfer assumptions on either parametrized distribution shift or apparent distribution similarities, e.g., identical conditionals or small distributional discrepancies. However, these assumptions may preclude the possibility of adaptation from intricately shifted and apparently very different distributions. To overcome this problem, we propose mechanism transfer, a metadistributional scenario in which a data generating mechanism is invariant among domains. This transfer assumption can accommodate nonparametric shifts resulting in apparently different distributions while providing a solid statistical basis for DA. We take the structural equations in causal modeling as an example and propose a novel DA method, which is shown to be useful both theoretically and experimentally. Our method can be seen as the first attempt to fully leverage the structural causal models for DA.
1 Introduction
Learning from a limited amount of data is a longstanding yet actively studied problem of machine learning. Domain adaptation (DA) (BenDavid et al., 2010) tackles this problem by leveraging auxiliary data sampled from related but different domains. In particular, we consider fewshot supervised DA for regression problems, where only a few labeled target domain data and many labeled source domain data are available.
A key component of DA methods is the transfer assumption (TA) to relate the source and the target distributions. Many of the previously explored TAs have relied on certain direct distributional similarities, e.g., identical conditionals (Shimodaira, 2000) or small distributional discrepancies (BenDavid et al., 2007). However, these TAs may preclude the possibility of adaptation from apparently very different distributions. Many others assume parametric forms of the distribution shift (Zhang et al., 2013) or the distribution family (Storkey and Sugiyama, 2007) which can highly limit the considered set of distributions. (we further review related work in Section 5.1).
To alleviate the intrinsic limitation of previous TAs due to relying on apparent distribution similarities or parametric assumptions, we focus on a metadistributional scenario where there exists a common generative mechanism behind the data distributions (Figures 1,2). Such a common mechanism may be more conceivable in applications involving structured table data such as medical records (Yadav et al., 2018). For example, in medical record analysis for disease risk prediction, it can be reasonable to assume that there is a pathological mechanism that is common across regions or generations, but the data distributions may vary due to the difference in cultures or lifestyles. Such a hidden structure (pathological mechanism, in this case), once estimated, may provide portable knowledge to enable DA, allowing one to obtain accurate predictors for underinvestigated regions or new generations.
Concretely, our assumption relies on the generative model of nonlinear independent component analysis (nonlinear ICA; Figure 1), where the observed labeled data are generated by first sampling latent independent components (ICs) and later transforming them by a nonlinear invertible mixing function denoted by (Hyvärinen et al., 2019). Under this generative model, our TA is that representing the mechanism is identical across domains (Figure 2). This TA allows us to formally relate the domain distributions and develop a novel DA method without assuming their apparent similarities or making parametric assumptions.
Our contributions.
Our key contributions can be summarized in three points as follows.

We formulate the flexible yet intuitively accessible TA of shared generative mechanism and develop a fewshot regression DA method (Section 3). The idea is as follows. First, from the source domain data, we estimate the mixing function by nonlinear ICA (Hyvärinen et al., 2019) because is the only assumed relation of the domains. Then, to transfer the knowledge, we perform data augmentation using the estimated on the target domain data using the independence of the IC distributions. In the end, the augmented data is used to fit a target predictor (Figure 3).

We theoretically justify the augmentation procedure by invoking the theory of generalized Ustatistics (Lee, 1990). The theory shows that the proposed data augmentation procedure yields the uniformly minimum variance unbiased risk estimator in an ideal case. We also provide an excess risk bound (Mohri et al., 2012) to cover a more realistic case (Section 4).
2 Problem Setup
In this section, we describe the problem setup and the notation. To summarize, our problem setup is homogeneous, multisource, and fewshot supervised domain adapting regression. That is, respectively, all data distributions are defined on the same data space, there are multiple source domains, and a limited number of labeled data is available from the target distribution (and we do not assume the availability of unlabeled data). In this paper, we use the terms domain and distribution interchangeably.
Notation.
Let us denote the set of real (resp. natural) numbers by (resp. ). For , we define . Throughout the paper, we fix and suppose that the input space is a subset of and the label space is a subset of . As a result, the overall data space is a subset of . We generally denote a labeled data point by . We denote by the set of independent distributions on with absolutely continuous marginals. For a distribution , we denote its induced expectation operator by . Table 3 in Supplementary Material provides a summary of notation.
Basic setup: Fewshot domain adaptation.
Let be a distribution (the target distribution) over , and let be a hypothesis class. Let be a loss function where is a constant. Our goal is to find a predictor which performs well for , i.e., the target risk is small. We denote . To this goal, we are given an independent and identically distributed (i.i.d.) sample . In a fully supervised setting where is large, a standard procedure is to select by empirical risk minimization (ERM), i.e., , where . However, when is not sufficiently large, may not accurately estimate , resulting in a high generalization error of . To compensate for the scarcity of data from the target distribution, let us assume that we have data from distinct source distributions over , that is, we have independent i.i.d. samples whose relations to are described shortly. We assume for simplicity.
Key assumption.
In this work, the key transfer assumption is that all domains follow nonlinear ICA models with identical mixing functions (Figure 2). To be precise, we assume that there exists a set of IC distributions , and a smooth invertible function (the transformation or mixing) such that is generated by first sampling and later transforming it by
(1) 
and similarly for . The above assumption allows us to formally relate and . It also allows us to estimate when sufficient identification conditions required by the theory of nonlinear ICA are met. Due to space limitation, we provide a brief review of the the nonlinear ICA method used in this paper and the known theoretical conditions in Supplementary Material A. The requirement for multiple source domains comes from the currently known identification condition of nonlinear ICA. Note that complex changes in are allowed, hence the assumption of invariant can accommodate intricate shifts in the apparent distribution . We discuss this further in Section 5.3 by taking a simple example.
Example: Structural equation models
A salient example of generative models expressed as Eq. (1) is structural equation models (SEMs; Pearl, 2009; Peters et al., 2017), more precisely, the reduced form (Reiss and Wolak, 2007) of Markovian SEMs (Pearl, 2009) under adequate assumptions such as invertibility. SEMs are used to describe the data generating mechanism involving the causality of random variables (Pearl, 2009). This interpretation of SEMs as Eq.(1) has been exploited in methods of causal structure discovery such as the linear nonGaussian additivenoise models and their successors (Kano and Shimizu, 2003; Shimizu et al., 2006; Monti et al., 2019). In the case of SEMs, the key assumption of this paper translates into the invariance of the structural equations among domains. Note that we do not require the estimation of the exact structural equations, but only the reduced form needs to be estimated. Since the reducedform equation can be derived given the structuralform equations by recursive imputation, we are only required to solve an easier estimation problem compared to causal discovery.
3 Proposed Method: Mechanism Transfer
In this section, we detail the proposed method, mechanism transfer (Algorithm 1). The method first estimates the common generative mechanism from the source domain data and then uses it to perform data augmentation of the target domain data to transfer the knowledge (Figure 3).
3.1 Step 1: Estimate using the source domain data
The first step estimates the common transformation by nonlinear ICA, namely via generalized contrastive learning (GCL; Hyvärinen et al., 2019). GCL uses auxiliary information for training a certain binary classification function, , equipped with a parametrized feature extractor . The trained feature extractor is used as an estimator of . The auxiliary information we use in our problem setup is the domain indices . The classification function to be trained in GCL is consisting of , and the classification task of GCL is logistic regression to classify as positive and as negative. This yields the following domaincontrastive learning criterion to estimate :
where and are sets of parametrized functions, denotes the expectation with respect to ( denotes the uniform distribution), and is the logistic loss . We use the solution as an estimator of . In experiments, is implemented by invertible neural networks (Kingma and Dhariwal, 2018), by multilayer perceptron, and is replaced by a random sampling renewed for every minibatch.
3.2 Step 2: Extract and inflate the target ICs using
The second step extracts and inflates the target domain ICs using the estimated . We first extract the ICs of the target domain data by applying the inverse of as
After the extraction, we inflate the set of IC values by taking all dimensionwise combinations of the estimated IC:
to obtain new plausible IC values . The intuitive motivation of this procedure stems from the independence of the IC distributions. Theoretical justifications are provided in Section 4. In our implementation, we use invertible neural networks (Kingma and Dhariwal, 2018) to model the function to enable the computation of the inverse .
3.3 Step 3: Synthesize target data from the inflated ICs
The third step estimates the target risk by the empirical distribution of the augmented data:
(2) 
and performs empirical risk minimization. In experiments, we use a regularization term to control the complexity of and select
The generated hypothesis is then used to make predictions in the target domain. In our experiments, we use , where and the norm is that of the reproducing kernel Hilbert space (RKHS) which we take the subset from. Note that we may well reduce the computing time by taking only a subset of combinations in Eq. (2).
4 Theoretical Insights
In this section, we state two theorems to investigate the statistical properties of the method proposed in Section 3 and provide plausibility beyond the intuition that we take advantage of the independence of the IC distributions.
4.1 Minimum variance property: Idealized case
The first theorem provides an insight into the statistical advantage of the proposed method: in the ideal case, the method attains the minimum variance among all possible unbiased risk estimators.
Theorem 1 (Minimum variance property of ).
Assume that . Then, for each , the proposed risk estimator is the uniformly minimum variance unbiased estimator of , i.e., for any unbiased estimator of ,
as well as holds.
The proof of Theorem 1 is immediate once we rewrite as a variate regular statistical functional and as its corresponding generalized Ustatistic (Lee, 1990). Details can be found in Supplementary Material D. Theorem 1 implies that the proposed risk estimator can have superior statistical efficiency in terms of the variance over the ordinary empirical risk.
4.2 Excess risk bound: More realistic case
In real situations, one has to estimate . The following theorem characterizes the statistical gain and loss arising from the estimation error . The intuition is that the increased number of points suppresses the possibility of overfitting because the hypothesis has to fit the majority of the inflated data, but the estimator has to be accurate so that fitting to the inflated data is meaningful. Note that the theorem is agnostic to how is obtained, hence it applies to more general problem setup as long as can be estimated.
Theorem 2 (Excess risk bound).
Let be the hypothesis generated by Eq. (2). Under appropriate assumptions (see Theorem 3 in Supplementary Material), for arbitrary , we have with probability at least ,
Here, is the Sobolev norm, and we define the effective Rademacher complexity by
(3) 
where are independent sign variables, is the expectation with respect to , the dummy variables are i.i.d. copies of , and is defined by using the degree symmetric group as
and and are higher order terms. The constants and depend only on and , respectively, while depends only on , and .
Details of the statement and the proof can be found in Supplementary Material C. The Sobolev norm (Adams and Fournier, 2003) emerges from the evaluation of the difference between the estimated IC distribution and the groundtruth IC distribution. In Theorem 2, the utility of the proposed method appears in the effective complexity measure. The complexity is defined by a set of functions which are marginalized over all but one argument, resulting in mitigated dependence on the input dimensionality from exponential to linear (Supplementary Material C, Remark 3).
5 Related Work and Discussion
In this section, we review some existing TAs for DA to clarify the relative position of the paper. We also clarify the relation to the literature of causalityrelated transfer learning.
5.1 Existing transfer assumptions
Here, we review some of the existing work and TAs. See Table 1 for a summary.
(1) Parametric assumptions.
Some TAs assume parametric distribution families, e.g., Gaussian mixture model in covariate shift (Storkey and Sugiyama, 2007). Some others assume parametric distribution shift, i.e., parametric representations of the target distribution given the source distributions. Examples include locationscale transform of class conditionals (Zhang et al., 2013; Gong et al., 2016), linearly dependent class conditionals (Zhang et al., 2015), and lowdimensional representation of the class conditionals after kernel embedding (Stojanov et al., 2019). In some applications, e.g., remote sensing, some parametric assumptions have proven useful (Zhang et al., 2013).
(2) Invariant conditionals and marginals.
Some methods assume invariance of certain conditionals or marginals (QuiñoneroCandela et al., 2009), e.g., in the covariate shift scenario (Shimodaira, 2000), for an appropriate feature transformation in transfer component analysis (Pan et al., 2011), for a feature selection in RojasCarulla et al. (2018), in the target shift (TarS) scenario (Zhang et al., 2013; Nguyen et al., 2016), and few components of regularvine copulas and marginals in Lopezpaz et al. (2012). For example, the covariate shift scenario has been shown to fit well to brain computer interface data (Sugiyama et al., 2007).
(3) Small discrepancy or integral probability metric.
Another line of work relies on certain distributional similarities, e.g., integral probability metric (Courty et al., 2017) or hypothesisclass dependent discrepancies (BenDavid et al., 2007; Blitzer et al., 2008; BenDavid et al., 2010; Kuroki et al., 2019; Zhang et al., 2019; Cortes et al., 2019). These methods assume the existence of the ideal joint hypothesis (BenDavid et al., 2010), corresponding to a relaxation of the covariate shift assumption. These TA are suited for unsupervised or semisupervised DA in computer vision applications (Courty et al., 2017).
(4) Transferable parameter.
Some others consider parameter transfer (Kumagai, 2016), where the TA is the existence of a parameterized feature extractor that performs well in the target domain for linearinparameter hypotheses and its learnability from the source domain data. For example, such a TA has been known to be useful in natural language processing or image recognition (Lee et al., 2009; Kumagai, 2016).
TA  AD  NP  Suited app. example 

(1) Parametric    Remote sensing  
(2) Invariant dist.    BCI  
(3) Disc. / IPM    Computer vision  
(4) Paramtransfer  Computer vision  
(Ours) Mechanism  Medical records 
5.2 Causality for transfer learning
Our method can be seen as the first attempt to fully leverage structural causal models for DA. Most of the causalityinspired DA methods express their assumptions in the level of graphical causal models (GCMs), which only has much coarser information than structural causal models (SCMs) (Peters et al., 2017, Table 1.1) exploited in this paper. Compared to previous work, our method takes one step further to assume and exploit the invariance of SCMs. Specifically, many studies assume the GCM (the anticausal scenario) following the seminal metaanalysis of Schölkopf et al. (2012) and use it to motivate their parametric distribution shift assumptions or the parameter estimation procedure (Zhang et al., 2013, 2015; Gong et al., 2016, 2018). Although such assumptions on the GCM have the virtue of being more robust to misspecification, they tend to require parametric assumptions to obtain theoretical justifications. On the other hand, our assumption enjoys a theoretical guarantee without relying on parametric assumptions.
5.3 Plausibility of the assumption
The invariance of causal mechanisms has been exploited in recent work of causal discovery such as Xu et al. (2014) and Monti et al. (2019), or under the name of the multienvironment setting in Ghassami et al. (2017). The SEMs are normally assumed to remain invariant unless explicitly intervened in (Hünermund and Bareinboim, 2019). As the first algorithm in the approach to fully exploit SCMs for DA, we consider the case where all variables are observable. Although it is often assumed in a causal inference problem that there are some unobserved confounding variables, we leave further extension to such a case for future work.
The relation between and can drastically change while is invariant. For example, even in a simple additive noise model , the conditional can shift drastically if the distribution of the independent noise changes in a complex manner, e.g., becoming multimodal from unimodal.
6 Experiment
In this section, we provide proofofconcept experiments to demonstrate the effectiveness of the proposed approach. Note that the primary purpose of the experiments is to confirm whether the proposed method can properly perform DA in realworld data, and it is not to determine which DA method and TA are the most suited for the specific dataset.
6.1 Implementation details of the proposed method
Estimation of (Step 1).
We model by an layer Glow neural network (Supplementary Material B.2). We model by a hiddenlayer neural network with a varied number of hidden units, output units, and the rectified linear unit activation (LeCun et al., 2015). We use its th output () as the value for . For training, we use the Adam optimizer (Kingma and Ba, 2017) with fixed parameters , fixed initial learning rate , and the maximum number of epochs . The other fixed hyperparameters of and its training process are described in Supplementary Material B.
Augmentation of target data (Step 3).
For each evaluation step, we take all combinations (with replacement) of the estimated ICs to synthesize target domain data. After we synthesize the data, we filter them by applying a novelty detection technique with respect to the union of source domain data. Namely, we use oneclass support vector machine (Schölkopf et al., 2000) with the fixed parameter and radial basis function (RBF) kernel with . This is because the estimated transform is not expected to be trained well outside the union of the supports of the source distributions.
Predictor hypothesis class .
As the predictor model, we use the kernel ridge regression (KRR) with RBF kernel. The bandwidth is chosen by the median heuristic similarly to Yamada et al. (2011) for simplicity. Note that the choice of the predictor model is for the sake of comparison with the other methods tailored for KRR (Cortes et al., 2019), and that an arbitrary predictor hypothesis class and learning algorithm can be easily combined with the proposed approach.
Hyperparameter selection.
We perform gridsearch for hyperparameter selection. The number of hidden units for is chosen from and the coefficient of weightdecay from . The regularization coefficient of KRR is chosen from following Cortes et al. (2019). To perform hyperparameter selection as well as earlystopping, we record the leaveoneout crossvalidation (LOOCV) meansquared error on the target training data every epochs and select its minimizer. The leaveoneout score is computed using the wellknown analytic formula instead of training the predictor for each split. Note that we only use the original target domain data and not the synthesized data as the heldout set.
Computation environment
6.2 Experiment using realworld data
Dataset.
We use the gasoline consumption data (Greene, 2012, p.284, Example 9.5), which is a panel data of gasoline usage in 18 of the OECD countries over 19 years. We consider each country as a domain, and we disregard the timeseries structure and consider the data as i.i.d. samples for each country in this proofofconcept experiment. The dataset contains four variables, all of which are logtransformed: motor gasoline consumption per car (the predicted variable), percapita income, motor gasoline price, and the stock of cars per capita (the predictor variables) (Baltagi and Griffin, 1983). For further details of the data, see Supplementary Material B. We used the dataset because there are very few public datasets for domain adapting regression tasks (Cortes and Mohri, 2014) especially for multisource DA, and also because the dataset has been used in econometric analyses involving SEMs (Baltagi, 2005), conforming to our approach.
Compared methods.
We compare the following DA methods, all of which apply to regression problems. Unless explicitly specified, the predictor class is chosen to be KRR with the same hyperparameter candidates as the proposed method (Section 6.1). Further details are described in Supplementary Material B.5.

Naive baselines (SrcOnly, TrgOnly, and S&TV): SrcOnly (resp. TrgOnly) trains a predictor on the source domain data (resp. target training data) without any device. SrcOnly can be effective if the source domains and the target domain have highly similar distributions. The S&TV baseline trains on both source and target domain data, but the LOOCV score is computed only from the target domain data.

TrAdaBoost: Twostage TrAdaBoost.R2; a boosting method tailored for fewshot regression transfer proposed in Pardoe and Stone (2010). It is an iterative method with earlystopping (Pardoe and Stone, 2010), for which we use the leaveoneout crossvalidation score on the target domain data as the criterion. As suggested in Pardoe and Stone (2010), we set the maximum number of outer loop iterations at . The base predictor is the decision tree regressor with the maximum depth (Hastie et al., 2009). Note that although TrAdaBoost does not have a clarified transfer assumption, we compare the performance for reference.

IW: Importance weighted KRR using RuLSIF (Yamada et al., 2011). The method directly estimates a relative joint density ratio function for , where is a hypothetical source distribution created by pooling all source domain data. Following Yamada et al. (2011), we experiment on and report the results separately. The regularization coefficient is selected from using importanceweighted crossvalidation (Sugiyama et al., 2007).

GDM: Generalized discrepancy minimization (Cortes et al., 2019). This method performs instanceweighted training on the source domain data with the weights that minimize the generalized discrepancy (via quadratic programming). We select the hyperparameters from as suggested by Cortes et al. (2019). The selection criterion is the performance of the trained predictor on the target training labels as the method trains on the source domain data and the target unlabeled data.

Copula: Nonparametric regularvine copula method (Lopezpaz et al., 2012). This method presumes using a specific joint density estimator called regularvine (Rvine) copulas. Adaptation is realized in two steps: the first step estimates which components of the constructed Rvine model are different by performing twosample tests based on maximum mean discrepancy (Lopezpaz et al., 2012), and the second step reestimates the components in which a change is detected using only the target domain data.

LOO (reference score): Leaveoneout crossvalidated error estimate is also calculated for reference. It is the average prediction error of predicting for a single heldout test point when the predictor is trained on the rest of the whole target domain data including those in the test set for the other algorithms.
Evaluation procedure.
The prediction accuracy was measured by the mean squared error (MSE). For each traintest split, we randomly select onethird (6 points) of the target domain dataset as the training set and use the rest as the test set. All experiments were repeated 10 times with different traintest splits of target domain data.
Results.
The results are reported in Table 2. We report the MSE scores normalized by that of LOO to facilitate the comparison, similarly to Cortes and Mohri (2014). In many of the target domain choices, the naive baselines (SrcOnly and S&TV) suffer from negative transfer, i.e., higher average MSE than TarOnly (in 12 out of 18 domains). On the other hand, the proposed method successfully performs better than TarOnly or is more resistant to negative transfer than the other compared methods. The performances of GDM, Copula, and IW are often inferior even compared to the baseline performance of SrcAndTarValid. For GDM and IW, this can be attributed to the fact that these methods presume the availability of abundant (unlabeled) target domain data, which is unavailable in the current problem setup. For Copula, the performance inferior to the naive baselines is possibly due to the restriction of the predictor model to its accompanied probability model (Lopezpaz et al., 2012). TrAdaBoost works reasonably well for many but not all domains. For some domains, it suffered from negative transfer similarly to others, possibly because of the very small number of training data points. Note that the transfer assumption of TrAdaBoost has not been stated (Pardoe and Stone, 2010), and it is not understood when the method is reliable.
Target  (LOO)  TrgOnly  Prop  SrcOnly  S&TV  TrAda  GDM  Copula  IW(.0)  IW(.5)  IW(.95) 

AUT  1  5.88 (1.60)  5.39 (1.86)  9.67 (0.57)  9.84 (0.62)  5.78 (2.15)  31.56 (1.39)  27.33 (0.77)  39.72 (0.74)  39.45 (0.72)  39.18 (0.76) 
BEL  1  10.70 (7.50)  7.94 (2.19)  8.19 (0.68)  9.48 (0.91)  8.10 (1.88)  89.10 (4.12)  119.86 (2.64)  105.15 (2.96)  105.28 (2.95)  104.30 (2.95) 
CAN  1  5.16 (1.36)  3.84 (0.98)  157.74 (8.83)  156.65 (10.69)  51.94 (30.06)  516.90 (4.45)  406.91 (1.59)  592.21 (1.87)  591.21 (1.84)  589.87 (1.91) 
DNK  1  3.26 (0.61)  3.23 (0.63)  30.79 (0.93)  28.12 (1.67)  25.60 (13.11)  16.84 (0.85)  14.46 (0.79)  22.15 (1.10)  22.11 (1.10)  21.72 (1.07) 
FRA  1  2.79 (1.10)  1.92 (0.66)  4.67 (0.41)  3.05 (0.11)  52.65 (25.83)  91.69 (1.34)  156.29 (1.96)  116.32 (1.27)  116.54 (1.25)  115.29 (1.28) 
DEU  1  16.99 (8.04)  6.71 (1.23)  229.65 (9.13)  210.59 (14.99)  341.03 (157.80)  739.29 (11.81)  929.03 (4.85)  817.50 (4.60)  818.13 (4.55)  812.60 (4.57) 
GRC  1  3.80 (2.21)  3.55 (1.79)  5.30 (0.90)  5.75 (0.68)  11.78 (2.36)  26.90 (1.89)  23.05 (0.53)  47.07 (1.92)  45.50 (1.82)  45.72 (2.00) 
IRL  1  3.05 (0.34)  4.35 (1.25)  135.57 (5.64)  12.34 (0.58)  23.40 (17.50)  3.84 (0.22)  26.60 (0.59)  6.38 (0.13)  6.31 (0.14)  6.16 (0.13) 
ITA  1  13.00 (4.15)  14.05 (4.81)  35.29 (1.83)  39.27 (2.52)  87.34 (24.05)  226.95 (11.14)  343.10 (10.04)  244.25 (8.50)  244.84 (8.58)  242.60 (8.46) 
JPN  1  10.55 (4.67)  12.32 (4.95)  8.10 (1.05)  8.38 (1.07)  18.81 (4.59)  95.58 (7.89)  71.02 (5.08)  135.24 (13.57)  134.89 (13.50)  134.16 (13.43) 
NLD  1  3.75 (0.80)  3.87 (0.79)  0.99 (0.06)  0.99 (0.05)  9.45 (1.43)  28.35 (1.62)  29.53 (1.58)  33.28 (1.78)  33.23 (1.77)  33.14 (1.77) 
NOR  1  2.70 (0.51)  2.82 (0.73)  1.86 (0.29)  1.63 (0.11)  24.25 (12.50)  23.36 (0.88)  31.37 (1.17)  27.86 (0.94)  27.86 (0.93)  27.52 (0.91) 
ESP  1  5.18 (1.05)  6.09 (1.53)  5.17 (1.14)  4.29 (0.72)  14.85 (4.20)  33.16 (6.99)  152.59 (6.19)  53.53 (2.47)  52.56 (2.42)  52.06 (2.40) 
SWE  1  6.44 (2.66)  5.47 (2.63)  2.48 (0.23)  2.02 (0.21)  2.18 (0.25)  15.53 (2.59)  2706.85 (17.91)  118.46 (1.64)  118.23 (1.64)  118.27 (1.64) 
CHE  1  3.51 (0.46)  2.90 (0.37)  43.59 (1.77)  7.48 (0.49)  38.32 (9.03)  8.43 (0.24)  29.71 (0.53)  9.72 (0.29)  9.71 (0.29)  9.79 (0.28) 
TUR  1  1.65 (0.47)  1.06 (0.15)  1.22 (0.18)  0.91 (0.09)  2.19 (0.34)  64.26 (5.71)  142.84 (2.04)  159.79 (2.63)  157.89 (2.63)  157.13 (2.69) 
GBR  1  5.95 (1.86)  2.66 (0.57)  15.92 (1.02)  10.05 (1.47)  7.57 (5.10)  50.04 (1.75)  68.70 (1.25)  70.98 (1.01)  70.87 (0.99)  69.72 (1.01) 
USA  1  4.98 (1.96)  1.60 (0.42)  21.53 (3.30)  12.28 (2.52)  2.06 (0.47)  308.69 (5.20)  244.90 (1.82)  462.51 (2.14)  464.75 (2.08)  465.88 (2.16) 
#Best    2  10  2  4  0  0  0  0  0  0 
7 Conclusion
In this paper, we proposed a novel fewshot supervised DA method for regression problems based on the assumption of shared generative mechanism. Through theoretical and experimental analysis, we demonstrated the effectiveness of the proposed approach. By considering the latent common structure behind the domain distributions, the proposed method successfully induces positive transfer even when a naive usage of the source domain data can suffer from negative transfer. Our future work includes making an experimental comparison with extensively more datasets and methods as well as an extension to the case where the underlying mechanism are not exactly identical but similar.
Acknowledgments
TT was supported by Masason Foundation. TT would like to thank Yuko Kuroki and Taira Tsuchiya for their feedback on the manuscript. MS was supported by JST CREST Grant Number JPMJCR18A2.
This is the Supplementary Material for “Fewshot Domain Adaptation by Causal Mechanism Transfer.” Table 3 summarizes the abbreviations and the symbols used in the paper.
Abbreviation / Symbol  Meaning 

DA  Domain adaptation 
TA  Transfer assumption 
SEM  Structural equation model 
GCM  Graphical causal model 
SCM  Structural causal model 
IC  Independent component 
ICA  Independent component analysis 
GCL  Generalized contrastive learning 
i.i.d.  Independent and identically distributed 
where  
The Sobolev norm  
The predictor random vector  
The predicted random variable  
The joint random variable  
The independent component vector  
The space of  
The space of  
The space of  
Predictor hypothesis class  
Loss function  
Target domain risk  
Minimizer of target domain risk  
The set of independent distributions  
Ground truth mixing function  
The target joint distribution  
The joint distribution of source domain  
The target independent component (IC) distribution  
The IC distribution of source domain  
The dimension of  
The number of source domains  
The size of the target labeled sample  
The size of the labeled sample from source domain  
Target labeled data set  
Source labeled data set of source domain  
The ordinary empirical risk estimator  
The proposed risk estimator (Eq. (2))  
The estimator of  
The penultimate layer functions composed with during GCL 
Appendix A Preliminary: Nonlinear ICA
Here, we use the same notation as the main text. The recently developed nonlinear ICA provides an algorithm to estimate the mixing function . For the case of nonlinear , the impossibility of identification (i.e., consistent estimation) of in the onesample i.i.d. case had been established more than two decades ago (Hyvärinen and Pajunen, 1999). However, recently, various conditions have been proposed under which can be identified with the help of auxiliary information (Hyvärinen and Morioka, 2016, 2017; Hyvärinen et al., 2019; Khemakhem et al., 2019).
The identification condition that is directly relevant to this paper is that of the generalized contrastive learning (GCL) proposed in Hyvärinen et al. (2019). Hyvärinen et al. (2019) assumes that an auxiliary variable from some measurable set is obtained for each data point as and that the ICs are conditionally independent given :
Under such conditions, GCL estimates by training a classification function
(4) 
parametrized by and with the logistic loss for classifying
where . The key condition for the identification of is the following.
Assumption 1 (Assumption of variability; Hyvärinen et al., 2019, Theorem 1).
For any , there exist distinct points in , denoted by , such that the set of dimensional vectors are linearly independent, where
Under Assumption 1 and some regularity conditions, Theorem 1 of Hyvärinen et al. (2019) states that the transformation in Eq. (4) trained by GCL is a consistent estimator of upto additional dimensionwise invertible transformations. Note that the assumption is intrinsically difficult to confirm based on data due to the unsupervised nature of the problem setting. In this paper, we use the source domain index as the auxiliary variable and employ GCL for domain adaptation. The present version of Assumption 1 requires that we have at least distinct source domains. Although this condition can be restrictive in highdimensional data, we conjecture that there is a possibility for this assumption to be made less stringent in the future because the identification condition is only known to be a sufficient condition, not a necessary condition. However, pursuing a refinement of the identification condition is out of the scope of this paper. Among the various methods for nonlinear ICA, we chose to use GCL (Hyvärinen et al., 2019) because it can operate under a nonparametric assumption on the IC distributions whereas other nonlinear ICA methods (Hyvärinen and Morioka, 2016, 2017; Khemakhem et al., 2019) may require parametric assumptions.
Appendix B Experiment Details
Here, we describe more implementation details of the experiment. We plan to publish the experiment code and the dataset if the manuscript is accepted for publication. A URL to the experiment code will appear here.
b.1 Dataset details
Gasoline consumption data.
The data was downloaded from http://bcs.wiley.com/hebcs/Books?action=resource&bcsId=4338&itemId=1118672321&resourceId=13452.
b.2 Model details: Invertible neural networks
Here, we describe the details of the Glow architecture (Kingma and Dhariwal, 2018) used in our experiments. Glow consists of three types of layers which are invertible by design, namely affine coupling layers, convolution layers, and actnorm layers. In our implementation, we do not include actnorm layers, and each layer of our Glow architecture consist of a convolution layer followed by an affine coupling layer.
Affine coupling layers.
The coefficients and for affine coupling layers in the notation of Kingma and Dhariwal (2018) are parametrized by two onehiddenlayer neural networks whose number of hidden units is the same and the first layer parameter is shared. The activation functions of the first layer, the second layer of , and the second layer of are the rectified linear unit (ReLU) activation (LeCun et al., 2015), the hyperbolic tangent function, and the linear activation function, respectively. A standard practice of affine coupling layers is to compose the coefficient with an exponential function so as to simplify the computation of the logdeterminant of the Jacobian (Kingma and Dhariwal, 2018). In our implementation, since we do not require the computation of the logdeterminant, we omit this device and instead compose . The addition of shifts the parameter space so that corresponds to the the identity map, where denotes the constant zero function. The split of the affine coupling layers is fixed at .
convolution layers.
We initialize the parameters of the neural networks by where is the number of parameters of each layer and is the normal distribution.
b.3 Model details: Penultimate layer networks
We initialize the parameter for each layer of by , where is the number of input features and is the uniform distribution.
b.4 Training details
During the training of GCL, we fix the batch size at 32.
b.5 Compared methods details
TrAdaBoost.
As suggested in Pardoe and Stone (2010), we use the linear loss function and set the maximum number of internal boosting iterations at .
Gdm.
We fix the number of sampling required for approximating the maximization in the generalization discrepancy at . This method presumes using hypothesis classes in a reproducing kernel Hilbert space (RKHS).
Copula.
For this model, the probabilistic model of nonparametric Rvine copula of depth is used following Lopezpaz et al. (2012). Kernel density estimators with RBF kernel are used for estimating the marginal distributions and the copulas. The bandwidths of the RBF kernels are determined using the ruleofthumb implemented as “normalreference” in the np package of R language (Hayfield and Racine, 2008). The predictions are made by numerically aggregating the estimated conditional distribution over the interval where denotes the square root of the unbiased variance of . The aggregation is performed by discretizing the interval into a grid of points. The level of the twosample test is fixed at for all combination of the twosample tests following the experiment code of Lopezpaz et al. (2012). This method is a singlesource domain adaptation method and we pool all source domain data for adaptation.
Appendix C Details and Proofs of Theorem 2
Here, we detail the assumptions, the statement, and the proof of Theorem 2.
c.1 Notation
To make the proof selfcontained, we first recall some general and problemspecific notation. In the notation here, we omit the domain identifiers from the distributions and the sample size, such as Tar or Src, because only the target domain data or their distributions appear in the proofs. The theorem holds regardless of how is estimated as long as is independent of the target domain data. In the proof, we extend the maximal discrepancy bound of Ustatistics previously proved for the case of degree in Rejchel (2012), to allow higher degrees.
General mathematical notation.
We denote the set of natural numbers (resp. real numbers) by (resp. ). For any , we define . We use to denote the number of combinations of elements. For a finite set , the notation denotes the operator to take an average over , i.e., . For a dimensional function , we denote its th dimension () by suffixing . For a vector , we denote its th element by . We denote the Jacobian determinant of a differentiable function at by . We denote the identity matrix by regardless of the size of the matrices when there is no ambiguity. For finite dimensional vectors, we denote the norm by and the norm by . For square matrices, we denote the operator norm by and the operator norm by . We use to denote the Sobolev space (on ) of order and define its associated norm by where is a multiindex and denotes the partial derivative (Adams and Fournier, 2003, Paragraph 3.1). We let be the degree symmetric group, be the set of grouping of indices in , and be the set of all size combinations (without replacement) of indices in .
Distributions and expectations.
We denote by the set of all factorized distributions on with absolutely continuous marginals. For a measure , we denote its product measure by (repeated times). We assume that all measures appearing in this proof are absolutely continuous with respect to the Lebesgue measure. The pushforward of a distribution by a function is denoted by . The expectation of a function with respect to measure is denoted by (if it exists) by abuse of notation. We also abuse the notation to use as the shorthand for where .
c.2 Problem setup
We denote the target domain distribution by .
We fix a hypothesis class , and our goal is to find a such that the risk functional
is small, where is a loss function. We denote by a minimizer of (assuming it exists). To this end, we are given the training data . Throughout, we assume . To complement the smallness of , we assume the existence of a generative mechanism. Concretely, we assume that there exists a diffeomorphism such that satisfies . With this transform, the original risk functional is also expressed as
As an estimator of , we are given another diffeomorphism such that . With this , the proposed method converts the dataset by . We can regard , where . We use (resp. ) to denote the probability measure corresponding to the density (resp. ). This conversion results in the relation:
As a candidate hypothesis , the proposed method selects a minimizer of the proposed risk estimator defined as
(5) 
In the proof, we evaluate its concentration around the expectation . We use to denote the expectation with respect to . Let denote a hypothesis which minimizes (assuming it exists).
In what follows, for notational simplicity, we define the variate symmetric function as
where indicates an averaging operation over all permutations (without replacement) of . We use to denote the sample average operator with respect to or , depending on the context.
c.3 Assumptions
Assumption 2 (The underlying density function is bounded and Lipschitz continuous).
Assume
Assumption 3 ( is Lipschitz continuous and Hölder continuous).
We assume where is the Hölder space (Adams and Fournier, 2003, Paragraph 1.29) and
Assumption 4 (Bounded derivatives of and ).
Assume that
where denotes the maximum absolute value of the elements of a matrix.
Assumption 5 (Loss function is bounded and uniformly Lipschitz continuous in ).
The considered loss function takes values in a bounded interval:
where . Also assume
Assumption 6 (Estimated feature extractor).
Assume is independent of and that for all .
Although and are assumed to be diffeomorphisms in the classical sense (implying that they are strongly differentiable), we introduce the Sobolev space because we want to measure their difference and their difference of derivatives in terms of integration.
Assumption 7 (Entropic condition: Euclidean class (Sherman, 1994)).
The function class is Euclidean for the envelope and constants and (Sherman, 1994), i.e., if is a measure for which , then
where is the pseudo metric defined by
for , and denotes the packing number of with respect to the pseudometric and radius . Without loss of generality, we take the envelope such that .
Assumption 8.
The hypothesis class is expressive enough so that the model approximation error does not expand due to , i.e.,
The following complexity measure of , which is a version of Rademacher complexity for our problem setting, is used to state the theorem.
Definition 1 (Effective Rademacher complexity).
Define
where are independent uniform sign variables and are independent of all other random variables.
We provide the definition of the ordinary Rademacher complexity in Section C.8 and make a comparison of the two complexity measures in terms of how they depend on the input dimensionality.
c.4 Theorem statement
Our goal is to prove the following theorem. This is a detailed version of the theorem appearing in the main body of the paper.
Theorem 3 (Excess risk bound).
Then for arbitrary , we have with probability at least ,
where
and are constants determined in Lemma 11.
Proof of Theorem 3.
By adding and subtracting terms, we have