X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=sidebyside;f=octave_packages%2Feconometrics-1.0.8%2Fgmm_obj.m;fp=octave_packages%2Feconometrics-1.0.8%2Fgmm_obj.m;h=0cbaa197a2cc20eae84e3fc231820439c46574d5;hb=f5f7a74bd8a4900f0b797da6783be80e11a68d86;hp=0000000000000000000000000000000000000000;hpb=1705066eceaaea976f010f669ce8e972f3734b05;p=CreaPhase.git diff --git a/octave_packages/econometrics-1.0.8/gmm_obj.m b/octave_packages/econometrics-1.0.8/gmm_obj.m new file mode 100644 index 0000000..0cbaa19 --- /dev/null +++ b/octave_packages/econometrics-1.0.8/gmm_obj.m @@ -0,0 +1,31 @@ +# Copyright (C) 2003,2004,2005 Michael Creel +# +# This program is free software; you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2 of the License, or +# (at your option) any later version. +# +# This program is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; If not, see . + + +# The GMM objective function, for internal use by gmm_estimate +# This is scaled so that it converges to a finite number. +# To get the chi-square specification +# test you need to multiply by n (the sample size) +function obj_value = gmm_obj(theta, data, weight, moments, momentargs) + + m = average_moments(theta, data, moments, momentargs); + + obj_value = m' * weight *m; + + if (((abs(obj_value) == Inf)) || (isnan(obj_value))) + obj_value = realmax; + endif + +endfunction