求助：请各位大牛看看，我用matlab做回归为什么会出现以下错误？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

clear all;
A=wk1read('D:\MATLAB701\work\基础设施\模型检验\lnall2.wk1',1,0);
W1=wk1read('D:\MATLAB701\work\基础设施\模型检验\W1.wk1');
T=16;
N=31;
W=normw(W1);
y=A(:,[3]);
x=A(:,[4,5,6,7,8,9,10,11,12]);
xconstant=ones(N*T,1);
[nobs K]=size(x);
>> results=ols(y,[xconstant x]);
vnames=strvcat('lninv','lncit','lnedu','lnune','lnroa','lnrai','lnktr','lnkqi','lnlab');
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
LMsarsem_panel(results,W,y,[xconstant x]);
Warning: Matrix is singular to working precision.
> In ols at 51
Wrong # of variable names in prt_reg -- check vnames argument
will use generic variable names

Ordinary Least-squares Estimates
R-squared      =       NaN
Rbar-squared   =       NaN
sigma^2        =       NaN
Durbin-Watson  =       NaN
Nobs, Nvars    =    496,    10
***************************************************************
??? Error using ==> betainc
X must be in the interval [0,1].

Error in ==> tdis_prb at 35
tmp = 1.0 - 0.5*betainc(x2,0.5*n,0.5);

Error in ==> prt_reg at 364
tout = tdis_prb(results.tstat,nobs-nvar); % find t-stat probabilities

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏4 回帖

关键词：MATLAB matla atlab Lab Atl matlab Matrix 模型

你应该挂上数据，让别人给你执行一下。

我也遇上同样情况。求高人指导[sweat]

kjfx2014 发表于 2015-4-26 21:48
我也遇上同样情况。求高人指导

同上回答

x=A(:,[4,5,6,7,8,9,10,11,12]);  是不是把这个改成x=A(:,[4,12]);

请问您的问题解决了吗

kjfx2014 发表于 2015-4-26 21:48
我也遇上同样情况。求高人指导

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

题主及各位同仁，我今天碰到同样的问题，一样的提示错。我已经觉得，特来共享。
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个数太多导致的！！！！

huangli 发表于 2017-2-18 21:04
补充：好像是X个数太多导致的！！！！

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