% Nikolas Tuchyna,318159,2.,BŠ,3ZEH20

% Vytvorenie úložiska
dt = datastore('airlinesmall.csv', 'TreatAsMissing', 'NA', 'MissingValue', 0);

% Načítanie dát
data = read(dt);

% Vytvorenie indexu pre riadky so zhodnými destináciami odletu 'DFW' a príletu 'AUS'
logicalIndex = strcmp(data.Origin, 'PHX') & strcmp(data.Dest, 'LAX');
data2 = data(logicalIndex, :); % Vyber iba riadky, ktoré spĺňajú podmienku zo vstupnej tabuľky

% Vstup od užívateľa pre výber cesty
disp('Vyberte možnosť:');
disp('1 pre "Meskanie"');
disp('2 pre "Meskanie pri vzlete"');
disp('3 pre "Cas letu"');
cesta = input('Zadajte číslo pre voľbu: ');

switch cesta
    case 1
        hodnota = 'ArrDelay';
        mes = data2.(hodnota)';
    case 2
        hodnota = 'DepDelay';
        mes = data2.(hodnota)';
    case 3
        hodnota = 'ArrTime';
        mes = data2.(hodnota)';
    otherwise
        disp('Neplatná voľba');
        return;
end

% Bodový odhad pre strednú hodnotu
mean_est = mean(mes);

% Bodový odhad pre rozptyl
var_est = var(mes);

% Bodový odhad pre smerodajnú odchylku
std_est = std(mes);

% Min a max hodnota vektora mes
min_mes = min(mes);
max_mes = max(mes);

% Počet prvkov v mes
n = length(mes);

% Statické hodnoty
priemer = mean_est; % Statická hodnota pre priemernú hodnotu
dlzka_vektora = length(mes); % Statická hodnota pre dĺžku vektora

% Intervalový odhad strednej hodnoty pre alpha = 5%
alpha = 0.05;
z = (1 - alpha / 2);
mean_conf = mean_est + z * std_est / sqrt(n) * [-1, 1];

% Intervalový odhad pre rozptyl
dolny_rozptyl = (n-1)*var_est/chi2inv(0.975,n-1);
horny_rozptyl = (n-1)*var_est/chi2inv(0.025,n-1);

% Intervalový odhad pre smerodajnú odchylku
dolny_odhad = sqrt(dolny_rozptyl)
horny_odhad = sqrt(horny_rozptyl);

% Výpis do konzoly
fprintf('Stĺpec: %s\n', hodnota);
fprintf('Priemer: %f\n', priemer);
fprintf('Dĺžka vektora: %d\n', dlzka_vektora);
fprintf('Min hodnota vektora mes: %f\n', min_mes);
fprintf('Max hodnota vektora mes: %f\n', max_mes);
fprintf('Bodový odhad pre strednú hodnotu: %f\n', mean_est);
fprintf('Bodový odhad pre rozptyl: %f\n', var_est);
fprintf('Bodový odhad pre smerodajnú odchylku: %f\n', std_est);
fprintf('Intervalový odhad strednej hodnoty pre alpha = 5%%: [%f, %f]\n', mean_conf);
fprintf('Intervalový odhad pre rozptyl: [%f,%f]\n', dolny_rozptyl,horny_rozptyl);
fprintf('Intervalový odhad pre smerodajnú odchylku: [%f,%f]\n', dolny_odhad,horny_odhad);









































































































































































