空间面板计量

空间面板计量的matlab代码；空间经济计量学主要研究存在空间效应的问题。空间效应主要包括空间相关和空间差异性。在研究中涉及空间相邻、空间相邻矩阵等概念。

function results =  acf（xkp）
% PURPOSE: Estimate the coefficients of the autocorrelation 
%          （covariance with its own lagged value）
%          function up to a particular lag {Xk = Cov[YtYt-k]}
%--------------------------------------------------------------
% USAGE: results = acf（xk）
% where:      x = a time-series vector
%             k = maximun lag considered 
%             p = Plot autocorrelations
%                 p =  1 Plot function 
%--------------------------------------------------------------
% RETURNS: a results structure
%         results.meth  = ‘acf‘
%         results.ac    = autocorrelation coefficients   
%         results.k     = dimension of the lag considered
%         results.lowb  = lower bound for 95% confidence interval
%         results.Topb  = Top bound for 95% confidence interval
%         results.qstat = Ljung-Box Q Statistic
%         results.prob  = area under de chi^（2） distribution
%                         used in testing the existance of 
%                         autocorrelation for the kth lag. 
%
%--------------------------------------------------------------
% SEE ALSO: 
%--------------------------------------------------------------
% References: 
% * Judge G. Hill C. Griffiths W. L黷kepohl H. Tsoung-Chao L
%   Introduction to the Theory and Practice of Econometrics 
%   1988 pg. 681-690.
% * Hill C. Computer Handbook to Accompany Introduction to the Theory 
%   and Practice of Econometrics1988 pg. 158-159.

% error checking on inputs
results.meth= ‘acf‘;

if （nargin > 3）
   error（‘Wrong # of arguments to acf‘）;
elseif（nargin == 3）
   p = 1;
elseif（nargin == 2） 
   p = 0;
else
   k = round（rows（x）/4）;
end;
 
if  （cols（x） > 1）
 error（‘acf cannot handle a matrix -- only vectors‘）;
end;

nobs = rows（x）;
results.k = k; 
%center data
xc= x-mean（x）;
%variance
v0= （xc‘*xc）/nobs; %not v0 = cov（x） see reference 

for i=1:k
   yt= xc（1:nobs-i1）;
   ytk= xc（1+i:nobs1）;
   ck = （yt‘*ytk）/nobs;
   xk（i1）=ck/v0;
   %ckbar = （yt‘*ytk）/（nobs-k）; 
   %xkbar（k1）=ck/v0; %asymptotically equivalent to xk
end

results.ac=xk;

%Significance of xk （95% confidence intervals）
   results.lowb= -2/sqrt（nobs）*ones（k1）;
   results.topb= 2/sqrt（nobs）*ones（k1）;
   
%Diagnostics （Ljung & Box 1978） Q-Statistic
qstat=zeros（k+11）;
for j=2:k+1
   qstat（j1） = （xk（j-11）^2）/（nobs-（j-1））+ qstat（j-11）;
   prob（j1）=  1-chis_prb（abs（（nobs*（nobs+2））*qstat（j1））j）;
end

qstat=trimr（（nobs*（nobs+2））*qstat10）;
prob =trimr（prob10）;
results.qstat=qstat;
results.prob=prob;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2631  2002-05-07 23:58  jplv7\coint\acf.m

     文件       2736  2005-03-18 10:31  jplv7\coint\adf.m

     文件        421  2009-01-24 10:48  jplv7\coint\adf_d.m

     文件       2861  2000-10-09 00:00  jplv7\coint\cadf.m

     文件        661  2001-04-20 10:04  jplv7\coint\cadf_d.m

     文件       1935  2003-11-11 12:00  jplv7\coint\contents.html

     文件       1497  2003-11-11 11:59  jplv7\coint\contents.m

     文件      16526  2003-10-23 22:02  jplv7\coint\crthegy.m

     文件       3042  2000-10-09 00:00  jplv7\coint\c_sja.m

     文件       3523  2000-10-09 00:00  jplv7\coint\c_sjt.m

     文件       1306  2000-10-09 00:00  jplv7\coint\detrend.m

     文件       9004  2003-10-27 13:58  jplv7\coint\FinalConsumption.wk1

     文件       9970  2003-10-27 13:48  jplv7\coint\hegy.m

     文件      10461  2003-11-11 06:41  jplv7\coint\hegy.zip

     文件       1754  2003-10-27 14:13  jplv7\coint\hegy_d.m

     文件       6381  2001-02-18 09:54  jplv7\coint\johansen.m

     文件       6382  2001-09-04 09:43  jplv7\coint\johansen.mbak

     文件        851  2001-04-20 10:05  jplv7\coint\johansen_d.m

     文件       4164  2000-10-09 00:00  jplv7\coint\make_html.m

     文件       4033  2000-10-09 00:00  jplv7\coint\phillips.m

     文件       1009  2001-04-20 10:05  jplv7\coint\phillips_d.m

     文件       6124  2000-10-09 00:00  jplv7\coint\prt_coint.m

     文件       1060  2000-10-09 00:00  jplv7\coint\ptrend.m

     文件      13194  2000-10-09 00:00  jplv7\coint\rztcrit.m

     文件       1003  2009-01-24 10:33  jplv7\coint\shenqiuadf.m

     文件      22317  2000-10-09 00:00  jplv7\coint\test.dat

     文件       6068  2000-10-09 00:00  jplv7\coint\ztcrit.m

     文件    9128430  2003-12-14 01:23  jplv7\data\agriculture.dat

     文件       2668  2003-12-14 01:59  jplv7\data\agriculture.txt

     文件       1669  2000-10-09 21:06  jplv7\data\anselin.dat

............此处省略1175个文件信息