data = dlmread('kroppstemperatur.txt');
 
%% Uppgift 1
%h0: myF = myM
%h1: myF != myM
%testar om datan är normalfördelad, H = 0 ==> a-alpha konfidens att datan
%kommer från normalfördelning
lillietest(data(:,1),'alpha',0.01);
meanTemp = mean(data(:,1));
h0 = 98.6;
[h, p0] = ttest(data(:,1),h0,'Tail','Both','alpha',0.05);
 
if h == 1
    disp('Forkastar h0 for h1')
end
 
%% Uppgift 3
close all
x = data(:,3); %freq
X = [ones(size(x)),x];
y = data(:,1); %temp
 
%y = betaX+alpha
[B,BINT,R,RINT,STATS] = regress(y,X);
yfit=B(1)+B(2)*x;
 
%Plot of observations >> 
hold off  
plot(x,y,'.'); 
%Plot of fit 
hold on 
plot(x,yfit,'r.')
%Plot of regression line  
xmesh=linspace(min(x)-1,max(x)+1,100); 
ymesh=B(1)+B(2)*xmesh; 
plot(xmesh,ymesh,'r-') 
legend({'Observed value','Fitted values','Regressin Line'},'Location','NorthWest') 
xlabel('Heart Freq');
ylabel('Body temp')
display(BINT)
disp('coefficient av korelation och p-värde')
disp(STATS([1,3]))
%% Uppgift 2, column 1 = temp, column 2 = freq
alpha = 0.05;
mData = data(data(:,2) == 1,[1 3]);
fData = data(data(:,2) == 2,[1 3]);
%testar om datan är normalfördelad, H = 0 ==> kan inte förkastas den kommer
%från normalfördelning
%kommer från normalfördelning
NormalantagandeMan=lillietest(mData(:,1),'alpha',alpha);
if NormalantagandeMan == 0
    disp('Cant reject that H_0 is from a Normaldist.' )
end
NormalantagandeFemale=lillietest(fData(:,1),'alpha',alpha);
if NormalantagandeFemale == 0
    disp('Cant reject that H_0 is from a Normaldist.' )
end

% testar om kvinnor och männens temperatur har samma varians med 
% 95 procent konfidensnivå
vt = vartest2(mData(:,1),fData(:,1));
% då de har samma varians
[h0,p0] = ttest2(mData(:,1),fData(:,1),'alpha',alpha,'tail','both','vartype','equal');
if h0 == 1
    disp('Forkastar h0 for h1')
end