indtest.m

计量工具箱

function [g2, bic] = indtest(d,n)
% PURPOSE: function called by raftery.m
% ------------------------------------------------
% SEE ALSO: coda(), prt()
% ------------------------------------------------
% REFERENCES: Geweke (1992), `Evaluating the accuracy of sampling-based
% approaches to the calculation of posterior moments', in J.O. Berger,
% J.M. Bernardo, A.P. Dawid, and A.F.M. Smith (eds.) Proceedings of
% the Fourth Valencia International Meeting on Bayesian Statistics,
% pp. 169-194, Oxford University Press
% Also: `Using simulation methods for Bayesian econometric models: 
% Inference, development and communication', at: www.econ.umn.edu/~bacc
% -----------------------------------------------------------------

% written by:
% James P. LeSage, Dept of Economics
% University of Toledo
% 2801 W. Bancroft St,
% Toledo, OH 43606
% [e-mail]
 
% NOTE: this code draws heavily on MATLAB programs written by
% Siddartha Chib available at: www.econ.umn.edu/~bacc
% I have repackaged it to make it easier to use.

 t = zeros(2,2);
 for i1 = 2:n;
  t(d(i1-1)+1,d(i1)+1)=t(d(i1-1)+1,d(i1)+1)+1;
 end;
 dcm1 = n - 1.0;  g2 = 0.0;
 for i1 = 1:2;
  for i2 = 1:2;
   if(t(i1,i2) ~= 0);
      t1 = t(i1,1) + t(i1,2); t2 = t(1,i2) + t(2,i2);
      fitted = (t1*t2)/dcm1; focus = t(i1,i2);
      g2 = g2 + log(focus/fitted)*focus;
      end;
     end; 
    end; 
    g2 = g2*2.0; bic = g2 - log(dcm1);