%% COMPLEMENTO ALLE LEZIONI DI TEORIA - METODI STATISTICI PER LA BIOINGEGNERIA
%% Canale B - A.A. 2026/2027
%% Docente: Enrico Longato

% Costruzione "manuale" del qqplot, per capire cosa ci va sugli assi

clc
clear all
close all

rng(42) % Setto un seed per riproducibilità

% Creo due campioni aleatori estratti da due distribuzioni diverse
N = 200;
gaussiana = 2*randn(200, 1) + 5; % Cambio la media e la SD per farvi vedere che viene uguale
uniforme = rand(200, 1);

% Calcolo i percentili campionari (andranno sull'asse y)
percentuali = 0:1:100;
perc_gaussiana = prctile(gaussiana, percentuali);
perc_uniforme = prctile(uniforme, percentuali);

% Calcolo i quantili teorici della distribuzione normale standard (asse x)
quantili_teorici_gaussiana = norminv(percentuali / 100);

% Faccio i plot con sovrapposta la retta di riferimento
figure
subplot(221)
plot(quantili_teorici_gaussiana, perc_gaussiana, '.')
hold on
plot(quantili_teorici_gaussiana, std(gaussiana)*quantili_teorici_gaussiana + mean(gaussiana))
title('QQ-plot GAUSSIANA')

subplot(222)
qqplot(gaussiana)
title('QQ-plot di MATLAB')

subplot(223)
plot(quantili_teorici_gaussiana, perc_uniforme, '.')
hold on
plot(quantili_teorici_gaussiana, std(uniforme)*quantili_teorici_gaussiana + mean(uniforme))
title('QQ-plot UNIFORME')

subplot(224)
qqplot(uniforme)
title('QQ-plot di MATLAB')