exponenta event banner

ecmnstd

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

Описание

пример

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

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

θ−E[θ]∼N (0, I 1 (start)),

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

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

пример

[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около-NUMSERIES матрица с NUMSAMPLES образцы NUMSERIES-мерный случайный вектор. Отсутствующие значения обозначаются NaNs.

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

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

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

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

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

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

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

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

свернуть все

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

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

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