r = rand(100)
a_hat = ones(100);
a_hat[find(r <= p_2)] = 2
a_hat[find(r > p_2 && r <= p_3)] = 3

votes = zeros(3)
for n=1:100
   votes[a_hat[n]] += 1
end
action = indmax(votes)