Understanding Knowledge Distillation in Non-autoregressive Machine Translation

Chunting Zhou, Graham Neubig, Jiatao Gu

Introduction

Traditional neural machine translation (NMT) systems (Bahdanau et al., 2015; Gehring et al., 2017; Vaswani et al., 2017) generate sequences in an autoregressive fashion; each target token is predicted step-by-step by conditioning on the previous generated tokens in a monotonic (e.g. left-to-right) order. While such autoregressive translation (AT) models have proven successful, the sequential dependence of decisions precludes taking full advantage of parallelism afforded by modern hardware (e.g. GPUs) at inference time. In contrast, non-autoregressive translation (NAT) models (Gu et al., 2018; Lee et al., 2018) predict the whole sequence or multi-token chunks of the sequence simultaneously, alleviating this problem by trading the model’s capacity for decoding efficiency. Such a non-autoregressive factorization assumes that the output tokens are independent from each other. However, this assumption obviously does not hold in reality and as a result NAT models generally perform worse than standard AT models.

One key ingredient in the training recipe for NAT models that is used in almost all existing works (Gu et al. (2018); Lee et al. (2018); Stern et al. (2019), inter alia) is creation of training data through knowledge distillation (Hinton et al., 2015). More precisely, sequence-level knowledge distillation (Kim & Rush, 2016) – a special variant of the original approach – is applied during NAT model training by replacing the target side of training samples with the outputs from a pre-trained AT model trained on the same corpus with a roughly equal number of parameters. It is usually assumed (Gu et al., 2018) that knowledge distillation’s reduction of the “modes” (alternative translations for an input) in the training data is the key reason why distillation benefits NAT training. However, this intuition has not been rigorously tested, leading to three important open questions:

Exactly how does distillation reduce the “modes”, and how we could we measure this reduction quantitatively? Why does this reduction consistently improve NAT models?

What is the relationship between the NAT model (student) and the AT model (teacher)? Are different varieties of distilled data better for different NAT models?

Due to distillation, the performance of NAT models is largely bounded by the choice of AT teacher. Is there a way to further close the performance gap with standard AT models?

In this paper, we aim to answer the three questions above, improving understanding of knowledge distillation through empirical analysis over a variety of AT and NAT models. Specifically, our contributions are as follows:

We first visualize explicitly on a synthetic dataset how modes are reduced by distillation (section 3.1). Inspired by the synthetic experiments, we further propose metrics for measuring complexity and faithfulness for a given training set. Specifically, our metrics are the conditional entropy and KL-divergence of word translation based on an external alignment tool, and we show that these metrics are correlated with NAT model performance (section 3.2).

We conduct a systematic analysis (section 4) over four AT teacher models and six NAT student models with various architectures on the standard WMT14 English-German translation benchmark. These experiments find a strong correlation between the capacity of an NAT model and the optimal dataset complexity that results in the best translation quality.

Inspired by these observations, we propose approaches to further adjust the complexity of the distilled data in order to match the model’s capacity (section 5). We also show that we can achieve the state-of-the-art performance for NAT models and largely match the performance of the AT model.

Background

In order to model the joint probability of the output sequence y{\bm{y}}, NMT models usually generate each output token conditioned on the previously generated ones p(y∣x)=∏t=1Tp(yt∣y<t,x)p({\bm{y}}|{\bm{x}})=\prod_{t=1}^{T}p(y_{t}|{\bm{y}}_{<t},{\bm{x}}). This is known as the autoregressive factorization. To generate a translation from this model, one could predict one token at a time from left to right and greedily take arg max⁡\operatorname*{arg\,max} over each output probability distribution, or use beam search to consider a fixed number of hypotheses. In this work, we study non-autoregressive translation (NAT), a special subset of NMT models with an additional restriction (the zeroth-order Markov assumption) upon the output predictions or a subset thereof. The simplest formulation of an NAT model independently factors the conditional distribution: p(y∣x)=∏t=1Tp(yt∣x)p({\bm{y}}|{\bm{x}})=\prod_{t=1}^{T}p(y_{t}|{\bm{x}}).

Standard NAT models (Gu et al., 2018) adopt an architecture similar to the Transformer (Vaswani et al., 2017) and make non-autoregressive predictions for the entire sequence with one forward pass of the decoder. However, because multiple translations are possible for a single input sentence (the so-called multi-modality problem; Gu et al. (2018)), vanilla NAT models can fail to capture the dependencies between output tokens. As a result, they tend to make egregious mistakes such as outputting tokens repeatedly. To improve the model’s ability to handle multi-modality, recent works have incorporated approaches including (1) relaxing the fully non-autoregressive restriction and adopting KK decoding passes (instead of just one) to iteratively refine the generated outputs (Lee et al., 2018; Ghazvininejad et al., 2019; Wang et al., 2018; Stern et al., 2018; 2019; Gu et al., 2019); (2) using latent variables (Kaiser et al., 2018; Ma et al., 2019; Shu et al., 2019) or structured information such as syntax trees (Akoury et al., 2019) to capture translation variation; (3) training NAT models with objectives other than maximum likelihood (Wang et al., 2019; Wei et al., 2019; Shao et al., 2019) which ameliorates the effects of multi-modality. However, to achieve competitive performance with the autoregressive model, almost all existing NAT models rely on training using data distilled from a pre-trained AT model instead of the real parallel training set, as described below.

2 Sequence-level Knowledge Distillation

Knowledge distillation (Liang et al., 2008; Hinton et al., 2015) was originally proposed for training a weaker student classifier on the targets predicted from a stronger teacher model. A typical approach is using the label probabilities produced by the teacher as “soft targets” q_{i}={\exp(z_{i}/\tau)}{\big{/}}{\sum_{j}\exp(z_{j}/\tau)} for training the student model, where qiq_{i} and ziz_{i} are the probability and the logit of class ii respectively and τ\tau is the temperature. Prior work has shown the effectiveness of adopting knowledge distillation in adversarial defense (Papernot et al., 2016), neural network compression (Howard et al., 2017), and fast inference for speech synthesis (Oord et al., 2018).

where t∈T{\bm{t}}\in\mathcal{T} is the space of possible target sequences. This can also be seen as a special case of standard distillation over the sentence space when the temperature τ\tau approaches , which is equivalent to taking the arg max⁡\operatorname*{arg\,max} over all feasible translations. While the “hard target” y^\hat{{\bm{y}}} is the most likely translation predicted by the teacher, in practice we use beam search as an approximation. As mentioned earlier, almost all the existing literature trains NAT models using sequence-level knowledge distillation from a pre-trained AT model to achieve competitive performance. Particularly, it is common to train the teacher model as a standard autoregressive Transformer (Vaswani et al., 2017) with a roughly equal number of trainable parameters as the desired NAT model on the real data. Next, we will first study how this knowledge distillation process affects the behavior of NAT models.

How does Distillation Improve NAT?

In this section, we start from an introductory example to illustrate how NAT models fail to capture the multi-modality of data. Then we propose a metric to assess the multi-modality of a data set and use it to test our hypothesis about how knowledge distillation affects NAT models.

Dataset. We start by investigating NAT’s difficulties in modeling multi-modality in output data using a synthetic setup where we explicitly include multiple modes in the training data. More specifically, we utilize three language pairs – English-German (En-De), English-French (En-Fr), and English-Spanish (En-Es) – from the Europarl parallel corpus.https://www.statmt.org/europarl/ We extract sentences that have aligned sentences for all languages, and create a multi-target En-De/Es/Fr corpus. In this case every English input sentence always corresponds to target sentences in three different languages, which forms three explicit output modes. Notably, this is similar to the one-to-many translation setting in Johnson et al. (2017) but in our case we do not have an explicit signal (e.g. target language tag) to tell the NMT model which target language to translate to.

Models. We train both the AT and NAT models on this concatenated data set, then compare the distributions of translations with each other. We use the standard Transformer(base) model (Vaswani et al., 2017) as the AT model, and a simplified version of Gu et al. (2018) as the NAT model where the decoder’s inputs are monotonically copied from the encoder embeddings and a length predictor is learned to predict the target sentence length. Both models are trained for 300,000300,000 steps using maximum likelihood. After training, we use both models to translate the English sentences in the validation and test sets.

where lil_{i} denotes the language class ii, and p(yt∣li)p(y_{t}|l_{i}) is the token frequency of yty_{t} in language lil_{i}. We assume p(li)p(l_{i}) follows a uniform distribution. As shown in Fig. 1 (a), points of the AT outputs are clustered closely to each vertex of the simplex, indicating that the AT model prefers to generate the whole sequence in one language. This phenomenon verifies our assumption that decoding with the AT model (distillation) is essentially selecting “modes” over the real data.

Visualization of NAT Outputs. We visualize outputs for the NAT model trained on the same data in Fig. 1 (b). In contrast to the AT results, the NAT points are scattered broadly inside the simplex, indicating that the NAT model fails to capture the mode of language types. Instead, it predicts tokens mixed with multiple languages, which corroborates our hypothesis that the NAT model has trouble consistently selecting a single mode when multiple modes exist.

Next, we create two datasets that have fewer modes than the original dataset. First, we randomly select a single target sentence from one of the three languages for each source sentence. Second, we perform distillation, decoding from the AT model trained on the combined training set. As noted in the AT results, distillation will also roughly be selecting a language mode, but we conjecture that this selection may be more systematic, selecting a particular language for a particular type of training sentence. As shown in Fig. 1(c) (d), NAT models trained on both of these datasets are more likely to choose one mode (language) when generating translations, showing that training with reduced modes is essential for NAT model. Furthermore, points in Fig. 1 (d) are clearly clustered better than (c) indicating that modes selected by AT models are indeed likely more systematic and easy to capture than those generated by randomly assigning a language for each sentence.

2 Quantitative Measures for Parallel Data

To better study why distillation is crucial for NAT models, in this section, we propose quantitative measures for analyzing the complexity and faithfulness of parallel data, two properties that we hypothesize are important for NAT training.

Measure of Complexity. Inspired by the observations in the synthetic experiments, we propose to use a measure of translation uncertainty, specifically operationalized as conditional entropy, as the measurement of complexity C(d)C(d) for any given dataset d={(x1,y1),...,(xN,yN)}d=\{({\bm{x}}_{1},{\bm{y}}_{1}),...,({\bm{x}}_{N},{\bm{y}}_{N})\}, where (x,y)({\bm{x}},{\bm{y}}) is sentence pair instantiation of (X,Y)(\mathbf{X},\mathbf{Y}) and X∈X,Y∈Y\mathbf{X}\in\mathcal{X},\mathbf{Y}\in\mathcal{Y}:

where we use xx and yy to denote a word in the source and target vocabulary respectively. TxT_{x} and TyT_{y} denote the length of the source and target sentences. To make the computation tractable, we make two additional assumptions on the conditional distribution p(y∣x)p({\bm{y}}|{\bm{x}}):

Assumption 1: We assume the target tokens are independent given the source sentence. Then the conditional entropy of a sentence can be converted into the sum of entropy of target words conditioned on the source sentence x{\bm{x}}.

The corpus level complexity C(d)C(d) is then calculated by adding up the conditional entropy H(Y∣X=x)\mathcal{H}(\mathbf{Y}|\mathbf{X}={\bm{x}}) of all sentences. To prevent C(d)C(d) from being dominated by frequent words, we calculate C(d)C(d) by averaging the entropy of target words conditioned on a source word, denoted C(d)=1∣Vx∣∑x∈VxH(y∣x)C(d)=\frac{1}{|\mathcal{V}_{x}|}\sum_{x\in\mathcal{V}_{x}}\mathcal{H}(y|x).

To illustrate that the proposed metric is a reasonable measure of complexity of a parallel corpus, in Tab. 1 we compute C(d)C(d) for parallel data from different language pairs, the concatenated data set, and the data distilled from the AT model described in section 3.1. We observe that the conditional entropy of the distilled data is much smaller than that of the concatenated or randomly selected data mentioned above. Additionally, we find that the conditional entropy of En-Es and En-Fr are similar but that of En-De is relatively larger, which can also explain why the student NAT model prefers to predict the modes of Es or Fr more often than De as shown in Fig. 1(d).

Measure of Faithfulness. C(d)C(d) reflects the level of multi-modality of a parallel corpus, and we have shown that a simpler data set is favorable to an NAT model. However, it is not fair to assess the data set only by its complexity; we can trivially construct a simple data set with no variations in the output, which obviously won’t be useful for training. The other important measurement of the data set is its faithfulness to the real data distribution. To measure the faithfulness of a parallel corpus dd, we use KL-divergence of the alignment distribution between the real parallel data set rr and an altered parallel data set dd, denoted F(d)F(d):

Empirical Study

In this section, we perform an extensive study over a variety of non-autoregressive (NAT) models trained from different autoregressive (AT) teacher models to assess how knowledge distillation affects the performance of NAT models.

Data. We use the data set commonly used by prior work as our evaluation benchmark: WMT14 English-German (En-De)http://www.statmt.org/wmt14/translation-task.html. We use newstest2013 as the validation set for selecting the best model, and newstest2014 as the test set. We learn a byte-pair encoding (BPE, Sennrich et al., 2016) vocabulary of 37,000 on the tokenized data.

AT Models. We set up four Transformer models with different parameter sizes: Transformer-tiny/small/base/big denoted as tiny, small, base, big respectively. We build base and big models following settings described in Vaswani et al. (2017), and reduce the model sizes for tiny, small to create weaker teacher models. Details of the model architectures can be found in Appendix A.

All the models are trained using the Adam optimizer (Kingma & Ba, 2014) with the maximum number of steps set to 300,000300,000. After training, we use the resulting AT models to decode the whole training set with beam size 55 and replace the real target sentences to create a new parallel corpus.

NAT Models. We consider the following NAT models, from vanilla to state-of-the-art. All the models are using the Transformer as the basic backbone and are (re-)implemented based on Fairseqhttps://github.com/pytorch/fairseq except for FlowSeq. We briefly outline the methods and parameters here, and describe detailed settings in the Appendix A.

Vanilla NAT (Gu et al., 2018): Similarly to section 3.1, we use a simplified version where the decoder’s inputs are directly copied from the encoder without considering latent variables.

FlowSeq (Ma et al., 2019): FlowSeq adopts normalizing flows (Kingma & Dhariwal, 2018) as the latent variables to model the mappings from source sentences to a latent space.

NAT with Iterative Refinement (iNAT, Lee et al., 2018): iNAT extends the vanilla NAT by iteratively reading and refining the translation. The number of iterations is set to 10 for decoding.

Insertion Transformer (InsT, Stern et al., 2019): InsT adopts a similar architecture as iNAT while generating the sequence by parallel insertion operations. Here, we only consider InsT trained with uniform loss as described in the original paper.

MaskPredict (MaskT, Ghazvininejad et al., 2019): MaskT adopts a masked language model (Devlin et al., 2018) to progressively generate the sequence from an entirely masked input. The number of iterations is set to be 10.

Levenshtein Transformer (LevT, Gu et al., 2019): LevT uses similar architectures as in InsT and MaskT while generating based on both insertion and deletion operations. We experiment with a base and big LevT model (LevT and LevT-big in Tab. 2).

We also summarize the parameter size, performance and relative decoding speed of the NAT models introduced in Tab. 2. We use the decoding time of vanilla NAT to represent one unit of time, and Iters ×\times Pass represents the relative time units used for each model.

As mentioned earlier, we analyze each model by training from both the real and 44 distilled targets. We train the NAT models for the same number of steps as the AT models. For a fair comparison of the actual ability of each NAT-based model, we test all the models based on greedy decoding without any advanced search algorithms (e.g. length beam (Ghazvininejad et al., 2019), noisy parallel decoding (Ma et al., 2019), or re-ranking from the teacher model (Gu et al., 2018)). Notably, the vanilla NAT and FlowSeq output translations with single forward pass, while the remaining models are based on the iterative refinement.

2 Analysis of the Distilled Data

We compare different dimensions of the data generated by the four AT models and the real data set in Fig. 3. First, Fig. 3 (a) shows that as the capacity of the AT model increases, the complexity C(d)C(d) of the distilled data increases, which indicates that the multi-modality increases as well. At the same time, we observe that F(d)F(d) defined in section 3.2 also decreases, showing that the distilled data more faithfully represents the word-level translation distribution of the original data. Second, we plot the BLEU score of the distilled data w.r.t to the real data set in (b) and we observe that the BLEU score of the distilled data from a higher-capacity teacher model is higher, which is both intuitive and in agreement with the results on KL divergence.

We also investigate how the relative ordering of words in the source and target sentences is changed during distillation. We use the fuzzy reordering score proposed in Talbot et al. (2011). A larger fuzzy reordering score indicates the more monotonic alignments. As shown in Fig 3 (c), the distilled data has significantly less reordering compared to the real parallel sentences, and the distilled data from a weaker AT teacher is more monotonic than a stronger AT teacher. We also show a randomly sampled example in Fig. 2 where compared to the real translation, the AT distilled target is much more monotonically aligned to the source sentence. This has potential benefits in that these simpler reordering patterns may be easier to learn for NAT models, but also disadvantages in that it may prevent NAT models from learning complex reordering patterns.

3 Analysis of Distillation Strategies

In section 4.2, we have shown that decoding with an AT model reduces the conditional entropy of the parallel data set, which mitigates multi-modality in the output data. But does the decoding method of the AT model affect this change in the data set? We also investigate different decoding strategies when creating distilled data, using the base Transformer model as the teacher and the vanilla NAT model as the student. In Tab. 3, four decoding methods are presented: sampling, sampling within the top-10 candidates, beam search, and greedy decoding. With the same AT model, the performance of the NAT model differs widely depending on the decoding approach, where distillation with beam search results in the best performance.

We can see that beam search or greedy decoding can reduce the complexity of the real data the most while maintaining high faithfulness. In contrast, sampling based decoding methods less aggressively reduce the modes in the output sequence. This finding is in concert with Ott et al. (2018), who demonstrate that because beam search approximately selects the most probable translation, it effectively reduces diversity in the output translations compared to sampling or the true distribution.

4 Distilled Data v.s. NAT Models

We next examine the relationship between the NAT students and distilled training data from different AT models. In Fig. 4, we demonstrate results for the NAT models listed in section 4.1. We use the test set performance on real data as a simple metric to measure the capacity of the NAT model and arrange the subfigures in an increasing order of the performance (left-to-right, top-to-bottom). The results in the figure demonstrate that, interestingly, weaker NAT students prefer distilled data with smaller complexity as measured above in section 4.2. The best performance of NAT models – from lower capacity ones to higher capacity ones – is achieved with distilled data of lower complexity to higher complexity, i.e. the vanilla NAT model performs best when using the distilled data from a small Transformer whereas LevT achieves the best performance when training with the distilled data from a big Transformer. Third, and notably, by simply changing the distilled data set upon which the models are trained, we are able to significantly improve the state-of-the-art results for models in a particular class. For example, FlowSeq increased to 2222, by simply changing from the distilled data of Transformer(base) to Transformer(small). Finally, we find that by distilling from a big AT model, LevT is able to close the gap with the Transformer (base) with a similar number of parameters. Both LevT and LevT-big achieve the state-of-the-art performance for NAT-based models.

Improvements to Knowledge Distillation

The previous section shows that the optimal complexity of the dataset is highly correlated with the capacity of the NAT model. In this section, we introduce three techniques that can be used to alter the distilled data to match the capacity of NAT model. Specifically, these techniques can be used to simplify the data further (BANs, MoE) for a lower-capacity student model or increase faithfulness of the data set (Interpolation) for a higher-capacity student model.

Born-Again Networks. We apply Born-Again neworks (BANs) to create a simplified dataset for NAT models. BANs were originally proposed as a self-distillation technique (Furlanello et al., 2018) that uses the output distribution of a trained model to train the original model. Starting from the real data, we repeatedly train new AT models with decoded sentences from the AT model at the previous iteration. This process is repeated for kk times and yields kk distilled data sets, upon which we perform NAT training and examine how the kk born-again teachers affect the performance of NAT students.

We conduct experiments using the vanilla NAT model (Gu et al., 2018) (which achieved the best performance with distilled data from a small Transformer in section 4.4) and the base Transformer as the AT model. As shown in Fig. 5, we can make the following observations: (i) The performance of the base AT model almost remains unchanged during the reborn iterations. (ii) The performance of the vanilla NAT model can be improved by 2 BLEU when using the distilled data from reborn iteration 6. (iii) As the reborn iterations continue, the complexity of the distilled data decreases and becomes constant eventually. Meanwhile, the quality of the distilled data compared to the real data decreases.

Mixture-of-Experts. The mixture-of-expert model (MoE; Shen et al. (2019)) learns different experts for diverse machine translation, and different mixture components were shown to capture consistent translation styles across examples. Inspired by this, we use one expert from the mixture model to translate the training data, which is supposed to generate a single style of translation and reduce the diversity in the original data set. Then we use the best single-expert translations as the distilled data to train the vanilla NAT model. Specifically, we follow Shen et al. (2019)’s setup, using the base Transformer model and uniform hard mixture model, varying the number of experts.

In Fig. 6, we observe that the performance of the best expert of MoE tends to decrease as the number of experts increases. However, the complexity (C(d)C(d)) and faithfulness (F(D)F(D)) of distilled data from different MoE models has a relatively large variance. Compared to using the distilled data from a plain base AT model, the performance of NAT model is improved by 1.21 BLEU when using the distilled data from the MoE model with the number of experts of 3 which produces the distilled data with the least complexity.

Sequence-Level Interpolation. section 4.4 shows stronger NAT models (e.g. MaskT, LevT) have the ability to learn from the dataset that is closer to the real data, and achieve better performance. We adopt the sequence-level interpolation proposed in Kim & Rush (2016) as a natural way to create a better dataset. Different from distillation, interpolation picks the sentence with the highest sentence-level BLEU score w.r.t. the ground truth from K−K-best beam search hypotheses. In our experiments, we first run beam search using the base Transformer model with a beam size of 55 then select the sentences with the highest BLEU score from the top-33 candidates.

Tab. 4 compares the performance of LevT trained with distilled data from the AT model with the standard distillation or interpolation. We observe that selection with BLEU score from the base AT model (base-inter) improves the performance of LevT ∼\sim 0.4 BLEU while the dataset complexity C(d)C(d) does not increase much.

Conclusion

In this paper, we first systematically examine why knowledge distillation improves the performance of NAT models. We conducted extensive experiments with autoregressive teacher models of different capacity and a wide range of NAT models. Furthermore, we defined metrics that can quantitatively measure the complexity of a parallel data set. Empirically, we find that a higher-capacity NAT model requires a more complex distilled data to achieve better performance. Accordingly, we propose several techniques that can adjust the complexity of a data set to match the capacity of an NAT model for better performance.

References

Appendix A Experimental Details

All the AT models are implemented based on the Transformer model using fairseq (Ott et al., 2019), and we basically follow the fairseq examples to train the transformershttps://github.com/pytorch/fairseq/blob/master/examples/translation.. Following the notation from Vaswani et al. (2017), we list the basic parameters of all the AT model we used:

Training

For all experiments, we adopt the Adam optimizer (Kingma & Ba, 2014) using β1=0.9,β2=0.98,ϵ=1e−8\beta_{1}=0.9,\beta_{2}=0.98,\epsilon=1e-8. The learning rate is scheduled using inverse_sqrt with a maximum learning rate 0.00050.0005 and 40004000 warmup steps. We set the label smoothing as 0.10.1. All the models are run on 88 GPUs for 300,000300,000 updates with an effective batch size of 32,00032,000 tokens. The best model is selected based on the validation loss except for FlowSeq which uses valid BLEU score.

Decoding

After training, we use beam-search with a fixed beam size 55 for all AT models to create the distilled dataset. We use length normalization without length penalty.

A.2 NAT Models

Tab. 2 also lists all the NAT models we test in this work. In general, all the NAT models except FlowSeq and LevT-big adopts a similar architecture and hyper-parameters as the Transformer-base (see Tab. 5). LevT-big is a naive extension of the original LevT model with a comparable parameter setting as Transformer-big (Tab. 5). For FlowSeq, we use the base model (FlowSeq-base) described in (Ma et al., 2019). We re-implemented the vanilla NAT as a simplified version of Gu et al. (2018) where instead of modeling fertility as described in the original paper, we monotonically copy the encoder embeddings to the input of the decoder. All the models except InsT require the additional module to predict the length of the output sequence, or the number of placeholders to be inserted, which is implemented as a standard softmax classifier over the lengths of [0, 256). For LevT, we also have a binary classifier to predict the deletion of the incorrect tokens.

Training

Similar to the AT models, all the NAT models are trained using the Adam optimizer with the same learning rate scheduler, in which the warmup steps are set to 10,00010,000. We train the FlowSeq model on 3232 GPUs with a batch size as 20482048 sentences, while all the other models are trained on 88 GPUs with an effective batch size of 64,00064,000 tokens. Note that, the batch sizes for training NAT is typically larger than the AT model, which improves final results. There are also specialized training settings for each models:

iNAT (Lee et al., 2018): following the original paper, we train the iNAT model jointly with 44 iterations of refinement during training. For each iteration, the model has the 50%50\% probability to learn as a denoising autoencoder, and the rest of the probability to learn from the model’s own prediction.

InsT (Stern et al., 2019): in this work, we only consider training the Insertion Transformer (InsT) using the slot-loss based on the uniform loss function (Stern et al., 2019). That is, we assign equal probabilities to all the insertable tokens inside each slot.

MaskT (Ghazvininejad et al., 2019): following the original paper, we train the model as a typical masked language model where the ratio of masked tokens is sampled from 0∼100%0\sim 100\%.

LevT (Gu et al., 2019): in this work, we only consider sequence generation tasks, which means the training of LevT is very similar to InsT. We use sentences with randomly deleted tokens to learn insertion, and learn deletion based on the model’s own prediction.

Decoding

For a fair comparison over all the NAT models, we use greedy decoding for all the models without considering any advanced decoding methods such as searching or re-ranking from a teacher model. For the vanilla NAT and FlowSeq, decoding is quite straight-forward and simply picks the arg max⁡\operatorname*{arg\,max} at every position. For iNAT and MaskT, we fix the decoding steps to 1010. Both InsT and LevT decode in an adaptive number of iterations, and we set the maximum iterations for both models to be 1010. A special EOS penalty that penalizes generating too short sequences is tuned based on the validation set for both InsT and LevT.

For all models, final results are calculated using tokenized BLEU score.

Appendix B Real Data Statistics

The detailed dataset split for WMT14 En-De is shown in Tab. 6. In Fig. 7, we also plot the histogram of the conditional entropy of each pair of sentences H(y∣x)\mathcal{H}({\bm{y}}|{\bm{x}}) in the real parallel data and different distilled data sets from the big-AT, base-AT, small-AT and tiny-AT respectively. It shows that the distribution of the sentence-level conditional entropy differs widely. The mode of H(y∣x)\mathcal{H}({\bm{y}}|{\bm{x}}) in the real data is the highest and follows by distilled data from the big-AT, base-AT, small-AT and tiny-AT. This observation aligns with the complexity value C(d)C(d) proposed in section 3.2.

Appendix C Additional Metrics

In Figure 8, we also showed results with different metrics together with BLEU scores considering that BLEU scores sometimes cannot fully capture the changes in the system. We considered 5 additional metrics in our experiments: METEOR (Banerjee & Lavie, 2005), RIBES (Isozaki et al., 2010), ChrF (Popović, 2015) TER (Snover et al., 2006), and BEER (Stanojevic & Sima’an, 2014). Not surprisingly, we find that all the metrics are correlated with the original BLEU scores quite well showing a similar trend as discussed earlier.

Appendix D Synthetic Data with access to the True Distribution

Bayesian decision theory is a fundamental statistical approach to the problem of pattern classification, which provides a principled rule of finding the optimal classification decision using probability and losses that accompany such decisions.

In the problem of structured prediction (Ma et al., 2017), let x{\bm{x}} denote the input sequence and y{\bm{y}} denote the output label sequence. Let H\mathcal{H} denote all the possible hypothesis functions from the input to the output space: H={h:X→Y}\mathcal{H}=\{h:\mathcal{X}\rightarrow\mathcal{Y}\}. Let r(y∣x)r({\bm{y}}|{\bm{x}}) denote the conditional risk on the input x{\bm{x}}, which is the expected loss of predicting y{\bm{y}} based on the posterior probabilities:

, where L(y,y′)L({\bm{y}},{\bm{y}}^{\prime}) is the loss function that penalizes predicting the true target y′{\bm{y}}^{\prime} as y{\bm{y}}. The classification task aims to find a hypothesis function hh that minimizes the overall risk RR given by

This is known as the Bayes risk. To minimize the overall risk, obviously we need to minimize the conditional risk for each input x{\bm{x}}. The Bayesian decision rule states that the global minimum of R(h)R(h) is achieved when the classifier make predictions that minimize each conditional risk given x{\bm{x}} and this gives the Bayes optimal classifier:

, which is the most probable output label sequence given the input sequence x{\bm{x}}.

This suggests that the Bayes classifier finds the most probable label at each time step given the input sequence.

D.2 Experimental Setups and Analysis

First, we use hseq∗(x)h_{seq}^{*}({\bm{x}}) to generate the distillation labels y′{\bm{y}}^{\prime} from the true HMM, which corresponds to applying the Viterbi decoding to each xi{\bm{x}}_{i} in DrealD_{real}. The training data set DseqD_{seq} is created with (xi{\bm{x}}_{i}, yi′{\bm{y}}^{\prime}_{i}). Next, we use htok∗(x)h_{tok}^{*}({\bm{x}}) to generate the distillation labels y^\hat{{\bm{y}}} and create the training data DtokD_{tok} of (xi,y^i)({\bm{x}}_{i},\hat{{\bm{y}}}_{i}). To generate y^\hat{{\bm{y}}}, we apply the forward-backward algorithm to each xi{\bm{x}}_{i} in DrealD_{real} and obtain P(yit∣xi)P(y_{i}^{t}|{\bm{x}}_{i}). We take arg max⁡\operatorname*{arg\,max} over the label space L\mathcal{L}: y^it=arg max⁡yit∈LP(yit∣xi)\hat{y}_{i}^{t}=\operatorname*{arg\,max}\limits_{y_{i}^{t}\in\mathcal{L}}P(y_{i}^{t}|{\bm{x}}_{i}).

We use these three training data (Dreal,Dtok,DseqD_{real},D_{tok},D_{seq}) to train the Bi-LSTM classifier respectively. We repeat the experiment for 50 times by constructing 50 HMM models with different random seeds as the data generator. We find that when evaluating with the token-level accuracy tacctacc, models trained with DtokD_{tok} yields the best performance (Bi-LSTM trained with DtokD_{tok} win 97.6% runs); when evaluating with the sequence-level accuracy saccsacc, models trained with DseqD_{seq} yields the best performance (Bi-LSTM trained with DseqD_{seq} win 98.5% runs). This is because the Bi-LSTM classifier has difficulty modeling the true data distribution defined by an HMM. On the other hand, it is easier for the Bi-LSTM classifier to model the distributions of DseqD_{seq} and DtokD_{tok}. Data sets DseqD_{seq} and DtokD_{tok} define deterministic conditional distributions over the input data, which are much simpler than the real data distribution. By definition, DtokD_{tok} is created by the optimal Bayes classifier htok∗(x)h_{tok}^{*}({\bm{x}}), this means that the Bi-LSTM classifier trained with DtokD_{tok} can better capture the distribution of P(yt∣x)=max⁡utP(ut∣x)P(y_{t}|{\bm{x}})=\max\limits_{u_{t}}P(u_{t}|{\bm{x}}), which can generalize better to the test data when evaluated with the token-level accuracy. Similarly, Bi-LSTM trained with DseqD_{seq} performs better on the test data with the sequence-level metric.

This corroborates our observation in machine translation task that NAT has difficulty in modeling the real conditional distribution of true sentence pairs. However, when using the distilled data translated from a pretrained autoregressive model with beam-search decoding, it performs better on the test set when evaluated with the BLEU score metric.