Title: LEARNING FACTORIZED MULTIMODAL REPRESENTATIONS

Abstract: Learning multimodal representations is a fundamentally complex research problem due to the presence of multiple heterogeneous sources of information. Although the presence of multiple modalities provides additional valuable information, there are two key challenges to address when learning from multimodal data: 1) models must learn the complex intra-modal and cross-modal interactions for prediction and 2) models must be robust to unexpected missing or noisy modalities during testing. In this paper, we propose to optimize for a joint generative-discriminative objective across multimodal data and labels. We introduce a model that factorizes representations into two sets of independent factors: multimodal discriminative and modality-specific generative factors. Multimodal discriminative factors are shared across all modalities and contain joint multimodal features required for discriminative tasks such as sentiment prediction. Modality-specific generative factors are unique for each modality and contain the information required for generating data. Experimental results show that our model is able to learn meaningful multimodal representations that achieve state-of-the-art or competitive performance on six multimodal datasets. Our model demonstrates flexible generative capabilities by conditioning on independent factors and can reconstruct missing modalities without significantly impacting performance. Lastly, we interpret our factorized representations to understand the interactions that influence multimodal learning.

1. INTRODUCTION:

Multimodal machine learning involves learning from data across multiple modalities (Baltrušaitis et al., 2017) . It is a challenging yet crucial research area with real-world applications in robotics (Liu et al., 2017) , dialogue systems (Pittermann et al., 2010) , intelligent tutoring systems (Petrovica et al., 2017) , and healthcare diagnosis (Frantzidis et al., 2010) . At the heart of many multimodal modeling tasks lies the challenge of learning rich representations from multiple modalities. For example, analyzing multimedia content requires learning multimodal representations across the language, visual, and acoustic modalities (Cho et al., 2015) . Although the presence of multiple modalities provides additional valuable information, there are two key challenges to address when learning from multimodal data: 1) models must learn the complex intra-modal and cross-modal interactions for prediction (Zadeh et al., 2017) , and 2) trained models must be robust to unexpected missing or noisy modalities during testing (Ngiam et al., 2011) .

In this paper, we propose to optimize for a joint generative-discriminative objective across multimodal data and labels. The discriminative objective ensures that the representations learned are rich in intra-modal and cross-modal features useful towards predicting the label, while the generative objective allows the model to infer missing modalities at test time and deal with the presence of noisy modalities. To this end, we introduce the Multimodal Factorization Model (MFM in Figure 1 ) that factorizes multimodal representations into multimodal discriminative factors and modality-specific generative factors. Multimodal discriminative factors are shared across all modalities and contain joint multimodal features required for discriminative tasks. Modality-specific generative factors are unique for each modality and contain the information required for generating each modality. We believe that factorizing multimodal representations into different explanatory factors can help each factor focus on learning from a subset of the joint information across multimodal data and labels. This method is in contrast to jointly learning a single factor that summarizes all generative and discriminative information (Srivastava & Salakhutdinov, 2012) . To sum up, MFM defines a joint distribution over multimodal data, and by the conditional independence assumptions in the assumed graphical model, both generative and discriminative aspects are taken into account. Our model design further provides interpretability of the factorized representations.

Through an extensive set of experiments, we show that MFM learns improved multimodal representations with these characteristics: 1) The multimodal discriminative factors achieve state-of-the-art or competitive performance on six multimodal time series datasets. We also demonstrate that MFM can generalize by integrating it with other existing multimodal discriminative models. 2) MFM allows flexible generation concerning multimodal discriminative factors (labels) and modality-specific generative factors (styles). We further show that we can perform reconstruction of missing modalities from observed modalities without significantly impacting discriminative performance. Finally, we interpret our learned representations using information-based and gradient-based methods, allowing us to understand the contributions of individual factors towards multimodal prediction and generation.

2. MULTIMODAL FACTORIZATION MODEL:

Multimodal Factorization Model (MFM) is a latent variable model (Figure 1 (a)) with conditional independence assumptions over multimodal discriminative factors and modality-specific generative factors. According to these assumptions, we propose a factorization over the joint distribution of multimodal data (Section 2.1). Since exact posterior inference on this factorized distribution can be intractable, we propose an approximate inference algorithm based on minimizing a joint-distribution Wasserstein distance over multimodal data (Section 2.2). Finally, we derive the MFM objective by approximating the joint-distribution Wasserstein distance via a generalized mean-field assumption.

Notation: We define X 1∶M as the multimodal data from M modalities and Y as the labels, with joint distribution P X 1∶M ,Y = P (X 1∶M , Y). Let X1∶M denote the generated multimodal data and Ŷ denote the generated labels, with joint distribution P X1∶M , Ŷ = P ( X1∶M , Ŷ).

2.1. FACTORIZED MULTIMODAL REPRESENTATIONS:

To factorize multimodal representations into multimodal discriminative factors and modality-specific generative factors, MFM assumes a Bayesian network structure as shown in Figure 1 (a) . In this graphical model, factors F y and F a{1∶M } are generated from mutually independent latent variables Z = [Z y , Z a{1∶M } ] with prior P Z . In particular, Z y generates the multimodal discriminative factor F y and Z a{1∶M } generate modality-specific generative factors F a{1∶M } . By construction, F y contributes to the generation of Ŷ while {F y , F ai } both contribute to the generation of Xi . As a result, the joint distribution P ( X1∶M , Ŷ) can be factorized as follows:

EQUATION

with dF = dF y ∏ M i=1 dF ai and dZ = dZ y ∏ M i=1 dZ ai . Exact posterior inference in Equation 1 may be analytically intractable due to the integration over Z. We therefore resort to using an approximate inference distribution Q(Z X 1∶M , Y) as detailed in the following subsection. As a result, MFM can be viewed as an autoencoding structure that consists of encoder (inference) and decoder (generative) modules (Figure 1 (c)). The encoder module for Q(⋅ ⋅) allows us to easily sample Z from an approximate posterior. The decoder modules are parametrized according to the factorization of P ( X1∶M , Ŷ Z) as given by Equation 1and Figure 1 (a).

2.2. MINIMIZING JOINT-DISTRIBUTION WASSERSTEIN DISTANCE OVER MULTIMODAL DATA:

Two common choices for approximate inference in autoencoding structures are Variational Autoencoders (VAEs) (Kingma & Welling, 2013) and Wasserstein Autoencoders (WAEs) (Zhao et al., 2017; Tolstikhin et al., 2017) . The former optimizes the evidence lower bound objective (ELBO), and the latter derives an approximation for the primal form of the Wasserstein distance. We consider the latter since it simultaneously results in better latent factor disentanglement (Zhao et al., 2017; Rubenstein et al., 2018) and better sample generation quality than its counterparts (Chen et al., 2016; Higgins et al., 2016; Kingma & Welling, 2013) . However, WAEs are designed for unimodal data and do not consider factorized distributions over latent variables that generate multimodal data. Therefore, we propose a variant for handling factorized joint distributions over multimodal data.

As suggested by Kingma & Welling (2013) , we adopt the design of nonlinear mappings (i.e. neural network architectures) in the encoder and decoder (Figure 1 Rubenstein et al., 2018; Tolstikhin et al., 2017) . For the decoder, we define the generation process from latent variables as

(c)). For the encoder Q(Z X 1∶M , Y), we learn a deterministic mapping Q enc ∶ X 1∶M , Y → Z (

G y ∶ Z y → F y , G a{1∶M } ∶ Z a{1∶M } → F a{1∶M } , D ∶ F y → Ŷ, and F 1∶M ∶ F y , F a{1∶M } → X1∶M

, where G y , G a{1∶M } , D and F 1∶M are deterministic functions parametrized by neural networks.

Let W c (P X 1∶M ,Y , P X1∶M , Ŷ ) denote the joint-distribution Wasserstein distance over multimodal data under cost function c X i and c Y . We choose the squared cost c(a, b) = a -b 2 2 , allowing us to minimize the 2-Wasserstein distance. The cost function can be defined not only on static data but also on time series data such as text, audio and videos. For example, given time series data

X = [X 1 , X 2 , ⋯, X T ] and X = [ X1 , X2 , ⋯, XT ], we define c(X, X) = ∑ T t=1 X t -Xt 2

2 . With conditional independence assumptions in Equation 1, we express W c (P X 1∶M ,Y , P X1∶M , Ŷ ) as:

Proposition 1. For any functions G y ∶ Z y → F y , G a{1∶M } ∶ Z a{1∶M } → F a{1∶M } , D ∶ F y → Ŷ, and F 1∶M ∶ F a{1∶M } , F y → X1∶M , we have W c (P X 1∶M ,Y , P X1∶M , Ŷ ) = inf Q Z =P Z E P X 1∶M ,Y E Q(Z X 1∶M ,Y) ⎡ ⎢ ⎢ ⎢ ⎢ ⎣ M i=1 c X i Xi, Fi Gai(Zai), Gy(Zy) + c Y Y, D Gy(Zy) ⎤ ⎥ ⎥ ⎥ ⎥ ⎦ , (2)

where P Z is the prior over Z = [Z y , Z a{1,M } ] and Q Z is the aggregated posterior of the proposed approximate inference distribution Q(Z X 1∶M , Y).

Proof: The proof is adapted from Tolstikhin et al. (Tolstikhin et al., 2017) . The two differences are: (1) we show that P ( X1∶M , Ŷ Z = z) are Dirac for all z ∈ Z, and (2) we use the fact that c

((X 1∶M , Y), ( X1∶M , Ŷ)) = ∑ M i=1 c X i (X i , Xi ) + c Y (Y, Ŷ).

Please refer to the supplementary material for proof details.

∎ The constraint on Q Z = P Z in Proposition 1 is hard to satisfy. To obtain a numerical solution, we first relax the constraint by performing a generalized mean field assumption on Q according to the conditional independence as shown in the inference network of Figure 1 (b) :

Q(Z X 1∶M , Y) ∶= Q(Z X 1∶M ) ∶= Q(Zy X 1∶M ) M i=1 Q(Za i Xi).

(3)

The intuition here is based on our design that Z y generates the multimodal discriminative factor F y and Z a{1∶M } generate modality-specific generative factors F a{1∶M } . Therefore, the inference for Z y should depend on all modalities X 1∶M and the inference for Z ai should depend only on the specific modality X i . Following this assumption, we define Q as a nonparametric set of all encoders that fulfill the factorization in Equation 3. A penalty term is added into our objective to find the Q(Z ⋅) ∈ Q that is the closest to prior P Z , thereby approximately enforcing the constraint where λ is a hyper-parameter and MMD is the Maximum Mean Discrepancy (Gretton et al., 2012) as a divergence measure between Q Z and P Z . The prior P Z is chosen as a centered isotropic Gaussian N (0, I), so that it implicitly enforces independence between the latent variables Z = [Z y , Z a{1,M } ] (Higgins et al., 2016; Kingma & Welling, 2013; Rubenstein et al., 2018) .

EQUATION

Equation 4 represents our hybrid generative-discriminative optimization objective over multimodal data: the first loss term

∑ M i=1 c Xi (X i , F (G ai (Z ai ), G y (Z y )))

is the generative objective based on reconstruction of multimodal data and the second term c Y (Y, D(G y (Z y ))) is the discriminative objective. In practice we compute the expectations in Equation 4 using empirical estimates over the training data. The neural architecture of MFM is illustrated in Figure 1(c) .

2.3. SURROGATE INFERENCE FOR MISSING MODALITIES:

A key challenge in multimodal learning involves dealing with missing modalities. A good multimodal model should be able to infer the missing modality conditioned on the observed modalities and perform predictions based only on the observed modalities. To achieve this objective, the inference process of MFM can be easily adapted using a surrogate inference network to reconstruct the missing modality given the observed modalities. Formally, let Φ denote the surrogate inference network. The generation of missing modality X1 given the observed modalities X 2∶M can be formulated as

Φ * = argmin Φ E P X 2∶M , X1 -log P Φ ( X1 X 2∶M ) with P Φ ( X1 X 2∶M ) ∶= P ( X1 Z a1 , Z y )Q Φ (Z a1 X 2∶M )Q Φ (Z y X 2∶M )dZ a1 dZ y .

(5)

Similar to Section 2.2, we use deterministic mappings in Q Φ (⋅ ⋅) and Q Φ (Z y ⋅) is also used for prediction P Φ ( Ŷ X 2∶M ) ∶= ∫ P ( Ŷ Z y )Q Φ (Z y X 2∶M )dZ y . Equation 5 suggests that in the presence of missing modalities, we only need to infer the latent codes rather than the entire modality.

2.4. ENCODER AND DECODER DESIGN:

We now discuss the implementation choices for the MFM neural architecture in Figure 1(c ). The encoder Q(Z y X 1∶M ) can be parametrized by any model that performs multimodal fusion (Morency et al., 2011; Zadeh et al., 2017) . For multimodal image datasets, we adopt Convolutional Neural Networks (CNNs) and Fully-Connected Neural Networks (FCNNs) with late fusion (Nojavanasghari et al., 2016) as our encoder Q(Z y X 1∶M ). The remaining functions in MFM are also parametrized by CNNs and FCNNs. For multimodal time series datasets, we choose the Memory Fusion Network (MFN) (Zadeh et al., 2018a) as our multimodal encoder Q(Z y X 1∶M ). We use Long Short-term Memory (LSTM) networks (Hochreiter & Schmidhuber, 1997) for functions Q(Z a{1∶M } X 1∶M ), decoder LSTM networks (Cho et al., 2014) for functions F 1∶M , and FCNNs for functions G y , G a{1∶M } and D. Details are provided in the appendix and the code is available at <anonymous>.

3. EXPERIMENTS:

In order to show that MFM learns multimodal representations that are discriminative, generative and interpretable, we design the following experiments. We begin with a multimodal synthetic image dataset that allows us to examine whether MFM displays discriminative and generative capabilities from factorized latent variables. Utilizing image datasets allows us to clearly visualize the generative capabilities of MFM. We then transition to six more challenging real-world multimodal video datasets 

∆ SOT A ↑ 1.5 ↑ 0.7 ↑ 0.8 ↓ 0.014 ↑ 0.030 ↑ 7.5 ↑ 6.1 ↑ 1.6 ↑ 0.8 ↑ 1.0 ↑ 0.

∆ SOT A ↑ 0.1 ↑ 0.5 ↑ 2.6 ↑ 3.3 ↑ 0.5 ↑ 0.7 ↑ 0.1 - ↑ 0.2 - ↑ 0.3 -

to 1) rigorously evaluate the discriminative capabilities of MFM in comparison with existing baselines, 2) analyze the importance of each design component through ablation studies, 3) assess the robustness of MFM's modality reconstruction and prediction capabilities to missing modalities, and 4) interpret the learned representations using information-based and gradient-based methods to understand the contributions of individual factors towards multimodal prediction and generation.

3.1. MULTIMODAL SYNTHETIC IMAGE DATASET:

In this section, we study MFM on a synthetic image dataset that considers SVHN (Netzer et al., 2011) and MNIST (Lecun et al., 1998) as the two modalities. SVHN and MNIST are images with different styles but the same labels (digits 0 ∼ 9). We randomly pair 100, 000 SVHN and MNIST images that have the same label, creating a multimodal dataset which we call SVHN+MNIST. 80, 000 pairs are used for training and the rest for testing. To justify that MFM is able to learn improved multimodal representations, we show both classification and generation results on SVHN+MNIST in Figure 2 .

Prediction: We perform experiments on both unimodal and multimodal classification tasks. UM denotes a unimodal baseline that performs prediction given only one modality as input and MM denotes a multimodal discriminative baseline that performs prediction given both images (Nojavanasghari et al., 2016) . We compare the results for UM(SVHN), UM(MNIST), MM and MFM on SVHN+MNIST in Figure 2 (b). We achieve better classification performance from unimodal to multimodal which is not surprising since more information is given. More importantly, MFM outperforms MM, which suggests that MFM learns improved factorized representations for discriminative tasks.

Generation: We generate images using the MFM generative network (Figure 2(a) ). We fix one variable out of Z = [Z a1 , Z a2 , and Z y ] and randomly sample the other two variables from prior P Z .

From Figure 2 (c), we observe that MFM shows flexible generation of SVHN and MNIST images based on labels and styles. This suggests that MFM is able to factorize multimodal representations into multimodal discriminative factors (labels) and modality-specific generative factors (styles).

3.2. MULTIMODAL TIME SERIES DATASETS:

In this section, we transition to more challenging multimodal time series datasets. All the datasets consist of monologue videos. Features are extracted from the language (GloVe word embeddings (Pen- Prediction: We first compare the performance of MFM with existing multimodal prediction methods.

M E M C M D M A M B MFM

For a detailed description of the baselines, please refer to the appendix. From Table 1 , we first observe that the best performing baseline results are achieved by different models across different datasets (most notably MFN, MARN, and TFN). On the other hand, MFM consistently achieves state-of-the-art or competitive results for all six multimodal datasets. We believe that the multimodal discriminative factor F y in MFM has successfully learned more meaningful representations by distilling discriminative features. This highlights the benefit of learning factorized multimodal representations towards discriminative tasks. Furthermore, MFM is model-agnostic and can be applied to other multimodal encoders Q(Z y X 1∶M ). We perform experiments to show consistent improvements in discriminative performance for several choices of the encoder: EF-LSTM (Morency et al., 2011) and TFN (Zadeh et al., 2017) . For Acc_2 on CMU-MOSI, our factorization framework improves the performance of EF-LSTM from 74.3 to 75.2 and TFN from 74.6 to 75.5.

Ablation Study: In Figure 3 , we present the models M {A,B,C,D,E} used for ablation studies. These models are designed to analyze the effects of using a multimodal discriminative factor, a hybrid generative-discriminative objective, factorized generative-discriminative factors and modality-specific generative factors towards both modality reconstruction and label prediction. The simplest variant is M A which represents a purely discriminative model without a joint multimodal discriminative factor (i.e. early fusion (Morency et al., 2011) ). M B models a joint multimodal discriminative factor which incorporates more general multimodal fusion encoders (Zadeh et al., 2018a) . M C extends From the table in Figure 3 , we observe the following general trends. For sentiment prediction, using 1) a multimodal discriminative factor outperforms modality-specific discriminative factors Missing Modalities: We now evaluate the performance of MFM in the presence of missing modalities using the surrogate inference model as described in Subsection 2.3. We compare with two baselines: 1) a purely generative Seq2Seq model (Cho et al., 2014 ) Φ G from observed modalities to missing modalities by optimizing E P X 1∶M (-log P Φ D (X 1 X 2∶M )), and 2) a purely discriminative model Φ D from observed modalities to the label by optimizing E P X 2∶M ,Y (-log P Φ G (Y X 2∶M )). Both models are modified from MFM by using only the two observed modalities as input and not explicitly accounting for missing modalities. We compare the reconstruction error of each modality (language, visual and acoustic) as well as the performance on sentiment prediction.

(M D > M C , M B > M A ),

Table 2 shows that MFM with missing modalities outperforms the generative (Φ G ) or discriminative baselines (Φ D ) in terms of modality reconstruction and sentiment prediction. Additionally, MFM with missing modalities performs close to MFM with all modalities observed. This fact indicates that MFM can learn representations that are relatively robust to missing modalities. In addition, discriminative performance is most affected when the language modality is missing, which is consistent with prior work which indicates that language is most informative in human multimodal language (Zadeh et al., 2017) . On the other hand, sentiment prediction is more robust to missing acoustic and visual features. Finally, we observe that reconstructing the low-level acoustic and visual features is easier as compared to the high-dimensional language features that contain high-level semantic meaning.

Interpretation of Multimodal Representations::

We devise two methods to study how individual factors in MFM influence the dynamics of multimodal prediction and generation. These interpretation methods represent both overall trends and fine-grained analysis that could be useful towards deeper understandings of multimodal representation learning. For more details, please refer to the appendix.

Firstly, an information-based interpretation method is chosen to summarize the contribution of each modality towards the multimodal representations. Since F y is a common cause of X1∶M , we can compare MI(F y , X1 ), ⋯, MI(F y , XM ), where MI(⋅, ⋅) denotes the mutual information measure between F y and generated modality Xi . Higher MI(F y , Xi ) indicates greater contribution from F y to Xi . Figure 4 reports the ratios r i = MI(F y , Xi ) MI(F ai , Xi ) which measure a normalized version of the mutual information between F ai and Xi . We observe that on CMU-MOSI, the language modality is most informative towards sentiment prediction, followed by the acoustic modality. We believe that this result represents a prior over the expression of sentiment in human multimodal language and is closely related to the connections between language and speech (Kuhl, 2000) .

Secondly, a gradient-based interpretation method to used analyze the contribution of each modality for every time step in multimodal time series data. We measure the gradient of the generated modality with respect to the target factors (e.g., F y ). Let {x 1 , x 2 , ⋯, x M } denote multimodal time series data where x i represents modality i, and xi = [x 1 i , ⋯, xt i , ⋯, xT i ] denote generated modality i across time steps t ∈ [1, T ]. The gradient ∇ fy (x i ) measures the extent to which changes in factor f y ∼ P (F y X 1∶M = x 1∶M ) influences the generation of sequence xi . Figure 4 plots ∇ fy (x i ) for a video in CMU-MOSI. We observe that multimodal communicative behaviors that are indicative of speaker sentiment such as positive words (e.g. "very profound and deep") and informative acoustic features (e.g. hesitant and emphasized tone of voice) indeed correspond to increases in ∇ fy (x i ).

4. RELATED WORK:

The two main pillars of research in multimodal representation learning have considered the discriminative and generative objectives individually. Discriminative representation learning (Liang et al., 2018; Chen et al., 2017; Chaplot et al., 2017; Frome et al., 2013; Socher et al., 2013; Tsai et al., 2017) models the conditional distribution P (Y X 1∶M ). Since these approaches are not concerned with modeling P (X 1∶M ) explicitly, they use parameters more efficiently to model P (Y X 1∶M ). For instance, recent works learn visual representations that are maximally dependent with linguistic attributes for improving one-shot image recognition (Tsai & Salakhutdinov, 2017) or introduce tensor product mechanisms to model interactions between the language, visual and acoustic modalities (Liu et al., 2018; Zadeh et al., 2017) . On the other hand, generative representation learning captures the interactions between modalities by modeling the joint distribution P (X 1 , ⋯, X M ) using either undirected graphical models (Srivastava & Salakhutdinov, 2012) , directed graphical models (Suzuki et al., 2016b) , or neural networks (Sohn et al., 2014) . Some generative approaches compress multimodal data into lower-dimensional feature vectors which can be used for discriminative tasks (Pham et al., 2018; Ngiam et al., 2011) . To unify the advantages of both approaches, MFM factorizes multimodal representations into generative and discriminative components and optimizes for a joint objective.

Factorized representation learning resembles learning disentangled data representations which have been shown to improve the performance on many tasks (Kulkarni et al., 2015; Lake et al., 2017; Higgins et al., 2016; Bengio et al., 2013) . Several methods involve specifying a fixed set of latent attributes that individually control particular variations of data and performing supervised training (Cheung et al., 2014; Karaletsos et al., 2015; Yang et al., 2015; Reed et al., 2014; Zhu et al., 2014) , assuming an isotropic Gaussian prior over latent variables to learn disentangled generative representations (Kingma & Welling, 2013; Rubenstein et al., 2018) and learning latent variables in charge of specific variations in the data by maximizing the mutual information between a subset of latent variables and the data (Chen et al., 2016) . However, these methods study factorization of a single modality. MFM factorizes multimodal representations and demonstrates the importance of modality-specific and multimodal factors towards generation and prediction. A concurrent and parallel work that factorizes latent factors in multimodal data was proposed by Hsu & Glass (2018) . They differ from us in the graphical model design, discriminative objective, prior matching criterion, and scale of experiments. We provide a detailed comparison with their model in the appendix.

5. CONCLUSION:

In this paper, we proposed the Multimodal Factorization Model (MFM) for multimodal representation learning. MFM factorizes the multimodal representations into two sets of independent factors: multimodal discriminative factors and modality-specific generative factors. The multimodal discriminative factor achieves state-of-the-art or competitive results on six multimodal datasets. The modalityspecific generative factors allow us to generate data based on factorized variables, account for missing modalities, and have a deeper understanding of the interactions involved in multimodal learning. Our future work will explore extensions of MFM for video generation, semi-supervised learning, and unsupervised learning. We believe that MFM sheds light on the advantages of learning factorizing multimodal representations and potentially opens up new horizons for multimodal machine learning. 

A PROOF OF PROPOSITION 1:

To simplify the proof, we first prove it for the unimodal case by considering the Wasserstein distance between P X,Y and P X, Ŷ .

A.1 UNIMODAL JOINT-DISTRIBUTION WASSERSTEIN DISTANCE Proposition 2. For any functions

EQUATION

)

where W c is the Wasserstein distance under cost function c X and c Y , P Z is the prior over Z = [Z a , Z y ] and Q Z is the aggregated posterior of the proposed inference distribution Q(Z X).

Proof: See the following.

To begin the proof, we abuse some notations as follows.

By definition, the Wasserstein distance under cost function c between P X,Y and P X, Ŷ is

EQUATION

where c (X, Y), ( X, Ŷ) ∶ (X , Y) × (X , Y) → R + is any measurable cost function. P (X, Y) ∼ P X,Y , ( X, Ŷ) ∼ P X, Ŷ is the set of all joint distributions of (X, Y), ( X, Ŷ) with marginals P X,Y and P X, Ŷ , respectively. Note that c (X, Y), (

X, Ŷ) = c X X, X + c Y Y, Ŷ .

Next, we denote the set of all joint distributions of (X, Y, X, Ŷ, Z) such that (X, Y) ∼ P X,Y , ( X, Ŷ, Z) ∼ P X, Ŷ,Z , and (X, Y) ⊥ ⊥ ( X, Ŷ) Z as P X,Y, X, Ŷ,Z . P X,Y, X, Ŷ and P X,Y,Z are the sets of the marginals (X, Y, X, Ŷ) and (X, Y, Z) induced by P X,Y, X, Ŷ,Z .

We now introduce two Lemmas to help the proof.

Lemma 1. P ( X, Ŷ Z = z) are Dirac for all z ∈ Z.

Proof: First, we have X = F (G a (Z a ), G y (Z y )) and Ŷ = D(G y (Z y )) with Z = {Z a , Z y }. Since the functions F, G a , G y , D are all deterministic, then P ( X, Ŷ Z) are Dirac measures. ◻ Lemma 2. P P X,Y , P X, Ŷ = P X,Y, X, Ŷ when P ( X, Ŷ Z = z) are Dirac for all z ∈ Z.

Proof: When X, Ŷ are deterministic functions of Z, for any A in the sigma-algebra induced by X, Ŷ, we have

E[I [ X, Ŷ∈A] X, Y, Z] = E[I [ X, Ŷ∈A] Z].

Therefore, this implies that (X, Y) ⊥ ⊥ ( X, Ŷ) Z which concludes the proof. A similar argument is made in Lemma 1 of (Tolstikhin et al., 2017) .

◻:

Now, we use the fact that P P X,Y , P X, Ŷ = P X,Y, X, Ŷ (Lemma 1 + Lemma 2),

c (X, Y), ( X, Ŷ) = c X X, X + c Y Y, Ŷ , X = F G a (Z a ), G y (Z y )

, and Ŷ = D G y (Z y ) , Eq. equation 7 becomes

EQUATION

)

Note that in Eq. equation 8, P X,Y,Z = P (X, Y) ∼ P X,Y , Z ∼ P Z and with a proposed Q(Z X), we can rewrite Eq. equation 8 as

EQUATION

) ∎ A.2 FROM UNIMODAL TO MULTIMODAL

The proof is similar to Proposition 2, and we present a sketch to it. We can first show P ( X1∶M , Ŷ Z = z) are Dirac for all z ∈ Z. Then we use the fact that c (X 1∶M , Y),

( X1∶M , Ŷ) = ∑ M i=1 c X i X i , Xi + c Y Y, Ŷ

. Finally, we follow the tower rule of expectation and the conditional independence property similar to the proof in Proposition 2 and this concludes the proof.

∎ B FULL BASELINE MODELS & RESULTS:

For a detailed description of the baselines, we point the reader to MFN (Zadeh et al., 2018a) , MARN (Zadeh et al., 2018b) , TFN (Zadeh et al., 2017) , BC-LSTM (Poria et al., 2017) , MV-LSTM (Rajagopalan et al., 2016) , EF-LSTM (Hochreiter & Schmidhuber, 1997; Graves et al., 2013; Schuster & Paliwal, 1997) , DF (Nojavanasghari et al., 2016) , MV-HCRF (Song et al., 2012; 2013) , EF-HCRF (Quattoni et al., 2007; Morency et al., 2007) , THMM (Morency et al., 2011) , SVM-MD (Zadeh et al., 2016) and RF (Breiman, 2001) .

We use the following extra notations for full descriptions of the baseline models described in Section 3.2, paragraph 3: Variants of EF-LSTM: EF-LSTM (Early Fusion LSTM) uses a single LSTM (Hochreiter & Schmidhuber, 1997) on concatenated multimodal inputs. We also implement the EF-SLSTM (stacked) (Graves et al., 2013) , EF-BLSTM (bidirectional) (Schuster & Paliwal, 1997) and EF-SBLSTM (stacked bidirectional) versions.

Variants of EF-HCRF: EF-HCRF: (Hidden Conditional Random Field) (Quattoni et al., 2007) uses a HCRF to learn a set of latent variables conditioned on the concatenated input at each time step.

EF-LDHCRF (Latent Discriminative HCRFs) (Morency et al., 2007) are a class of models that learn hidden states in a CRF using a latent code between observed concatenated input and hidden output. EF-HSSHCRF: (Hierarchical Sequence Summarization HCRF) (Song et al., 2013 ) is a layered model that uses HCRFs with latent variables to learn hidden spatio-temporal dynamics.

Variants of MV-HCRF: MV-HCRF: Multi-view HCRF (Song et al., 2012) is an extension of the HCRF for Multi-view data, explicitly capturing view-shared and view specific sub-structures. MV-LDHCRF: (Morency et al., 2007) is a variation of the MV-HCRF model that uses LDHCRF instead of HCRF. MV-HSSHCRF: (Song et al., 2013) further extends EF-HSSHCRF by performing Multiview hierarchical sequence summary representation.

In the following, we provide the full results for all baselines models described in Section 3.2, paragraph 3. Table 3 contains results for multimodal speaker traits recognition on the POM dataset. Table 4 contains results for the multimodal sentiment analysis on the CMU-MOSI, ICT-MMMO, YouTube, and MOUD datasets. Table 5 contains results for multimodal emotion recognition on the IEMOCAP dataset. MFM consistently achieves state-of-the-art or competitive results for all six multimodal datasets. We believe that by our MFM design, the multimodal discriminative factor F y has successfully learned more meaningful representations by distilling discriminative features. This highlights the benefit of learning factorized multimodal representations towards discriminative tasks.

C MULTIMODAL FEATURES:

For each of the multimodal time series datasets as mentioned in Section 3.2, paragraph 3, we extracted the following multimodal features: Language: We use pre-trained word embeddings (glove.840B.300d) (Pennington et al., 2014) to convert the video transcripts into a sequence of 300 dimensional word vectors. Visual: We use Facet (iMotions, 2017) to extract a set of features including per-frame basic and advanced emotions and facial action units as indicators of facial muscle movement (Ekman et al., 1980; Ekman, 1992) . Acoustic: We use COVAREP (Degottex et al., 2014) to extract low level acoustic features including 12 Mel-frequency cepstral coefficients (MFCCs), pitch tracking and voiced/unvoiced segmenting features, glottal source parameters, peak slope parameters and maxima dispersion quotients. To reach the same time alignment between different modalities we choose the granularity of the input to be at the level of words. The words are aligned with audio using P2FA (Yuan & Liberman, 2008) to get their exact utterance times. We use expected feature values across the entire word for visual and acoustic features since they are extracted at a higher frequencies.

We make a note that the features for some of these datasets are constantly being updated. The authors of Zadeh et al. (2018a) notified us of a discrepancy in the sampling rate for acoustic feature extraction in the ICT-MMMO, YouTube and MOUD datasets which led to inaccurate word-level alignment between the three modalities. They publicly released the updated multimodal features. We performed all experiments on the latest versions of these datasets which can be accessed from https: //github.com/A2Zadeh/CMU-MultimodalSDK. All baseline models were retrained with extensive hyperparameter search for fair comparison.

D INFORMATION AND GRADIENT-BASED INTERPRETATION:

Information-Based Interpretation: We choose the normalized Hilbert-Schmidt Independence Criterion (Gretton et al., 2005; Wu et al., 2018) as the approximation (see Sugiyama & Yamada (2012) ; Wu et al. (2018) ) of our MI measure:

EQUATION

where ⋅ represents y or a i , n is the number of {F ⋅ , Xi } pairs,

H = I -1 n 11 ⊺ , K F⋅ ∈ R n×n is the Gram matrix of F ⋅ with K F⋅ ij = k 1 (F ⋅i , F ⋅j ), K Xi ∈ R n×n is the Gram matrix of Xi with K Xijk = k 2 ( Xij , Xik ). k 1 (⋅, ⋅) and k 2 (⋅, ⋅) are predefined kernel functions.

The most common choice for the kernel is the RBF kernel. However, if we consider time series data with various time steps, we need to either perform data augmentation or choose another kernel choice. For example, we can adopt the Global Alignment Kernel (Cuturi et al., 2007) which considers the 

↑ 0.1 ↑ 0.5 ↑ 2.6 ↑ 3.3 ↑ 0.5 ↑ 0.7 ↑ 0.1 - ↑ 0.2 - ↑ 0.

3 alignment between two varying-length time series when computing the kernel score between them.

To simplify our analysis, we choose to augment data before we calculate the kernel score with the RBF kernel. More specifically, we perform averaging over time series data:

EQUATION

The bandwidth of the RBF kernel is set as 1.0 throughout the experiments. , i ∈ {( )anguage, (v)isual, (a)coustic} for the POM dataset for personality traits prediction. Here, we provide an additional interpretation result for the POM dataset in Table 6 . We observe that the language modality is also the most informative while the visual and acoustic modalities are almost equally informative. This result is in agreement with behavioral studies which have observed that non-verbal behaviors are particularly informative of personality traits (Guimond & Massrieh, 2012; Levine et al., 2009; Mohammadi et al., 2010) . For example, the same sentence "this movie was great" can convey significantly different messages on speaker confidence depending on whether it was said in a loud and exciting voice, with eye contact, or powerful gesticulation.

Gradient-Based Interpretation: MFM reconstructs x i as follows:

EQUATION

)

Equation equation 12 also explains how we obtain f y ∼ P (F y X 1∶M = x 1∶M ). The gradient flow through time is defined as:

EQUATION

E ENCODER AND DECODER DESIGN FOR MULTIMODAL SYNTHETIC IMAGE DATASET:

For experiments on the multimodal synthetic image dataset, we use convolutional+fully-connected layers for the encoder and deconvolutional+fully-connected layers for the decoder (Zeiler et al., 2010) . Different convolutional layers are each applied on the input SVHN and MNIST images to learn modality-specific generative factors. Next, we concatenate the features from two more convolutional layers on SVHN and MNIST to learn the multimodal-discriminative factor. The multimodal discriminative factor is passed through fully-connected layers to predict the label. For generation, we concatenate the multimodal discriminative factors and the modality-specific generative factor together and use a deconvolutional layer to generate digits.

F ENCODER AND DECODER DESIGN FOR MULTIMODAL TIME SERIES DATASETS:

Figure 5 illustrates how MFM operates on multimodal time series data. The encoder Q(Z y X 1∶M ) can be parametrized by any model that performs multimodal fusion (Nojavanasghari et al., 2016; Zadeh et al., 2018a) . We choose the Memory Fusion Network (MFN) (Zadeh et al., 2018a) as our encoder Q(Z y X 1∶M ). We use encoder LSTM networks and decoder LSTM networks (Cho et al., 2014) to parametrize functions Q(Z a1∶M X 1∶M ) and F 1∶M respectively, and FCNNs to parametrize functions G y , G a{1∶M } and D. ) can be parametrized by any model that performs multimodal fusion (Nojavanasghari et al., 2016; Zadeh et al., 2018a) . We use encoder LSTM networks and decoder LSTM networks (Cho et al., 2014) to parametrize functions Q(Z a1∶M X 1∶M ) and F 1∶M respectively, and FCNNs to parametrize functions G y , G a{1∶M } and D.

G SURROGATE INFERENCE GRAPHICAL MODEL:

We illustrate the surrogate inference for addressing the missing modalities issue in Figure 6 . The surrogate inference model infers the latent codes given the present modalities. These inferred latent codes can then be used for reconstructing the missing modalities or label prediction in the presence of missing modalities. 

H COMPARISON WITH HSU & GLASS (2018):

A similar approach for factorizing the latent factors was recently proposed by Hsu & Glass (2018) in work that was performed independently and in parallel. In comparison with MFM, there are several major differences that can be categorized into the (1) prior matching discrepancy, (2) inference network, (3) discriminative objective, (4) multimodal fusion, (5) scale of experiments.

1. MFM uses MMD(Q Z , P Z ) (see Equation 4) as the prior matching discrepancy while Hsu & Glass (2018) use KL(Q Z X , P Z )) (see Section 2.2.1 in Hsu & Glass (2018) ). 2. MFM considers multimodal and unimodal inference in a single network (see Figure 1 They further propose to match the coherence between these two networks using an additional loss term (see Equation 7in Hsu & Glass (2018) ). 3. MFM learns to predict the labels using a generative framework (see Figure 1 2018)). 5. We evaluate the performance of MFM over a much larger scale of datasets. We perform experiments on six multimodal time-series datasets that take on the form of videos with the language, visual, and acoustic modalities. These datasets span three core research areas of multimodal personality traits recognition, multimodal sentiment analysis, and multimodal emotion recognition. On the other hand, Hsu & Glass (2018) evaluates their model on a spoken digit dataset which randomly combines a digit image with a spoken digit (see Section 4 in Hsu & Glass (2018) ). MFM further considers experiments to evaluate reconstruction and prediction in the presence of missing modalities (see Section 3.2) which Hsu & Glass (2018) do not. Lastly, we compares to over 20 baseline models in our experiments (see Section 3.2) and explore the choice of various multimodal encoders in MFM. Hsu & Glass (2018) only compares to the JMVAE baseline model (Suzuki et al., 2016a) which resembles the M D model in our ablation study (see Section 4.4 in Hsu & Glass (2018) ).

In Table 7 , we provide a comparison on the CMU-MOSI, ICT-MMMO, YouTube and MOUD datasets to test the disentanglement and prediction performance for the model described in Hsu & Glass (2018) . These experimental results show that across these datasets and metrics, MFM performs better than the model proposed in Hsu & Glass (2018) . We would like to highlight that at the time of submission, the code for (Hsu & Glass, 2018) had not been made public and we reimplemented their model to experiment on our datasets.

I COMPARISON WITH β-VAE:

Although β-VAE (Higgins et al., 2017) was designed to handle unimodal data, we provide an extension to multimodal data. To achieve this, we set the choice of prior matching discrepancy as the KL-divergence KL(Q Z X , P Z ) and set β large (i.e. β ∈ {10, 50, 100, 200}) to encourage disentanglement of latent variables. We train a β-VAE to model multimodal data using the same 

:



Figures and Tables:

Figure 1: Illustration of the proposed Multimodal Factorization Model (MFM) with three modalities. MFM factorizes multimodal representations into multimodal discriminative factors Fy and modality-specific generative factors Fa {1∶M } . (a) MFM Generative Network with latent variables {Zy, Za {1∶M } }, factors {Fy, Fa {1∶M } }, generated multimodal data X1∶3 and labels Ŷ. (b) MFM Inference Network. (c) MFM Neural Architecture. Best viewed zoomed in and in color.

(a) MFM generative network for multimodal image dataset SVHN+MNIST, (b) unimodal and multimodal classification accuracies, and (c) conditional generation for SVHN and MNIST digits. MFM shows improved capabilities in digit prediction as well as flexible generation of both images based on labels and styles.

Figure 3: Models used in the ablation studies of MFM. Each model removes a design component from our model. Modality reconstruction and sentiment prediction results are reported on CMU-MOSI with best results in bold. Factorizing multimodal representations into multimodal discriminative factors and modality-specific generative factors are crucial for improved performance. nington et al., 2014)), visual (Facet (iMotions, 2017)), and acoustic (COVAREP (Degottex et al., 2014)) modalities. For a detailed description of feature extraction, please refer to the appendix. We consider the following six datasets across three domains: 1) Multimodal Personality Trait Recognition: POM (Park et al., 2014) contains 903 movie review videos annotated for the following personality traits: confident (con), passionate (pas), voice pleasant (voi), dominant (dom), credible (cre), vivid (viv), expertise (exp), entertaining (ent), reserved (res), trusting (tru), relaxed (rel), outgoing (out), thorough (tho), nervous (ner), persuasive (per) and humorous (hum). The short form is indicated in parenthesis. 2) Multimodal Sentiment Analysis: CMU-MOSI (Zadeh et al., 2016) is a collection of 2199 monologue opinion video clips annotated with sentiment. ICT-MMMO (Wöllmer et al., 2013) consists of 340 online social review videos annotated for sentiment. YouTube (Morency et al., 2011) contains 269 product review and opinion video segments from YouTube each annotated for sentiment. MOUD (Perez-Rosas et al., 2013) consists of 79 product review videos in Spanish. Each video consists of multiple segments labeled as either positive, negative or neutral sentiment.3) Multimodal Emotion Recognition: IEMOCAP(Busso et al., 2008) consists of 302 videos of recorded dyadic dialogues. The videos are divided into multiple segments each annotated for the presence of 6 discrete emotions (happy, sad, angry, frustrated, excited and neutral), resulting in a total of 7318 segments in the dataset. We report results using the following metrics: Acc_C = multiclass accuracy across C classes, F1 = F1 score, MAE = Mean Absolute Error, r = Pearson's correlation.

and 2) adding generative capabilities to the model improves performance (M C > M A , M E > M B ). For both sentiment prediction and modality reconstruction, 3) factorizing into separate generative and discriminative factors improves performance (M E > M D ), and 4) using modality-specific generative factors outperforms multimodal generative factors (MFM > M E ). These observations support our design decisions of factorizing multimodal representations into multimodal discriminative factors and modality-specific generative factors.

Figure 4: Analyzing the multimodal representations learnt in MFM via information-based (entire dataset) and gradient-based interpretation methods (single video) on CMU-MOSI.

Figure 5: Recurrent neural architecture for MFM. The encoder Q(Z y X 1∶M) can be parametrized by any model that performs multimodal fusion(Nojavanasghari et al., 2016; Zadeh et al., 2018a). We use encoder LSTM networks and decoder LSTM networks(Cho et al., 2014) to parametrize functions Q(Z a1∶M X 1∶M ) and F 1∶M respectively, and FCNNs to parametrize functions G y , G a{1∶M } and D.

Figure 6: The surrogate inference graphical model to deal with missing modalities in MFM. Red lines denote original inference in MFM and green lines denote surrogate inference to infer latent codes given present modalities.

(b)), whileHsu & Glass (2018) considers separate networks (see Figure1inHsu & Glass (2018)).

(a)), whileHsu & Glass (2018) use an additional hinge loss to separate the latent factors from different labels (see Equation 9 in Hsu & Glass (2018)). 4. MFM is a flexible framework that can be combined with any multimodal fusion encoder (see Section 2.4), while Hsu & Glass (2018) considers a fixed multimodal encoder (similar to early fusion) (see Section 4.1 in Hsu & Glass (

Results for multimodal speaker traits recognition on POM, multimodal sentiment analysis on CMU-MOSI, ICT-MMMO, YouTube, MOUD, and multimodal emotion recognition on IEMOCAP. SOTA1 and SOTA2 refer to the previous best and second best state-of-the-art respectively, and ∆ SOT A shows improvement over SOTA1. Symbols depict the baseline giving the result: # MFN, ‡ MARN, * TFN, † BC-LSTM, ◇ MV-LSTM, § EF-LSTM, ♭ DF, ♡ SVM, • RF. For detailed tables with results for all models, please refer to the appendix. SOTA2 0.359 † 0.425 † 0.166 ‡ 0.235 ‡ 0.358 † 0.417 † 0.450 † 0.378 ‡ 0.295 ◇ 0.237 ◇ 0.215 ‡ 0.238 ◇ 0.363 † 0.258 ◇ 0.344 † 0.319 † SOTA1 0.395 # 0.428 # 0.193 # 0.313 # 0.367 # 0.431 # 0.452 # 0.395 # 0.333 # 0.296 # 0.255 # 0.259 # 0.381 # 0.318 # 0.377 # 0.386 #

Table content: <table><tr><td>Dataset</td><td/><td/><td/><td/><td/><td/><td colspan="4">POM Personality Traits</td><td/><td/><td/><td/><td/><td/></tr><tr><td>Task</td><td>Con</td><td>Pas</td><td>Voi</td><td>Dom</td><td>Cre</td><td>Viv</td><td>Exp</td><td>Ent</td><td>Res</td><td>Tru</td><td>Rel</td><td>Out</td><td>Tho</td><td>Ner</td><td>Per</td><td>Hum</td></tr><tr><td>Metric</td><td/><td/><td/><td/><td/><td/><td/><td>r</td><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>MFM</td><td colspan="16">0.431 0.450 0.197 0.411 0.380 0.448 0.467 0.452 0.368 0.212 0.309 0.333 0.404 0.333 0.334 0.408</td></tr><tr><td colspan="10">∆ SOT A ↑ 0.036 ↑ 0.022 ↑ 0.004 ↑ 0.097 ↑ 0.013 ↑ 0.017 ↑ 0.015 ↑ 0.057 ↑ 0.035</td><td>-</td><td colspan="4">↑ 0.054 ↑ 0.074 ↑ 0.023 ↑ 0.015</td><td>-</td><td>↑ 0.022</td></tr><tr><td>Dataset</td><td/><td/><td/><td colspan="2">CMU-MOSI</td><td/><td/><td colspan="3">ICT-MMMO</td><td colspan="3">YouTube</td><td colspan="3">MOUD</td></tr><tr><td>Task</td><td/><td/><td/><td>Sentiment</td><td/><td/><td/><td colspan="3">Sentiment</td><td colspan="3">Sentiment</td><td colspan="3">Sentiment</td></tr><tr><td>Metric</td><td>Acc_7</td><td/><td>Acc_2</td><td>F1</td><td colspan="2">MAE</td><td>r</td><td colspan="2">Acc_2</td><td>F1</td><td colspan="2">Acc_3</td><td>F1</td><td>Acc_2</td><td/><td>F1</td></tr><tr><td>SOTA2</td><td>34.1 #</td><td/><td>77.1  ‡</td><td>77.0  ‡</td><td colspan="2">0.968  ‡</td><td>0.625  ‡</td><td colspan="2">72.5  *</td><td>72.6  *</td><td>48.3  ‡</td><td/><td>45.1  †</td><td>81.1 #</td><td/><td>80.4 #</td></tr><tr><td>SOTA1</td><td>34.7  ‡</td><td/><td>77.4 #</td><td>77.3 #</td><td colspan="2">0.965 #</td><td>0.632 #</td><td colspan="2">73.8 #</td><td>73.1 #</td><td colspan="2">51.7 #</td><td>51.6 #</td><td>81.1  ‡</td><td/><td>81.2  ‡</td></tr><tr><td>MFM</td><td>36.2</td><td/><td>78.1</td><td>78.1</td><td colspan="2">0.951</td><td>0.662</td><td>81.3</td><td/><td>79.2</td><td>53.3</td><td/><td>52.4</td><td>82.1</td><td/><td>81.7</td></tr></table>

The effect of missing modalities on multimodal data reconstruction and sentiment prediction on CMU-MOSI. MFM with surrogate inference is able to better handle missing modalities during test time as compared to the purely generative (Seq2Seq) or purely discriminative baselines. A by optimizing a hybrid generative-discriminative objective over modality-specific factors. M D extends M B by optimizing a hybrid generative-discriminative objective over a joint multimodal factor (resemblingSrivastava & Salakhutdinov (2012)). M E factorizes the representation into separate generative and discriminative factors. Finally, MFM is obtained from M E by using modality-specific generative factors instead of a joint multimodal generative factor.

Table content: <table><tr><td>Task</td><td/><td colspan="2">X⋅ Reconstruction</td><td/><td colspan="2">Ŷ Prediction</td></tr><tr><td>Metric</td><td colspan="3">MSE ( ) MSE (a) MSE (v)</td><td colspan="2">Acc_7 Acc_2</td><td>F1</td><td>MAE</td><td>r</td></tr><tr><td/><td colspan="5">Purely Generative and Discriminative Baselines</td><td/></tr><tr><td>(anguage) missing</td><td>0.0411</td><td>-</td><td>-</td><td>19.4</td><td>59.6</td><td colspan="2">59.7 1.386 0.225</td></tr><tr><td>a(udio) missing</td><td>-</td><td>0.0533</td><td>-</td><td>34.0</td><td>73.5</td><td colspan="2">73.4 1.024 0.615</td></tr><tr><td>v(isual) missing</td><td>-</td><td>-</td><td>0.0220</td><td>33.7</td><td>75.4</td><td colspan="2">75.4 0.996 0.634</td></tr><tr><td/><td/><td colspan="3">Multimodal Factorization Model (MFM)</td><td/><td/></tr><tr><td>(anguage) missing</td><td>0.0403</td><td>-</td><td>-</td><td>21.7</td><td>62.0</td><td colspan="2">61.7 1.313 0.236</td></tr><tr><td>a(udio) missing</td><td>-</td><td>0.0468</td><td>-</td><td>35.4</td><td>74.3</td><td colspan="2">74.3 1.011 0.603</td></tr><tr><td>v(isual) missing</td><td>-</td><td>-</td><td>0.0215</td><td>35.0</td><td>76.4</td><td colspan="2">76.3 0.990 0.635</td></tr><tr><td>all present</td><td>0.0391</td><td>0.0384</td><td>0.0182</td><td>36.2</td><td>78.1</td><td colspan="2">78.1 0.951 0.662</td></tr></table>

Shengjia Zhao, Jiaming Song, and Stefano Ermon. Infovae: Information maximizing variational autoencoders. arXiv preprint arXiv:1706.02262, 2017. Zhenyao Zhu, Ping Luo, Xiaogang Wang, and Xiaoou Tang. Multi-view perceptron: a deep model for learning face identity and view representations. In NIPS, 2014.

Table content: <table/>

Results for personality trait recognition on the POM dataset. The best results are highlighted in bold and ∆ SOT A shows the change in performance over previous state of the art. Improvements are highlighted in green. MFM achieves state-of-the-art or competitive performance on all datasets and metrics. ↑ 0.036 ↑ 0.022 ↑ 0.004 ↑ 0.097 ↑ 0.013 ↑ 0.017 ↑ 0.015 ↑ 0.057 ↑ 0.035 -↑ 0.054 ↑ 0.074 ↑ 0.023 ↑ 0.015 -↑ 0.022

Table content: <table><tr><td>Dataset</td><td/><td/><td/><td/><td/><td colspan="4">POM Speaker Personality Traits</td><td/><td/><td/><td/><td/><td/></tr><tr><td>Task</td><td>Con</td><td>Pas</td><td>Voi</td><td>Dom Cre</td><td>Viv</td><td>Exp</td><td>Ent</td><td>Res</td><td>Tru</td><td>Rel</td><td>Out</td><td>Tho</td><td>Ner</td><td>Per</td><td>Hum</td></tr><tr><td>Metric</td><td/><td/><td/><td/><td/><td/><td>r</td><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>Majority</td><td colspan="15">-0.041 -0.029 -0.104 -0.031 -0.122 -0.044 -0.065 -0.105 0.006 -0.077 -0.024 -0.085 -0.130 0.097 -0.127 -0.069</td></tr><tr><td>SVM</td><td colspan="15">0.063 0.086 -0.004 0.141 0.113 0.076 0.134 0.141 0.166 0.168 0.104 0.066 0.134 0.068 0.064 0.147</td></tr><tr><td>DF</td><td colspan="15">0.240 0.273 0.017 0.139 0.112 0.173 0.118 0.217 0.148 0.143 0.019 0.093 0.041 0.136 0.168 0.259</td></tr><tr><td>EF-LSTM</td><td colspan="15">0.200 0.302 0.031 0.079 0.170 0.244 0.265 0.240 0.142 0.062 0.083 0.152 0.260 0.105 0.217 0.227</td></tr><tr><td colspan="16">EF-SLSTM 0.221 0.327 0.042 0.151 0.177 0.239 0.268 0.248 0.204 0.069 0.092 0.215 0.252 0.159 0.218 0.196</td></tr><tr><td colspan="16">EF-BLSTM 0.162 0.289 -0.034 0.135 0.191 0.279 0.274 0.231 0.184 0.154 0.093 0.147 0.245 0.166 0.243 0.272</td></tr><tr><td colspan="16">EF-SBLSTM 0.174 0.310 0.021 0.088 0.170 0.224 0.261 0.241 0.155 0.163 0.097 0.120 0.215 0.121 0.216 0.171</td></tr><tr><td colspan="16">MV-LSTM 0.358 0.416 0.131 0.146 0.280 0.347 0.323 0.326 0.295 0.237 0.119 0.238 0.284 0.258 0.239 0.317</td></tr><tr><td>BC-LSTM</td><td colspan="15">0.359 0.425 0.081 0.234 0.358 0.417 0.450 0.361 0.293 0.109 0.075 0.078 0.363 0.184 0.344 0.319</td></tr><tr><td>TFN</td><td colspan="15">0.089 0.201 0.030 0.020 0.124 0.204 0.171 0.223 -0.051 -0.064 0.114 0.060 0.048 -0.002 0.106 0.213</td></tr><tr><td>MARN</td><td colspan="15">0.340 0.410 0.166 0.235 0.340 0.374 0.406 0.378 0.282 0.147 0.215 0.204 0.348 0.235 0.303 0.287</td></tr><tr><td>MFN</td><td colspan="15">0.395 0.428 0.193 0.313 0.367 0.431 0.452 0.395 0.333 0.296 0.255 0.259 0.381 0.318 0.377 0.386</td></tr><tr><td>MFM</td><td colspan="15">0.431 0.450 0.197 0.411 0.380 0.448 0.467 0.452 0.368 0.212 0.309 0.333 0.404 0.333 0.334 0.408</td></tr><tr><td>∆ SOT A</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr></table>

Sentiment prediction results on CMU-MOSI, ICT-MMMO, YouTube and MOUD. The best results are highlighted in bold and ∆ SOT A shows the change in performance over previous state of the art (SOTA). Improvements are highlighted in green. MFM achieves state-of-the-art or competitive performance on all datasets and metrics.

Table content: <table/>

Emotion recognition results on IEMOCAP test set. The best results are highlighted in bold and ∆ SOT A shows the change in performance over previous SOTA. Improvements are highlighted in green. MFM achieves state-of-the-art or competitive performance on all datasets and metrics.

Table content: <table><tr><td>Dataset</td><td/><td/><td/><td/><td/><td colspan="2">IEMOCAP Emotions</td><td/><td/><td/><td/><td/></tr><tr><td>Task</td><td colspan="2">Happy</td><td>Sad</td><td/><td colspan="2">Angry</td><td colspan="2">Frustrated</td><td colspan="2">Excited</td><td colspan="2">Neutral</td></tr><tr><td>Metric</td><td>Acc_2</td><td>F1</td><td>Acc_2</td><td>F1</td><td>Acc_2</td><td>F1</td><td>Acc_2</td><td>F1</td><td>Acc_2</td><td>F1</td><td>Acc_2</td><td>F1</td></tr><tr><td>Majority</td><td>85.6</td><td>79.0</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>SVM</td><td>86.1</td><td>81.5</td><td>81.1</td><td>78.8</td><td>82.5</td><td>82.4</td><td>77.3</td><td>71.1</td><td>86.4</td><td>86.0</td><td>65.2</td><td>64.9</td></tr><tr><td>RF</td><td>85.5</td><td>80.7</td><td>80.1</td><td>76.5</td><td>81.9</td><td>82.0</td><td>78.6</td><td>75.3</td><td>88.9</td><td>85.1</td><td>63.2</td><td>57.3</td></tr><tr><td>THMM</td><td>85.6</td><td>79.2</td><td>79.5</td><td>79.8</td><td>79.3</td><td>73.0</td><td>71.6</td><td>69.6</td><td>86.0</td><td>84.6</td><td>58.6</td><td>46.4</td></tr><tr><td>EF-HCRF</td><td>85.7</td><td>79.2</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>EF-LDHCRF</td><td>85.8</td><td>79.5</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>MV-HCRF</td><td>15.0</td><td>4.9</td><td>79.4</td><td>70.3</td><td>24.2</td><td>9.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>MV-LDHCRF</td><td>85.7</td><td>79.2</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>CMV-HCRF</td><td>14.4</td><td>3.6</td><td>79.4</td><td>70.3</td><td>24.2</td><td>9.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td colspan="2">CMV-LDHCRF 85.8</td><td>79.5</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>EF-HSSHCRF</td><td>85.8</td><td>79.5</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td colspan="2">MV-HSSHCRF 85.8</td><td>79.5</td><td>79.4</td><td>70.3</td><td>75.8</td><td>65.4</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>DF</td><td>86.0</td><td>81.0</td><td>81.8</td><td>81.2</td><td>75.8</td><td>65.4</td><td>78.4</td><td>76.8</td><td>89.6</td><td>84.7</td><td>59.1</td><td>44.0</td></tr><tr><td>EF-LSTM</td><td>85.2</td><td>83.3</td><td>82.1</td><td>81.1</td><td>84.5</td><td>84.3</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>68.2</td><td>67.1</td></tr><tr><td>EF-SLSTM</td><td>85.6</td><td>79.0</td><td>80.7</td><td>80.2</td><td>82.8</td><td>82.2</td><td>77.5</td><td>69.7</td><td>89.3</td><td>86.2</td><td>68.8</td><td>68.5</td></tr><tr><td>EF-BLSTM</td><td>85.0</td><td>83.7</td><td>81.8</td><td>81.6</td><td>84.2</td><td>83.3</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>67.1</td><td>66.6</td></tr><tr><td>EF-SBLSTM</td><td>86.0</td><td>84.2</td><td>80.2</td><td>80.5</td><td>85.2</td><td>84.5</td><td>79.5</td><td>70.4</td><td>89.6</td><td>84.7</td><td>67.8</td><td>67.1</td></tr><tr><td>MV-LSTM</td><td>85.9</td><td>81.3</td><td>80.4</td><td>74.0</td><td>85.1</td><td>84.3</td><td>79.5</td><td>73.8</td><td>88.9</td><td>85.8</td><td>67.0</td><td>66.7</td></tr><tr><td>BC-LSTM</td><td>84.9</td><td>81.7</td><td>83.2</td><td>81.7</td><td>83.5</td><td>84.2</td><td>80.0</td><td>76.1</td><td>86.9</td><td>85.4</td><td>67.5</td><td>64.1</td></tr><tr><td>TFN</td><td>84.8</td><td>83.6</td><td>83.4</td><td>82.8</td><td>83.4</td><td>84.2</td><td>74.1</td><td>74.3</td><td>75.6</td><td>78.0</td><td>67.5</td><td>65.4</td></tr><tr><td>MARN</td><td>86.7</td><td>83.6</td><td>82.0</td><td>81.2</td><td>84.6</td><td>84.2</td><td>79.5</td><td>76.6</td><td>89.6</td><td>87.1</td><td>66.8</td><td>65.9</td></tr><tr><td>MFN</td><td>90.1</td><td>85.3</td><td>85.8</td><td>79.2</td><td>87.0</td><td>86.0</td><td>80.3</td><td>76.9</td><td>89.8</td><td>86.3</td><td>71.8</td><td>61.7</td></tr><tr><td>MFM</td><td>90.2</td><td>85.8</td><td>88.4</td><td>86.1</td><td>87.5</td><td>86.7</td><td>80.4</td><td>74.5</td><td>90.0</td><td>87.1</td><td>72.1</td><td>68.1</td></tr><tr><td>∆ SOT A</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr></table>

Information-Based interpretation results showing ratios r i =

Table content: <table/>

Comparison withHsu & Glass (2018) for sentiment analysis on CMU-MOSI, ICT-MMMO, YouTube, and MOUD. MFM outperforms the baselines across these datasets and metrics.

Table content: <table><tr><td>Dataset</td><td/><td colspan="3">CMU-MOSI</td><td/><td colspan="2">ICT-MMMO</td><td colspan="2">YouTube</td><td>MOUD</td></tr><tr><td>Task</td><td/><td/><td>Sentiment</td><td/><td/><td colspan="2">Sentiment</td><td colspan="2">Sentiment</td><td>Sentiment</td></tr><tr><td>Metric</td><td colspan="3">Acc_7 Acc_2 F1</td><td>MAE</td><td>r</td><td colspan="2">Acc_2 F1</td><td colspan="2">Acc_3 F1</td><td>Acc_2 F1</td></tr><tr><td>without factorization (EF)</td><td>32.4</td><td>74.3</td><td colspan="3">74.3 1.023 0.622</td><td>72.5</td><td>70.9</td><td>44.1</td><td>43.6</td><td>67.0</td><td>64.3</td></tr><tr><td colspan="2">without factorization (MFN) 34.1</td><td>77.4</td><td colspan="3">77.3 0.965 0.632</td><td>73.8</td><td>73.1</td><td>51.7</td><td>51.6</td><td>81.1</td><td>80.4</td></tr><tr><td>Hsu &amp; Glass (2018)</td><td>33.8</td><td>75.2</td><td colspan="3">75.2 1.049 0.584</td><td>77.5</td><td>75.0</td><td>51.7</td><td>48.6</td><td>66.0</td><td>62.9</td></tr><tr><td>MFM</td><td>36.2</td><td>78.1</td><td colspan="3">78.1 0.951 0.662</td><td>81.3</td><td>79.2</td><td>53.3</td><td>52.4</td><td>82.1</td><td>81.7</td></tr></table>

Comparison with β-VAE for multimodal sentiment analysis on CMU-MOSI, ICT-MMMO, YouTube, and MOUD. MFM outperforms β-VAE across these datasets and metrics. factorization as proposed in our model (i.e. modality-specific generative factors Z a{1∶M } and a multimodal discriminative factor Z y ). To provide a fair comparison to our discriminative model, we fine tune by training a classifier on top of the multimodal discriminative factor Z y to the label Y. We provide experimental results in Table8on the CMU-MOSI, ICT-MMMO, YouTube and MOUD datasets. MFM outperforms β-VAE across these datasets and metrics.

Table content: <table><tr><td>Dataset</td><td/><td/><td>CMU-MOSI</td><td/><td/><td colspan="2">ICT-MMMO</td><td colspan="2">YouTube</td><td colspan="2">MOUD</td></tr><tr><td>Task</td><td/><td/><td>Sentiment</td><td/><td/><td colspan="2">Sentiment</td><td colspan="2">Sentiment</td><td colspan="2">Sentiment</td></tr><tr><td colspan="3">Metric Acc_7 Acc_2</td><td>F1</td><td>MAE</td><td>r</td><td>Acc_2</td><td>F1</td><td>Acc_3</td><td>F1</td><td>Acc_2</td><td>F1</td></tr><tr><td>β-VAE</td><td>29.7</td><td>71.3</td><td>71.3</td><td>1.094</td><td>0.552</td><td>65.0</td><td>59.8</td><td>46.7</td><td>31.3</td><td>60.4</td><td>54.2</td></tr><tr><td>MFM</td><td>36.2</td><td>78.1</td><td>78.1</td><td>0.951</td><td>0.662</td><td>81.3</td><td>79.2</td><td>53.3</td><td>52.4</td><td>82.1</td><td>81.7</td></tr></table>

Peer review (iclr2019):
 Multimodality learning is an important topic in multimedia and human computer interaction. How to efficiently leverage the additional information cross multimodality is the key to the task. Authors proposed the Bayesian latent variable model to factorize the multimodality representation into multimodal discriminative factors and modality-specific generating factors, which is interesting. Approximate inference is also proposed to learn this model via a generalised mean-field assumption. The technical quality of the paper is sound and significant, The problem to solve in this paper is also well motivated and important. In general, this is a well-written paper, I have a few minor questions which requires authors for further elaboration. 1. If I understand it correctly, in the current work, the feature Xs are continuous. Does the approach apply to categorical or binary features? 2. In equation(4), MMD is used. How to solve the computation complexity problem since the complexity of MMD is O(n^2)? It is true that the batch size should be small? How to select the hyper-parameters of kernels?