load('../Data/StudData2012.txt')

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

indexKvinne=find(kjonn==2);
hogdeKvinne=hogde(indexKvinne);


mean(hogdeKvinne)
var(hogdeKvinne)
std(hogdeKvinne)

alfa=0.05;
n1=40;
n2=5;


%Antar kjent varians/standaravvik:  sigma = 6
sigma=6;
z_alfahalve=norminv(1-alfa/2, 0,1);

%Case 1
nedre1=mean(hogdeKvinne)-z_alfahalve*sigma/sqrt(n1);
ovre1=mean(hogdeKvinne)+z_alfahalve*sigma/sqrt(n1);

KI1Z=[nedre1,ovre1]

%Case 2
nedre2=mean(hogdeKvinne)-z_alfahalve*sigma/sqrt(n2);
ovre2=mean(hogdeKvinne)+z_alfahalve*sigma/sqrt(n2);

KI2Z=[nedre2,ovre2]




%T-test: 'Estimerer varians', antar at estimatet blir 36, for å 
%få samanliknbare resultat
sigma=6;
%Case 1
nedre1=mean(hogdeKvinne)-tinv(1-alfa/2,n1-1)*sigma/sqrt(n1);
ovre1=mean(hogdeKvinne)+tinv(1-alfa/2,n1-1)*sigma/sqrt(n1);

KI1T=[nedre1, ovre1]

%Case 2
nedre2=mean(hogdeKvinne)-tinv(1-alfa/2,n2-1)*sigma/sqrt(n2);
ovre2=mean(hogdeKvinne)+tinv(1-alfa/2,n2-1)*sigma/sqrt(n2);

KI2T=[nedre2, ovre2]

%For menn: Antar kjent varians/standaravvik:  sigma = 6
sigma=6;
z_alfahalve=norminv(1-alfa/2, 0,1);
n3=90;
%Case 1
nedre1=mean(hogdeMann)-z_alfahalve*sigma/sqrt(n1);
ovre1=mean(hogdeMann)+z_alfahalve*sigma/sqrt(n1);
KI1ZMann=[nedre1, ovre1]
