estimate
максимизирует функцию логарифмической правдоподобности использование fmincon
от Optimization Toolbox™. fmincon
имеет много опций оптимизации, таких как выбор алгоритма оптимизации и допуска нарушения ограничений. Выберите опции оптимизации с помощью optimoptions
.
estimate
использует fmincon
опции оптимизации по умолчанию, за этими исключениями. Для получения дополнительной информации смотрите fmincon
и optimoptions
в Optimization Toolbox.
Свойства optimoptions | Описание | оцените Настройки |
---|---|---|
Algorithm | Алгоритм для минимизации отрицательной функции логарифмической правдоподобности | 'sqp' |
Display | Level of display для прогресса оптимизации | 'off' |
Diagnostics | Отобразитесь для диагностической информации о функции, которая будет минимизирована | 'off' |
ConstraintTolerance | Допуск завершения на нарушениях ограничений | 1e-7 |
Если вы хотите использовать опции оптимизации, которые отличаются от значения по умолчанию, то установленный ваше собственное использование optimoptions
.
Например, предположите, что вы хотите estimate
отобразить диагностику оптимизации. Лучшая практика состоит в том, чтобы установить аргумент пары "имя-значение" 'Display','diagnostics'
в estimate
. В качестве альтернативы можно направить оптимизатор, чтобы отобразить диагностику оптимизации.
Задайте модель регрессии с AR (1) ошибки (Mdl
) и симулируйте данные из него.
Mdl0 = regARIMA('AR',0.5,'Intercept',0,'Variance',1); rng(1); % For reproducibility y = simulate(Mdl0,25);
Mdl
не имеет компонента регрессии. По умолчанию, fmincon
не отображает диагностику оптимизации. Используйте optimoptions
установить его отображать диагностику оптимизации и устанавливать другой fmincon
свойства к настройкам по умолчанию estimate
перечисленный в предыдущей таблице.
options = optimoptions(@fmincon,'Diagnostics','on','Algorithm',... 'sqp','Display','off','ConstraintTolerance',1e-7)
options = fmincon options: Options used by current Algorithm ('sqp'): (Other available algorithms: 'active-set', 'interior-point', 'sqp-legacy', 'trust-region-reflective') Set properties: Algorithm: 'sqp' ConstraintTolerance: 1.0000e-07 Display: 'off' Default properties: CheckGradients: 0 FiniteDifferenceStepSize: 'sqrt(eps)' FiniteDifferenceType: 'forward' MaxFunctionEvaluations: '100*numberOfVariables' MaxIterations: 400 ObjectiveLimit: -1.0000e+20 OptimalityTolerance: 1.0000e-06 OutputFcn: [] PlotFcn: [] ScaleProblem: 0 SpecifyConstraintGradient: 0 SpecifyObjectiveGradient: 0 StepTolerance: 1.0000e-06 TypicalX: 'ones(numberOfVariables,1)' UseParallel: 0 Show options not used by current Algorithm ('sqp')
% @fmincon is the function handle for fmincon
Опции, которые вы устанавливаете, появляются под Set by user:
заголовок. Свойства под Default:
заголовок является другими опциями, которые можно установить.
Подходящий Mdl
к y
использование новых опций оптимизации.
Mdl = regARIMA(1,0,0);
EstMdl = estimate(Mdl,y,'Options',options);
____________________________________________________________ Diagnostic Information Number of variables: 3 Functions Objective: @(X)nLogLike(X,YData,XData,E,U,Mdl,AR.Lags,MA.Lags,maxPQ,T,isDistributionT,userSpecifiedU0,trapValue) Gradient: finite-differencing Hessian: finite-differencing (or Quasi-Newton) Nonlinear constraints: @(x)internal.econ.arimaNonLinearConstraints(x,LagsAR,LagsSAR,LagsMA,LagsSMA,tolerance) Nonlinear constraints gradient: finite-differencing Constraints Number of nonlinear inequality constraints: 1 Number of nonlinear equality constraints: 0 Number of linear inequality constraints: 0 Number of linear equality constraints: 0 Number of lower bound constraints: 3 Number of upper bound constraints: 3 Algorithm selected sqp ____________________________________________________________ End diagnostic information ARMA(1,0) Error Model (Gaussian Distribution): Value StandardError TStatistic PValue ________ _____________ __________ _________ Intercept -0.12097 0.44747 -0.27034 0.7869 AR{1} 0.46386 0.15781 2.9393 0.0032895 Variance 1.2308 0.47275 2.6035 0.0092266
Примечание
estimate
численно максимизирует функцию логарифмической правдоподобности, потенциально с помощью равенства, неравенства и ограничений нижней и верхней границы. Если вы устанавливаете Algorithm
к чему-либо кроме sqp
, убедитесь, что алгоритм поддерживает подобные ограничения, такие как interior-point
. Например, trust-region-reflective
не поддерживает ограничения неравенства.
estimate
устанавливает ограничительный уровень ConstraintTolerance
таким образом, ограничения не нарушены. Оценка с активным ограничением имеет ненадежные стандартные погрешности, потому что оценка ковариации отклонения принимает, что функция правдоподобия локально квадратична вокруг оценки наибольшего правдоподобия.
Программное обеспечение осуществляет эти ограничения при оценке модели регрессии с ошибками ARIMA:
Устойчивость несезонных и сезонных полиномов оператора AR
Обратимость несезонных и сезонных полиномов оператора MA
Инновационное отклонение, строго больше, чем нуль
Степени свободы, строго больше, чем два для инновационного распределения t
estimate
| fmincon
| optimoptions
| regARIMA