用matlab程序运行空间面板过程中出现出现如下问题,要怎么解决呢?请各位大牛指点-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 会计>>

会计库

>>

用matlab程序运行空间面板过程中出现出现如下问题,要怎么解决呢?请各位大牛指点

用matlab程序运行空间面板过程中出现出现如下问题,要怎么解决呢?请各位大牛指点

发布:zwt991232076 | 分类:会计库

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

经管之家新媒体交易平台

提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

提供微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

程序:%dimensionsoftheproblemA=xlsread('C:\Users\lenovo\Documents\MATLAB\Zandberg\cigarette.xls');W1=xlsread('C:\Users\lenovo\Documents\MATLAB\Zandberg\Spat-Sym-US.xls');T=29;%numberoftimeperiodsN=46; ...
免费学术公开课,扫码加入


程序:
% dimensions of the problem
A=xlsread('C:\Users\lenovo\Documents\MATLAB\Zandberg\cigarette.xls');
W1=xlsread('C:\Users\lenovo\Documents\MATLAB\Zandberg\Spat-Sym-US.xls');
T=29; % number of time periods
N=46; % number of regions
% row-normalize W
W=normw(W1); % function of LeSage
y=A(:,3); % column number in the data matrix that corresponds to the dependent variable
x=A(:,[4,6]); % column numbers in the data matrix that correspond to the independent variables
for t=1:T+1
t1=(t-1)*N+1;t2=t*N;
wx(t1:t2,:)=W*x(t1:t2,:);
Wy(t1:t2,1)=W*y(t1:t2,1);
end
[nobs K]=size(x);
xconstant=ones(nobs,1);
info.lflag=0;
info.tl=1;
info.stl=1;
info.ted=1;
info.dyn=1;
info.model=3;
info.fe=0;
info.bc=0;
results=sar_panel_FE(y(N+1:end),[y(1:end-N) Wy(1:end-N) x(N+1:end,:) wx(N+1:end,:)],W,T,info);
vnames=strvcat('logcit','timelag logcit','spacetimelag logcit','logp','logy','W*logp','W*logy');
results1=sar_jihai_time(y(1:nobs),[x(N+1:nobs,:) wx(N+1:nobs,:)],W,info);
results.beta=results1.theta1(1:end-2);
results.rho=results1.theta1(end-1);
results.tstat=results1.tstat(1:end-1);
prt_sp(results,vnames,1);
btemp=results1.theta1;
varcov=results1.varcov;
%
% Direct and indirect effects estimates
% st=short term, lt=long term
% c=(long term) convergence effect of dependent variable
%
% Direct and indirect effects estimates
% st=short term, lt=long term
% c=(long term) convergence effect of dependent variable
%
NSIM=1000;
px=2;
[npar dummy]=size(btemp);
simresults=zeros(npar-1,NSIM);
simdirst=zeros(px,NSIM);
simindst=zeros(px,NSIM);
simtotst=zeros(px,NSIM);
simdirlt=zeros(px,NSIM);
simindlt=zeros(px,NSIM);
simtotlt=zeros(px,NSIM);
simdirc=zeros(1,NSIM);
simindc=zeros(1,NSIM);
simtotc=zeros(1,NSIM);
for sim=1:NSIM
parms = chol(varcov)'*randn(size(btemp)) + btemp;
deltasim = parms(npar-1,1); % coef WY(t)
betasim = parms(3:npar-2,1);
tausim = parms(1,1); % Coef Y(t-1)
etasim = parms(2,1); % Coef WY(t-1)
simresults(:,sim)=[tausim;etasim;betasim;deltasim];
SS=(eye(N)-deltasim*W)\eye(N);
SC=SS*((tausim-1)*eye(N)+(deltasim+etasim)*W);
simdirc(1,sim)=sum(diag(SC))/N; % average direct effect
simindc(1,sim)=sum(sum(SC,2)-diag(SC))/N; % average indirect effect
simtotc(1,sim)=simdirc(1,sim)+simindc(1,sim);
for p=1:px
C=zeros(N,N);
for i=1:N
for j=1:N
if (i==j) C(i,j)=betasim(p);
else C(i,j)=betasim(p+2)*W(i,j);
end
end
end
SC=SS*C;
simdirst(p,sim)=sum(diag(SC))/N; % average direct effect
simindst(p,sim)=sum(sum(SC,2)-diag(SC))/N; % average indirect effect
simtotst(p,sim)=simdirst(p,sim)+simindst(p,sim);
SC=((1-tausim)*eye(N)-(deltasim+etasim)*W)\C;
simdirlt(p,sim)=sum(diag(SC))/N; % average direct effect
simindlt(p,sim)=sum(sum(SC,2)-diag(SC))/N; % average indirect effect
simtotlt(p,sim)=simdirlt(p,sim)+simindlt(p,sim);
end
end
fprintf(1,'Convergence effect \n');
[mean(simdirc,2) mean(simdirc,2)./std(simdirc,0,2) mean(simindc,2) mean(simindc,2)./std(simindc,0,2)...
mean(simtotc,2) mean(simtotc,2)./std(simtotc,0,2)]
fprintf(1,'Short term effects \n');
[mean(simdirst,2) mean(simdirst,2)./std(simdirst,0,2) mean(simindst,2) mean(simindst,2)./std(simindst,0,2)...
mean(simtotst,2) mean(simtotst,2)./std(simtotst,0,2)]
fprintf(1,'Long term effects \n');
[mean(simdirlt,2) mean(simdirlt,2)./std(simdirlt,0,2) mean(simindlt,2) mean(simindlt,2)./std(simindlt,0,2)...
mean(simtotlt,2) mean(simtotlt,2)./std(simtotlt,0,2)]
clear simresults simdirc simindc simtotc simdirst simindst simtotst simdirlt simindlt simtotlt;
运行结果:
Pooled model with spatially lagged dependent variable, spatial and time period fixed effects
Dependent Variable = logcit
R-squared = 0.9772
corr-squared = 0.8615
sigma^2 = 0.0011
Nobs,Nvar,#FE = 1334, 7, 80
log-likelihood = 2623.2622
# of iterations = 1
min and max rho = -1.3924, 1.0000
total time in secs = 0.4510
time for optimiz = 0.1110
time for lndet = 0.0130
time for eigs = 0.1290
time for t-stats = 0.0080
No lndet approximation used
***************************************************************
Variable CoefficientAsymptot t-stat z-probability
timelag logcit 0.864568 65.041225 0.000000
spacetimelag logcit -0.015347 -0.287520 0.773714
logp -0.265979 -13.192203 0.000000
logy 0.099796 4.165387 0.000031
W*logp 0.170413 3.664979 0.000247
W*logy -0.022369 -0.868904 0.384900
W*dep.var. 0.076402 2.025253 0.042841
??? Reference to non-existent field 'varcov'.
Error in ==> demo at 33
varcov=results1.varcov;
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-4840848-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
经管之家 人大经济论坛 大学 专业 手机版