Title: LEARNING STATE REPRESENTATIONS IN COMPLEX SYSTEMS WITH MULTIMODAL DATA

Abstract: Representation learning becomes especially important for complex systems with multimodal data sources such as cameras or sensors. Recent advances in reinforcement learning and optimal control make it possible to design control algorithms on these latent representations, but the field still lacks a large-scale standard dataset for unified comparison. In this work, we present a large-scale dataset and evaluation framework for representation learning for the complex task of landing an airplane. We implement and compare several approaches to representation learning on this dataset in terms of the quality of simple supervised learning tasks and disentanglement scores. The resulting representations can be used for further tasks such as anomaly detection, optimal control, model-based reinforcement learning, and other applications.

1. INTRODUCTION:

In order to act in real world scenarios and control a complex system such as an airplane, car, or industrial facility, an automated agent needs to process complex high-dimensional data coming from different domains: video feeds from cameras, LIDAR sensors on a car, altitude and speed sensors on an airplane, sensors related to the internal state of the system, and so on. An important problem in this regard would be to map this rich stream of multimodal information into a lower-dimensional space that compresses all modalities into a uniform latent representation (embedding); the agent could then use this embedding to learn or otherwise construct control algorithms. Thus, representation learning lies at the heart of optimal control for complex systems with multimodal features.

Over the last decade, deep neural networks have surpassed other methods in processing nearly all modalities of high-dimensional unstructured data, including images, natural language texts, sounds, and time series. One of the most important properties of neural networks that has made the deep learning revolution possible is their ability to extract meaningful low-dimensional representations of raw unstructured input data. Representation learning with deep neural networks is a large and wellestablished area of research (Kingma & Welling, 2013; Hinton & Salakhutdinov, 2006; Radford et al., 2015) . Latent features learned by deep neural networks find applications in various domains, including reinforcement learning for complex systems (Zhang et al., 2018; Watter et al., 2015; van Hoof et al., 2016) . In these works, the authors often propose special techniques and architectures to design the latent space in different ways suitable for further use: make the latent space locally linear, capture the dynamics, and so on. Designing such feature extractors is a complex task, usually done by hand. Moreover, it is hard to compare different architectures, to a large extent because it is far from obvious how to measure the quality of resulting representations. One reason for that is that the best metric for the quality of learned representations would be the quality of the final task in question, which is hard to obtain in reinforcement learning due to the sheer scale of this final task.

A common way to deal with such problems is to use a unified dataset and a unified set of benchmarks, such as, for example, ImageNet and the ILSVRC benchmark in computer vision (Russakovsky et al., 2015) . Such unified datasets might also prove useful for transfer learning tasks. However, the field of reinforcement learning for complex systems is yet to agree on a common representative large-scale dataset. In this work, we present such a multimodal dataset for the representation learning problem together with a unified benchmark framework for feature extractors suitable for a comprehensive comparative evaluation of different feature extractors. This dataset has been gathered with largescale computer simulations based on the X-Plane simulator and consists of data streams from various sensors along with images taken from the frontal camera of the plane. We propose and implement several different metrics for comparison between extracted features. We also survey, implement, and compare different neural architectures for learning multimodal state representations.

The paper is organized as follows. In Section 2 we survey related work on representation learning and evaluation of representations. Section 3 presents the X-Plane dataset and explains its characteristic features. In Section 4, we present the various representation learning models that we have implemented and compared on this dataset. Section 5 shows the evaluation metrics and presents a large-scale comparison across different models, and Section 6 concludes the paper.

2. RELATED WORK:

Modern agent control methods commonly use techniques based on deep learning as feature extractors to deal with complex multimodal data, either explicitly (Zhang et al., 2018; Watter et al., 2015) or implicitly (Mnih et al., 2015) . Explicit techniques separate representations of the learning environment from the agents operating in this enviroment; the environment representation can be learned either independently (Zhang et al., 2018) or in a common end-to-end architecture (Parisi et al., 2017) . There are several major approaches to constructing such models (Lesort et al., 2018) :

(1) autoencoders that reconstruct the input observation data, producing the latent representation between encoder and decoder; (2) forward models that predict the next state either in the latent space or in the raw data, basing the prediction on the current latent representation; (3) inverse models that use two states to predict the actions between them; this approach can be combined with forward models; (4) models with prior knowledge of the system that impose additional constraints on the latent space according to some fundamental system properties such as causality, temporal continuity, or some more specific knowledge about the system (Jonschkowski & Brock, 2015); (5) multimodal temporal fusion models such as the Memory Fusion Network (Zadeh et al., 2018) .

If there is no additional data labeling available such as, e.g., a list of factors, the most direct way to measure the quality of representations would be by evaluating the final quality of solving the main task that a model can achieve based on this representation. However, in real world cases, in particular in reinforcement learning, the main task is often hard to solve and unstable to train, so it cannot be consistently used to evaluate latent representations. Therefore, many indirect ways have been proposed to measure the quality of representations that introduce proxies that can be expected to lead to better solutions of the final control problem. The most common indirect approaches include (see also a comprehensive survey by Lesort et al. (2018) ): (1) task performance, the most intuitive metric, where representation quality is measured by the quality of performing some other relatively simple task, e.g., by predicting some available target variables with simple models that take the latent representation as input (Higgins et al., 2016; van Hoof et al., 2016) ; (2) KNN-MSE, proposed by Lesort et al. (2017) , measures the degree of preservation of the same neighbors between the latent space and the ground truth: KNN-MSE(I) = 1 k I ∈KNN(I,k) φ(I) -φ(I ) , where I is the initial raw input, I is a neighbour of I, and φ is the feature extractor; KNN-MSE is a good metric in situations where the distance in the original input space is well-defined but becomes hard to apply for highly variable multimodal data; (3) a similar approach with humans in the loop is to evaluate whether similar input states according to human evaluation do indeed map to close representations in the latent space (Sermanet et al., 2018) ; however, for complex multimodal inputs this is again inapplicable since a human would not have an intuitive notion of similarity between two sets of several hundred sensor readings; (4) disentanglement scores (Eastwood & Williams, 2018; Higgins et al., 2016) measure the disentanglement (mutual independence) of extracted features; if there are some known generative factors, these metrics assess whether individual elements of the latent representation capture individual generative factors independently; we will consider such metrics in detail in Section 5.1.

As a representative multimodal dataset for a complex system, we have used the X-Plane flight simulator (Laminar Research, 2018) , well-known for its faithful simulation of all systems of an aircraft. It has already been used to solve optimal control problems for aircraft; e.g., Bittar et al. (2014) develop separate control laws for stable and maneuvering flight, and Garcia & Barnes (2009) use X-Plane to simulate a system of several unmanned aerial vehicles. However, to the best of our knowledge this is the first attempt to produce a large-scale dataset for representation learning from X-Plane or any other flight simulator.

3. THE X-PLANE DATASET:

We present the X-Plane Dataset (Dataset, 2018; Code, 2018) as a benchmark for evaluating the quality of state representations (embeddings), where the main problem is to encode the state of some complex system, learn a mapping from high-dimensional multimodal data to a low-dimensional latent representation. We have used the X-Plane simulation environment because it is extremely accurate and can produce a lot of useful and different sensor data. In total, we have recorded 8011 landings, with mean duration of a landing of 115 seconds. For every landing, the dataset contains the readings of 1090 sensors arranged in time series recorded with with a frequency of 5 frames per second, together with the corresponding image taken from the camera located at the front of the airplane; the images are recorded at the same frequency. Table 1 summarizes the various groups of sensors recorded in the dataset, showing the number of sensors in a group together with a brief description. The total dataset size is 93GB uncompressed; it contains sensor readings and about 7M 256 × 256 images, joined into per-flight videos for better compression. To ensure diversity in the dataset, we recorded landings with random perturbations and in different environments: with different airports and runways, time of day, weather conditions etc. We have used at most 12 perturbations (abrupt changes in the environment such as, e.g., wind gusts, malfunction of various systems in the aircraft, and so on) applied during each flight; Table 2 summarizes various failures and perturbations. We have used a custom control unit based on the Boeing 737 guide (Brady, 2014) in order to make automatic landings stable.

During dataset collection, for every landing we chose a random airport and a random runway and spawned a plane in a random position, usually at a distance of 3-5 miles from the runway. There are, in total, 114 different runways in 70 airports used in the dataset. Then we applied different landing conditions-landing speed, flaps position, time of day, visibility, precipitation and so onand enabled the autopilot that lowered the plane on the glide slope path. To achieve a better landing, we disable the autopilot before touchdown, increase the pitch of the plane, and after touchdown we decrease the pitch and enable reverse thrust. At the end of the flight, there are three possible situations: successful landing, aircraft crash, or time limit reached. We set a strict time limit of 160 seconds for every flight. We have also done an airport-stratified split of the dataset into 4 parts: training set for feature extractors, validation set used for early stopping, training set for benchmark supervised learning models, and test set for scoring the benchmark models. We have made the dataset available at (Dataset, 2018) .

4. MODELS:

For experimental evaluation of different representation learning models, we used a wide variety of different models, from basic autoencoders to dynamic actions-aware encoders. In this section, we present these models, from simple to more complex ones. Our models were constructed with four core building blocks: standard recurrent LSTM cell (Hochreiter & Schmidhuber, 1997) , onedimensional convolutions (Conv1D), ConvLSTM (Shi et al., 2015) , and an attention cell (Bahdanau et al., 2014) . The attention cell A(x) operates as follows (see Fig. 2 ):

v i = ReLU(Wx i + b), s = softmax(c V + b), a = Vs, Out = (x, a),

where x is the input, W is a d × N matrix for embedding size d and input dimension N , V is the matrix with columns v i , c is a vector of size d, s are attention weights. We considered the following specific models (Table 4 lists all models with brief descriptions and numerical results).

One-dimensional autoencoders. A simple autoencoder that trains to reconstruct a given set of timesteps (Fig. 2a ). We used four types of autoencoders (AE) in the comparison: autoencoder with 1-layer LSTM cells for encoder and decoder, with 2-layer LSTM cells, with a 6-layer convolutional autoencoder with kernel size 7. In LSTM autoencoders, after the encoder block we used simple averaging of vectors from all timestamps in the time series, using the result as embedding. In Conv1D, after the encoder block we used max-pooling to get a single vector from the time series; in the decoder, it is copied the necessary number of times and fed to another Conv1D block. In our experiments, we used 6 layers of Conv1D with kernel size 7.

Image Autoencoders. We trained two image autoencoders for reconstructing individual pictures from the flight time series. The first autoencoder uses a PCA encoder and a PCA decoder. To obtain Autoregression Models. These models are trained to predict a state at time t + k given some subset of states up to time t. In our experiments, we trained autoregression models to predict with k = 1 (next time step) and k = 30. In this task, as an encoder we have compared LSTM, Conv1D, ConvLSTM, and attention-based architectures, with a simple linear layer after each.

Context Models. These models employ ideas similar to the word2vec CBOW model (Mikolov et al., 2013) . Given a short sequence of timesteps, the model trains to predict several timesteps before and after the sequence; in our experiments, we used window size of 10 before and 10 after given states. To train context models, we have used LSTM, attention cells, and Conv1D blocks.

After an encoder block, we repeated an embedding 2C times, where C is the context size, and then applied a simple linear layer to each vector with different weights for each of the C timestamps.

Multimodal Temporal Encoder. The multimodal temporal encoder (MTE), introduced by Yang et al. ( 2017), is a state of the art model intended to fuse inputs from different modalities (see Fig. 3 ). MTE uses an LSTM unit for each modality and shares hidden states across these units. By doing this, MTE forces the model to learn a fused representation across modalities; formally, it adds a correlation loss that computes the correlation between projections of differ-

ent modalities L corr (H 1 t , H 2 t ) = M i=1 (h 1 ti -H 1 t )(h 2 ti -H 2 t ) M i=1 (h 1 ti -H 1 t ) 2 M i=1 (h 2 ti -H 2 t ) 2

. The final loss is calculated as

L = M i=1 L recon,i -λL corr ,

where M is the batch size, h j t is the representation of modality j at timestep t, H j t = {h j ti } M i=1 , H j t = 1 M M i h j ti , λ is a constant hyperparameter, and h t is the final embedding (see Fig. 3 ). The sensors data is encoded with a three-layer fully connected neural network with dropout and ReLU activations, and a pretrained ResNet18 architecture for the images.

The Dynamics Model. This approach, proposed by Zhang et al. (2018) , introduces a way to decouple the training process (for a reinforcement learning task) into learning the dynamics model and the reward model. The dynamics model (shown in Fig. 4 ) is trained using a combination of four loss functions: reconstruction loss L t,recon (θ enc , θ dec ) = (s t -ŝt ) 2 , state loss L t,state (θ for , θ dec ) = (s t+1 -ŝt+1 ) 2 , forward loss L t,for (θ for , θ enc ) = (z t+1 -ẑt+1 ) 2 , and inverse loss (with a trainable LSTM unit) L t,inv (θ inv ) = (a t -ât ) 2 . The final loss is computed as L(θ dynamics ) = T t=0 (λ dec (L t,recon + L t,state ) + λ for L t,for + λ inv L t,inv ), where ẑt+1 , h t = f for (z t , a t , h t-1 ; θ for ), ât = f inv (z t , z t+1 ; θ inv ), ŝt+1 = f dec (ẑ t+1 ; θ dec ), s t is the state at time t, a t is the action at time t, h t is the hidden state, and z t is the latent representation at time t; λ dec , λ for , and λ inv are (constant) hyperparameters. If the dynamics change, we need only to re-train the LSTM unit and can keep the encoder and decoder unchanged, with the assumption that already learned representation contains all the necessary information about the new dynamics. Since the original approach was presented only for inputs with a single modality, we expanded this idea for the multimodal case. To encode data from the sensors, we used a three-layer fully connected neural network with dropouts and ReLU activations. For the images, we used a pretrained ResNet18 architecture. After that, a linear layer was used to concatenate outputs of encoders for each modality into the embeddings, and the final embedding results by averaging over the time steps.

5. EXPERIMENTAL EVALUATION:



5.1. EVALUATION FRAMEWORK:

As part of the dataset package, we have implemented and made available the framework which is designed to make a comprehensive evaluation of learned representation based on a number of fixed predefined tests. The main idea behind our framework is to combine the two main approaches to measuring representation quality. First, we measure quality of representations by using them as features for a number of simple tasks. This approach was used, in particular, in the Black Box Learning Challenge (Goodfellow et al., 2013) and Unsupervised and Transfer Learning Challenge (Guyon et al., 2011) , whose main objective was to learn a good representation from rich unlabeled data, and representations were evaluated on supervised learning tasks that were not known to the participants.

The second approach is to evaluate the disentanglement in representations. For this we use a QEDR (Quantitative Evaluation of Disentangled Representations) framework proposed by Eastwood & Williams (2018) . The idea is that the ideal representation of data is a vector of separated and independent generative factors for the data (perhaps scaled and permuted). The matrix R, where R ij is the relative importance of code (representation) variable c i in predicting generative factor z j , is used to compute three evaluation metrics for the quality of a representation. Disentanglement is a measure of the degree to which a representation factorizes the factors of variation in original data; for a code variable c i it is defined as D i = 1 -H(P i. ) where H(P i. ) is the entropy of the pseudodistribution P i. ; P ij = R ij / k R ik denotes the "probability" of c i being important for predicting generative factor z j , and total disentanglement is computed as the weighted average i ρ i D i , where ρ i = j R ij / kj R kj is the relative code variable importance. Completeness is a measure of the degree to which a factor of variation in the original data is captured by a single code variable; for a factor z j it is defined as C j = 1-H( P.j ) where H( P.j ) is the entropy of the pseudo-distribution P.j ; total completeness is the average of C j . Informativeness measures the amount of information that a representation captures about the underlying factors of variation; following Eastwood & Williams (2018) , we use normalised root-mean-square (N RM SE) as the informativeness metric. We use failure scores as generative factors since we know the ground truth for them and they are mutually independent. We define R ij = |W ij |, where W is the weight matrix of lasso regression learned on the representation vector to predict the vector of factors.

Our evaluation framework is shown in Fig. 5 . We assume that a feature extractor constructs a single vector representation for a time series of sensor readings. We have implemented disentanglement, completeness, and informativeness scores computed on 12 failures scores as generative factors and 5 supervised learning targets, including: wind speed, turbulence power, vertical acceleration at the time of landing (a quality of touchdown), autoregression (given t states, predict state t + 1), failure classification (in dataset generation we randomly applied various failures to the airplane, and the task is to find out whether a failure is present in a given time series). To evaluate supervised learning, we first apply the provided feature extractors for every training example, then train simple models (namely, a three-layer fully connected neural network with ReLU activations) on these features for each target, and finally compute the accuracy on the hold-out set; these accuracies are the final performance metrics. We have made the code for the evaluation pipeline available at (Code, 2018).

5.2. EVALUATION RESULTS:

Table 4 shows the main results of our evaluation study; in this section we interpret these quantitative results. The models were trained on the X-Plane dataset with regression tasks evaluated on a holdout set; time series lengths varied from 25 to 75 during training. To establish a simple baseline for our models we use the mean target value over the training set as a prediction. In Table 4 , the "MEAN baseline" row shows the absolute MSE values for MEAN and the other rows are normalized to the MEAN baseline. First, we see that different regression tasks have very different models doing well on them; e.g., most models cannot outperform even the mean baseline for wind regression, and models that show good results on wind regression perform poorly on auto regression and vice versa. The effect of images on the results is contradictory: we have trained TS Regression and Dynamics models on three types of data (sensors, images, and both) each, and for TS Regression adding images improves scores on benchmarks, while for the Dynamics model the results deteriorate. The results show that embedding size d needs to be tuned for each model separately; some models even "explode" and give unstable results (very large errors) for some values of d; there is no general correlation between d and benchmark scores. It appears that the Dynamics model trained on sensors only is the best tradeoff between efficiency and quality since it is by far the fastest model. Both disentanglement and completeness scores are small, perhaps because we used only a small subset of generative factors to compute them. The disentanglement score is better for smaller d, which shows that information is packed more efficiently in this case. Performance of the models that use images strongly depends on the quality of the images themselves. For example, if the horizon and runway are clearly visible (as in Table 5a ), models with images outperform models that do not use them (Table 5 ), while on night-time noisy images (Table 5b ) models that use images lose.

6. CONCLUSION:

In this work, we have presented the dataset for learning state representation and a unified evaluation framework to measure the quality of multimodal joint representations produced by different encoders, through both secondary supervised learning tasks and quantitative metrics of disentanglement and informativeness quality. This work represents an attempt to help advance the research in model-based reinforcement learning and state representation learning by providing a unified. At the same time, the large-scale comparison between baseline and state of the art models that we perform in this work has produced some interesting results by itself. In general, we believe that this dataset can become the new standard for evaluation in representation learning for complex systems, and we hope it will inspire many novel techniques to advance the state of the art that we have attempted to establish and quantify in our practical evaluation. 

:



Figures and Tables:

Figure 1: Attention cell architecture.

Figure 2: Three main ideas used to construct baseline models: (a) autoencoder for the state vectors, (b) context prediction model, (c) autoregression (forward) model.

Figure 3: The Multimodal Temporal Encoder.Figure4: The Dynamics Model.

Figure 5: Evaluation pipeline.

Groups of sensors represented in the X-Plane dataset.

Table content: <table><tr><td>Group</td><td>Description</td></tr></table>



Table content: <table/>

Sample images from the dataset.

Table content: <table/>

Evaluation results; left to right: model name and description, whether it uses images, latent dimension d, mean inference time t (s), MSE for regression tasks normalized to the MEAN baseline (lower is better), QEDR scores (Section 5.1).

Table content: <table><tr><td>Model</td><td/><td colspan="2">Img d</td><td>MSE for regression tasks</td><td colspan="2">QEDR Scores</td></tr><tr><td/><td/><td/><td colspan="4">Wind Turb. Land Auto Fail Overall Inform. Disent. Compl.</td></tr><tr><td colspan="2">MEAN Baseline</td><td/><td colspan="2">-0.969 1.036 2.936 2.039 0.961 7.940</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="2">LSTM Autoencoder</td><td/><td colspan="2">32 1.014 0.985 0.976 0.953 0.788 0.953</td><td>0.912</td><td>0.095</td><td>0.107</td></tr><tr><td colspan="2">1 layer encoder,</td><td/><td colspan="2">64 1.017 0.970 0.986 0.795 0.763 0.912</td><td>0.893</td><td>0.099</td><td>0.126</td></tr><tr><td colspan="2">1 layer decoder</td><td/><td colspan="2">128 1.096 0.954 0.934 0.786 0.724 0.893</td><td>0.885</td><td>0.067</td><td>0.121</td></tr><tr><td/><td/><td/><td colspan="2">256 1.076 1.012 8.770 0.777 0.941 3.820</td><td>0.897</td><td>0.054</td><td>0.112</td></tr><tr><td colspan="2">LSTM Autoencoder</td><td/><td colspan="2">32 1.013 0.956 0.925 0.933 0.949 0.945</td><td>0.892</td><td>0.135</td><td>0.153</td></tr><tr><td colspan="2">2 layers bidir. encoder</td><td/><td colspan="2">64 1.024 0.967 1.075 0.807 0.783 0.951</td><td>0.892</td><td>0.089</td><td>0.131</td></tr><tr><td colspan="2">2 layers bidir. decoder</td><td/><td colspan="2">128 1.062 0.978 0.958 0.832 0.797 0.921</td><td>0.889</td><td>0.070</td><td>0.124</td></tr><tr><td/><td/><td/><td colspan="2">256 1.521 4.236 0.958 0.846 1.797 1.527</td><td>0.941</td><td>0.059</td><td>0.122</td></tr><tr><td colspan="2">Conv1D AutoEncoder</td><td/><td colspan="2">32 1.008 0.994 1.068 0.893 0.992 0.997</td><td>0.924</td><td>0.086</td><td>0.111</td></tr><tr><td colspan="2">6-layer 1D convolutional</td><td/><td colspan="2">64 1.002 0.982 0.960 0.862 0.894 0.935</td><td>0.924</td><td>0.078</td><td>0.126</td></tr><tr><td colspan="2">autoencoder with</td><td/><td colspan="2">128 1.002 0.975 0.979 0.867 0.901 0.943</td><td>0.912</td><td>0.070</td><td>0.128</td></tr><tr><td colspan="2">kernel size 7 × 7</td><td/><td colspan="2">256 1.000 0.978 1.009 0.893 0.955 0.968</td><td>0.922</td><td>0.059</td><td>0.123</td></tr><tr><td colspan="2">PCA Autoencoder</td><td/><td colspan="2">32 0.999 0.961 0.901 0.862 0.861 0.906</td><td>0.921</td><td>0.084</td><td>0.111</td></tr><tr><td colspan="2">1 layer decoder</td><td/><td colspan="2">64 1.000 0.997 1.000 0.997 0.999 0.999</td><td>0.917</td><td>0.075</td><td>0.119</td></tr><tr><td colspan="2">256 × 256 × 3 images</td><td/><td colspan="2">128 1.000 1.000 0.999 0.997 0.998 0.999</td><td>0.917</td><td>0.064</td><td>0.119</td></tr><tr><td colspan="2">mean features over time series</td><td/><td colspan="2">256 1.000 0.998 1.000 0.999 0.999 0.999</td><td>1.000</td><td>0.067</td><td>0.132</td></tr><tr><td colspan="2">ResNet34 Autoencoder</td><td/><td colspan="2">32 0.988 0.960 0.853 0.927 0.942 0.913</td><td>0.924</td><td>0.137</td><td>0.157</td></tr><tr><td colspan="3">ResNet34 encoder with 1 FC layer</td><td colspan="2">64 1.006 0.958 0.898 0.913 0.916 0.925</td><td>0.931</td><td>0.080</td><td>0.136</td></tr><tr><td colspan="2">15 layer conv2d decoder</td><td/><td colspan="2">128 0.994 0.972 0.910 0.908 0.927 0.930</td><td>0.940</td><td>0.068</td><td>0.127</td></tr><tr><td colspan="2">images only</td><td/><td colspan="2">256 0.992 0.963 0.901 0.916 0.935 0.928</td><td>0.978</td><td>0.071</td><td>0.130</td></tr><tr><td colspan="2">TS Regression LSTM</td><td/><td colspan="2">32 1.010 0.950 1.037 0.816 0.948 0.955</td><td>0.925</td><td>0.081</td><td>0.108</td></tr><tr><td colspan="2">2-layer bidir. LSTM</td><td/><td colspan="2">64 1.002 1.211 1.122 1.436 5.566 1.738</td><td>0.920</td><td>0.066</td><td>0.106</td></tr><tr><td colspan="2">trained to predict</td><td/><td colspan="2">128 1.035 1.119 0.872 0.886 1.558 1.011</td><td>0.921</td><td>0.062</td><td>0.116</td></tr><tr><td colspan="2">state at t + 30</td><td/><td colspan="2">256 1.074 1.084 0.952 0.888 1.034 0.978</td><td>0.925</td><td>0.058</td><td>0.109</td></tr><tr><td colspan="2">TS Regression Attention</td><td/><td colspan="2">32 1.020 0.985 0.999 0.940 0.948 0.978</td><td>0.920</td><td>0.084</td><td>0.115</td></tr><tr><td colspan="2">3-layer attention encoder</td><td/><td colspan="2">64 1.004 1.036 1.039 1.464 1.005 1.139</td><td>0.923</td><td>0.072</td><td>0.115</td></tr><tr><td colspan="2">trained to predict</td><td/><td colspan="2">128 1.136 0.962 0.869 1.241 0.788 0.999</td><td>0.922</td><td>0.059</td><td>0.108</td></tr><tr><td colspan="2">state at t + 1</td><td/><td colspan="2">256 1.039 1.174 0.991 1.363 1.887 1.225</td><td>0.923</td><td>0.054</td><td>0.114</td></tr><tr><td colspan="3">TS Regr. ConvLSTM+LSTM</td><td colspan="2">32 1.028 0.969 0.961 0.846 0.870 0.930</td><td>0.922</td><td>0.083</td><td>0.115</td></tr><tr><td colspan="3">5-layer conv. LSTM with 3 × 3 kernel</td><td colspan="2">64 1.034 1.004 0.949 0.813 0.929 0.929</td><td>0.923</td><td>0.077</td><td>0.116</td></tr><tr><td colspan="3">for images and 1-layer LSTM for sensors,</td><td colspan="2">128 1.007 1.009 1.009 1.043 1.070 1.025</td><td>0.922</td><td>0.064</td><td>0.117</td></tr><tr><td colspan="3">trained to predict state at t + 30</td><td colspan="2">256 1.018 0.977 0.970 0.791 0.857 0.917</td><td>0.915</td><td>0.055</td><td>0.114</td></tr><tr><td colspan="2">TS Regr. ConvLSTM</td><td/><td colspan="2">32 1.000 0.998 0.999 0.996 0.999 0.998</td><td>0.917</td><td>0.084</td><td>0.109</td></tr><tr><td colspan="2">5-layer conv. LSTM with</td><td/><td colspan="2">64 1.000 0.999 0.999 0.996 0.999 0.998</td><td>0.921</td><td>0.072</td><td>0.113</td></tr><tr><td colspan="2">3 × 3 kernel for images,</td><td/><td colspan="2">128 1.000 0.998 1.000 0.999 0.998 0.999</td><td>0.921</td><td>0.062</td><td>0.112</td></tr><tr><td colspan="3">trained to predict state at t + 30</td><td colspan="2">256 1.000 0.998 1.000 1.000 0.999 1.000</td><td>0.923</td><td>0.059</td><td>0.120</td></tr><tr><td colspan="2">Context LSTM Regressor</td><td/><td colspan="2">32 1.011 1.079 1.206 0.931 1.005 1.071</td><td>0.923</td><td>0.089</td><td>0.122</td></tr><tr><td colspan="2">2-layer bidir. LSTM</td><td/><td colspan="2">64 1.003 0.993 1.023 0.834 0.855 0.948</td><td>0.924</td><td>0.073</td><td>0.116</td></tr><tr><td colspan="2">trained to predict states</td><td/><td colspan="2">128 1.271 1.059 1.262 1.045 0.906 1.138</td><td>0.922</td><td>0.059</td><td>0.110</td></tr><tr><td/><td/><td/><td/><td>170</td><td>0.916</td><td>0.052</td><td>0.110</td></tr><tr><td colspan="2">Context Attention Regressor</td><td/><td colspan="2">32 1.148 0.950 0.931 1.103 0.820 0.991</td><td>0.924</td><td>0.079</td><td>0.107</td></tr><tr><td colspan="2">3-layer attention net,</td><td/><td colspan="2">64 1.041 0.922 0.868 1.012 0.691 0.912</td><td>0.921</td><td>0.069</td><td>0.109</td></tr><tr><td colspan="2">trained to predict states</td><td/><td colspan="2">128 1.021 1.026 1.243 1.153 0.762 1.106</td><td>0.922</td><td>0.059</td><td>0.110</td></tr><tr><td/><td/><td/><td/><td>086</td><td>0.924</td><td>0.053</td><td>0.112</td></tr><tr><td colspan="2">Context Conv1D Regressor</td><td/><td colspan="2">32 1.001 0.979 0.977 0.858 0.893 0.940</td><td>0.921</td><td>0.087</td><td>0.121</td></tr><tr><td colspan="2">6-layer 1D convolutional</td><td/><td colspan="2">64 1.002 0.978 1.000 0.842 0.863 0.940</td><td>0.923</td><td>0.071</td><td>0.116</td></tr><tr><td colspan="3">network trained to predict states</td><td colspan="2">128 1.002 0.970 0.978 0.829 0.883 0.930</td><td>0.923</td><td>0.062</td><td>0.115</td></tr><tr><td/><td/><td/><td/><td>934</td><td>0.923</td><td>0.059</td><td>0.122</td></tr><tr><td colspan="3">Multimodal Temporal Encoder</td><td colspan="2">32 1.002 0.976 0.999 0.814 0.913 0.938</td><td>0.923</td><td>0.088</td><td>0.115</td></tr><tr><td colspan="2">LSTM with shared weights</td><td/><td colspan="2">64 1.001 0.962 0.999 0.804 0.884 0.931</td><td>0.922</td><td>0.067</td><td>0.105</td></tr><tr><td colspan="2">between inputs,</td><td/><td colspan="2">128 1.010 0.982 1.000 0.816 0.903 0.940</td><td>0.920</td><td>0.060</td><td>0.108</td></tr><tr><td>ResNet18</td><td/><td/><td colspan="2">256 1.003 0.988 0.998 0.847 0.949 0.953</td><td>0.920</td><td>0.055</td><td>0.114</td></tr><tr><td colspan="2">The Dynamics Model (X)</td><td/><td colspan="2">32 1.000 0.998 0.999 0.950 0.981 0.984</td><td>0.922</td><td>0.095</td><td>0.124</td></tr><tr><td colspan="3">Decouple dynamics with ResNet18,</td><td colspan="2">64 1.000 0.984 1.000 0.951 0.968 0.981</td><td>0.910</td><td>0.073</td><td>0.115</td></tr><tr><td colspan="2">embedding is the mean of</td><td/><td colspan="2">128 1.001 0.971 0.925 0.890 0.836 0.920</td><td>0.923</td><td>0.063</td><td>0.114</td></tr><tr><td colspan="2">embs. over time series</td><td/><td colspan="2">256 1.000 1.000 1.000 1.000 0.998 1.000</td><td>0.923</td><td>0.061</td><td>0.127</td></tr><tr><td colspan="2">The Dynamics Model (Img)</td><td/><td colspan="2">32 1.000 0.998 0.999 0.929 0.998 0.981</td><td>0.924</td><td>0.097</td><td>0.126</td></tr><tr><td colspan="3">Decouple dynamics with ResNet18,</td><td colspan="2">64 1.000 0.998 0.999 0.996 0.998 0.998</td><td>0.923</td><td>0.074</td><td>0.119</td></tr><tr><td colspan="2">embedding is the mean of</td><td/><td colspan="2">128 1.000 0.999 1.000 0.986 0.999 0.996</td><td>0.924</td><td>0.068</td><td>0.123</td></tr><tr><td colspan="2">embs. over time series</td><td/><td colspan="2">256 1.000 0.998 0.949 0.897 0.881 0.940</td><td>0.925</td><td>0.060</td><td>0.124</td></tr><tr><td colspan="3">The Dynamics Model (Img+X)</td><td colspan="2">32 1.000 0.998 0.998 0.996 0.998 0.998</td><td>0.926</td><td>0.082</td><td>0.108</td></tr><tr><td colspan="3">Decouple dynamics with ResNet18,</td><td colspan="2">64 1.000 0.998 1.000 0.998 0.999 0.999</td><td>0.923</td><td>0.075</td><td>0.122</td></tr><tr><td colspan="2">embedding is the mean of</td><td/><td colspan="2">128 1.000 0.998 1.000 0.929 0.998 0.981</td><td>0.921</td><td>0.067</td><td>0.126</td></tr><tr><td colspan="2">embs. over time series</td><td/><td colspan="2">256 1.000 0.999 1.000 0.943 0.998 0.985</td><td>0.912</td><td>0.060</td><td>0.126</td></tr><tr><td>(a)</td><td/><td>(b)</td><td/><td>(c)</td><td/></tr><tr><td>(a)</td><td>0.163</td><td>0.123</td><td>0.207</td><td>4.488</td><td/><td>0.121</td></tr><tr><td>(b)</td><td>1.698</td><td>1.812</td><td>1.274</td><td>1.060</td><td/><td>1.805</td></tr><tr><td>(c)</td><td>1.106</td><td>1.080</td><td>0.978</td><td>1.187</td><td/><td>1.076</td></tr></table>



Table content: <table/>

Peer review (iclr2019):
 Overview and contributions: The authors present a newly collected dataset and evaluation framework for learning representations for landing an airplane. The dataset is collected from the X-Plane simulation environment and consists of 8011 landings, each landing consists of time series data from 1090 sensors. Their evaluation metric is a combination of disentanglement score, regression tasks, and failure classification. The authors test a combination of baseline models from basic autoencoders to dynamic actions-aware encoders. The writing is generally clear but I have doubts about the proposed evaluation metrics, experiments, and significance of the dataset (details below). Strengths: 1. The task seems to be novel and complex. The authors have done a good job of collecting the dataset and ensuring that the data is clean and comprehensive. 2. The authors have performed a comprehensive job of evaluating the combinations of baseline models for their proposed task. Weaknesses: 1. Table 4 on evaluation results, while comprehensive, lacks some explanation. The issue with MSE is that it is hard to interpret what these values mean. Specifically, how difficult is this task? How well can a human perform on this task? How well are the baselines doing relative to human-level performance, and is there room for improvement? The answers to these questions are important towards whether this new dataset will be a strong benchmark for representation learning. 2. There is less novelty in terms of the models presented for evaluation since they are composed of existing models. What are some state-of-the-art models for similar tasks, and do they constitute fair comparison? Questions to authors: 1. Refer to weakness points 1 and 2. 2. What biases do you think might exist in the dataset during the collection process? How might these biases affect what the models learn, and how can they be mitigated? 3. How do you ensure that all sensors are active at all times and that all sensors provide useful information for predicting the label? Are there cases where the multisensor data is noisy in certain modalities or missing in other modalities? If so, what are some models that can remain robustness to noisy or missing modalities? 4. Why do you think disentangled representations will help? Sure, they have been generally shown to help learn more interpretable representations, and help in flexible generation from disentangled factors. But in terms of discriminative or generative performance on your newly proposed dataset, does learning disentangled representations help? What are some models that can learn effectively learn such disentangled representations? Presentation improvements, typos, edits, style, missing references: Section 3, line 7, *with with a frequency* -> *with* I would suggest referring to some recent work on multimodal temporal fusion, such as "Memory Fusion Network for Multi-view Sequential Learning. Amir Zadeh, Paul Pu Liang, Navonil Mazumder, Soujanya Poria, Erik Cambria, Louis-Philippe Morency, AAAI 2018"