
%% 
%Eksempel p? ordningsvariable
%Ser paa hogde for n menn. 

n=25;

% Antar at hogde er normalfordelt med
mu=180
sigma=6;

% Lager vektor vi skal plotte for
del=8*sigma/1000;
xvec=[mu-4*sigma:del:mu+4*sigma];

cumvec=normcdf(xvec,mu,sigma);
pdfvec=normpdf(xvec,mu,sigma);
% Fordelinga til hogde
figure;
clf;
plot(xvec,pdfvec,'k', 'LineWidth',2)
hold
%%
% Minste (minimum)
f1=n*pdfvec.*((1-cumvec).^(n-1));
% Tredje minste 
f3=n*nchoosek(n-1,3-1)*((pdfvec.*cumvec.^(3-1)).*(1-cumvec).^(n-3));
% Median
m=(n+1)/2;
fm=n*nchoosek(n-1,m-1)*((pdfvec.*cumvec.^(m-1)).*(1-cumvec).^(n-m));
% Tredje storste
f3top=n*nchoosek(n-1,3-1)*((pdfvec.*(cumvec.^(n-3))).*(1-cumvec).^(3-1));
% Storste (maksimum)
ftop=n*pdfvec.*(cumvec.^(n-1));

%figure(2);
%clf;
plot(xvec,f1,xvec,f3,xvec,fm,xvec,f3top,xvec,ftop,'LineWidth',2);
h_legend=legend('Hogde','Minste','3. minste','Median','3. storste','Storste');
set(h_legend,'FontSize',16);

