%%zad 1
close all;
age = load('3_age.txt');
height = load('3_height.txt');
speed_city = load('3_speed_city.txt');
speed_out = load('3_speed_out.txt');

%
% mean_age = mean(age)
% mean_height = mean(height)
% mean_speed_city = mean(speed_city)
%
%
% mean_age = std(age)
% mean_height = std(height)
% mean_speed_city = std(speed_city)
%
% boxplot(age)%   nic 92
% figure
% boxplot(height)%    nic 194
% figure
% boxplot(speed_city) %28.7244 71.95555
% figure
%
% age = age(age <= 92);
% height = height(height <= 194);
% speed_city = speed_city(speed_city >= 28.7244 & speed_city <= 71.95555);
%
% mean_age = mean(age)
% mean_height = mean(height)
% mean_speed_city = mean(speed_city)
%
% mean_age = std(age)
% mean_height = std(height)
% mean_speed_city = std(speed_city)

%% zad 2

% % speed_city = sort(speed_city);
% % speed_out = sort(speed_out);
% plot(speed_city, speed_out, '.');
%
% wsp_kor = corrcoef(speed_city, speed_out)


%% zad 3
% tab = [
%     1.520     4.020
%       0.710    -0.790
%       1.700     4.390
%       0.700    -0.880
%       0.720    -0.760
%       1.700     4.460
%       0.510    -0.090
%       0.810    -0.790
%       0.970     0.030
%       1.100     2.040
%       1.520     4.140
%       0.110     1.000
%       1.520     4.090
%       0.110     1.030
%       0.310     1.580
%       0.970     0.110
%       0.820    -0.780
%       0.510    -0.080
%       0.810    -0.770
%       1.520     4.060
%       0.320     1.610
%       0.710    -0.780
%       0.960     0.010
%       0.700    -0.770
%       1.100     2.050
%       2.020     5.150
%       1.710     4.340
%       0.800    -0.810
%       0.960     0.050
%       0.320     1.630
%       ]';
% x = tab(1,:);
% y = tab(2,:);
% X=0:0.01:2.5;
% plot(x, y, '.r')
% hold on
%
% a = polyfit(x, y, 1)
% y5 = a(1)*x+a(2);
% plot(x, y5, 'g')
%
%
% a = polyfit(x, y, 2)
% y = a(1)*X.*X + a(2)*X + a(3);
% plot(X, y, 'y')

%% zad 4
% Sr1 =2 od1=3
% Sr2 =5 od2=1
%
% Sr1 =2 od1=2
% Sr2 =5 od2=1

x = -10:0.01:10;

y1 = normpdf(x, 2, 3);
y2 = normpdf(x, 5, 1);
plot(x, y1);
hold on
plot(x,y2)
figure

y3 = normpdf(x, 2, 2);
y4 = normpdf(x, 5, 1);
plot(x, y3);
hold on
plot(x,y4)
figure

[X,Y,T,AUC] = perfcurve(y1,y2,'virginica');
% plot(X, Y)