# Text Generation Tutorial

In this tutorial, we will use recurrent neural networks to **generate** 
sequences. 
Generating sequences is more involved comparing to making predictions about
sequences. However, it is a very interesting task. This tutorial will help
prepare you for project 4.

Much of today's content is an adaptation of the "Practical PyTorch" github 
repository, found here:
https://github.com/spro/practical-pytorch/blob/master/char-rnn-generation/char-rnn-generation.ipynb

```
import torch
import torch.nn as nn
import torch.nn.functional as F
import torch.optim as optim
```

## Reccurent Neural Networks Revisited

To use a recurrent neural network, our input sequence needs to be
broken down into tokens. We could choose whether to
tokenize based on words, or based on characters. The representation
of each token (GloVe or one-hot) is processed by the RNN one step at a time
to update the hidden (or context) state.

In a predictive RNN, the value of the
hidden states  is a representation of **all the text that was processed thus far**.
Similarly, in a generative RNN, The value of the hidden
state will be a representation of **all the text that still needs to be generated**.
We will use this hidden state to produce the sequence, one token at a time.

Similar to last class we will break up the problem of generating text
to generating one token at a time.

We will do so with the help of two functions:

1. We need to be able to generate the *next* token, given the current 
   hidden state. In practice, we get a probability distribution over 
   the next token, and sample from that probability distribution.
2. We need to be able to update the hidden state somehow. To do so,
   we need two piece of information: the old hidden state, and the actual
   token that was generated in the previous step. The actual token generated
   will inform the subsequent tokens.

We will repeat both functions until a special "END OF SEQUENCE" token is
generated. 

Note that there are several tricky things that we will have to figure out.
For example, how do we actually sample the actual token from the probability
distribution over tokens? What would we do during training, and how might 
that be different from during testing/evaluation? We will answer those
questions as we implement the RNN.

For now, let's start with our training data.

## Data: Happy Tweets from 2018

The training set we use is a collection of "happy" tweets from the Sentiment140
data set.
We will only use tweets that are 140 characters or shorter, and tweets
that contains more than just a URL.
Since tweets often contain creative spelling and numbers, and upper vs lower
case characters are read very differently, we will use a character-level RNN.

```
import csv
tweets = list(line[0] for line in csv.reader(open('happy.csv')))
len(tweets)
```

Let's look at a few of these tweets, just to get a sense of the kind of text
we're dealing with:

```
print(tweets[100])
print(tweets[1000])
print(tweets[10000])
```

## Generating One Tweet

Our generative RNN is a **generative model**. At a high level, generative
models are trained to maximize the probability of generating our training data.
We'll need to figure out what that means, and to manipulate our data 
Our RNN model generates text one character at a time. In each iteration,
we'll train our data to generate a different tweet.

Normally, when we build a new machine learn model, we want to make sure
that our model can overfit. To that end, we will first build a neural network
that can generate _one_ tweet really well. We can choose any tweet (or any other text)
we want.  Let's choose to build an RNN that generates `tweet[200]`.

```
tweet = tweets[200]
tweet
```

Since tweets are often poorly spelled, we will treat a tweet as a **sequence of characters**.
Our RNN will generate one **character** as a time.

Since PyTorch works with numbers instead of strings, we will need to convert
characters into integers. This of this integer as a sparse representation of a
one-hot encoding.  We'll build dictionary mappings
from the character to the index of that character (a unique integer identifier),
and from the index to the character. We'll use the same naming scheme that `torchtext`
uses (`stoi` and `itos`, which means "string to index" and "index to string").

For simplicity, we'll work with a limited vocabulary containing
just the characters in `tweet[200]`, plus two special tokens:

- `<EOS>` represents "End of String", which we'll append to the end of our tweet.
  Since tweets are variable-length, this is a way for the RNN to signal
  that the entire sequence has been generated.
- `<BOS>` represents "Beginning of String", which we'll prepend to the beginning of 
  our tweet. This is the first token that we will feed into the RNN.

The way we use these special tokens will become more clear as we build the model.

### Part (a)

Build two Python dictionary mappings
from the character to the index of that character (a unique integer identifier),
and from the index to the character. We'll use the same naming scheme that `torchtext`
uses (`stoi` and `itos`).

```
vocab = list(set(tweet)) + ["<BOS>", "<EOS>"]
vocab_stoi = {} # TODO: build a dictionary mapping of word to unique index
vocab_itos = {} # TODO: build a dictionary mapping of a unique index to a word (string)
vocab_size = len(vocab)
```

### Part (b)

Our model have three parts:

1. An **embedding** layer that takes a character index (sparse representation of a one-hot vector
   representing the character), and returns an **embedding** ${\bf x}^{(t)}$ for the character.
   We could use a one-hot embedding here, but using a low-dimensional but dense embedding
   is more powerful, and is more in line with what we will do in Project 4.
2. A **recurrent neural network** layer. We will use a `nn.GRU` unit. At each time step,
   this layer will take the previous hidden state ${\bf h}^{(t-1)}$ and the embedding of
   a new generated token ${\bf x}^{(t)}$ and compute the new hidden state ${\bf h}^{(t)}$ representing
   *tokens not yet generated*.
3. The **projection MLP** (a fully-connected layer) that takes the hidden state ${\bf h}^{(t)}$
   and computes teh distribution of the next character to generate.

Filling in the missing numbers in the `__init__` method using
the parameters `vocab_size`, `emb_size`, and `hidden_size`.

```
class TextGenerator(nn.Module):
    def __init__(self,
                 vocab_size,     # number of unique characters in our vocabulary
                 embedding_size, # size of the word embeddings ${\bf x}^{(t)}$
                 hidden_size):   # size of the hidden state in the RNN
        super(TextGenerator, self).__init__()

        # Embedding 
        self.embed = nn.Embedding(num_embeddings=None, # TODO
                                  embedding_dim=None)  # TODO

        # recurrent neural network
        self.rnn = nn.GRU(input_size=None, #TODO
                          hidden_size=None, #TODO
                          batch_first=True)

        # a fully-connect layer that outputs a distribution over
        # the next token, given the RNN output
        self.proj = nn.Linear(in_features=None, # TODO
                              out_features=None) # TODO

    def forward(self, inp, hidden=None):
        emb = self.embed(inp)                  # generate one-hot vectors of input
        output, hidden = self.rnn(emb, hidden) # get the next output and hidden state
        output = self.proj(output)             # predict distribution over next tokens
        return output, hidden

model = TextGenerator(vocab_size, 128, 128)
```

## Training with Teacher Forcing

At a very high level, we want our RNN model to have a high probability
of generating the tweet. An RNN model generates text
one character at a time based on the hidden state value.
At each time step, we will check whether the mdoel generated the
correct character. That is, at each time step,
we are trying to select the correct next character out of all the 
characters in our vocabulary. Recall that this problem is a multi-class
classification problem, and we can use Cross-Entropy loss to train our
network to become better at this type of problem.

```
criterion = nn.CrossEntropyLoss()
```

However, we don't just have a single multi-class classification problem.
Instead, we have **one classification problem per time-step** (per token)!
So, how do we predict the first token in the sequence? 
How do we predict the second token in the sequence? 

To help you understand what happens durign RNN training, we'll start with a
inefficient training code that shows you what happens step-by-step. We'll
start with computing the loss for the first token generated, then the second token,
and so on.
Later on, we'll switch to a simpler and more performant version of the code.

So, let's start with the first classification problem: the problem of generating
the **first** token (`tweet[0]`).

To generate the first token, we'll feed the RNN network (with an initial, empty
hidden state) the "<BOS>" token. Then, the output 

```
bos_input = torch.Tensor([vocab_stoi["<BOS>"]]).long().unsqueeze(0)
output, hidden = model(bos_input, hidden=None)
output # distribution over the first token
```

We can compute the loss using `criterion`. Since the model is untrained,
the loss is expected to be high. (For now, we won't do anything
with this loss, and omit the backward pass.)

```
target = torch.Tensor([vocab_stoi[tweet[0]]]).long().unsqueeze(0)
criterion(output.reshape(-1, vocab_size), # reshape to 2D tensor
          target.reshape(-1))             # reshape to 1D tensor
```

Now, we need to update the hidden state and generate a prediction
for the next token. To do so, we need to provide the current token to
the RNN. We already said that during test time, we'll need to sample
from the predicted probabilty over tokens that the neural network
just generated. 

Right now, we can do something better: we can **use the ground-truth,
actual target token**. This technique is called **teacher-forcing**, 
and generally speeds up training. The reason is that right now, 
since our model does not perform well, the predicted probability
distribution is pretty far from the ground truth. So, it is very,
very difficult for the neural network to get back on track given bad
input data.

```
# Use teacher-forcing: we pass in the ground truth `target`,
# rather than using the NN predicted distribution
output, hidden = model(target, hidden)
output # distribution over the second token
```

Similar to the first step, we can compute the loss, quantifying the
difference between the predicted distribution and the actual next
token. This loss can be used to adjust the weights of the neural
network (which we are not doing yet).

```
target = torch.Tensor([vocab_stoi[tweet[1]]]).long().unsqueeze(0)
criterion(output.reshape(-1, vocab_size), # reshape to 2D tensor
          target.reshape(-1))             # reshape to 1D tensor
```

We can continue this process of:

- feeding the previous ground-truth token to the RNN,
- obtaining the prediction distribution over the next token, and
- computing the loss,

for as many steps as there are tokens in the ground-truth tweet.

```
for i in range(2, len(tweet)):
    output, hidden = model(target, hidden)
    target = torch.Tensor([vocab_stoi[tweet[1]]]).long().unsqueeze(0)
    loss = criterion(output.reshape(-1, vocab_size), # reshape to 2D tensor
                     target.reshape(-1))             # reshape to 1D tensor
    print(i, output, loss)
```

Finally, with our final token, we should expect to output the "<EOS>"
token, so that our RNN learns when to stop generating characters.

```
output, hidden = model(target, hidden)
target = torch.Tensor([vocab_stoi["<EOS>"]]).long().unsqueeze(0)
loss = criterion(output.reshape(-1, vocab_size), # reshape to 2D tensor
                 target.reshape(-1))             # reshape to 1D tensor
print(i, output, loss)
```

In practice, we don't really need a loop. Recall that in a predictive RNN,
the `nn.RNN` module can take an entire sequence as input. We can do the
same thing here:

```
tweet_ch = ["<BOS>"] + list(tweet) + ["<EOS>"]
tweet_indices = [vocab_stoi[ch] for ch in tweet_ch]
tweet_tensor = torch.Tensor(tweet_indices).long().unsqueeze(0)

print(tweet_tensor.shape)

output, hidden = model(tweet_tensor[:,:-1]) # <EOS> is never an input token
target = tweet_tensor[:,1:]                 # <BOS> is never a target token
loss = criterion(output.reshape(-1, vocab_size), # reshape to 2D tensor
                 target.reshape(-1))             # reshape to 1D tensor
```

Here, the input to our neural network model is the *entire*
sequence of input tokens (everything from "<BOS>" to the
last character of the tweet). The neural network generates a prediction distribution
of the next token at each step. We can compare each of these  with the ground-truth
`target`.

### Part (c)

Complete the code for the training loop to generate a single tweet.

```
optimizer = torch.optim.Adam(model.parameters(), lr=0.001)
criterion = nn.CrossEntropyLoss()
for it in range(1000): # You might need to increase this
    optimizer.zero_grad()
    output, _ = None # TODO
    loss = None      # TODO
    loss.backward()
    optimizer.step()

    if (it+1) % 100 == 0:
        print("[Iter %d] Loss %f" % (it+1, float(loss)))
```

Once this code runs, your training code should decrease with training,
which is what we expect.

## Generating a Token

At this point, we want to see whether our model is actually learning
something. So, we need to talk about how to
actually use the RNN model to generate text. If we can 
generate text, we can make a qualitative asssessment of how well
our RNN is performing.

The main difference between training and test-time (generation time)
is that we don't have the ground-truth tokens to feed as inputs
to the RNN. Instead, we need to actually **sample** a token based
on the neural network's prediction distribution.

But how can we sample a token from a distribution?

On one extreme, we can always take
the token with the largest probability (argmax). This has been our
go-to technique in other classification tasks. However, this idea
will fail here. The reason is that in practice, 
**we want to be able to generate a variety of different sequences from
the same model**. An RNN that can only generate a single new tweet
is fairly useless.

In short, we want some randomness. We can do so by using the logit
outputs from our model to construct a multinomial distribution over
the tokens, then and sample a random token from that multinomial distribution.

One natural multinomial distribution we can choose is the 
distribution we get after applying the softmax on the outputs.
However, we will do one more thing: we will add a **temperature**
parameter to manipulate the softmax outputs. We can set a
**higher temperature** to make the probability of each token
**more even** (more random), or a **lower temperature** to assighn
more probability to the tokens with a higher logit (output).
A **higher temperature** means that we will get a more diverse sample,
with potentially more mistakes. A **lower temperature** means that we
may see repetitions of the same high probability sequence.

```
def sample_sequence(model, max_len=100, temperature=0.8):
    generated_sequence = ""
   
    inp = torch.Tensor([vocab_stoi["<BOS>"]]).long()
    hidden = None
    for p in range(max_len):
        output, hidden = model(inp.unsqueeze(0), hidden)
        # Sample from the network as a multinomial distribution
        output_dist = output.data.view(-1).div(temperature).exp()
        top_i = int(torch.multinomial(output_dist, 1)[0])
        # Add predicted character to string and use as next input
        predicted_char = vocab_itos[top_i]
        
        if predicted_char == "<EOS>":
            break
        generated_sequence += predicted_char       
        inp = torch.Tensor([top_i]).long()
    return generated_sequence
```

## Part (d)

Try sampling from your model with different temperature settings
(e.g. from 0.8 to 5). Make sure to sample multiple times per setting.
Since we only trained the model on a single sequence, we won't see
the full effect of the temperature parameter yet. 

```python
print(sample_sequence(model, temperature=1.0))
```

For now, the output of the calls to the `sample_sequence` function
assures us that our training code looks reasonable, and we can
proceed to training on our full dataset!

## Training the Tweet Generator

For the actual training, let's use `torchtext` so that we can use
the `BucketIterator` to make batches. Like in Lab 5, we'll create a 
`torchtext.data.Field` to use `torchtext` to read the CSV file, and convert
characters into indices. The object has convient parameters to specify
the BOS and EOS tokens.

```
import torchtext

text_field = torchtext.data.Field(sequential=True,      # text sequence
                                  tokenize=lambda x: x, # because are building a character-RNN
                                  include_lengths=True, # to track the length of sequences, for batching
                                  batch_first=True,
                                  use_vocab=True,       # to turn each character into an integer index
                                  init_token="<BOS>",   # BOS token
                                  eos_token="<EOS>")    # EOS token

fields = [('text', text_field)]
tweets = torchtext.data.TabularDataset("happy.csv", "csv", fields)
len(tweets)
```

We will use the version of `vocab_stoi` and `vocab_itos` that torchtext provides.

```
text_field.build_vocab(tweets)
vocab_stoi = text_field.vocab.stoi # so we don't have to rewrite sample_sequence
vocab_itos = text_field.vocab.itos # so we don't have to rewrite sample_sequence
vocab_size = len(text_field.vocab.itos)
vocab_size
```

Let's just verify that the `BucketIterator` works as expected, but start with batch_size of 1.

```
data_iter = torchtext.data.BucketIterator(tweets, 
                                          batch_size=1,
                                          sort_key=lambda x: len(x.text),
                                          sort_within_batch=True)
for (tweet, lengths), label in data_iter:
    print(label)   # should be None
    print(lengths) # contains the length of the tweet(s) in batch
    print(tweet.shape) # should be [1, max(length)]
    break
```

## Part (e)

Modify the training code to account for batching.

```
def train(model, data, batch_size=1, num_epochs=1, lr=0.001, print_every=100):
    optimizer = torch.optim.Adam(model.parameters(), lr=lr)
    criterion = nn.CrossEntropyLoss()
    it = 0
    
    data_iter = torchtext.data.BucketIterator(data,
                                              batch_size=batch_size,
                                              sort_key=lambda x: len(x.text),
                                              sort_within_batch=True)
    for e in range(num_epochs):
        # get training set
        avg_loss = 0
        for (tweet, lengths), label in data_iter:
            target = None # TODO
            inp = None # TODO
            # cleanup
            optimizer.zero_grad()
            # forward pass
            output, _ = None # TODO
            loss = None # TODO
            # backward pass
            loss.backward()
            optimizer.step()

            avg_loss += loss
            it += 1 # increment iteration count
            if it % print_every == 0:
                print("[Iter %d] Loss %f" % (it+1, float(avg_loss/print_every)))
                print("    " + sample_sequence(model, 140, 0.8))
                avg_loss = 0

model = TextGenerator(vocab_size, 128, 128)
```

## Part (f)

Train your model. This model might take a looong time too train. However, the model
learns very quickly that many tweets begin with `@`.

```
#train(model, tweets, batch_size=1, num_epochs=1, lr=0.004, print_every=100)
#train(model, tweets, batch_size=32, num_epochs=1, lr=0.004, print_every=100)
```


## Part (g)

Try generating some sequences using your model. Vary your temperature settings.

```python
#print(sample_sequence(model, temperature=0.8))
#print(sample_sequence(model, temperature=0.8))
#print(sample_sequence(model, temperature=1.0))
#print(sample_sequence(model, temperature=1.0))
#print(sample_sequence(model, temperature=1.5))
#print(sample_sequence(model, temperature=1.5))
#print(sample_sequence(model, temperature=2.0))
#print(sample_sequence(model, temperature=2.0))
#print(sample_sequence(model, temperature=5.0))
#print(sample_sequence(model, temperature=5.0))
```
