Eu estou usando a função plotgauss2d do BNT para visualizar como a resposta de um nó Gaussiano 2D muda quando a evidência é observada em outras partes da rede.

eng = jtree_inf_engine(bnet);
evidence = cell(1, 2)
eng = enter_evidence(eng, evidence);
marginals = marginal_nodes(eng, 1); p_1 = marginals.T
marginals = marginal_nodes(eng, 2); p_2 = marginals.T
marginals
clf; plotgauss2d(marginals.mu, marginals.Sigma);
hold all;
evidence{1} = 1;
marginals = marginal_nodes(enter_evidence(eng, evidence), 2);
p = plotgauss2d(marginals.mu, marginals.Sigma);
set(p, 'Color', 'green');
evidence{1} = 2;
marginals = marginal_nodes(enter_evidence(eng, evidence), 2);
p = plotgauss2d(marginals.mu, marginals.Sigma);
set(p, 'Color', 'red');
legend({'Unknown', 'Class 1', 'Class 2'});
hold off;

Como você pode ver, a legenda não está pegando as cores de plotagem alteradas que precisei definir manualmente. (Infelizmente, plotgauss2d não percorre as cores automaticamente, como você pode desejar.)

Existe uma maneira de definir as cores da linha usadas na legenda também?