%% data for vo2max kan ikke legges ut offentlig, s? her er koden med simulerte data

data=normrnd(0,1,1,1000)
%data=exprnd(1,1,1000)
%% lage boksplott
boxplot(data);
%% histogram - med areal 1, og legge til normalkurve
figure
histogram(data,'Normalization','pdf');
hold on
mu=mean(data);
sigma=std(data);
del=8*sigma/1000;
xvec=[mu-3*sigma:del:mu+3*sigma];
pdfpopulasjon=normpdf(xvec,mu,sigma);
plot(xvec, pdfpopulasjon, 'b','LineWidth',3)
set(gca,'FontSize',20)

%% empirisk F, og legge til normal cdf
figure
cdfplot(data)
hold on
plot(xvec,normcdf(xvec,mu,sigma),'color','r','LineWidth',2)

%% Normal probability plott - der observasjoner p?? x-akse
figure
normplot(data);

%% QQplott er aksene er byttet om og fokus er p?? kvantiler
figure
qqplot(data);