clc; clear all
load data_test
x = X1
obem = length(x)
sum_elem = sum(x)
min_elem = min(x)
max_elem = max(x)
razmah = max_elem - min_elem
stepen_svobody = obem-1
sred_arifm = mean(x)
Med = median(x)
Disp = var(x)
Koef = std(x)
Assim = skewness(x)
Eccess = kurtosis(x)

p = sred_arifm;
t = tinv(1-p/2,v);
Mxd = [Mx-Sx*t/sqrt(n), Mx+Sx*t/sqrt(n)]';
disp('Доверительный интервал математического ожидания');
fprintf('%8.4f < MX < %8.4f\n',Mxd)


p = Disp;
chi2l = chi2inv(1-p/2,v);
chi2r = chi2inv(p/2,v);
Dxd = [v*Dx./chi2l,v*Dx./chi2r]';
disp('Доверительный интервал для генеральной дисперсии');
fprintf('%8.4f < DX <%8.4f\n',Dxd)


Da = 6*(obem-1)/(obem+1)/(obem+3);
De = 24*obem*(obem-2)*(obem-3)/(obem+1)^2/(obem+3)/(obem+5);
Axd = [Ax-sqrt(Da./p),Ax+sqrt(Da./p)]';
disp('Доверительный интервал для генеральной асимметрии')
fprintf('%8.4f <Ax< %8.4f\n',Axd)
Exd =[Ex-sqrt(De./p),Ex+sqrt(De./p)]';
disp('Доверительный интервал для генерального эксцесса')
fprintf('%8.4f <Ex< %8.4f\n',Exd


Q1 = prctile(x,25)
Q2 = prctile(x,50)
Q3 = prctile(x,75)
h = boxplot(x,'notch','on');
hold on
plot(1,Mx,'ok','LineWidth',2,'MarkerSize',8)
set(h(1:7),'LineWidth',2)%7 возм. нужно исправить
set(gca,'Box','on','LineWidth',2) title('Блочная диаграмма','FontSize',14)
ylabel('{\bf Значения элементов выборки}','FontSize',12)
text(1.1,Q1,'{\bf Q{1}-25%}','FontSize',12,'Color','b')
text(1.1,Q2,'{\bf Q{2}-50% медиана}','FontSize',12,'Color','r')
text(1.1,Q3,'{\bf Q_{3}-75%}','FontSize',12,'Color','b')
text(0.75,Mx,'{\bf Среднее}','FontSize',12,'Color','k')
grid on