Effacer les filtres
Effacer les filtres

sticky information model, fmincon function

3 vues (au cours des 30 derniers jours)
Houda
Houda le 19 Jan 2018
Commenté : Houda le 24 Jan 2018
Hey, I am using Reis (2009) model, and I just replicate his results using his data. Starting with the SIBE folder, to first get the mode of the posterior distribution with the SIGE_mode: I get the following errors (matlab 2015a):
Warning: R-square and the F statistic are not well-defined unless X has a column of ones.
Type "help regress" for more information.
> In regress (line 164)
In SIGE_mode (line 33)
Subscripted assignment dimension mismatch.
Error in SIGEv (line 32)
Omega(:,1:5)=PHI;
Error in SIGEloglik (line 44)
V=SIGEv(ma_m,ma_a,ma_g,ma_nu,ma_gam,Sigma);
Error in SIGEpost (line 25)
L = SIGEloglik(pars)'; % Minus Log-likelihood
Error in fmincon (line 534)
initVals.f = feval(funfcn{3},X,varargin{:});
Error in SIGE_mode (line 67)
[Mode fval a b c d hess] = ...
Caused by:
Failure in initial user-supplied objective function evaluation. FMINCON cannot continue.
How should I fix this?? thanks for your help Houda
  4 commentaires
Afficher 2 commentaires plus anciens
Walter Roberson
Walter Roberson le 20 Jan 2018
I had to go to the dataverse site to find some of the files, but even that site did not have SIGE_mode . SIGE_mode is mentioned in one of the .doc files but it does not appear to be available.
Houda
Houda le 20 Jan 2018
Modifié(e) : Houda le 20 Jan 2018
the code SIGE_mode is included in the zipped file attached above, I attach it again (file name 2008-0077_data), the folder that contains this code is SIBE folder. I am looking also to the codes in dataverse, I didn't know about it before you told me. Thanks

Connectez-vous pour commenter.

Connectez-vous pour répondre à cette question.

Réponse acceptée

Walter Roberson
Walter Roberson le 20 Jan 2018
I see it now; I am not sure why I could not find it earlier.
The code in SIGEv has a bug unless N (the number of periods, configured in SIGE_mode) is less than or equal to T (the number of rows in the input data, read in SIGE_mode). You will need to reduce N (currently 100) to something <= 75
  5 commentaires
Afficher 3 commentaires plus anciens
Walter Roberson
Walter Roberson le 20 Jan 2018
The code in SIGEv.m fails if N > T. I traced the source of N and T back to SIGE_mode where N is assigned in the line you indicate and T is assigned a few lines earlier as the number of rows of data read in.
Houda
Houda le 24 Jan 2018
It works for N=74 Thanks dear Walter Houda

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by