%%
load('../Data/StudDataV16.txt')
%%
hogde=StudDataV16(:,1);
hogdeFar=StudDataV16(:,2);
hogdeMor=StudDataV16(:,3);
kjonn=StudDataV16(:,4);
har=StudDataV16(:,5);
auge=StudDataV16(:,6);

%% Make vectors with hight of men and women

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

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



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


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))

%% Styrke for mu_1=170
mu1=170
normcdf(cutoff,mu1,sigma/sqrt(N))


%% Som z-observator

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

%% Styrkefunksjon hogde kvinner

mu1 = 167:0.1:174;
y = normcdf(cutoff,mu1,sigma/sqrt(N));
%y5=y
h1 = plot(mu1,y,'LineWidth',3);
set(gca,'FontSize',20);
hold
 plot([167,174],[0.05,0.05],'r-','LineWidth',3)
 
 
 xlabel('\mu_1','FontSize',24)
 ylabel('Styrke: Sannsyn for a forkaste H_0 nar \mu = \mu_1','FontSize',24)
 
 
 
%%
plot(mu1,y5,'g-','LineWidth',3)
h_legend=legend('Styrke Case 1','Signifikansniva \alpha','Styrke Case 2');
set(h_legend,'FontSize',24)

%% Hypotese Menn H0: mu=179.8 H1: mu>179.8

N=antmenn
mu0=179.8;
sigma=6.5;
alpha=0.05;


ourObs=mean(hogdeMann)
meanSigma=sigma/sqrt(N)

pverdi = 1- normcdf(ourObs,mu0,meanSigma)

%% Forkastningsomraade

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

%% Som Z-observator

ourZ=(mean(hogdeMann)-mu0)/(meanSigma)

cutoffZ = norminv(1-alpha, 0, 1)

