您好 请问您的这个问题解决了吗

题主及各位同仁，我今天碰到同样的问题，一样的提示错。我已经觉得，特来共享。
1.这是软件问题，不是数据问题。
2.具体出问题原因不清楚，但知道和格式有关系，似乎和回车有关系。
解决办法：
回到demo文件，从新复制出来到txt文本里面修改，然后再帖进去，不要用word。
下面是我成功解决的代码：
T=11;
N=67;
W=normw(W);
y=inc(:,1);
x=inc(:,[2:7]);
xconstant=ones(N*T,1);
[nobs K]=size(x);
%最小二乘法
results=ols(y,[xconstant x]);
vnames=char('y','intercept','x1','x2','x3','x4','x5','x6','x7');
prt_reg(results,vnames,1);
sige=results.sige*((nobs-K)/nobs);
loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid;

%空间依赖性的LM检验

LMsarsem_panel(results,W,y,[xconstant x]);

lm1=lmlag_panel(y,[xconstant x],W);
prt_tests(lm1);

lm2=lmerror_panel(y,[xconstant x],W);
prt_tests(lm2);

lm3=lmlag_robust_panel(y,[xconstant x],W);
prt_tests(lm3);

lm4=lmerror_robust_panel(y,[xconstant x],W);
prt_tests(lm4);

%空间固定效应的LM检验
model=1;
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
results=ols(ywith,xwith);
vnames=strvcat('y', 'x1','x2','x3','x4','x5','x6','x7');
prt_reg(results,vnames);
sfe=meanny-meannx*results.beta; % including the constant term
yme = y - mean(y);
et=ones(T,1);
error=y-kron(et,sfe)-x*results.beta;
rsqr1 = error'*error;
rsqr2 = yme'*yme;
FE_rsqr2 = 1.0 - rsqr1/rsqr2
sige=results.sige*((nobs-K)/nobs);
logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
LMsarsem_panel(results,W,ywith,xwith); % (Robust) LM tests

lm1=lmlag_panel(ywith,xwith,W);
prt_tests(lm1);

lm2=lmerror_panel(ywith,xwith,W);
prt_tests(lm2);

lm3=lmlag_robust_panel(ywith,xwith,W);
prt_tests(lm3);

lm4=lmerror_robust_panel(ywith,xwith,W);
prt_tests(lm4);
%----------------------------------------------------------------
%时间固定的LM检验
model=2;
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
results=ols(ywith,xwith);
vnames=strvcat('y', 'x1','x2','x3','x4','x5','x6','x7'); % should be changed if x is changed
prt_reg(results,vnames);
tfe=meanty-meantx*results.beta; % including the constant term
yme = y - mean(y);
en=ones(N,1);
error=y-kron(tfe,en)-x*results.beta;
rsqr1 = error'*error;
rsqr2 = yme'*yme;
FE_rsqr2 = 1.0 - rsqr1/rsqr2 % r-squared including fixed effects
sige=results.sige*((nobs-K)/nobs);
logliktfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
LMsarsem_panel(results,W,ywith,xwith); % (Robust) LM tests

lm1=lmlag_panel(ywith,xwith,W);
prt_tests(lm1);

lm2=lmerror_panel(ywith,xwith,W);
prt_tests(lm2);

lm3=lmlag_robust_panel(ywith,xwith,W);
prt_tests(lm3);

lm4=lmerror_robust_panel(ywith,xwith,W);
prt_tests(lm4);
%时空双固定的LM检验
model=3;
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
results=ols(ywith,xwith);
vnames=strvcat('y', 'x1','x2','x3','x4','x5','x6','x7'); % should be changed if x is changed
prt_reg(results,vnames);
intercept=mean(y)-mean(x)*results.beta;
sfe=meanny-meannx*results.beta-kron(en,intercept);
tfe=meanty-meantx*results.beta-kron(et,intercept);
yme = y - mean(y);
ent=ones(N*T,1);
error=y-kron(tfe,en)-kron(et,sfe)-x*results.beta-kron(ent,intercept);
rsqr1 = error'*error;
rsqr2 = yme'*yme;
FE_rsqr2 = 1.0 - rsqr1/rsqr2 % r-squared including fixed effects
sige=results.sige*((nobs-K)/nobs);
loglikstfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid

LMsarsem_panel(results,W,ywith,xwith); % (Robust) LM tests

lm1=lmlag_panel(ywith,xwith,W);
prt_tests(lm1);

lm2=lmerror_panel(ywith,xwith,W);
prt_tests(lm2);

lm3=lmlag_robust_panel(ywith,xwith,W);
prt_tests(lm3);

lm4=lmerror_robust_panel(ywith,xwith,W);
prt_tests(lm4);

% -----------------------------------------
%固定效应的LR显著性检验
LR=-2*(logliktfe-loglikstfe);
dof=N;
probability=1-chis_prb(LR,dof);
% Note: probability > 0.05 implies rejection of spatial fixed effects
fprintf(1,'LR-test joint significance spatial fixed effects, degrees of freedom and probability = %9.4f,%6d,%9.4f \n',LR,dof,probability);
LR=-2*(logliksfe-loglikstfe);
dof=T;
probability=1-chis_prb(LR,dof);
% Note: probability > 0.05 implies rejection of spatial fixed effects
fprintf(1,'LR-test joint significance time-periode fixed effects, degrees of freedom and probability = %9.4f,%6d,%9.4f \n',LR,dof,probability);

补充：好像是X个数太多导致的！！！！

我的demo就是从text中copy的，也是出错了，具体要怎么修改呢？我的是空间计量，和你的编程还不太一样，可否告知一下，谢谢

仔细检查数据的共线性！！！！

O(∩_∩)O谢谢

我粘过去你的代码 改完之后还是出现之前的问题 但是自变量选择三个就没有问题

出现同样的问题，解释变量选择两个的话就没有问题，多加一个解释变量就出现了问题，例如我的三个解释变量会出现：
代码：T=10;
N=16;
W=normw(W1);
y=A(:,3);
x=A(:,[4,6]);
xconstant=ones(N*T,1);
[nobs K]=size(x);
results=ols(y,[xconstant x]);
vnames=strvcat('y','intercept','x1','x2','x3');
prt_reg(results,vnames,1);
sige=results.sige*((nobs-K)/nobs);
loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
% The (robust)LM tests developed by Elhorst
LMsarsem_panel(results,W,y,[xconstant x]); % (Robust) LM tests
结果出现了：
Wrong # of variable names in prt_reg -- check vnames argument
will use generic variable names

Ordinary Least-squares Estimates
R-squared      =    0.2004
Rbar-squared   =    0.1902
sigma^2        =    0.0017
Durbin-Watson  =    0.4011
Nobs, Nvars    =    160,     3
***************************************************************
Variable        Coefficient      t-statistic    t-probability
variable 1         0.078386        16.458904         0.000000
variable 2         0.000038         4.919880         0.000002
variable 3        -0.000347        -4.332922         0.000026


loglikols =

   286.1946

LM test no spatial lag, probability          =  247.2800,    0.000
robust LM test no spatial lag, probability   =  181.9879,    0.000
LM test no spatial error, probability        =  169.5590,    0.000
robust LM test no spatial error, probability =  104.2669,    0.000
本文来自: 人大经济论坛 MATLAB等数学软件专版 版，详细出处参考： https://bbs.pinggu.org/forum.php?mod=viewthread&tid=1345696&page=2&from^^uid=8219544

应该是多重共线性问题导致的，我去除了导致多重共线性的变量后就不报错了

你好  能分享一下空间计量  matlab 里里面这个 prt_tests函数吗 我运行提示找不到该函数

你好  能分享一下空间计量  matlab 里里面这个 prt_tests函数吗 我运行提示找不到该函数