Добавить маркер в гистограмму с накоплением - PullRequest
Новая
       17

Добавить маркер в гистограмму с накоплением

0 голосов
/ 13 мая 2019

с некоторой помощью я смог реализовать сложенную гистограмму так, как мне нужно.Отсутствует только одна часть: маркер, который представляет медиану.К сожалению, я не мог придумать какие-либо идеи о том, как добавить один и не мог найти какую-либо помощь в Интернете.Эта картинка показывает мою диаграмму и то, как я хочу добавить маркер так, как он показан.Это не должно быть так, swaure или что-то еще сделает эту работу тоже.

without marker

with marker

код для создания графика: (breite - это лист Excel, который я импортировал. Вы можете использовать testw и удалить testw = breite для проверки кода) 

% testw = [0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
%     0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ];
testw = breite;

%Datenvorbereitung
gerundet = round(testw,1);

maximum = max(gerundet),1;
anzahl = ceil(maximum / 0.1);

werte = zeros([1 anzahl]);
sz = size(werte);
haeufigkeit = zeros(sz);

for i= 1:anzahl
    werte(i) = 0.1*i;
end

wertebereich = zeros(sz);

for i = 1:anzahl
    wertebereich(i) = 0.1;
end

%Häufigkeitsverteilung
for i = 1:10000
    if gerundet(i) <= 0.2
        haeufigkeit(1) = haeufigkeit(1) + 1
    else
        for j = 2:anzahl
            if gerundet(i) >= werte(j-1) && gerundet(i) < werte(j)
                haeufigkeit(j) = haeufigkeit(j) + 1
            end
        end
    end
end

%Grauwerte
prozentual = zeros(sz);

for i = 1:anzahl
    prozentual(i) = haeufigkeit(i)/sum(haeufigkeit)
end

grauwerte = zeros(sz);

for i = 1:anzahl
    grauwerte(i) = round(1-1*prozentual(i),3);
end

colormap gray
b = bar([wertebereich; nan(size(wertebereich))],'stacked');
% colorbar('Direction','reverse')
colorbar('Ticks',[0, 0.2, 0.4, 0.6 ,0.8, 1],...
         'TickLabels',{'100%','80%','60%','40%','20%', '0%'})

for i = 1:anzahl
     b(1,i).FaceColor = [grauwerte(i) grauwerte(i) grauwerte(i)];
     b(1,i).EdgeColor = 'none';
end

1 Ответ

2 голосов
/ 13 мая 2019

Если я правильно понимаю, просто добавьте что-то вроде этого после вашего кода: 

hold on
med = median(gerundet);
plot([.5 1.5], [med med], '-', 'color', [.2 .6 1], 'linewidth', 2)

enter image description here

...