楼主: shijun19880808
3979 8

[问答] matlab做空间出现问题Wrong # of variable names in prt_reg [推广有奖]

  • 0关注
  • 7粉丝

等待验证会员

硕士生

7%

还不是VIP/贵宾

-

威望
0
论坛币
8834 个
通用积分
0.1800
学术水平
1 点
热心指数
2 点
信用等级
2 点
经验
1813 点
帖子
38
精华
0
在线时间
216 小时
注册时间
2012-11-27
最后登录
2019-4-29

10论坛币
T=13;N=11;y=log(ss(:,[1]));x=log(ss(:,2:5));xconstant=ones(N*T,1);[nobs K]=size(x);
results=ols(y, [xconstant x]);
vnames=strvcat('cons','x1','x2','x3','x4');
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
I=eye(T);W1=kron(I,dd);
res=moran(y,[ xconstant x],W1);
prt(res);
Wrong # of variable names in prt_reg -- check vnames argument
will use generic variable names

Ordinary Least-squares Estimates
R-squared      =    0.4483
Rbar-squared   =    0.4323
sigma^2        =    0.0338
Durbin-Watson  =    2.2329
Nobs, Nvars    =    143,     5
***************************************************************
Variable        Coefficient      t-statistic    t-probability
variable 1         3.458572         6.686364         0.000000
variable 2        -0.056446        -3.452232         0.000738
variable 3         0.141382         8.495703         0.000000
variable 4        -0.392556        -6.554149         0.000000
variable 5         0.209887         3.988067         0.000108

loglikols =
   41.8904
Moran I-test for spatial correlation in residuals
                     
Moran I                    0.32644387
Moran I-statistic          5.33393633
Marginal Probability       0.00000010
mean                      -0.01051809
standard deviation         0.06317323

关键词:Variable MATLAB Names wrong atlab matlab names moran 空间
沙发
黑白电视机 发表于 2014-4-16 12:21:37 |只看作者 |坛友微信交流群
试试将“vnames=strvcat('cons','x1','x2','x3','x4'); ”改成“vnames=strvcat(‘y’,'cons','x1','x2','x3','x4'); ”
难得糊涂

使用道具

藤椅
shijun19880808 发表于 2014-4-16 12:33:35 |只看作者 |坛友微信交流群
黑白电视机 发表于 2014-4-16 12:21
试试将“vnames=strvcat('cons','x1','x2','x3','x4'); ”改成“vnames=strvcat(‘y’,'cons','x1','x2','x ...
哦,好了,谢谢啊

使用道具

板凳
gotosky2007 发表于 2018-5-18 22:02:55 |只看作者 |坛友微信交流群
感谢感谢

使用道具

报纸
珊珊来吃^^ 发表于 2018-6-16 13:33:09 |只看作者 |坛友微信交流群
黑白电视机 发表于 2014-4-16 12:21
试试将“vnames=strvcat('cons','x1','x2','x3','x4'); ”改成“vnames=strvcat(‘y’,'cons','x1','x2','x ...
请问大牛,初学空间计量出现这个错误,怎么解决呀?谢谢!

>> results=ols(y,[xconstant x]);
vnames=strvcat('farmer','mig','out farm','incomefarm','income','gov','land','farmlabor','cooperation','zheng','jiufenfarmer','jiufencoop','zhiliang');
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
未定义函数或变量 'prt_reg'

使用道具

地板
徐子墨 发表于 2019-5-23 08:28:26 |只看作者 |坛友微信交流群
珊珊来吃^^ 发表于 2018-6-16 13:33
请问大牛,初学空间计量出现这个错误,怎么解决呀?谢谢!

>> results=ols(y,[xconstant x]);
...
文件路径没有设置好,用cd"路径名"

使用道具

7
誉满天下. 发表于 2019-5-27 02:01:08 |只看作者 |坛友微信交流群
黑白电视机 发表于 2014-4-16 12:21
试试将“vnames=strvcat('cons','x1','x2','x3','x4'); ”改成“vnames=strvcat(‘y’,'cons','x1','x2','x ...
您好,我也是这个问题Wrong # of variable names in prt_reg -- check vnames argument
will use generic variable names

使用道具

8
誉满天下. 发表于 2019-5-27 02:02:44 |只看作者 |坛友微信交流群
黑白电视机 发表于 2014-4-16 12:21
试试将“vnames=strvcat('cons','x1','x2','x3','x4'); ”改成“vnames=strvcat(‘y’,'cons','x1','x2','x ...
T=11; % number of time periods
N=31; % number of regions
% row-normalize W
W=normw(W1); % function of LeSage
y=A(:,[4]); % column number in the data matrix that corresponds to the dependent variable
x=A(:,[5,11]); % column numbers in the data matrix that correspond to the independent variables
xconstant=ones(N*T,1);
[nobs K]=size(x);
% ----------------------------------------------------------------------------------------
% ols estimation
results=ols(y,[xconstant x]);
vnames=char('CXST','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

% The (robust)LM tests developed by Elhorst

LMsarsem_panel(results,W,y,[xconstant x]); % (Robust) LM tests

% The lm tests developed by Donald Lacombe
% see http://www.rri.wvu.edu/lacombe/~lacombe.htm

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);


结果:
Wrong # of variable names in prt_reg -- check vnames argument
will use generic variable names

Ordinary Least-squares Estimates
R-squared      =    0.0440
Rbar-squared   =    0.0383
sigma^2        =  119.5661
Durbin-Watson  =    0.3193
Nobs, Nvars    =    341,     3
***************************************************************
Variable        Coefficient      t-statistic    t-probability
variable 1         9.906080        11.754521         0.000000
variable 2         0.175004         1.236923         0.216974
variable 3         0.045317         3.524113         0.000483


loglikols =

  -1.2980e+03

LM test no spatial lag, probability          =  143.5858,    0.000
robust LM test no spatial lag, probability   =    0.7348,    0.391
LM test no spatial error, probability        =  144.8356,    0.000
robust LM test no spatial error, probability =    1.9846,    0.159
LM lag test for omitted spatial lag in panel data
                     
LM value                 143.58581315
Marginal Probability       0.00000000
chi(1) .01 value           6.64000000

LM error test for spatial errors in panel data
                     
LM value                 144.83556667
Marginal Probability       0.00000000
chi(1) .01 value           6.64000000

Robust LM lag test for omitted spatial lag in panel data
                     
LM value                   0.73480738
Marginal Probability       0.39132935
chi(1) .01 value           6.64000000

Robust LM error test for spatial errors in panel data
                     
LM value                   1.98456090
Marginal Probability       0.15891075
chi(1) .01 value           6.64000000

使用道具

9
17547972789 学生认证  发表于 2020-7-9 14:40:12 |只看作者 |坛友微信交流群
誉满天下. 发表于 2019-5-27 02:02
T=11; % number of time periods
N=31; % number of regions
% row-normalize W
请问您解决这个问题了吗

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-20 03:27