1 function p = plot_mu_frequency(mu_values,params)
2 %
function p = plot_mu_frequency(mu_values,params)
4 %
function plotting the mu_values (columnwise 2D/3D vectors) as a
5 % point distribution, the colors and size indicating the frequency
8 % required fields of params:
9 % mu_names: cell array of the names of the components of mu_values
11 % determine frequency vectors in a list
13 % marker sizes minimum and maximum
17 % eliminate possible NaN entries in list:
19 while ~isempty(find(isnan(mu_values(:,1))))
20 mu_values = [mu_values(:,2:end)];
24 if (length(params.mu_names) == 2)
25 mu_values = [mu_values; zeros(1,size(mu_values,2))];
26 params.mu_names = [params.mu_names, {
''}];
30 [mu_single, mu_frequency] = count_column_frequency(mu_values);
32 %l = plot3(mu_values(1,:),mu_values(2,:),mu_values(3,:),
'.');
34 cmax = max(mu_frequency);
43 j = find(mu_frequency==i);
45 plot_args = [plot_args, ...
46 {mu_single(1,j)},{mu_single(2,j)},{mu_single(3, ...
51 % first plot all black
52 %p = plot3(plot_args{:}), hold on;
54 % wi = (mmax*((pfreq(i)-1)/(cmax-1)) + mmin*(1-(pfreq(i)-1)/(cmax-1)))*1.1;
55 % set(p(i),
'Markersize',wi,
'Color',cm(1,:));
57 % then plot all with gray shades
58 p = plot3(plot_args{:});
61 wi = mmax*((pfreq(i)-1)/(cmax-1)) ...
62 + mmin*(1-(pfreq(i)-1)/(cmax- 1));
66 set(p(i),'Markersize',wi,'MarkerFaceColor',cm(pfreq(i),:),'Marker','o',...
67 'MarkerEdgeColor',[0 0 0]);
73 set(c,'Limits',[0.5,cmax+0.5]);
76 xlabel(params.mu_names{1});
77 ylabel(params.mu_names{2});
78 zlabel(params.mu_names{3});
84 % TO BE ADJUSTED TO NEW SYNTAX