感谢楼组

运行检验程序结果显示为     未定义函数或变量 'lmlag_panel'   怎么处理

运行这个文件demoLMsarsem_panel.m 结果怎么会不行能

'lmlag_panel'是一个函数，在你的文件包中可能没有包括这个函数。文件demoLMsarsem_panel.m 并不包括所有函数，有很多函数要在外部调用，所以你先检查下你存放程序的目录里是否包含哪个函数。

我用这个代码得不出结果   怎么办

运行结果总是显示：错误使用 mprint (line 231)
Wrong # rnames in mprint
我确实认真核实了每个模块的变量代码，没发现写错，请问大神还有什么原因呢？急求

原程序是两个解释变量，自己换成四个解释变量的话，除了解释变量以外，还有什么需要改动的地方吗？

楼主，就是还需要改哪些子程序呢？

你把mprint函数找出来，仔细对照下，你的函数变量名有没有传递过来。这个错误是你修改了原函数，但是修改后的结果没有传递到打印函数中，变量数量出错导致的。

在我前面的帖子里有介绍，除了修改主程序的变量名和数量以外，你还需要修改相应的子程序中的变量名和数量。

谢谢您，我找出来mprint函数，很长，看不明白，但是里面没有找到函数名字的地方，能不能问您的qq，我把mprint函数传给您，您帮我指出来，可以吗？

谢谢您，我找出来mprint函数，很长，看不明白，但是里面没有找到函数名字的地方，能不能问您的qq，我把mprint函数传给您，您帮我指出来，可以吗？

% setup defaults
fid = 1; rflag = 0; cflag = 0; rnum = 0; nfmts = 1; cwidth = 80;
[nobs nvars] = size(y);
begr = 1; endr = nobs; begc = 1; endc = nvars; fmt = '%10.4f';
if nargin == 1
% rely on defaults
elseif nargin == 2
  if ~isstruct(info)
    error('mprint: you must supply the options as a structure variable');
  end;
fields = fieldnames(info);
nf = length(fields);
for i=1:nf
    if strcmp(fields{i},'fmt')
        fmts = info.fmt;
  [nfmts junk] = size(fmts);
  if nfmts == nvars
   fmt = fmts;
  elseif nfmts == 1
   fmt = fmts;
  else
   error('mprint: wrong # of formats in string -- need nvar');
  end;
    elseif strcmp(fields{i},'fid')
        fid = info.fid;
    elseif strcmp(fields{i},'begc');
        begc = info.begc;
    elseif strcmp(fields{i},'begr');
        begr = info.begr;
    elseif strcmp(fields{i},'endc');
        endc = info.endc;
    elseif strcmp(fields{i},'endr');
        endr = info.endr;
    elseif strcmp(fields{i},'width');
      cwidth = info.width;
    elseif strcmp(fields{i},'cnames');
        cnames = info.cnames;
        cflag = 1;
    elseif strcmp(fields{i},'rnames');
        rnames = info.rnames;
        rflag = 1;
    elseif strcmp(fields{i},'rflag');
        rnum = info.rflag;
    end;
end;

else
error('Wrong # of arguments to mprint');
   
end; % end of if-elseif input checking


% see if the user supplied row names and set rnum
% correct her mistake if she did this
if rflag == 1
rnum = 0;
end;

% parse formats
if nfmts == 1
   f1 = strtok(fmt,'%');
   f2 = strtok(f1,'.');
    if strcmp(f1,f2)
     f2 = strtok(f2,'d');
     dflag = 1;
     fflag = 0;
    else
     tmp1 = strtok(fmt,'f');
     tmp2 = strtok(tmp1,'.');
     tmp1 = tmp1(2:length(tmp1));
     tmp2 = tmp2(2:length(tmp2));
     opoint = num2str(str2num(tmp1) - str2num(tmp2));
     decimal = opoint(1,length(opoint));
     f2 = strtok(f2,'f');
     fflag = 1;
     dflag = 0;
    end;
   f2 = str2num(f2);
   nwide = floor(cwidth/f2); % 80 columns divided by format
   nvar = endc-begc+1;
   nsets = ceil(nvar/nwide);
else %  wrapping in this case is based on widest format in the list
nwidev = zeros(nfmts,1);
nsetsv = zeros(nfmts,1);
f2v = zeros(nfmts,1);
dflagv = zeros(nfmts,1);
fflagv = zeros(nfmts,1);
decimalv = zeros(nfmts,1);
   for ii=1:nfmts;
   f1 = strtok(fmt(ii,:),'%');
   f2 = strtok(f1,'.');
    if strcmp(f1,f2)
     f2 = strtok(f2,'d');
     dflagv(ii,1) = 1;
     fflagv(ii,1) = 0;     
    else
     tmp1 = strtok(fmt(ii,:),'f');
     tmp2 = strtok(tmp1,'.');
     tmp1 = tmp1(2:length(tmp1));
     tmp2 = tmp2(2:length(tmp2));
     opoint = num2str(str2num(tmp1) - str2num(tmp2));
     decimalv(ii,1) = opoint(1,length(opoint));     
     f2 = strtok(f2,'f');
     fflagv(ii,1) = 1;
     dflagv(ii,1) = 0;     
    end;
   f2v(ii,1) = str2num(f2);
   nwidev(ii,1) = floor(cwidth/f2v(ii,1)); % cwidth columns divided by format
   nvar = endc-begc+1;
   nsetsv(ii,1) = ceil(nvar/nwidev(ii,1));   
end;
nsets = min(nsetsv);
nwide = max(nwidev);
end;

% if we have row and column labels
% adjust variable labels and column heading strings
% to match the width of the printing format

if rnum == 1
dstr = 'Obs#';
end;

if cflag == 1 % we have column headings
[vsize nsize] = size(cnames); % error check cnames argument
if vsize ~= nvars; error('Wrong # cnames in mprint'); end;   
if nfmts == 1 % case of only 1 format string
  nmax = max(f2,nsize); % build format strings
                        % based on widest format              
  sfmt = ['%', num2str(nmax)];
  sfmt = [sfmt,'s '];
  ffmt = ['%', num2str(nmax)];
   if dflag == 1
   ffmt = [ffmt,'d '];
   elseif fflag == 1
   ffmt = [ffmt,'.'];
   ffmt = [ffmt,decimal];
   ffmt = [ffmt,'f '];
   end;
else % we have multiple format strings, process each
sfmtv = []; fmtv = [];
  for ii=1:nfmts % find and parse multiple formats
  nmax = max(f2v(ii,:),nsize); % build format strings
                        % based on widest format              
  sfmtv{ii} = ['%', num2str(nmax)];
  sfmtv{ii} = [sfmtv{ii},'s '];
  ffmtv{ii} = ['%', num2str(nmax)];
   if dflagv(ii,1) == 1
   ffmtv{ii} = [ffmtv{ii},'d '];
   elseif fflagv(ii,1) == 1
   ffmtv{ii} = [ffmtv{ii},'.'];
   ffmtv{ii} = [ffmtv{ii},decimalv(ii,1)];   
   ffmtv{ii} = [ffmtv{ii},'f '];
   end;
  end; % end of for ii loop
end; % end of if-else
elseif cflag == 0 % we have no column headings
if nfmts == 1 % case of only 1 format string
  nmax = f2; % augment format string with a space (the hard way)
  ffmt = ['%', num2str(nmax)];
   if dflag == 1
   ffmt = [ffmt,'d '];
   elseif fflag == 1
   ffmt = [ffmt,'.'];
   ffmt = [ffmt,decimal];
   ffmt = [ffmt,'f '];
   end;
else % we have multiple format strings, process each
sfmtv = []; fmtv = [];
  for ii=1:nfmts % find and parse multiple formats
  nmax = f2v(ii,:); % augment format strings with a space
  ffmtv{ii} = ['%', num2str(nmax)];
   if dflagv(ii,1) == 1
   ffmtv{ii} = [ffmtv{ii},'d '];
   elseif fflagv(ii,1) == 1
   ffmtv{ii} = [ffmtv{ii},'.'];
   ffmtv{ii} = [ffmtv{ii},decimalv(ii,1)];   
   ffmtv{ii} = [ffmtv{ii},'f '];
   end;
  end; % end of for ii loop
end; % end of if-else   
end; % end of if-elseif cflag == 0,1
   
if rflag == 1 % we have row labels
[vsize nsize] = size(rnames); % error check cnames argument
if vsize ~= nobs+1; error('Wrong # rnames in mprint'); end;  
rfmt = ['%', num2str(nsize)];
rfmt = [rfmt,'s '];
end; % end of if rflag == 1

if (rflag == 0 & cflag == 0)
    ffmt = fmt;
end;

% print matrix
for j=1:nsets;
if nfmts == 1 % print row header and column headers
if rnum == 1;fprintf(fid,'%5s ',dstr);     
     elseif rflag == 1   
  fprintf(fid,rfmt,rnames(1,:));
     end;  
     if cflag == 1
    for i = (j-1)*nwide+begc:j*nwide+begc-1
  if i <= endc
% find version #;
  %[version,junk] = version; vers = str2num(version);
   %if vers == 5.2
   fprintf(fid,sfmt,strjust(cnames(i,:),'right'));
   %else
   %fprintf(fid,sfmt,strjust(cnames(i,:)));
   %end;
  end;
end;
     end;
  fprintf(fid,'\n');
else % we have multiple formats
if rnum == 1;fprintf(fid,'%5s ',dstr);     
    elseif rflag == 1   
fprintf(fid,rfmt,rnames(1,:));
    end;
    if cflag == 1
   for i = (j-1)*nwide+begc:j*nwide+begc-1
  if i <= endc
% find version #;
  %[version,junk] = version; vers = str2num(version);
   %if vers == 5.2
   fprintf(fid,sfmtv{i},strjust(cnames(i,:),'right'));
   %else
   %fprintf(fid,sfmtv{i},strjust(cnames(i,:)));
   %end;
  end;
   end;
    end;
fprintf(fid,'\n');
end; % end of if-else nfmts
for k = begr:endr; % print row labels and numbers in matrix
  if rnum == 1; fprintf(fid,'%5d ',k);
        elseif rflag == 1        
  fprintf(fid,rfmt,rnames(k+1,:));
        end;
  for l = (j-1)*nwide+begc:j*nwide+begc-1
   if l <= endc
    if nfmts == 1
    fprintf(fid,ffmt,y(k,l));
    else
    fprintf(fid,ffmtv{l},y(k,l));
    end;
   end;
  end; % end of for l
  fprintf(fid,'\n');
end; % end of for k
fprintf(fid,'\n');
end; % end of for j

麻烦您帮看看，我需要修改什么地方可以吗？实在是麻烦了

T=12; % number of time periods
N=3; % number of regions
% row-normalize W
W=normw(W1); % function of LeSage
y=A(:,[1]); % column number in the data matrix that corresponds to the dependent variable
x=A(:,[2:4]); % column numbers in the data matrix that correspond to the independent variables
for t=1:T
    t1=(t-1)*N+1;t2=t*N;
    wx(t1:t2,:)=W*x(t1:t2,:);
end
xconstant=ones(N*T,1);
[nobs K]=size(x);
% ----------------------------------------------------------------------------------------
% No fixed effects + spatially lagged dependent variable
info.lflag=0; % required for exact results
info.model=0;
info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on
% New routines to calculate effects estimates
results=sar_panel_FE(y,[xconstant x],W,T,info);
vnames=strvcat('logcit','intercept','logp','logy','logk');

我把变量改成了3个就，后面的Vnames我都改了啊，求楼主解答
提示错误Subscripted assignment dimension mismatch.

你好，请问你这个问题解决了么，我也一直这样报错，不知道怎么回事，还请大神指导

请问sdm的随机效应模型是把FE改成RE么
怎么显示错误呢
??? Reference to non-existent field 'parm'.

Error in ==> direct_indirect_effects_estimates at 21
parm=results.parm;
做固定效应的sdm时都好好的

只有这段代码，不清楚整个状况。你这个错误从经验上判断是处在参数传递有误上，你好好检查下参数传递的情况吧。

你的变量设置有误，你检查下，变量的数量和程序中的数量是不是一致。

Wald_spatial_lag=(Rafg*btemp)'*inv(Rafg*varcov*Rafg')*Rafg*btemp
prob_spatial_lag=1-chis_cdf (Wald_spatial_lag, K)
Undefined function or variable 'k'.
楼主您好，这是做warld 空间滞后检验出的问题，您能帮我看看为什么吗、。谢谢了

楼主您好，这个 bias correction 做出来是这样的，能麻烦您告知一下是为什么吗？
>> Spatial and time period fixed effects + spatially lagged dependent variable + spatially
% independent variables
% No bias correction
info.bc=0;
info.lflag=0; % required for exact results
info.model=3;
info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on
% New routines to calculate effects estimates
Undefined function 'Spatial' for input arguments of type 'char'.

你对变量“K”没有事先定义，你要检查前面的定义情况。看你现在贴出的语句，k是滞后期吗？

“char”是一种字符变量，前面的语句中对'Spatial'的变量符号的定义没有传递或者是在嵌套的程序里面没有顺利的传递，你可以试试用步进方式检查你的程序语言是否正确。

请问楼主有木有看过demodynfeerrorw文件，里面的m=63  % Value that indicates the period of time between the first observation in the sample and the moment the process being modeled has started in the past 这是什么意思啊

请问您的第6点：“不要光改了主程序中的变量数量和名称，在子程序中的也要修改哦”
是指主程序中引用的函数（子程序）里也要修改嘛？可以问下例如哪个要改嘛

我在运行到“ random effects estimator by ML ”，总是提醒如下错误，能帮忙看看是哪里的问题嘛？

错误使用  *
内部矩阵维度必须一致。

出错 direct_indirect_effects_estimates (line 74)
    parms = chol(results.cov)'*randn(size(parm))+ parm;

出错 Untitled3 (line 259)
direct_indirect_effects_estimates(results,W,spat_model);

请问楼主，为什么我LM检验的结果没有log-l的结果，其他都有。是不是数据有问题

请问我在空间计量工具包中没有发现有杜宾模型的代码，可否请楼主指点一下，在哪里能够运行

楼主你好，我在进行SDEM检验时出现这个错误：
Error using  *
Inner matrix dimensions must agree.

Error in direct_indirect_effects_estimates (line 86)
     parms = chol(results.cov)'*randn(size(parm)) + parm;

Error in model1 (line 34)
direct_indirect_effects_estimates(results,W,spat_model);
请问是什么原因造成的。

同问，请问你解决了吗

可能是标点出错了，或者字母错了，不知道还有没有其他的原因