fgls

Выполнимые обобщенные наименьшие квадраты

Описание

пример

coeff = fgls(X,y) возвращает содействующие оценки модели y многофакторной линейной регрессии = Xβ + ε с помощью выполнимых обобщенных наименьших квадратов (FGLS) первой оценкой ковариации инновационного процесса ε.

NaNs в данных указывают на отсутствующие значения, который fgls удаляет использующее мудрое списком удаление. fgls наборы Data = [X y], затем это удаляет любую строку в Data содержа по крайней мере один NaN. Мудрое списком удаление уменьшает эффективный объем выборки и изменяет основу времени ряда.

пример

coeff = fgls(Tbl) возвращает содействующие оценки FGLS с помощью данных о предикторе в первом numPreds столбцы таблицы Tbl и данные об ответе в последнем столбце.

fgls удаляет все отсутствующие значения в Tbl, обозначенный NaNs, с помощью мудрого списком удаления. Другими словами, fgls удаляет все строки в Tbl содержа по крайней мере один NaN. Мудрое списком удаление уменьшает эффективный объем выборки и изменяет основу времени ряда.

пример

coeff = fgls(___,Name,Value) задает опции с помощью одного или нескольких аргументов пары "имя-значение" в дополнение к входным параметрам в предыдущих синтаксисах. Например, можно выбрать инновационную модель ковариации, задать количество итераций и построить оценки после каждой итерации.

пример

[coeff,se,EstCoeffCov] = fgls(___) дополнительно возвращает вектор из содействующих стандартных погрешностей FGLS, se = sqrt(diag(EstCov)), и FGLS оценил содействующую ковариационную матрицу (EstCoeffCov).

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

[coeff,se,EstCoeffCov,iterPlots] = fgls(___) возвращает указатели на нанесенные на график графические объекты. Используйте элементы iterPlots изменить свойства графиков после того, как вы создаете их.

Примеры

свернуть все

Предположим, что чувствительность американского Индекса потребительских цен (CPI) к изменениям в заплаченной компенсации сотрудников (COE) представляет интерес.

Загрузите США макроэкономический набор данных. Постройте серию CPI и COE.

load Data_USEconModel

figure;
subplot(2,1,1)
plot(dates,DataTable.CPIAUCSL);
title '{\bf Consumer Price Index, Q1 in 1947 to Q1 in 2009}';
datetick;
axis tight;
subplot(2,1,2);
plot(dates,DataTable.COE);
title '{\bf Compensation Paid to Employees, Q1 in 1947 to Q1 in 2009}';
datetick;
axis tight;

Figure contains 2 axes objects. Axes object 1 with title blank C o n s u m e r blank P r i c e blank I n d e x , blank Q 1 blank i n blank 1 9 4 7 blank t o blank Q 1 blank i n blank 2 0 0 9 contains an object of type line. Axes object 2 with title blank C o m p e n s a t i o n blank P a i d blank t o blank E m p l o y e e s , blank Q 1 blank i n blank 1 9 4 7 blank t o blank Q 1 blank i n blank 2 0 0 9 contains an object of type line.

Ряды являются неустановившимися. Стабилизируйте их путем применения журнала, и затем первого различия.

CPI = diff(log(DataTable.CPIAUCSL));
COE = diff(log(DataTable.COE));

Регресс CPI на COE включая точку пересечения, чтобы получить оценки обычных наименьших квадратов (OLS). Сгенерируйте изолированный остаточный график.

Mdl = fitlm(COE,CPI);

figure;
plotResiduals(Mdl,'lagged')

Figure contains an axes object. The axes object with title Plot of residuals vs. lagged residuals contains 3 objects of type line.

В остаточном графике существует восходящий тренд, который предлагает, чтобы инновации включили авторегрессивный процесс. Это нарушает одно из классических линейных предположений модели. Следовательно, тесты гипотезы на основе коэффициентов регрессии являются неправильными, даже асимптотически.

Оцените коэффициенты регрессии с помощью FGLS. По умолчанию, fgls включает точку пересечения в модель регрессии и налагает модель AR (1) на инновации. Опционально, отобразите OLS и оценки FGLS путем определения 'final' для 'display' аргумент пары "имя-значение".

coeff = fgls(CPI,COE,'display','final');
OLS Estimates:

       |  Coeff    SE   
------------------------
 Const | 0.0122  0.0009 
 x1    | 0.4915  0.0686 

FGLS Estimates:

       |  Coeff    SE   
------------------------
 Const | 0.0148  0.0012 
 x1    | 0.1961  0.0685 

Если серия COE является внешней относительно CPI, то FGLS оценивает (coeff) сопоставимы и асимптотически более эффективны, чем оценки OLS.

Предположим, что чувствительность американского Индекса потребительских цен (CPI) к изменениям в заплаченной компенсации сотрудников (COE) представляет интерес. Этот пример улучшает анализ, обрисованный в общих чертах в Оценке в качестве примера Коэффициенты FGLS Используя Опции по умолчанию.

Загрузите американский макроэкономический набор данных.

load Data_USEconModel

Ряды являются неустановившимися. Стабилизируйте их путем применения журнала, и затем первого различия.

CPI = diff(log(DataTable.CPIAUCSL));
COE = diff(log(DataTable.COE));

Регресс CPI на COE включая точку пересечения, чтобы получить оценки OLS. Постройте коррелограммы для остаточных значений.

Mdl = fitlm(COE,CPI);
u = Mdl.Residuals.Raw;

figure;
subplot(2,1,1)
autocorr(u);
subplot(2,1,2);
parcorr(u);

Figure contains 2 axes objects. Axes object 1 with title Sample Autocorrelation Function contains 4 objects of type stem, line. Axes object 2 with title Sample Partial Autocorrelation Function contains 4 objects of type stem, line.

Коррелограммы предлагают, чтобы инновации оказали значительные влияния AR. Согласно Методологии Поля-Jenkins, инновации, кажется, включают серию AR (3).

Оцените коэффициенты регрессии с помощью FGLS. По умолчанию, fgls принимает, что инновации авторегрессивны. Укажите, что инновациями является AR (3) использование 'arLags' аргумент пары "имя-значение".

[coeff,se] = fgls(CPI,COE,'arLags',3,'display','final');
OLS Estimates:

       |  Coeff    SE   
------------------------
 Const | 0.0122  0.0009 
 x1    | 0.4915  0.0686 

FGLS Estimates:

       |  Coeff    SE   
------------------------
 Const | 0.0148  0.0012 
 x1    | 0.1972  0.0684 

Если серия COE является внешней относительно CPI, то FGLS оценивает (coeff) сопоставимы и асимптотически более эффективны, чем оценки OLS.

Смоделируйте номинальный GNP (GNPN) темп роста, составляющий эффекты темпов роста индекса потребительских цен (CPI), действительная заработная плата (WR), и денежный запас (MS). Объясните классические линейные отъезды модели.

Загрузите набор данных Нельсона Плоссера.

load Data_NelsonPlosser
varIdx = [8,10,11,2];               % Variable indices
idx = ~any(ismissing(DataTable),2); % Identify nonmissing values 
Tbl = DataTable(idx,varIdx);        % Tabular array of variables
T = sum(idx);                       % Sample size

Постройте ряд.

figure;
for j = 1:4;
    subplot(2,2,j);
    plot(dates(idx),Tbl{:,j});
    title(Tbl.Properties.VariableNames{j});
    axis tight;
end;

Figure contains 4 axes objects. Axes object 1 with title CPI contains an object of type line. Axes object 2 with title WR contains an object of type line. Axes object 3 with title MS contains an object of type line. Axes object 4 with title GNPN contains an object of type line.

Все ряды кажутся неустановившимися.

Примените журнал, и затем первое различие для каждого ряда.

dLogTbl = array2table(diff(log(Tbl{:,:})),...
    'VariableNames',strcat(Tbl.Properties.VariableNames,'Rate'));

Регресс GNPNRate на другие переменные в dLogTbl. Исследуйте график рассеивания и коррелограммы остаточных значений.

Mdl = fitlm(dLogTbl);

figure;
plotResiduals(Mdl,'caseorder');
axis tight;

Figure contains an axes object. The axes object with title Case order plot of residuals contains 2 objects of type line.

figure;
subplot(2,1,1);
autocorr(Mdl.Residuals.Raw);
subplot(2,1,2);
parcorr(Mdl.Residuals.Raw);

Figure contains 2 axes objects. Axes object 1 with title Sample Autocorrelation Function contains 4 objects of type stem, line. Axes object 2 with title Sample Partial Autocorrelation Function contains 4 objects of type stem, line.

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

Оцените коэффициенты FGLS путем объяснения heteroscedasticity остаточных значений. Укажите, что предполагаемая инновационная ковариация является диагональной с квадратами остатков как веса.

fgls(dLogTbl,'innovMdl','HC0','display','final');
OLS Estimates:

         |  Coeff     SE   
---------------------------
 Const   | -0.0076  0.0085 
 CPIRate |  0.9037  0.1544 
 WRRate  |  0.9036  0.1906 
 MSRate  |  0.4285  0.1379 

FGLS Estimates:

         |  Coeff     SE   
---------------------------
 Const   | -0.0102  0.0017 
 CPIRate |  0.8853  0.0169 
 WRRate  |  0.8897  0.0294 
 MSRate  |  0.4874  0.0291 

Создайте эту модель регрессии с ARMA (1,2) ошибки, где εt является Гауссовым со средним значением 0 и отклонением 1.

yt=1+xt[23]+utut=0.6ut-1+εt-0.3εt-1+0.1εt-1.

beta = [2 3];
phi = 0.2;
theta = [-0.3 0.1];
Mdl = regARIMA('AR',phi,'MA',theta,'Intercept',1,'Beta',beta,'Variance',1);

Mdl regARIMA модель. Можно получить доступ к его свойствам с помощью записи через точку.

Симулируйте 500 периодов 2D стандартных Гауссовых значений для xt, и затем симулируйте ответы с помощью Mdl.

numObs = 500;
rng(1); % For reproducibility
X = randn(numObs,2);
y = simulate(Mdl,numObs,'X',X);

fgls AR поддержек (p) инновационные модели. Можно преобразовать полином модели ARMA в полином модели AR бесконечной задержки использование arma2ar. По умолчанию, arma2ar возвращает коэффициенты для первых 10 терминов. После преобразования определите, сколько задержек получившейся модели AR является практически значительным путем проверки длины возвращенного вектора коэффициентов. Выберите количество терминов, которые превышают 0.00001.

format long
arParams = arma2ar(phi,theta)
arParams = 1×3

   -0.1000    0.0700    0.0310

arLags = sum(abs(arParams) > 0.00001);
format short

Некоторые параметры имеют маленькую величину. Вы можете хотеть сократить количество задержек, чтобы включать в инновационную модель для fgls.

Оцените коэффициенты и их стандартные погрешности с помощью FGLS и симулированных данных. Укажите, что инновации включают AR (arLags) процесс.

[coeff,~,EstCoeffCov] = fgls(X,y,'innovMdl','AR','arLags',arLags)
coeff = 3×1

    1.0372
    2.0366
    2.9918

EstCoeffCov = 3×3

    0.0026   -0.0000    0.0001
   -0.0000    0.0022    0.0000
    0.0001    0.0000    0.0024

Предполагаемые коэффициенты близко к их истинным значениям.

Этот пример подробно останавливается на анализе в Оценке Коэффициенты FGLS Моделей, Содержащих Ошибки ARMA. Создайте эту модель регрессии с ARMA (1,2) ошибки, где εt является Гауссовым со средним значением 0 и отклонением 1.

yt=1+xt[23]+utut=0.6ut-1+εt-0.3εt-1+0.1εt-1.

beta = [2 3];
phi = 0.2;
theta = [-0.3 0.1];
Mdl = regARIMA('AR',phi,'MA',theta,'Intercept',1,'Beta',beta,'Variance',1);

Симулируйте 500 периодов 2D стандартных Гауссовых значений для xt, и затем симулируйте ответы с помощью Mdl.

numObs = 500;
rng(1); % For reproducibility
X = randn(numObs,2);
y = simulate(Mdl,numObs,'X',X);

Преобразуйте полином модели ARMA в полином модели AR бесконечной задержки использование arma2ar. По умолчанию, arma2ar возвращает коэффициенты для первых 10 терминов. Найдите количество терминов, которые превышают 0.00001.

arParams = arma2ar(phi,theta);
arLags = sum(abs(arParams) > 0.00001);

Оцените коэффициенты регрессии с помощью трех итераций FGLS и задайте количество задержек в инновационной модели AR (arLags). Кроме того, задайте, чтобы построить содействующие оценки и их стандартные погрешности для каждой итерации, и отобразить итоговые оценки и оценки OLS в табличной форме.

[coeff,~,EstCoeffCov] = fgls(X,y,'innovMdl','AR','arLags',arLags,...
    'numIter',3,'plot',{'coeff','se'},'display','final');
OLS Estimates:

       |  Coeff    SE   
------------------------
 Const | 1.0375  0.0480 
 x1    | 2.0409  0.0473 
 x2    | 2.9860  0.0488 

FGLS Estimates:

       |  Coeff    SE   
------------------------
 Const | 1.0372  0.0514 
 x1    | 2.0366  0.0470 
 x2    | 2.9919  0.0486 

Figure contains an axes object. The axes object with title blank C o e f f i c i e n t s contains 9 objects of type line. These objects represent Const, x1, x2.

Figure contains an axes object. The axes object with title blank S t a n d a r d blank E r r o r s contains 9 objects of type line. These objects represent Const, x1, x2.

Алгоритм, кажется, сходится после того, как первая итерация и оценки близко к оценкам OLS со стандартными погрешностями, являющимися немного меньшим.

Свойства итеративных оценок FGLS в конечных выборках затрудняют, чтобы установить. Для асимптотических свойств одна итерация FGLS достаточна. fgls поддерживает итеративный FGLS для экспериментирования.

Если оценки или стандартные погрешности показывают нестабильность после последовательных итераций, то предполагаемая инновационная ковариация может быть плохо обусловлена. Рассмотрите масштабирование остаточных значений с помощью 'resCond' аргумент пары "имя-значение", чтобы улучшить создание условий предполагаемой инновационной ковариации.

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

свернуть все

Данные о предикторе для модели многофакторной линейной регрессии в виде numObs- numPreds числовая матрица.

numObs количество наблюдений и numPreds количество переменных предикторов.

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

Данные об ответе для модели многофакторной линейной регрессии в виде numObs- 1 вектор с числовыми или логическими записями.

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

Предиктор и данные об ответе для модели многофакторной линейной регрессии в виде numObs- numPreds + 1 табличный массив.

Первый numPreds переменные Tbl данные о предикторе, и последняя переменная является данными об ответе.

Данные о предикторе должны быть числовыми, и данные об ответе должны быть числовыми или логическими.

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

Оси, на которых можно построить в виде вектора из Axes объекты с длиной равняются количеству графиков, заданных plot аргумент пары "имя-значение".

По умолчанию, fgls создает отдельную фигуру для каждого графика.

Примечание

NaNs в XY, или Tbl укажите на отсутствующие значения и fgls удаляет наблюдения, содержащие по крайней мере один NaN. Таким образом, чтобы удалить NaNs в X или y, программное обеспечение объединяет их ([X y]), и затем использует мудрое списком удаление, чтобы удалить любую строку, которая содержит по крайней мере один NaN. Программное обеспечение также удаляет любую строку Tbl содержа по крайней мере один NaN. Удаление NaNs в данных уменьшает объем выборки и может также создать неправильные временные ряды.

Аргументы name-value

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

Пример: 'innovMdl','HC0','numIter',10,'plot','coeff' задает устойчивую инновационную модель ковариации Белого, 10 итераций FGLS, и построить содействующие оценки после каждой итерации.

Имена переменных используются в отображениях и графиках результатов в виде разделенной запятой пары, состоящей из 'varNames' и длина numCoeffs вектор ячейки из векторов символов. Программное обеспечение обрезает все имена переменных до первых пяти символов.

varNames должен включать имена переменных для всех переменных в модели, таких как термин точки пересечения (например, 'Const') или термины высшего порядка (например, 'x1^2' или 'x1:x2'). Если значение 'intercept' true, затем первым элементом является имя точки пересечения. Порядок всех других элементов соответствует порядку столбцов X или переменные предикторы в Tbl.

Если 'Intercept' true, затем его именем по умолчанию является 'Const'. Имена переменных по умолчанию для:

  • Переменные предикторы в X вектор ячейки из векторов символов {'x1','x2',...}

  • Табличный массив Tbl свойство Tbl.Properties.VariableNames

Пример: 'varNames',{'Const','AGE','BBD'}

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

Укажите, включать ли точку пересечения модели когда fgls подбирает модель в виде разделенной запятой пары, состоящей из 'intercept' и true или false. Количество коэффициентов модели, numCoeffs, numPreds + intercept.

ЗначениеОписание
trueВключайте точку пересечения в модель.
falseИсключите точку пересечения из модели.

Пример: 'intercept',false

Типы данных: логический

Модель для инновационной ковариации оценивает в виде разделенной запятой пары, состоящей из 'innovMdl' и вектор символов.

Установите 'innovMdl' задавать структуру инновационного средства оценки ковариации Ω^.

  • Для диагональных инновационных моделей ковариации (i.e., модели с heteroscedasticity), Ω^=diag(ω), где ω = {ωi; i = 1..., T} является вектором из инновационных оценок отклонения для наблюдений и T = numObs.

    fgls оценивает управляемый данными векторный ω с помощью соответствующих остаточных значений модели (ε), их рычаги hi=xi(XX)1xi, и степени свободы dfe.

    Используйте эту таблицу, чтобы выбрать 'innovMdl'.

    ЗначениеВесСсылка
    'CLM'

    ωi=1dfei=1Tεi2

    [4]
    'HC0'

    ωi=εi2

    [6]
    'HC1'

    ωi=Tdfeεi2

    [5]
    'HC2'

    ωi=εi21hi

    [5]
    'HC3'

    ωi=εi2(1hi)2

    [5]
    'HC4'

    ωi=εi2(1hi)di

    где di=min(4,hih¯),

    [1]

  • Для полных инновационных моделей ковариации (т.е. моделей, имеющих heteroscedasticity и автокорреляции), задайте 'AR'. Программное обеспечение налагает модель AR (p) на инновации и построения Ω^ использование количества задержек, p, задано аргументом пары "имя-значение" arLags и уравнения Уокера Рождества.

Если numIter 1 и вы задаете InnovCov0, затем fgls игнорирует InnovMdl.

Пример: 'innovMdl',HC0

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

Количество задержек, чтобы включать в инновационную модель AR в виде разделенной запятой пары, состоящей из 'arLags' и положительное целое число.

Если innovMdl не 'AR' (т.е. для диагональных моделей), затем программное обеспечение игнорирует значение 'arLags'.

Для общих инновационных моделей ARMA преобразуйте в эквивалентную форму AR:

  • Построение инновационной модели ARMA изолирует полином оператора использование LagOp. Затем разделите полином AR на использование полинома MA, например, mrdivide. Результатом является бесконечный порядок, представление AR модели ARMA.

  • Используя arma2ar, который возвращает коэффициенты бесконечного порядка, представление AR модели ARMA.

Пример: 'arLags',4

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

Начальная инновационная ковариация в виде заданной запятой пары, состоящей из 'InnovCov0' и вектор из положительных скалярных величин, положительной полуопределенной матрицы или положительной определенной матрицы.

InnovCov0 заменяет управляемую данными оценку инновационной ковариации (Ω^) в первой итерации GLS.

  • Для диагональных инновационных моделей ковариации (т.е. моделей с heteroscedasticity), задайте numObs- 1 вектор. InnovCov0 (j) отклонение инноваций j.

  • Для полных инновационных моделей ковариации (т.е. моделей, имеющих heteroscedasticity и автокорреляции), задайте numObs- numObs матрица. InnovCov0 (jK) ковариация инноваций j и k.

  • По умолчанию, fgls использует управляемое данными Ω^ (см. innovMdl).

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

Количество итераций, чтобы реализовать для алгоритма FGLS в виде разделенной запятой пары, состоящей из 'numIter' и положительное целое число.

fgls оценивает инновационную ковариацию (Ω^) в каждой итерации от остаточного ряда согласно инновационной модели ковариации (innovMdl). Затем программное обеспечение вычисляет оценки GLS коэффициентов модели.

Пример: 'numIter',10

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

Отметьте указание, чтобы масштабировать остаточные значения в каждой итерации FGLS в виде разделенной запятой пары, состоящей из 'resCond' и true или false.

ЗначениеОписание
truefgls масштабирует остаточные значения в каждой итерации.
falsefgls не масштабирует остаточные значения в каждой итерации.

Масштабирование остаточных значений в каждой итерации FGLS имеет тенденцию улучшать создание условий оценки инновационной ковариации (Ω^).

Типы данных: логический

Отображение Командного окна управляет в виде разделенной запятой пары, состоящей из 'display' и значение в этой таблице.

ЗначениеОписание
'final'Отобразите итоговые оценки.
'iter'Отобразите оценки после каждой итерации.
'off'Подавите отображение Командного окна.

fgls показывает результаты оценки в табличной форме.

Пример: 'display','iter'

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

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

Чтобы исследовать сходимость алгоритма FGLS, это - хорошая практика, чтобы задать графический вывод оценок для каждой итерации. Эта таблица содержит доступные значения.

ЗначениеОписание
'allПостройте предполагаемые коэффициенты, их стандартные погрешности и остаточную среднеквадратическую ошибку (MSE) на отдельных графиках.
'coeff'Постройте предполагаемые коэффициенты.
'mse'Постройте MSE.
'off'Не стройте результаты.
'se'Постройте предполагаемый коэффициент.

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

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

свернуть все

Содействующие оценки FGLS, возвращенные как numCoeffs- 1 числовой вектор.

Порядок оценок соответствует порядку столбцов матрицы предиктора или Tbl.VariableNames. Например, в модели с точкой пересечения, значением β^1 (соответствие предиктору x 1), находится в положении 2 coeff.

Содействующие оценки стандартной погрешности, возвращенные как numCoeffs- 1 числовое. Элементы se sqrt(diag(EstCoeffCov)).

Порядок оценок соответствует порядку коэффициентов в coeff. Например, в модели с точкой пересечения, предполагаемой стандартной погрешностью β^1 (соответствие предиктору x 1), находится в положении 2 se, и квадратный корень из значения в положении (2,2) EstCoeffCov.

Содействующая оценка ковариации, возвращенная как numCoeffs- numCoeffs числовая матрица.

Порядок строк и столбцов EstCoeffCov соответствует порядку коэффициентов в coeff. Например, в модели с точкой пересечения, предполагаемой ковариацией β^1 (соответствие предиктору x 1) и β^2 (соответствие предиктору x 2), находятся в положениях (2,3) и (3,2) EstCoeffCov, соответственно.

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

iterPlots не доступно если значение plot аргументом пары "имя-значение" является 'off'.

Больше о

свернуть все

Выполнимые обобщенные наименьшие квадраты

Feasible generalized least squares (FGLS) оценивает коэффициенты модели линейной регрессии кратного и их ковариационной матрицы в присутствии несферических инноваций с неизвестной ковариационной матрицей.

Позвольте yt = Xt β + εt быть кратным модель линейной регрессии, где инновационный процесс, εt является Гауссовым со средним значением 0, но с истинной, несферической ковариационной матрицей Ω (например, инновациями является heteroscedastic или автокоррелируемый). Кроме того, предположите, что объемом выборки является T и существуют предикторы p (включая точку пересечения). Затем средство оценки FGLS β

β^FGLS=(XΩ^1X)1XΩ^1y,

где Ω^ инновационная оценка ковариации на основе модели (e.g., инновационный процесс формирует модель AR (1)). Предполагаемая содействующая ковариационная матрица

Σ^FGLS=σ^FGLS2(XΩ^1X)1,

где

σ^FGLS2=y[Ω^1Ω^1X(XΩ^1X)1XΩ^1]y/(Tp).

Оценки FGLS вычисляются можно следующим образом:

  1. OLS применяется к данным, и затем остаточным значениям (ε^t) вычисляются.

  2. Ω^ оценивается на основе модели для инновационной ковариации.

  3. β^FGLS оценивается, наряду с его ковариационной матрицей Σ^FGLS.

  4. Дополнительный: Этот процесс может быть выполнен с помощью итераций путем выполнения следующих шагов до β^FGLS сходится.

    1. Вычислите остаточные значения подобранной модели с помощью оценок FGLS.

    2. Примените шаги 2-3.

Если Ω^ сопоставимое средство оценки Ω и предикторы, которые включают X, являются внешними, затем средства оценки FGLS сопоставимы и эффективны.

Асимптотические распределения средств оценки FGLS неизменны повторной итерацией. Однако итерации могут изменить конечные демонстрационные распределения.

Обобщенные наименьшие квадраты

Generalized least squares (GLS) оценивает коэффициенты модели линейной регрессии кратного и их ковариационной матрицы в присутствии несферических инноваций с известной ковариационной матрицей.

Настройка и процесс для получения оценок GLS эквивалентны в FGLS, но замене Ω^ с известной инновационной ковариационной матрицей Ω.

В присутствии несферических инноваций, и с известной инновационной ковариацией, средства оценки GLS являются несмещенными, эффективными, и сопоставимыми, и тесты гипотезы на основе оценок допустимы.

Метод взвешенных наименьших квадратов

Weighted least squares (WLS) оценивает коэффициенты модели линейной регрессии кратного и их ковариационной матрицы в присутствии некоррелированых, но heteroscedastic инноваций с известной, диагональной ковариационной матрицей.

Настройка и процесс, чтобы получить оценки WLS эквивалентны в FGLS, но замене Ω^ с известным, диагональной матрицей весов, обычно диагональными элементами является инверсия отклонений инноваций.

В присутствии heteroscedastic инноваций, и когда отклонения инноваций известны, средства оценки WLS являются несмещенными, эффективными, и сопоставимыми, и тесты гипотезы на основе оценок допустимы.

Советы

  • Получить стандартные оценки обобщенных наименьших квадратов (GLS):

    • Установите InnovCov0 аргумент пары "имя-значение" известной инновационной ковариации.

    • Установите numIter аргумент пары "имя-значение" 1.

  • Чтобы получить оценки WLS, установите InnovCov0 аргумент пары "имя-значение" вектору из обратных весов (e.g., инновационные оценки отклонения).

  • В определенных моделях и с повторными итерациями, различия в шкале в остаточных значениях могут произвести плохо обусловленную предполагаемую инновационную ковариацию и вызвать числовую нестабильность. Если вы устанавливаете 'resCond',true, затем создание условий улучшается.

Алгоритмы

  • В присутствии несферических инноваций GLS производит эффективные оценки относительно OLS, и сопоставимые содействующие ковариации, условное выражение на инновационной ковариации. Степень, к который fgls обеспечивает эти свойства, зависит от точности и модели и оценки инновационной ковариации.

  • Вместо того, чтобы оценить FGLS оценивает обычный путь, fgls методы использования, которые быстрее и более устойчивы, и применимы к случаям неполного ранга.

  • Традиционные методы FGLS, такие как процедура Кокрейна-Оркатта, используют младший разряд, авторегрессивные модели. Эти методы, однако, оценивают параметры в инновационной ковариационной матрице с помощью OLS, где fgls оценка наибольшего правдоподобия (MLE) использования [2].

Ссылки

[1] Cribari-Neto, F. "Асимптотический Вывод Под Heteroskedasticity Неизвестной Формы". Computational Statistics & Data Analysis. Издание 45, 2004, стр 215–233.

[2] Гамильтон, J. D. Анализ Временных Рядов. Принстон, NJ: Издательство Принстонского университета, 1994.

[3] Судья, Г. Г., В. Э. Гриффитс, Р. К. Хилл, Х. Латкеполь и Т. К. Ли. Теория и практика эконометрики. Нью-Йорк, Нью-Йорк: John Wiley & Sons, Inc., 1985.

[4] Kutner, M. H. К. Дж. Нахцхайм, Дж. Нетер и В. Ли. Прикладные Линейные Статистические модели. 5-й редактор Нью-Йорк: Макгроу-Хилл/ирвин, 2005.

[5] Маккиннон, J. G. и H. Белый. "Некоторые Heteroskedasticity-сопоставимые Средства оценки Ковариационной матрицы с Улучшенными Конечными Демонстрационными Свойствами". Журнал Эконометрики. Издание 29, 1985, стр 305–325.

[6] Белый, H. "Heteroskedasticity-сопоставимая Ковариационная матрица и Прямой Тест для Heteroskedasticity". Econometrica. Издание 48, 1980, стр 817–838.

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