我有以下问题:我需要构建数据的散点图。一切都很好,但那里有一些重复的数据:
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),或者可能位于该点内。
知道如何实现这个目标吗?
答案 0 :(得分:8)
你可以很容易地做到这一点,让我们把它分成两部分:
首先,您需要确定唯一的2d点并计算它们。这就是我们拥有unique和accumarray功能的原因。如果您没有立即了解他们正在做什么以及他们有什么输出,请仔细阅读文档:
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(Auniq(:,1), Auniq(:,2), [], cnt);
colormap(jet(max(cnt))); % just for the looks of it