%% Load data and make variables

load('../Data/KjonnOgHoyde.mat')

kjonn=KjonnOgHoyde(:,1);
hogde=KjonnOgHoyde(:,2);


% Make vectors with hight of men and women

indexMann=find(kjonn==0);
indexKvinne=find(kjonn==1);
hogdeMann=hogde(indexMann);
hogdeKvinne=hogde(indexKvinne);

antkvinner=length(indexKvinne)
antmenn=length(indexMann)


%% Hypotesetest H0 mu=172 mot H1 mu<172, antar kjent varians, N=antkvinner
N=antkvinner


mu0=172;
sigma=6;
alpha=0.05;


ourObs=mean(hogdeKvinne)

meanSigma=sigma/sqrt(N)

pvalue = normcdf(ourObs,mu0,meanSigma)



%% Som forkastningsomraade

cutoff = norminv(alpha, mu0, sigma/sqrt(N))

%% Som z-observator

ourZ=(ourObs-mu0)/(meanSigma)
cutoffZ = -norminv(1-alpha, 0, 1)

%% Hypotesetest H0 mu=172 mot H1 mu<172, antar ukjent varians, N=antkvinner

%Estimert varians til obs og gj.snitt
s=std(hogdeKvinne)
meanSigmaEst=s/sqrt(N)

%% Som t-observator
ourT=(ourObs-mu0)/(meanSigmaEst)
cutoffT=-tinv(1-alpha, N-1)



