function obslik = mk_hmm_obs_lik_matrix(engine, evidence)T  = size(evidence,2);Q = length(engine.startprob);obslik = ones(Q, T);bnet = bnet_from_engine(engine);% P(o1,o2| Q1,Q2) = P(o1|Q1,Q2) * P(o2|Q1,Q2)onodes = bnet.observed;for i=1:length(onodes)  data = cell2num(evidence(onodes(i),:));  if bnet.auto_regressive(onodes(i))    params = engine.obsprob{i};    obslik_i = mk_arhmm_obs_lik(data, params.big_mu, params.big_Sigma, params.big_W, ...				params.big_mu0, params.big_Sigma0);  elseif myismember(onodes(i), bnet.dnodes)    obslik_i = mk_dhmm_obs_lik(data, engine.obsprob{i}.big_CPT);  else    obslik_i = mk_ghmm_obs_lik(data, engine.obsprob{i}.big_mu, engine.obsprob{i}.big_Sigma);  end  obslik = obslik .* obslik_i;end