function B = mk_ghmm_obs_lik(data, mu, Sigma)
% MK_GHMM_OBS_LIK Make the observation likelihood vector for a Gaussian HMM.
% B = mk_ghmm_obs_lik(data, mu, Sigma)
%
% Inputs:
% data(:,t) = y(t) = observation vector at time t
% mu(:,j) = E[Y(t) | Q(t)=j]
% Sigma(:,:,j) = Cov[Y(t) | Q(t)=j]
%
% Output:
% B(i,t) = Pr(y(t) | Q(t)=i)

Q = size(mu, 2);
O = size(data, 1);
T = size(data, 2);
B = zeros(Q,T);

for j=1:Q
  B(j,:) = gaussian_prob(data, mu(:,j), Sigma(:,:,j))';
end
