Boxplot/Box_24.m at main · DeepanJayaraman/Boxplot · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
function F = Box_24(X1,X2,X3,X4)
% X1:4 - data; Name1 - name of the distribution

group10 = [zeros(length(X1(:,1)),1)',ones(length(X1(:,2)),1)',...
    zeros(length(X1(:,1)),1)'+2,zeros(length(X1(:,2)),1)'+3,...
    zeros(length(X1(:,1)),1)'+4,zeros(length(X1(:,2)),1)'+5];
group50 = [zeros(length(X2(:,1)),1)'+6,zeros(length(X2(:,2)),1)'+7,...
    zeros(length(X2(:,1)),1)'+8,zeros(length(X2(:,2)),1)'+9,...
    zeros(length(X2(:,1)),1)'+10,zeros(length(X2(:,2)),1)'+11];
group90 = [zeros(length(X3(:,1)),1)'+12,zeros(length(X3(:,2)),1)'+13,...
    zeros(length(X3(:,1)),1)'+14,zeros(length(X3(:,2)),1)'+15,...
    zeros(length(X3(:,1)),1)'+16,zeros(length(X3(:,2)),1)'+17];
group99 = [zeros(length(X4(:,1)),1)'+18,zeros(length(X4(:,2)),1)'+19,...
    zeros(length(X4(:,1)),1)'+20,zeros(length(X4(:,2)),1)'+21,...
    zeros(length(X4(:,1)),1)'+22,zeros(length(X4(:,2)),1)'+23];


group =[group10 group50 group90 group99];

D = 0.25;
positions = [D*5 D*6 D*7 D*8 D*9 D*10 D*12 D*13 D*14 D*15 D*16 D*17 D*19 D*20 D*21 D*22 D*23 D*24 D*26 D*27 D*28 D*29 D*30 D*31];

X1 = reshape(X1,[],1);X2 = reshape(X2,[],1);X3 = reshape(X3,[],1); X4 = reshape(X4,[],1);


X = [X1' X2' X3' X4'];

% difference
figure;

hB = boxplot(X,group, 'positions', positions,'Color','kkkkkkkkkkkkkkkkkk','symbol','');
h = findobj(gca,'Tag','Box');

Hatch_fill2(h,[1,7,13,19],'b','single',90,'w',1,1)
Hatch_fill2(h,[2,8,14,20],'b','cross',90,'w',1,0.5)
Hatch_fill2(h,[3,9,15,21],'m','single',45,'w',1,1)
Hatch_fill2(h,[4,10,16,22],'m','cross',45,'w',1,1)
Hatch_fill2(h,[5,11,17,23],'r','single',135,'w',1,1)
Hatch_fill2(h,[6,12,18,24],'r','fill',0,'w',1,1)

hB = boxplot(X,group, 'positions', positions,'Color','kkkkkkkkkkkkkkkkkkkkkkkk','symbol','');
lw = 1.5;
set(gca,'fontsize',14,'FontWeight','bold','FontName', 'Times');
set(hB,'linew',lw);

set(gca,'xTick',[1.875,3.625,5.375,7.125])
% Tick = cellstr(num2str(Percentile));
set(gca,'xTickLabel',{'10','15','20','25'})


ylabel('\bf\it{R_P}','fontsize',14,'FontWeight',...
    'bold','FontName', 'Times');

% title(strcat(Name1));
set(gca,'fontsize',14,'FontWeight','bold','FontName', 'Times')
xlabel('Sample size','fontsize',14,'FontWeight','bold','FontName', 'Times')
% xtickangle(45)

[hLegend,objh,~] = legend(findall(gca,'Tag','Box'),...
    {'C-Moment: Load - 2.5 kN','L-Moment: Load - 2.5 kN',...
    'C-Moment: Load - 4 kN','L-Moment: Load - 4 kN',...
    'C-Moment: Load - 7 kN','L-Moment: Load - 7 kN','Ideal'},'numcolumns',2);

% set(objh,'linewidth',1.5);
% hChildren = findall(get(hLegend,'Children'), 'Type','Line');
set(gca,'fontsize',14,'FontWeight','bold','FontName', 'Times')

for j = 1:24
    lines = findobj(gcf, 'type', 'line', 'Tag', 'Median');
    set(lines(j), 'Color', 'w','LineStyle','-','LineWidth',1.5);
end
ylim([0.9 1.6])

line([2.75,2.75],ylim,'Color','k','LineWidth',0.5,'LineStyle','--')
line([4.5,4.5],ylim,'Color','k','LineWidth',0.5,'LineStyle','--')
line([6.25,6.25],ylim,'Color','k','LineWidth',0.5,'LineStyle','--')
line(xlim,[1,1],'Color','k','LineWidth',1.5,'LineStyle',':')