% x = 0:0.1:1;
% y = 0:0.1:2.5;
% [X,Y] = meshgrid(x,y);
a = [errorxsub_initial(letter,:)' entropyxsub_initial(letter,:)']
%a = [errorxsub_final(letter,:)' entropyxsub_final(letter,:)']
y = [zeros(9,1);ones(6,1)];
options = statset('Display','final');
%h = gscatter(a(:,1),a(:,2),y);
obj1 = gmdistribution.fit(a(1:9,:),1,'Options',options);
obj2 = gmdistribution.fit(a(10:15,:),1,'Options',options);
%obj = gmdistribution.fit(a,2,'Options',options);
%contour(X,Y,a(1:9,:),[1,1]);
hold on
ezcontour(@(x,y)pdf(obj2,[x y]),get(gca,{'XLim','YLim'}));
ezcontour(@(x,y)pdf(obj1,[x y]),get(gca,{'XLim','YLim'}));
%ezcontour(@(x,y)pdf(obj,[x y]),[0 0.2], [0 1.5])%get(gca,{'XLim','YLim'}));
title('{\bf Scatter Plot and Fitted Gaussian Mixture Contours}')
%legend(h,'Model 0','Model1')
hold off