parcorr

Частичная автокорреляция выборки

Описание

пример

parcorr(y) строит графики выборки функции частичной автокорреляции ( PACF) одномерных стохастических временных рядов y с доверительными границами.

пример

parcorr(y,Name,Value) использует дополнительные опции, заданные одним или несколькими аргументами пары "имя-значение". Для примера, parcorr(y,'NumLags',10,'NumSTD',2) строит график выборки PACF y для 10 отстает и отображает доверительные границы, состоящие из 2 стандартные ошибки.

пример

pacf = parcorr(___) возвращает выборку PACF y использование любого из входных параметров в предыдущих синтаксисах.

пример

[pacf,lags,bounds] = parcorr(___) дополнительно возвращает числа задержек, которые MATLAB® использует для вычисления PACF, а также возвращает приблизительные верхние и более низкие доверительные границы.

parcorr(ax,___) графики на осях, заданных ax вместо текущих систем координат (gca). ax может предшествовать любой комбинации входных аргументов в предыдущих синтаксисах.

[pacf,lags,bounds,h] = parcorr(___) строит график выборки PACF y и дополнительно возвращает указатели на графические объекты. Используйте элементы h для изменения свойств графика после его создания.

Примеры

свернуть все

Задайте модель AR (2):

yt=0.6yt-1-0.5yt-2+εt,

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

rng(1); % For reproducibility
Mdl = arima('AR',{0.6 -0.5},'Constant',0,'Variance',1)
Mdl = 
  arima with properties:

     Description: "ARIMA(2,0,0) Model (Gaussian Distribution)"
    Distribution: Name = "Gaussian"
               P: 2
               D: 0
               Q: 0
        Constant: 0
              AR: {0.6 -0.5} at lags [1 2]
             SAR: {}
              MA: {}
             SMA: {}
     Seasonality: 0
            Beta: [1×0]
        Variance: 1

Симулируйте 1000 наблюдений за Mdl.

y = simulate(Mdl,1000);

Вычислите PACF.

[partialACF,lags,bounds] = parcorr(y,'NumAR',2);
bounds
bounds = 2×1

    0.0632
   -0.0632

bounds отображения (-0.0633, 0.0633), которые являются верхней и нижней доверительными границами.

Постройте график PACF.

parcorr(y)

Figure contains an axes. The axes with title Sample Partial Autocorrelation Function contains 4 objects of type stem, line.

PACF отрезает после второй задержки. Это поведение указывает на процесс AR (2).

Задайте мультипликативную сезонную ARMA (2,0,1)×(3,0,0)12 модель:

(1-0.75L-0.15L2)(1-0.9L12+0.75L24-0.5L36)yt=2+εt-0.5εt-1,

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

Mdl = arima('AR',{0.75,0.15},'SAR',{0.9,-0.75,0.5},...
    'SARLags',[12,24,36],'MA',-0.5,'Constant',2,...
    'Variance',1);

Моделируйте данные из Mdl.

rng(1);
y = simulate(Mdl,1000); 

Постройте график функции частичной автокорреляции (PACF) по умолчанию.

figure
parcorr(y)

Figure contains an axes. The axes with title Sample Partial Autocorrelation Function contains 4 objects of type stem, line.

Коррелограмма по умолчанию не отображает структуру зависимости для более высоких лагов.

Постройте график PACF для 40 лагов.

figure
parcorr(y,'NumLags',40)

Figure contains an axes. The axes with title Sample Partial Autocorrelation Function contains 4 objects of type stem, line.

Коррелограмма показывает большие корреляции в лагах 12, 24 и 36.

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

свернуть все

Наблюдались одномерные временные ряды, для которых программное обеспечение вычисляет или строит график PACF, заданный как вектор. Последний элемент y содержит последнее наблюдение.

Задайте отсутствующие наблюдения используя NaN. parcorr функция обрабатывает отсутствующие значения как отсутствующие полностью случайным образом.

Типы данных: double

Оси, на которых нужно построить график, заданные как Axes объект.

По умолчанию, parcorr графики для текущей системы координат (gca).

Аргументы в виде пар имя-значение

Задайте необязательные разделенные разделенными запятой парами Name,Value аргументы. Name - имя аргумента и Value - соответствующее значение. Name должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке Name1,Value1,...,NameN,ValueN.

Пример: parcorr(y,'NumLags',10,'NumSTD',2) строит график выборки PACF y для 10 отстает и отображает доверительные границы, состоящие из 2 стандартные ошибки.

Количество лагов в выборке PACF, заданное как разделенная разделенными запятой парами, состоящая из 'NumLags' и положительное целое число. parcorr использует лаги 0:NumLags для оценки PACF.

Значением по умолчанию является min ([20, T – 1]), где T является эффективным размером выборки y.

Пример: parcorr(y,'Numlags',10) строит график выборки PACF y для лагов 0 через 10.

Типы данных: double

Количество лагов в теоретической AR- модели y, заданная как разделенная разделенными запятой парами, состоящая из 'NumAR' и неотрицательное целое число, меньше NumLags.

parcorr использует NumAR для оценки доверительных границ. Для лагов > NumAR, parcorr принимает, что y - гауссов процесс белого шума длины n. Следовательно, стандартная ошибка приблизительно 1/T, где T - эффективный размер выборки y.

Пример: parcorr(y,'NumAR',10) задает, что y является AR (10) обрабатывает, и строит доверительные границы для всех лагов, больше 10.

Типы данных: double

Количество стандартных ошибок в доверительных границах, заданное как разделенная разделенными запятой парами, состоящая из 'NumSTD' и неотрицательный скаляр. Для всех лагов > NumARграницы доверия равны 0 ± NumSTD*σ^, где σ^ - оцененная стандартная ошибка частичной автокорреляции выборки.

Значение по умолчанию дает приблизительно 95% доверительных границ.

Пример: parcorr(y,'NumSTD',1.5) строит графики PACF y с доверительными границами 1.5 стандартные ошибки на расстоянии от 0.

Типы данных: double

Метод оценки PACF, заданный как разделенная разделенными запятой парами, состоящая из 'Method' и значение в этой таблице.

ЗначениеОписаниеОграничения
'ols'Обычные наименьшие квадраты (OLS)y должен быть полностью наблюдаемым рядом (то есть не содержать никаких NaN значения)
'yule-walker'Уравнения Юла-УокераНичего

Если y - это полностью наблюдаемый ряд, тогда по умолчанию это 'ols'. В противном случае значение по умолчанию является 'yule-walker'.

Пример: parcorr(y,'Method','yule-walker') оценивает PACF y используя уравнения Юла-Уокера, а затем строит график PACF.

Типы данных: char | string

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

свернуть все

Выборка PACF одномерных временных рядов y, возвращенный как числовой вектор длины NumLags + 1.

Элементы pacf соответствуют лагам 0,1,2,..., NumLags (то есть элементы lags). Для всех временных рядов y, частичная автокорреляция с задержкой 0 pacf(1) = 1.

Числа задержек, используемые для оценки PACF, возвращенные как числовой вектор длины NumLags + 1.

Аппроксимируйте верхние и более низкие частичные доверительные границы автокорреляции, принимая y является AR (NumAR) обрабатывает, возвращается как двухэлементный числовой вектор.

Указатели на графические объекты, возвращенные как графический массив. h содержит уникальные идентификаторы графика, которые можно использовать для запроса или изменения свойств графика.

Подробнее о

свернуть все

Функция частичной автокорреляции

partial autocorrelation function измеряет корреляцию между yt и y t + k после корректировки на линейные эффекты y t + 1,..., y t + k - 1.

Оценка PACF включает решение уравнений Юла-Уокера относительно автокорреляций. Однако, если временные ряды полностью соблюдаются, то PACF может быть оценен путем подбора последовательных авторегрессивных моделей порядков 1 , 2,... использование обычных наименьших квадратов. Для получения дополнительной информации см. [1], Глава 3.

Отсутствующие полностью случайным образом

Наблюдения случайной переменной missing completely at random, если тенденция наблюдения к отсутствию независима как от случайной переменной, так и от тенденции всех других наблюдений к отсутствию.

Совет

Чтобы построить график PACF без доверительных границ, установите 'NumSTD',0.

Алгоритмы

parcorr строит график PACF, когда вы не запрашиваете никаких выходов или когда вы запрашиваете четвертый выход.

Ссылки

[1] Box, G. E. P., G. M. Jenkins, and G. C. Reinsel. Анализ временных рядов: прогнозирование и управление. 3-й эд. Englewood Cliffs, Нью-Джерси: Prentice Hall, 1994.

[2] Гамильтон, Дж. Д. Анализ временных рядов. Princeton, NJ: Princeton University Press, 1994.

Представлено до R2006a