ecmnstd

Стандартные погрешности для среднего значения и ковариации неполных данных

Синтаксис

[StdMean,StdCovariance] = ecmnstd(Data,Mean,Covariance,Method)

Аргументы

Data

NUMSAMPLES-by-NUMSERIES матрица с выборками NUMSAMPLES NUMSERIES - размерный случайный вектор. Отсутствующие значения обозначаются NaN s.

Mean

NUMSERIES-by-1 вектор-столбец параметра наибольшего правдоподобия оценивает для среднего значения Data с помощью алгоритма условной максимизации ожидания (ECM)

Covariance

NUMSERIES-by-NUMSERIES матрица ковариации наибольшего правдоподобия оценивает для ковариации Data с помощью алгоритма ECM

Method

(Необязательно) метод указания Вектора символов оценки для вычислений стандартной погрешности. Методы:

  • hessian — Гессиан (по умолчанию) наблюдаемой отрицательной логарифмической функции правдоподобия.

  • fisher — Матрица информации о Фишере.

Описание

[StdMean, StdCovariance] = ecmnstd(Data,Mean,Covariance,Method) вычисляет стандартные погрешности для среднего значения и ковариации неполных данных.

StdMean является NUMSERIES-by-1 вектор-столбец стандартных погрешностей оценок для каждого элемента среднего векторного Mean.

StdCovariance является NUMSERIES-by-NUMSERIES матрица стандартных погрешностей оценок для каждого элемента ковариационной матрицы Covariance.

Используйте эту стандартную программу после оценки среднего значения и ковариации Data с ecmnmle. Если средние и отличные элементы ковариации обработаны как параметр θ по оценке наибольшего правдоподобия полных данных, то как количество демонстрационных увеличений, θ достигает асимптотической нормальности, таким образом что

θE[θ]N(0,I1(θ)),

где E [θ] является средним значением, и I (θ) является матрицей информации о Фишере.

С недостающими данными Гессиан H (θ) является хорошим приближением для получения информации о Фишере (который может только быть аппроксимирован, когда данные отсутствуют).

Обычно желательно использовать Method по умолчанию, поскольку результирующие стандартные погрешности включают увеличенную неуверенность из-за недостающих данных. В частности, стандартные погрешности, вычисленные с Гессианом, обычно больше, чем стандартные погрешности, вычисленные с матрицей информации о Фишере.

Примечание

Эта стандартная программа является медленной для NUMSERIES > 10 или NUMSAMPLES > 1000.

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