code: def main():
      
      model = load_model('../data/hmmdefs')
      hmm = model.hmms['r-We']
      for state_name in hmm.state_names:
          print(state_name)
          state = model.states[state_name]
          print(state.means_)
      print(model)
      model2 = load_model('../data/prior.hmm1mixSI.rate32')
      print(model2)
docs: Test code called from commandline