ecmnstd

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

Описание

пример

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

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

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

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

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

пример

[StdMean,StdCovar] = ecmnstd(___,Method) добавляет необязательный аргумент для Method.

Примеры

свернуть все

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

load ecmtechdemo.mat

Период времени для этих данных простирается с 19 апреля 2000 года до 18 апреля 2005 года. Шестой акции в Assets - Google (GOOG), которая начала торговаться 19 августа 2004 года. Так, все возвраты до 20 августа 2004 года отсутствуют и представлены как NaNs. Также у Amazon (AMZN) было несколько дней с отсутствующими значениями, рассеянными в течение последних пяти лет.

[ECMMean, ECMCovar] = ecmnmle(Data)
ECMMean = 12×1

    0.0008
    0.0008
   -0.0005
    0.0002
    0.0011
    0.0038
   -0.0003
   -0.0000
   -0.0003
   -0.0000
      ⋮

ECMCovar = 12×12

    0.0012    0.0005    0.0006    0.0005    0.0005    0.0003    0.0005    0.0003    0.0006    0.0003    0.0005    0.0006
    0.0005    0.0024    0.0007    0.0006    0.0010    0.0004    0.0005    0.0003    0.0006    0.0004    0.0006    0.0012
    0.0006    0.0007    0.0013    0.0007    0.0007    0.0003    0.0006    0.0004    0.0008    0.0005    0.0008    0.0008
    0.0005    0.0006    0.0007    0.0009    0.0006    0.0002    0.0005    0.0003    0.0007    0.0004    0.0005    0.0007
    0.0005    0.0010    0.0007    0.0006    0.0016    0.0006    0.0005    0.0003    0.0006    0.0004    0.0007    0.0011
    0.0003    0.0004    0.0003    0.0002    0.0006    0.0022    0.0001    0.0002    0.0002    0.0001    0.0003    0.0016
    0.0005    0.0005    0.0006    0.0005    0.0005    0.0001    0.0009    0.0003    0.0005    0.0004    0.0005    0.0006
    0.0003    0.0003    0.0004    0.0003    0.0003    0.0002    0.0003    0.0005    0.0004    0.0003    0.0004    0.0004
    0.0006    0.0006    0.0008    0.0007    0.0006    0.0002    0.0005    0.0004    0.0011    0.0005    0.0007    0.0007
    0.0003    0.0004    0.0005    0.0004    0.0004    0.0001    0.0004    0.0003    0.0005    0.0006    0.0004    0.0005
      ⋮

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

StdMeanF = ecmnstd(Data,ECMMean,ECMCovar,'fisher')
StdMeanF = 12×1

    0.0010
    0.0014
    0.0010
    0.0009
    0.0011
    0.0013
    0.0009
    0.0006
    0.0009
    0.0007
      ⋮

Вычислите стандартные ошибки, которые используют сгенерированную данными матрицу Гессиана (которая учитывает возможную потерю информации из-за недостающих данных) с помощью опции 'hessian'.

StdMeanH = ecmnstd(Data,ECMMean,ECMCovar,'hessian')
StdMeanH = 12×1

    0.0010
    0.0014
    0.0010
    0.0009
    0.0011
    0.0021
    0.0009
    0.0006
    0.0009
    0.0007
      ⋮

В различие стандартных ошибок показано увеличение неопределенности оценки ожидаемых возвратов активов из-за недостающих данных. Чтобы просмотреть различия:

Assets
Assets = 1x12 cell
  Columns 1 through 6

    {'AAPL'}    {'AMZN'}    {'CSCO'}    {'DELL'}    {'EBAY'}    {'GOOG'}

  Columns 7 through 12

    {'HPQ'}    {'IBM'}    {'INTC'}    {'MSFT'}    {'ORCL'}    {'YHOO'}

StdMeanH'
ans = 1×12

    0.0010    0.0014    0.0010    0.0009    0.0011    0.0021    0.0009    0.0006    0.0009    0.0007    0.0010    0.0012

StdMeanF'
ans = 1×12

    0.0010    0.0014    0.0010    0.0009    0.0011    0.0013    0.0009    0.0006    0.0009    0.0007    0.0010    0.0012

StdMeanH' - StdMeanF'
ans = 1×12
10-3 ×

   -0.0000    0.0021   -0.0000   -0.0000   -0.0000    0.7742   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000

Два актива с отсутствующими данными, AMZN и GOOG, являются единственными активами, которые имеют различия из-за недостающей информации.

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

свернуть все

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

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

Максимальные оценки параметра правдоподобия для среднего значения Data использование алгоритма ECM, заданного как NUMSERIES-by- 1 Вектор-столбец.

Максимальные оценки параметра правдоподобия для ковариации Data использование алгоритма ECM, заданного как NUMSERIES-by- NUMSERIES матрица.

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

  • 'hessian' - Гессиан наблюдаемой отрицательной логарифмической функции логарифмической правдоподобности. Этот метод рекомендуется, поскольку результирующие стандартные ошибки включают неопределенности увеличения из-за недостающих данных. В частности, стандартные ошибки, вычисленные с помощью Hessian, обычно больше, чем стандартные ошибки, вычисленные с помощью информационной матрицы Фишера.

  • 'fisher' - Информационная матрица Фишера.

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

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

свернуть все

Стандартные ошибки оценок для каждого элемента Mean вектор, возвращенный как NUMSERIES-by- 1 Вектор-столбец.

Стандартные ошибки оценок для каждого элемента Covariance матрица, возвращенная как NUMSERIES-by- NUMSERIES матрица.

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