Hidden Markov Models Part 2: Algorithms

Transcription

1 Hidden Markov Models Part 2: Algorithms CSE 6363 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1

2 Hidden Markov Model An HMM consists of: A set of states s 1,, s K. An initial state probability function π k = p(z 1 = s k ). A state transition matrix A, of size K K, where A k,j = p z n = s j z n 1 = s k ) Observation probability functions, also called emission probabilities, defined as: φ k x = p x n = x z n = s k ) 2

3 The Basic HMM Problems We will next see algorithms for the four problems that we usually want to solve when using HMMs. Problem 1: learn an HMM using training data. This involves: Learning the initial state probabilities π k. Learning the state transition matrix A. Learning the observation probability functions P k x. Problems 2, 3, 4 assume we have already learned an HMM. Problem 2: given a sequence of observations X = x 1,, x N, compute p X HMM): the probability of X given the model. Problem 3: given an observation sequence X = x 1,, x N, find the hidden state sequence Z = z 1,, z N maximizing p Z X, HMM) Problem 4: given an observation sequence X = x 1,, x N, compute, for any n, k, the probability p z n = s k X, HMM) 3

4 The Basic HMM Problems Problem 1: learn an HMM using training data. Problem 2: given a sequence of observations X = x 1,, x N, compute p X HMM): the probability of X given the model. Problem 3: given an observation sequence X = x 1,, x N, find the hidden state sequence Z = z 1,, z N maximizing p Z X, HMM) Problem 4: given an observation sequence X = x 1,, x N, compute, for any n, k, the probability p z n = s k X, HMM) We will first look at algorithms for problems 2, 3, 4. Last, we will look at the standard algorithm for learning an HMM using training data. 4

5 Probability of Observations Problem 2: given an HMM model θ, and a sequence of observations X = x 1,, x N, compute p X θ): the probability of X given the model. Inputs: θ = (π, A, φ k ): a trained HMM, with probability functions π, A, φ k specified. A sequence of observations X = x 1,, x N. Output: p X θ). 5

6 Probability of Observations Problem 2: given a sequence of observations X = x 1,, x N, compute p X θ). Why do we care about this problem? What would be an example application? 6

7 Probability of Observations Problem 2: given a sequence of observations X = x 1,, x N, compute p X θ). Why do we care about this problem? We need to compute p X θ) to classify X. Suppose we have multiple models θ 1,, θ C. For example, we can have one model for digit 2, one model for digit 3, one model for digit 4 A Bayesian classifier classifies X by finding the θ c that maximizes p(θ c X. Using Bayes rule, we must maximize p(x θ c p(θ c ) p(x). Therefore, we need to compute p(x θ c for each c. 7

8 The Sum Rule The sum rule, which we saw earlier in the course, states that: p X = p X, Y = y y Y We can use the sum rule, to compute p X θ as follows: p X θ = p X, Z θ Z According to this formula, we can compute p X θ by summing p X, Z θ over all possible state sequences Z. 8

9 The Sum Rule p X θ = p X, Z θ Z According to this formula, we can compute p X θ by summing p X, Z θ over all possible state sequences Z. An HMM with parameters θ specifies a joint distribution function p X, Z θ as: p X, Z θ = p z 1 θ p z n z n 1, θ) Therefore: N n=2 p X θ = p z 1 θ p z n z n 1, θ) Z N n=2 N n=1 N n=1 p x n z n, θ) p x n z n, θ) 9

10 The Sum Rule p X θ = p z 1 θ p z n z n 1, θ) N N p x n z n, θ) Z n=2 n=1 According to this formula, we can compute p X θ by summing p X, Z θ over all possible state sequences Z. What would be the time complexity of doing this computation directly? It would be linear to the number of all possible state sequences Z, which is typically exponential to N (the length of X). Luckily, there is a polynomial time algorithm for computing p X θ, that uses dynamic programming. 10

11 Dynamic Programming for p X θ p X θ = p z 1 θ p z n z n 1, θ) N N p x n z n, θ) Z n=2 n=1 To compute p X θ, we use a dynamic programming algorithm that is called the forward algorithm. We define a 2D array of problems, of size N K. N: length of observation sequence X. K: number of states in the HMM. Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ 11

12 The Forward Algorithm - Initialization p X θ = p z 1 θ p z n z n 1, θ) N N p x n z n, θ) Z n=2 n=1 Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Problem (1, k):??? 12

13 The Forward Algorithm - Initialization p X θ = p z 1 θ p z n z n 1, θ) N N p x n z n, θ) Z n=2 n=1 Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Problem (1, k): Compute α 1, k = p x 1, z 1 = s k θ Solution:??? 13

14 The Forward Algorithm - Initialization p X θ = p z 1 θ p z n z n 1, θ) N N p x n z n, θ) Z n=2 n=1 Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Problem (1, k): Compute α 1, k = p x 1, z 1 = s k θ Solution: α 1, k = p x 1, z 1 = k θ = π k φ k (x 1 ) 14

15 The Forward Algorithm - Main Loop Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Solution: 15

16 The Forward Algorithm - Main Loop Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Solution: K p x 1,, x n, z n = s k θ = p x 1,, x n, z n 1 = s j, z n = s k θ j=1 K = φ k (x n ) p x 1,, x n 1, z n 1 = s j, z n = s k θ j=1 K = φ k x n p x 1,, x n 1, z n 1 = s j θ A j,k j=1 16

17 The Forward Algorithm - Main Loop Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Solution: p x 1,, x n, z n = k θ K = φ k x n p x 1,, x n 1, z n 1 = s j θ A j,k j=1 K = φ k (x n ) α[n 1, j]a j,k j=1 17

18 The Forward Algorithm - Main Loop Problem (n, k): Compute α n, k = p x 1,, x n, z n = s k θ Solution: p x 1,, x n, z n = k θ = φ k (x n ) α[n 1, j]a j,k K K j=1 Thus, α n, k = φ k (x n ) j=1 α[n 1, j]a j,k. α n, k is easy to compute once all α[n 1, j] values have been computed. 18

19 The Forward Algorithm Our aim was to compute p X θ. Using the dynamic programming algorithm we just described, we can compute α N, k = p X, z N = s k θ. How can we use those α N, k values to compute p X θ? We use the sum rule: K p X θ = p X, z N = s k θ k=1 K p X θ = α[n, k] k=1 19

20 Problem 3: Find the Best Z Problem 3: Inputs: a trained HMM θ, and an observation sequence X = x 1,, x N Output: the state sequence Z = z 1,, z N maximizing p Z θ, X). Why do we care? Some times, we want to use HMMs for classification. In those cases, we do not really care about maximizing p Z θ, X), we just want to find the HMM parameters θ that maximize p θ X). Some times, we want to use HMMs to figure out the most likely values of the hidden states given the observations. In the tree ring example, our goal was to figure out the average temperature for each year. Those average temperatures were the hidden states. Solving problem 3 provides the most likely sequence of average temperatures. 20

21 Problem 3: Find the Best Z Problem 3: Input: an observation sequence X = x 1,, x N Output: the state sequence Z = z 1,, z N maximizing p Z θ, X). We want to maximize p Z θ, X). Using the definition of conditional probabilities: p Z θ, X) = p X, Z θ) p X θ) Since X is known, p X θ) will be the same over all possible state sequences Z. Therefore, to find the Z that maximizes p Z θ, X), it suffices to find the Z that maximizes p Z, X θ). 21

22 Problem 3: Find the Best Z Problem 3: Input: an observation sequence X = x 1,, x N Output: the state sequence Z = z 1,, z N maximizing p Z θ, X), which is the same as maximizing p Z, X θ). Solution: (again) dynamic programming. Problem (n, k): Compute G n, k and H n, k, where: G n, k = argmax z 1,,z n 1 p x 1,, x n, z 1,, z n 1, z n = s k θ Note that: We optimize over all possible values z 1,, z n 1. However, z n is constrained to be equal to s k. Values x 1,, x n are known, and thus they are fixed. 22

23 Problem 3: Find the Best Z Problem 3: Input: an observation sequence X = x 1,, x N Output: the state sequence Z = z 1,, z N maximizing p Z θ, X), which is the same as maximizing p Z, X θ). Solution: (again) dynamic programming. Problem (n, k): Compute G n, k and H n, k, where: H n, k = p x 1,, x n, G n, k θ H n, k is the joint probability of: the first n observations x 1,, x n the sequence we stored in G n, k 23

24 Problem 3: Find the Best Z Problem 3: Input: an observation sequence X = x 1,, x N Output: the state sequence Z = z 1,, z N maximizing p Z θ, X), which is the same as maximizing p Z, X θ). Solution: (again) dynamic programming. Problem (n, k): Compute G n, k and H n, k, where: H n, k = p x 1,, x n, G n, k θ H n, k is the joint probability of: the first n observations x 1,, x n the sequence we stored in G n, k G n, k stores a sequence of state values. H n, k stores a number (a probability). 24

25 The Viterbi Algorithm Problem (n, k): Compute G n, k and H n, k, where: G n, k = argmax z 1,,z n 1 p x 1,, x n, z 1,, z n 1, z n = s k θ H n, k = p x 1,, x n, G n, k θ We use a dynamic programming algorithm to compute G n, k and H n, k for all n, k such that 1 n N, 1 k K. This dynamic programming algorithm is called the Viterbi algorithm. 25

26 The Viterbi Algorithm - Initialization Problem (1, k): Compute G 1, k and H 1, k, where: G 1, k = argmax z 1,,z n 1 p x 1, z 1 = s k θ H 1, k = p x 1, z 1 = s k θ What is G 1, k? It is the empty sequence. We must maximize over the previous states z 1,, z n 1, but n = 1, so there are no previous states. What is H 1, k? H 1, k = p x 1, z 1 = s k θ = π k φ k x 1. 26

27 The Viterbi Algorithm Main Loop Problem (n, k): Compute G n, k and H n, k, where: G n, k = argmax z 1,,z n 1 p x 1,, x n, z 1,, z n 1, z n = s k θ H n, k = p x 1,, x n, G n, k θ How do we find G n, k? G n, k = G n, k s j for some G n, k and some j. The last element of G n, k is some s j. Therefore, G n, k is formed by appending some s j to some sequence G n, k 27

28 The Viterbi Algorithm Main Loop Problem (n, k): Compute G n, k and H n, k, where: G n, k = argmax z 1,,z n 1 p x 1,, x n, z 1,, z n 1, z n = s k θ H n, k = p x 1,, x n, G n, k θ G n, k = G n, k s j for some G n, k and some j. We can prove that G n, k = G n 1, j, for that j. The proof is very similar to the proof that DTW finds the optimal alignment. If G n 1, j is better than G n, k, then G n 1, j s j is better than G n, k, which is a contradiction. 28

29 The Viterbi Algorithm Main Loop Problem (n, k): Compute G n, k and H n, k, where: G n, k = argmax z 1,,z n 1 p x 1,, x n, z 1,, z n 1, z n = s k θ H n, k = p x 1,, x n, G n, k θ G n, k = G n 1, j s j for some j. Let's define j as: j = argmax j=1,,k Then: G n, k = G n 1, j s j. H n 1, j A j,k φ k x n H n, k = H n 1, j A j,kφ k x n 29

30 The Viterbi Algorithm Output Our goal is to find the Z maximizing p Z θ, X), which is the same as finding the Z maximizing p Z, X θ). The Viterbi algorithm computes G n, k and H n, k, where: G n, k = argmax z 1,,z n 1 p x 1,, x n, z 1,, z n 1, z n = s k θ H n, k = p x 1,, x n, G n, k θ Let's define k as: k = argmax H N, k k=1,,k Then, the Z maximizing p Z, X θ) is G N, k s k 30

31 State Probabilities at Specific Times Problem 4: Inputs: a trained HMM θ, and an observation sequence X = x 1,, x N Output: an array γ, of size N K, where γ n, k = p z n = s k X, θ). In words: given an observation sequence X, we want to compute, for any moment in time n and any state s k, the probability that the hidden state at that moment is s k. 31

32 State Probabilities at Specific Times Problem 4: Inputs: a trained HMM θ, and an observation sequence X = x 1,, x N Output: an array γ, of size N K, where γ n, k = p z n = s k X, θ). We have seen, in our solution for problem 2, that the forward algorithm computes an array α, of size N K, where α n, k = p x 1,, x n, z n = s k θ. We can also define another array β, of size N K, where β n, k = p x n+1,, x N θ, z n = s k. In words, β n, k is the probability of all observations after moment n, given that the state at moment n is s k. 32

33 State Probabilities at Specific Times Problem 4: Inputs: a trained HMM θ, and an observation sequence X = x 1,, x N Output: an array γ, of size N K, where γ n, k = p z n = s k X, θ). We have seen, in our solution for problem 2, that the forward algorithm computes an array α, of size N K, where α n, k = p x 1,, x n, z n = s k θ. We can also define another array β, of size N K, where β n, k = p x n+1,, x N θ, z n = s k. Then, α n, k β n, k = p x 1,, x n, z n = s k θ p x n+1,, x N θ, z n = s k = p(x, z n = s k θ) Therefore: γ n, k = p z n = s k X, θ) = α n,k β n,k P X θ) 33

34 State Probabilities at Specific Times The forward algorithm computes an array α, of size N K, where α n, k = p x 1,, x n, z n = s k θ. We define another array β, of size N K, where β n, k = p x n+1,, x N θ, z n = s k. Then, γ n, k = p z n = s k X, θ) = In the above equation: α n,k β n,k P X θ) α n, k and P X θ) are computed by the forward algorithm. We need to compute β n, k. We compute β n, k in a way very similar to the forward algorithm, but going backwards this time. The resulting combination of the forward and backward algorithms is called (not surprisingly) the forward-backward algorithm. 34

35 The Backward Algorithm We want to compute the values of an array β, of size N K, where β n, k = p x n+1,, x N θ, z n = s k. Again, we will use dynamic programming. Problem (n, k) is to compute value β n, k. However, this time we start from the end of the observations. First, compute β N, k = p {} θ, z N = s k. In this case, the observation sequence x n+1,, x N is empty, because n = N. What is the probability of an empty sequence? 35

36 Backward Algorithm - Initialization We want to compute the values of an array β, of size N K, where β n, k = p x n+1,, x N θ, z n = s k. Again, we will use dynamic programming. Problem (n, k) is to compute value β n, k. However, this time we start from the end of the observations. First, compute β N, k = p {} θ, z N = s k. In this case, the observation sequence x n+1,, x N is empty, because n = N. What is the probability of an empty sequence? β N, k = 1 36

37 Backward Algorithm - Initialization Next: compute β N 1, k. β N 1, k = p x N θ, z N 1 = s k = K j=1 p x N, z N = s j θ, z N 1 = s k = K j=1 p x N z N = s j p z N = s j θ, z N 1 = s k 37

38 Backward Algorithm - Initialization Next: compute β N 1, k. β N 1, k = p x N θ, z N 1 = s k K = p x N z N = s j p(z N = s j θ, z N 1 = s k j=1 K = φ j x N A k,j j=1 38

39 Backward Algorithm Main Loop Next: compute β n, k, for n < N 1. β n, k = p x n+1,, x N θ, z n = s k = K j=1 p x n+1,, x N, z n+1 = s j θ, z n = s k = K j=1 p x n+1 z n+1 = s j p x n+2,, x N, z n+1 = s j θ, z N 1 = s k 39

40 Backward Algorithm Main Loop Next: compute β n, k, for n < N 1. β n, k = p x n+1,, x N θ, z n = s k = K j=1 p x n+1 z n+1 = s j p x n+2,, x N, z n+1 = s j θ, z N 1 = s k K j=1 φ j x n+1 p x n+2,, x N θ, z n+1 = s j A k,j We will take a closer look at the last step 40

41 Backward Algorithm Main Loop K j=1 p x n+1 z n+1 = s j p x n+2,, x N, z n+1 = s j θ, z n = s k = K j=1 φ j x n+1 p x n+2,, x N θ, z n+1 = s j A k,j p x n+1 z n+1 = s j = φ j x n+1, by definition of φ j. 41

42 Backward Algorithm Main Loop K p x n+1 z n+1 = s j p x n+2,, x N, z n+1 = s j θ, z n = s k = j=1 K φ j x n+1 p x n+2,, x N θ, z n+1 = s j A k,j j=1 p x n+2,, x N, z n+1 = s j θ, z n = s k = p x n+2,, x N θ, z n = s k, z n+1 = s j p z n+1 = s j z n = s k by application of chain rule: P A, B = P A B P(B), which can also be written as P A, B C = P A B, C P B C) 42

43 Backward Algorithm Main Loop K j=1 p x n+1 z n+1 = s j p x n+2,, x N, z n+1 = s j θ, z n = s k = K j=1 φ j x n+1 p x n+2,, x N θ, z n+1 = s j A k,j p x n+2,, x N, z n+1 = s j θ, z n = s k = p x n+2,, x N θ, z n = s k, z n+1 = s j p z n+1 = s j z n = s k = p x n+2,, x N θ, z n = s k, z n+1 = s j A k,j (by definition of A k,j ) 43

44 Backward Algorithm Main Loop K j=1 p x n+1 z n+1 = s j p x n+2,, x N, z n+1 = s j θ, z n = s k = K j=1 φ j x n+1 p x n+2,, x N θ, z n+1 = s j A k,j p x n+2,, x N, z n+1 = s j θ, z n = s k = p x n+2,, x N θ, z n = s k, z n+1 = s j p z n+1 = s j z n = s k = p x n+2,, x N θ, z n = s k, z n+1 = s j A k,j = p x n+2,, x N θ, z n+1 = s j A k,j (x n+2,, x N are conditionally independent of z n given z n+1 ) 44

45 Backward Algorithm Main Loop In the previous slides, we showed that β n, k = p x n+1,, x N θ, z n = s k = K j=1 φ j x n+1 p x n+2,, x N θ, z n+1 = s j A k,j Note that p x n+2,, x N θ, z n+1 = s j is β n + 1, j. Consequently, we obtain the recurrence: K β n, k = φ j x n+1 β n + 1, j A k,j j=1 45

46 Backward Algorithm Main Loop K β n, k = φ j x n+1 β n + 1, j A k,j j=1 Thus, we can easily compute all values β n, k using this recurrence. The key thing is to proceed in decreasing order of n, since values β n, depend on values β n + 1,. 46

47 The Forward-Backward Algorithm Inputs: A trained HMM θ. An observation sequence X = x 1,, x N Output: An array γ, of size N K, where γ n, k = p z n = s k X, θ). Algorithm: Use the forward algorithm to compute α n, k. Use the backward algorithm to compute β n, k. For n = 1 to N: For k = 1 to K: γ n, k = p z n = s k X, θ) = α n,k β n,k P X θ) 47

48 Problem 1: Training an HMM Goal: Estimate parameters θ = (π k, A k,j, φ k ). The training data can consist of multiple observation sequences X 1, X 2, X 3,, X M. We denote the length of training sequence X i as N m. We denote the elements of each observation sequence X m as: X m = (x m,1, x m,2,, x m,nm ) Before we start training, we need to decide on: The number of states. The transitions that we will allow. 48

49 Problem 1: Training an HMM While we are given X 1, X 2, X 3,, X M, we are not given the corresponding hidden state sequences Z 1, Z 2, Z 3,, Z M. We denote the elements of each hidden state sequence Z m as: Z m = (z m,1, z m,2,, z m,nm ) The training algorithm is called Baum-Welch algorithm. It is an Expectation-Maximization algorithm, similar to the algorithm we saw for learning Gaussian mixtures. 49

50 Expectation-Maximization When we wanted to learn a mixture of Gaussians, we had the following problem: If we knew the probability of each object belonging to each Gaussian, we could estimate the parameters of each Gaussian. If we knew the parameters of each Gaussian, we could estimate the probability of each object belonging to each Gaussian. However, we know neither of these pieces of information. The EM algorithm resolved this problem using: An initialization of Gaussian parameters to some random or non-random values. A main loop where: The current values of the Gaussian parameters are used to estimate new weights of membership of every training object to every Gaussian. The current estimated membership weights are used to estimate new parameters (mean and covariance matrix) for each Gaussian. 50

51 Expectation-Maximization When we want to learn an HMM model using observation sequences as training data, we have the following problem: If we knew, for each observation sequence, the probabilities of the hidden state values, we could estimate the parameters θ. If we knew the parameters θ, we could estimate the probabilities of the hidden state values. However, we know neither of these pieces of information. The Baum-Welch algorithm resolves this problem using EM: At initialization, parameters θ are given (mostly) random values. In the main loop, these two steps are performed repeatedly: The current values of parameters θ are used to estimate new probabilities for the hidden state values. The current probabilities for the hidden state values are used to estimate new values for parameters θ. 51

52 Baum-Welch: Initialization As we said before, before we start training, we need to decide on: The number of states. The transitions that we will allow. When we start training, we initialize θ = (π k, A k,j, φ k ) to random values, with these exceptions: For any s k that we do not want to allow to ever be an initial state, we set the corresponding π k to 0. The rest of the training will keep these π k values always equal to 0. For any s k, s j such that we do not want to allow transitions from s k to s j, we set the corresponding A k,j to 0. The rest of the training will keep these A k,j values always equal to 0. 52

53 Baum-Welch: Initialization When we start training, we initialize θ = (π k, A k,j, φ k ) to random values, with these exceptions: For any s k that we do not want to allow to ever be an initial state, we set the corresponding π k to 0. The rest of the training will keep these π k values always equal to 0. For any s k, s j such that we do not want to allow transitions from s k to s j, we set the corresponding A k,j to 0. The rest of the training will keep these A k,j values always equal to 0. These initial choices constrain the topology of the resulting HMM model. For example, they can force the model to be fully connected, to be a forward model, or to be some other variation. 53

54 Baum-Welch: Expectation Step Before we start the expectation step, we have some current values for the parameters of θ = (π k, A k,j, φ k ). The goal of the expectation step is to use those values to estimate two arrays: γ m, n, k = p x m,1,, x m,n, z m,n = s k θ, X m This is the same γ n, k that we computed earlier, using the forwardbackward algorithm. Here we just need to make the array three-dimensional, because we have multiple observation sequences X m. We can compute γ m,, by running the forward-backward algorithm with θ and X m as inputs. ξ m, n, j, k = p z m,n 1 = s j, z m,n = s k θ, X m 54

55 Baum-Welch: Expectation Step To facilitate our computations, we will also extend arrays α and β to three dimensions, in the same way that we extended array γ. α m, n, k = p x m,1,, x m,n, z m,n = s k θ β m, n, k = p x m,n+1,, x m,nm θ, z m,n = s k Where needed, values α[m,, ], β[m,, ], γ[m,, ] can be computed by running the forward-backward algorithm with inputs θ and X m. 55

56 Computing the ξ Values Using Bayes rule: ξ m, n, j, k = p z m,n 1 = s j, z m,n = s k θ, X m = p X m θ, z m,n 1 = s j, z m,n = s k p z m,n 1 = s j, z m,n = s k θ p X m θ p X m θ can be computed using the forward algorithm. We will see how to simplify and compute: p X m θ, z m,n 1 = s j, z m,n = s k p z m,n 1 = s j, z m,n = s k θ 56

57 Computing the ξ Values ξ m, n, j, k = p z m,n 1 = s j, z m,n = s k θ, X m = p X m θ, z m,n 1 = s j, z m,n = s k p z m,n 1 = s j, z m,n = s k θ p X m θ p X m θ, z m,n 1 = s j, z m,n = s k = p x m,1,, x m,n 1 θ, z m,n 1 = s j p x m,n,, x m,nm θ, z m,n = s k Why? Earlier observations x m,1,, x m,n 1 are conditionally independent of state z m,n given state z m,n 1. Later observations x m,n,, x m,nm are conditionally independent of state z m,n 1 given state z m,n. 57

58 Computing the ξ Values ξ m, n, j, k = p z m,n 1 = s j, z m,n = s k θ, X m = p X m θ, z m,n 1 = s j, z m,n = s k p z m,n 1 = s j, z m,n = s k θ p X m θ p X m θ, z m,n 1 = s j, z m,n = s k = p x m,1,, x m,n 1 θ, z m,n 1 = s j p x m,n,, x m,nm θ, z m,n = s k p x m,n,, x m,nm θ, z m,n = s k = p x m,n θ, z m,n = s k p x m,n+1,, x m,nm θ, z m,n = s k = φ k x m,n β m, n, k 58

59 Computing the ξ Values ξ m, n, j, k = p z m,n 1 = s j, z m,n = s k θ, X m = p X m θ, z m,n 1 = s j, z m,n = s k p z m,n 1 = s j, z m,n = s k θ p X m θ p X m θ, z m,n 1 = s j, z m,n = s k = p x m,1,, x m,n 1 θ, z m,n 1 = s j φ k x m,n β m, n, k = p x m,1,, x m,n 1, z m,n 1 = s j θ p z m,n 1 = s j θ φ k x m,n β m, n, k = α m, n 1, j γ m, n 1, j φ k x m,n β m, n, k 59

60 Computing the ξ Values ξ m, n, j, k = α m, n 1, j γ m, n 1, j φ k x m,n β m, n, k p z m,n 1 = s j, z m,n = s k θ p X m θ We can further process p z m,n 1 = s j, z m,n = s k θ : p z m,n 1 = s j, z m,n = s k θ = p z m,n = s k θ, z m,n 1 = s j p z m,n 1 = s j θ = A j,k γ m, n 1, j 60

61 Computing the ξ Values So, we get: ξ m, n, j, k = α m, n 1, j γ m, n 1, j φ k x m,n β m, n, k A j,k γ m, n 1, j p X m θ = α m, n 1, j φ k x m,n β m, n, k A j,k p X m θ 61

62 Computing the ξ Values Based on the previous slides, the final formula for ξ is: ξ m, n, j, k = α m, n 1, j φ k x m,n β m, n, k A j,k p X m θ This formula can be computed using the parameters θ and the algorithms we have covered: α m, n 1, j is computed with the forward algorithm. β m, n, k is computed by the backward algorithm. p X m θ is computed with the forward algorithm. φ k and A j,k are specified by θ. 62

63 Baum-Welch: Summary of E Step Inputs: Training observation sequences X 1, X 2, X 3,, X M. Current values of parameters θ = (π k, A k,j, φ k ). Algorithm: For m = 1 to M: Run the forward-backward algorithm to compute: α m, n, k for n and k such that 1 n N m, 1 k K. β m, n, k for n and k such that 1 n N m, 1 k K. γ m, n, k for n and k such that 1 n N m, 1 k K. p X m θ For n, j, k such that 1 n N m, 1 j, k K: ξ m, n, j, k = α m, n 1, j φ k x m,n β m, n, k A j,k p X m θ 63

64 Baum-Welch: Maximization Step Before we start the maximization step, we have some current values for: γ m, n, k = p z m,n = s k θ, X m ξ m, n, j, k = p z m,n 1 = s j, z m,n = s k θ, X m The goal of the maximization step is to use those values to estimate new values for θ = (π k, A k,j, φ k ). In other words, we want to estimate new values for: Initial state probabilities π k. State transition probabilities A k,j. Observation probabilities φ k. 64

65 Baum-Welch: Maximization Step Initial state probabilities π k are computed as: π k = M m=1 M m=1 γ m, 1, k K j=1 γ m, 1, j In words, we compute for each k the ratio of these two quantities: The sum of probabilities, over all observation sequences X m, that state s k was the initial state for X m. The sum of probabilities over all observation sequences X m and all states s j that states j was the initial state for X m. 65

66 Baum-Welch: Maximization Step State transition probabilities A k,j are computed as: A k,j = M m=1 M m=1 N m n=2 N m n=2 ξ m, n, k, j K i=1 ξ m, n, k, i In words, we compute, for each k, j the ratio of these two quantities: The sum of probabilities, over all state transitions of all observation sequences X m, that state s k was followed by s j. The sum of probabilities, over all state transitions of all observation sequences X m and all states s i, that state s k was followed by state s i. 66

67 Baum-Welch: Maximization Step Computing the observation probabilities φ k depends on how we model those probabilities. For example, we could choose: Discrete probabilities. Histograms. Gaussians. Mixtures of Gaussians. We will show how to compute distributions φ k for the cases of: Discrete probabilities. Gaussians. 67

68 Baum-Welch: Maximization Step Computing the observation probabilities φ k depends on how we model those probabilities. In all cases, the key idea is that we treat observation x m,n as partially assigned to distribution φ k. γ m, n, k = p z m,n = s k θ, X m is the weight of the assignment of x m,n to φ k. In other words: We don't know what hidden state x m,n corresponds to, but we have computed, for each state s k, the probability γ m, n, k that x m,n corresponds to s k. Thus, x m,n influences our estimate of distribution φ k with weight proportional to γ m, n, k. 68

69 Baum-Welch: Maximization Step Suppose that the observations are discrete, and come from a finite set Y = y 1,, y R. In that case, each x m,n is an element of Y. Define an auxiliary function Eq x, y as: Eq x, y = 1 if x = y 0 if x y Then, φ k y r = p x m,n = y r z m,n = s k, and it can be computed as: M N m m=1 n=1 γ m, n, k Eq x m,n, y r φ k y r = M N m γ m, n, k m=1 n=1 69

70 Baum-Welch: Maximization Step If observations come from a finite set Y = y 1,, y R : φ k y r = M m=1 N m n=1 M m=1 γ m, n, k Eq x m,n, y r N m n=1 γ m, n, k The above formula can be seen as a weighted average of the times x m,n is equal to y r, where the weight of each x m,n is the probability that x m,n corresponds to hidden state s k. 70

71 Baum-Welch: Maximization Step Suppose that observations are vectors in R D, and that we model φ k as a Gaussian distribution. In that case, for each φ k we need to estimate a mean μ k and a covariance matrix S k. μ k y r = M m=1 M m=1 N m n=1 γ m, n, k x m,n N m n=1 γ m, n, k Thus, μ k is a weighted average of values x m,n, where the weight of each x m,n is the probability that x m,n corresponds to hidden state s k. 71

72 Baum-Welch: Maximization Step Suppose that observations are vectors in R D, and that we model φ k as a Gaussian distribution. In that case, for each φ k we need to estimate a mean μ k and a covariance matrix S k. μ k = M m=1 M m=1 N m n=1 γ m, n, k x m,n N m n=1 γ m, n, k S k = M m=1 N m n=1 γ m, n, k x m,n μ k x m,n μ k T M m=1 N m n=1 γ m, n, k 72

73 Baum-Welch: Summary of M Step π k = M m=1 M m=1 γ m, 1, k K j=1 γ m, 1, j A k,j = M m=1 M m=1 N m n=2 N m n=2 ξ m, n, k, j K i=1 ξ m, n, k, i For discrete observations: φ k y r = M m=1 N m n=1 M m=1 γ m, n, k Eq x m,n, y r N m n=1 γ m, n, k 73

74 Baum-Welch: Summary of M Step π k = M m=1 M m=1 γ m, 1, k K j=1 γ m, 1, j A k,j = M m=1 M m=1 N m n=2 N m n=2 ξ m, n, k, j K i=1 ξ m, n, k, i For Gaussian observation distributions: φ k is a Gaussian with mean μ k and covariance matrix S k, where: μ k = M m=1 M m=1 N m n=1 γ m, n, k x m,n N m n=1 γ m, n, k S k = M m=1 N m n=1 γ m, n, k x m,n μ k x m,n μ k T N m M m=1 n=1 γ m, n, k 74

75 Baum-Welch Summary Initialization: Initialize θ = (π k, A k,j, φ k ) with random values, but set to 0 values π k and A k,j to specify the network topology. Main loop: E-step: Compute all values p X m θ, α m, n, k, β m, n, k, γ m, n, k using the forward-backward algorithm. Compute all values ξ m, n, j, k = α m,n 1,j φ k x m,n β m,n,k A j,k p X m θ M-step: Update θ = (π k, A k,j, φ k ), using the values γ m, n, k and ξ m, n, j, k computed at the E-step.. 75

76 Hidden Markov Models - Recap HMMs are widely used to model temporal sequences. In HMMs, an observation sequence corresponds to a sequence of hidden states. An HMM is defined by parameters θ = (π k, A k,j, φ k ), and defines a joint probability p X, Z θ. An HMM can be used as a generative model, to produce synthetic samples from distribution p X, Z θ. HMMs can be used for Bayesian classification: One HMM θ c per class. Find the class c that maximizes p θ c X. 76

77 Hidden Markov Models - Recap HMMs can be used for Bayesian classification: One HMM θ c per class. Find the class c that maximizes p θ c X, using probabilities p X θ c calculated by the forward algorithm. HMMs can be used to find the most likely hidden state sequence Z for observation sequence X. This is done using the Viterbi algorithm. HMMs can be used to find the most likely hidden state z n for a specific value x n of an observation sequence X. This is done using the forward-backward algorithm. HMMs are trained using the Baum-Welch algorithm, which is an Expectation-Maximization algorithm. 77