подвести итог

Класс: regARIMA

Отобразите результаты оценки модели регрессии с ошибками ARIMA

Синтаксис

summarize(Mdl)
results = summarize(Mdl)

Описание

пример

summarize(Mdl) отображает сводные данные модели регрессии с ошибками ARIMA Mdl.

  • Если Mdl является предполагаемой моделью, возвращенной estimate, то summarize распечатывает результаты оценки к Командному окну MATLAB®. Отображение включает сводные данные оценки и таблицу оценок параметра с соответствующими стандартными погрешностями, статистикой t и p - значения. Сводные данные оценки включают подходящую статистику, такую как Критерий информации о Akaike (AIC) и предполагаемое инновационное отклонение.

  • Если Mdl является непредполагаемой моделью, возвращенной regARIMA, то summarize распечатывает отображение стандартного объекта (то же отображение, которое regARIMA распечатывает во время образцового создания).

пример

results = summarize(Mdl) возвращает одну из следующих переменных и не распечатывает к Командному окну.

  • Если Mdl является предполагаемой моделью, то results является структурой, содержащей результаты оценки.

  • Если Mdl является непредполагаемой моделью, то results является объектом модели regARIMA, который равен Mdl.

Входные параметры

развернуть все

Модель Regression с ошибками ARIMA, заданными как объект модели regARIMA, возвращенный estimate или regARIMA.

Выходные аргументы

развернуть все

Образцовые сводные данные, возвращенные как массив структур или объект модели regARIMA.

  • Если Mdl является предполагаемой моделью, то results является массивом структур, содержащим поля в этой таблице.

    Поле Описание
    DescriptionОбразцовое итоговое описание (строка)
    SampleSizeЭффективный объем выборки (числовой скаляр)
    NumEstimatedParametersКоличество предполагаемых параметров (числовой скаляр)
    LogLikelihoodОптимизированное loglikelihood значение (числовой скаляр)
    AICКритерий информации о Akaike (числовой скаляр)
    BICБайесов информационный Критерий (числовой скаляр)
    TableОценки наибольшего правдоподобия параметров модели с соответствующими стандартными погрешностями, статистика t (оценка, разделенная на стандартную погрешность), и p - значения (принимающий нормальность); таблица со строками, соответствующими параметрам модели

  • Если Mdl является непредполагаемой моделью, то results является объектом модели regARIMA, который равен Mdl.

Примеры

развернуть все

Регресс валовой внутренний продукт (ВВП) США на индекс потребительских цен (CPI) США с помощью модели регрессии с ARMA (1,1) ошибки, и обобщает результаты.

Загрузите США Макроэкономический набор данных и предварительно обработайте данные.

load Data_USEconModel;
logGDP = log(DataTable.GDP);
dlogGDP = diff(logGDP);
dCPI = diff(DataTable.CPIAUCSL);

Соответствуйте модели к данным.

ToEstMdl = regARIMA('ARLags',1,'MALags',1);
EstMdl = estimate(ToEstMdl,dlogGDP,'X',dCPI,'Display','off');

Отобразите оценки.

summarize(EstMdl)
 
   ARMA(1,1) Error Model (Gaussian Distribution)
 
    Effective Sample Size: 248
    Number of Estimated Parameters: 5
    LogLikelihood: 798.406
    AIC: -1586.81
    BIC: -1569.24
 
                   Value       StandardError    TStatistic      PValue  
                 __________    _____________    __________    __________

    Intercept      0.014776      0.0014627        10.102      5.4236e-24
    AR{1}           0.60527        0.08929        6.7787      1.2125e-11
    MA{1}          -0.16165        0.10956       -1.4755         0.14009
    Beta(1)        0.002044     0.00070616        2.8946       0.0037969
    Variance     9.3578e-05     6.0314e-06        15.515      2.7338e-54

 

Оцените несколько моделей путем передачи данных estimate. Отличайтесь степени авторегрессивного и скользящего среднего значения p и q, соответственно. Результаты оценки содержат AIC, который можно извлечь и затем сравнить среди моделей.

Моделируйте ответ и данные о предикторе для модели регрессии с ошибками ARMA:

yt=2+Xt[-21.5]+utut=0.75ut-1-0.5ut-2+εt+0.7εt-1,

где εt является Гауссовым со средним значением 0 и отклонением 1.

Mdl = regARIMA('Intercept',2,'Beta',[-2; 1.5],...
    'AR',{0.75, -0.5},'MA',0.7,'Variance',1);

rng(2);            % For reproducibility
X = randn(1000,2); % Predictors
y = simulate(Mdl,1000,'X',X);

Чтобы определить количество AR и задержек MA, создайте и оцените модели множественной регрессии с ARMA (p, q) ошибки. Отличайтесь p = 1.., 3 и q = 1..., 3 среди моделей. Подавите все отображения оценки. Извлеките AIC от структуры результатов оценки. Поле AIC хранит AIC.

pMax = 3;
qMax = 3;
AIC = zeros(pMax,qMax); % Preallocation

for p = 1:pMax
    for q = 1:qMax
        ToEstMdl = regARIMA(p,0,q);
        EstMdl = estimate(ToEstMdl,y,'X',X,'Display','off');
        results = summarize(EstMdl);
        AIC(p,q) = results.AIC;
    end
end

Сравните значения AIC среди моделей.

minAIC = min(min(AIC))
minAIC = 2.9280e+03
[bestP,bestQ] = find(AIC == minAIC)
bestP = 2
bestQ = 1

Модель оптимальной подгонки является моделью регрессии с AR (2,1) ошибки, потому что его соответствующий AIC является самым низким. Эта модель также имеет структуру модели, используемой, чтобы моделировать данные.

Смотрите также

Объекты

Функции

Введенный в R2018a

Для просмотра документации необходимо авторизоваться на сайте