6

我有以下问题:我需要构建数据的散点图。一切都很好,但那里有一些重复的数据:

x = [11, 10, 3, 8, 2, 6, 2, 3, 3, 2, 3, 2, 3, 2, 2, 2, 3, 3, 2, 2];
y = [29, 14, 28, 19, 25, 21, 27, 15, 24, 23, 23, 18, 0, 26, 11, 27, 23, 30, 30, 25];

可以看到有两个元素(2, 25); (2,27); (3,24); 因此如果使用常规构建此数据,scatter(x,y)我将丢失此信息: 在此处输入图像描述

我发现的解决方法是使用未记录的'jitter'参数

scatter(x,y, 'jitter','on', 'jitterAmount', 0.06);

但我不喜欢这样的前景: 在此处输入图像描述

我试图实现的是:

在此处输入图像描述

重复的数量在点旁边(如果数量大于 1),或者可能在点内。

知道如何实现这一目标吗?

4

1 回答 1

8

你可以很容易地做到这一点,让我们把它分成两部分:

首先,您需要识别唯一的 2d 点并计算它们。这就是我们拥有uniqueaccumarray函数的目的。如果您不立即了解他们在做什么以及他们有什么输出,请通读文档:

x = [11 10 3  8  2  6  2  3  3  2  3  2  3  2  2  2  3  3  2  2];
y = [29 14 28 19 25 21 27 15 24 23 23 18 0  26 11 27 23 30 30 25];
A=[x' y'];

[Auniq,~,IC] = unique(A,'rows');
cnt = accumarray(IC,1);

现在每一行都Auniq包含唯一的 2d 点,同时cnt包含每个点的出现次数:

>> [cnt Auniq]

ans =

     1     2    11
     1     2    18
     1     2    23
     2     2    25
     1     2    26
     ...etc

对于显示出现次数,有很多可能性。正如您所提到的,您可以将数字放在散布标记的内部/旁边,其他选项是颜色编码、标记的大小......让我们做所有这些,您当然也可以组合!

标记旁边的数字

scatter(Auniq(:,1), Auniq(:,2));
for ii=1:numel(cnt)
    if cnt(ii)>1
        text(Auniq(ii,1)+0.2,Auniq(ii,2),num2str(cnt(ii)), ...
            'HorizontalAlignment','left', ...
            'VerticalAlignment','middle', ...
            'FontSize', 6);
    end
end
xlim([1 11]);ylim([0 30]);

在此处输入图像描述

标记内的数字

scatter(Auniq(:,1), Auniq(:,2), (6+2*(cnt>1)).^2); % make the ones where we'll put a number inside a bit bigger

for ii=1:numel(cnt)
    if cnt(ii)>1
        text(Auniq(ii,1),Auniq(ii,2),num2str(cnt(ii)), ...
            'HorizontalAlignment','center', ...
            'VerticalAlignment','middle', ...
            'FontSize', 6);
    end
end

正如你所看到的,我使用 scatter 函数本身非常简单地放大了标记的大小。

在此处输入图像描述

颜色编码

scatter(Auniq(:,1), Auniq(:,2), [], cnt);
colormap(jet(max(cnt))); % just for the looks of it

在此处输入图像描述

之后,您可以添加颜色条或图例以指示每种颜色的出现次数。

于 2012-12-08T16:06:32.593 回答