File:Papier weibull 3 cycles log dense.svg

clear;
fig = scf(1);
fig.figure_size = [630, 891];
clf;

function [y]=transfo(F)
    y = log(-log(1-F));
endfunction

xbase = [1:0.1:4, 4.2:0.2:8, 8.5:0.5:9.5];
xtraits = [xbase, 10*xbase, 100*xbase, 1000];
xetiq = [1:10, 20:10:100, 200:100:1000];

lnxtraits = log(xtraits);
lnxetiq = log(xetiq);
lnxmin = min(lnxtraits);
lnxmax = max(lnxtraits);
nxtraits = size(xtraits, 'c');
nxetiq = size(xetiq, 'c');

Ftraits = [0.01:0.001:0.03, 0.032:0.002:0.06, 0.065:0.005:0.1,...
    0.11:0.01:0.4, 0.4:0.02:0.8, 0.81:0.01:0.9, 0.91:0.01:0.99];
ntraits = size(Ftraits, 'c');
Fetiq = [0.01, 0.02, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.99]
netiq = size(Fetiq, 'c');

ytraits = transfo(Ftraits);
yetiq = transfo(Fetiq);

ymin = min(ytraits);
ymax = max(ytraits);

for i = 1:ntraits
    yy = ytraits(i); FF = Ftraits(i);
    xpoly([lnxmin, lnxmax],[yy, yy]);
    h = gce();
    h.thickness = 0.5;
end

for i = 1:netiq
    yy = yetiq(i); FF = Fetiq(i);
    xpoly([lnxmin, lnxmax],[yy, yy]);
    xstring(-0.6, yy-0.08, string(FF));
    h = gce();
    h.font_size = 0.3;
end

FF = 1-exp(-1);
yy = log(log(1-FF));
xpoly([lnxmin, lnxmax],[yy, yy])
h = gce();
h.thickness = 0.5;
h.line_style=3;
xstring(-0.6, yy-0.08, "0.632");
h = gce();
h.font_size = 0.3;

for i = 1:nxtraits
    xx = lnxtraits(i);
    xpoly([xx, xx], [ymin, ymax])
    h1 = gce();
    h1.thickness = 0.5;
end

for i = 1:nxetiq
    xx = lnxetiq(i);
    xstring(xx-0.15, ymin-0.05, string(xetiq(i)))
    h = gce();
    h.font_size = 0.3;
    h.font_angle = 90;
end

a = gca();
a.data_bounds=[0, ymin ; 10, ymax];