clc
clear all
close all

% Matrice dei dati (campioni appaiati)
X = [134 140
     122 130
     132 135
     130 126
     128 134
     140 138
     118 124
     127 126
     125 132
     142 144];

x_pre = X(:,1);
x_post = X(:,2);

%livello di significatività
alpha = 0.05;

% t test per dati appaiati

% Differenze tra i campioni
w = x_pre-x_post;

% Sistema di ipotesi
% H0: mu_w = 0
% H1: mu_w diversa da 0

% Valore osservato per la statistica del test
t_oss = mean(w)/std(w)*sqrt(length(w))

% Valore critico
t_c = tinv(1-alpha/2,length(w)-1)

% Esito finale: poichè |t_oss|> t_c rifiutiamo H0

% P-value
pvalue = 2*(1-tcdf(abs(t_oss),length(w)-1))

% Soluzione con ttest (3 chiamate della function equivalenti nel nostro caso)
[h,p,ci,stat]=ttest(w)
[h,p,ci,stat]=ttest(w,0,'Alpha',alpha,'Tail','both')
[h,p,ci,stat]=ttest(x_pre,x_post)