fgls

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

Синтаксис

coeff = fgls(X,y)
coeff = fgls(Tbl)
coeff = fgls(___,Name,Value)
[coeff,se,EstCoeffCov] = fgls(___)
coeff = fgls(ax,___)
[coeff,se,EstCoeffCov,iterPlots] = fgls(___)

Описание

пример

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

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

пример

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

fgls удаляет все отсутствующие значения в Tbl, обозначенном NaN s, с помощью мудрого списком удаления. Другими словами, 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;

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

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

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

Mdl = fitlm(COE,CPI);

figure;
plotResiduals(Mdl,'lagged')

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

Оцените коэффициенты регрессии с помощью 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);

Коррелограммы предлагают, чтобы инновации имели значительные эффекты 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;

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

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

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

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

Mdl = fitlm(dLogTbl);

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

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

Невязки, кажется, вспыхивают в, и таким образом, они показывают 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.100000000000000   0.070000000000000   0.031000000000000

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 

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

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

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

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

свернуть все

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

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

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

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

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

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

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

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

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

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

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

Примечание

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

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

Укажите необязательные аргументы в виде пар ""имя, значение"", разделенных запятыми. Имя (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' задавать структуру инновационного средства оценки ковариации Ω^.

  • Для диагональных инновационных моделей ковариации (т.е. моделей с 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-by-1 вектор. InnovCov0(j) является отклонением инноваций j.

  • Для полных инновационных моделей ковариации (т.е. моделей, имеющих heteroscedasticity и автокорреляции), задайте numObs-by-numObs матрица. InnovCov0(j,k) является ковариацией инноваций 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-by-1 числовой вектор.

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

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

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

Содействующая оценка ковариации, возвращенная как numCoeffs-by-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,

где Ω^ инновационная оценка ковариации на основе модели (например, инновационный процесс формирует модель 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 вектору обратных весов (например, инновационные оценки отклонения).

  • В определенных моделях и с повторными итерациями, различия в шкале в невязках могут произвести плохо обусловленную предполагаемую инновационную ковариацию и вызвать числовую нестабильность. Если вы устанавливаете '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. Анализ timeseries. Принстон, 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