CIR MODEL оценка логарифмического правдоподобия исходных параметров

Question

для этого

% Optimization using fminsearch 
if ~isfield(Model, 'MatlabDisp'), Model.MatlabDisp = 'off'; end; 

options = optimset('LargeScale', 'off', 'MaxIter', 300, 'MaxFunEvals',300,'Display',Model.MatlabDisp, 'TolFun', 1e-4, 'TolX', 1e-4, 'TolCon', 1e-4); 
if ~isfield(Model, 'Method'), Model.Method = 'besseli'; end; 
if strcmp(Model.Method, 'ncx2pdf') 
    [Params, Fval, Exitflag] = fminsearch(@(Params) CIRobjective2(Params, Model), InitialParams, options); 
else 
    [Params, Fval, Exitflag] = fminsearch(@(Params) CIRobjective1(Params, Model), InitialParams, options); 
end 
%tol = 1e-6; 
%lb = [ tol, tol, tol]; 
%ub = [ inf, inf, inf];  
%[Params, Fval, Exitflag] = fmincon(@(Params) CIRobjective1(Params, Model),                 
InitialParams, [], [], [], [], lb, ub, [], options); 

Results.Params = Params; 
Results.Fval = -Fval/Nobs; 
Results.Exitflag = Exitflag; 

if strcmp(Model.Disp, 'y') 
    fprintf('\n alpha = %+3.6f\n mu = %+3.6f\n sigma = %+3.6f\n', Params(1), Params(2), Params(3)); 
    fprintf(' log-likelihood = %+3.6f\n', -Fval/Nobs); 
end 

Я получил

неопределенной функции или метод «mtimes» для входных аргументов типа «структура».

Error in ==> CIRobjective2 at 26 
u = c*exp(-alpha*TimeStep)*DataL; 

Error in ==> @(Params)CIRobjective2(Params,Model) 


Error in ==> fminsearch at 205 
fv(:,1) = funfcn(x,varargin{:}); 

Это реальный исходный файл,

http://www.mathworks.com/matlabcentral/fileexchange/37297-maximum-likelihood-estimation-of-the-cox-ingersoll-ross-process-the-matlab-implementation/content/CIRestimation.m

Не могли бы вы помочь мне, почему я получаю эти ошибки?

С уважением, Масуд

Answer 1

Вы уверены, что вы положили все свои функции права и определить его хорошо? Возможно, ваша функция не в правильном направлении.

Сначала попробуйте выполнить начальный параметр, определите InitialParams, чтобы он отображался в ваших переменных. Затем сделайте CIRobjective1, затем 2, затем yo может выполнить оптимизацию с использованием fminsearch. Надеюсь, это поможет. Я новый пользователь Matlab, но должен использовать его для своих тезисов.

С уважением, Лита.