%Регресионный анализ, динамические ряды данных
clc;
cla reset;
clear;
t=[1 2 3 4 5 6 7 8 9 10	11 12 13 14 15 16 17 18];%вариант 1
Yt=[10 14 21 24 33 41 44 47 49 43 47 50 48 54 57 61 59 65];
%параметры регрессии
m=1;
n=size(t);
format bank;
disp('Вариант #1 коэф корреляции')
Koeff_regr=corrcoef(t, Yt)
Whichstats={'tstat','fstat','dwstat'};

stats=regstats(Yt',t','linear',Whichstats);
disp('Дисперсионный анализ');
fstat=stats.fstat %расчетная эф статистика
F_tabl=finv(0.95,m,n(2)-m-1) %табличное значение f
tstat=stats.tstat;%расчетная эф статистика
T_tabl=tinv(0.95,n(2)-m-1) %табличное значение tстьюдента
CoeffTable=dataset({tstat.beta,'Coef'},{tstat.se,'StdEr'},{tstat.t,'tStat'})
kritDW=0;
kritF=0;
dwstat=stats.dwstat;
%Стьюдента
if abs(tstat.t)>T_tabl  
    disp('Модель адекватна, критерий стьюдента выполняется');
    kritDW=1;
else
    disp('Модель не адекватна, критерий стьюдента не выполняется');
    kritDW=0;
end
%КН1
sumznam=0;
sumchis=0;
for j=1:1:n(2)
    sumznam=(Yt(i)^2)+sumznam;
    sumchis=((Yt(i)-Yrash(i))^2)+sumchis;
end
KH1=sqrt(sumchis/sumznam)

    
