matlab/zero_prob_transitions.m

   1 function [zero_transitions] = zero_prob_transitions(model, varmap)\r
   2 \r
   3 % Load the traces from the file\r
   4 % Skip over the first line\r
   5 % For each subsequent line, load the observation and the name of the obs\r
   6 % Store this in X(i).obs and X(i).name\r
   7 \r
   8 % Add up all of the transitions matrices, and then return the entries that\r
   9 % are still 0.\r
  10 \r
  11 transitions = zeros(model.num_states, model.num_states + 1);\r
  12 for i=1:model.num_mcs\r
  13     transitions = model.chain(i).T + transitions;\r
  14 end\r
  15 [row col] = find(transitions == 0);\r
  16 \r
  17 % Convert row col into a vector of strings\r
  18 for i=1:size(row, 1)\r
  19     from = char(varmap{row(i), 1});\r
  20     to = char(varmap{col(i), 1});\r
  21     zero_transitions{i, 1} = char(sprintf('%s:%s', from, to));\r
  22 end