hmmTest.jl

states = ["fair"; "unfair"]; 
coins = ["head", "tail"]; 

h = jHMM.HMM()
h = set_states(h,(:state,states))
h = set_transitions(h, f(state,statep) ~ P_tr(state,statep)  )

h = set_observation_space(h,(:coin,coins))
h = set_emission(h, f(state,coin) ~ P_em(state,coin) )

#flipping coins
cflip = (rand(30) .< 0.5 ) ;
cflip = [cflip; (rand(30) .< 0.5+bias ) ];
cflip = [cflip; (rand(30) .< 0.5 ) ];

obs = fill("head",length(cflip));
obs[cflip] = "tail";

h = set_observations(h,(:coin,obs))

 forward(h)