全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Gauss专版
2011-12-23 12:28:49
CDFWALD是对应每个地区都产生了N个值,N表示模拟次数
对应1%,5%,10%的应该分别只有一个临界值的,如何计算
二维码

扫码加我 拉你入群

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

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

2011-12-24 10:29:41
ywh19860616 发表于 2011-12-23 12:28
CDFWALD是对应每个地区都产生了N个值,N表示模拟次数
对应1%,5%,10%的应该分别只有一个临界值的,如何计 ...
FRML 有两种表示法
tsp50rm.pdf page 173/466
  normalized form
  unnormalized (implicit) form
你说的是属于第二种
Hypothesis testing 都有用到
#######
底下例子取自MSGMm0mx.tsp
PARAM S11 S12 S13 B11 B1 LAMBDA1 d11 d12 d13 d14 d15;
PARAM     S22 S23 B22 B2 LAMBDA2 d21 d22 d23 d24 d25;
PARAM         S33 B33 B3 LAMBDA3 d31 d32 d33 d34 d35;
PARAM             B44 B4 LAMBDA4 d41 d42 d43 d44 d45;
param beta d1 d2;

? Specify the SGM system: 41 parameters
FRML E1 Q1=(S11*Pd1+S12*Pd2+S13*Pd3
          -THETA1*(S11*P11+S12*P12+S13*P13
                          +S22*P22+S23*P23
                                  +S33*P33))*(Y1+BETA*Y2)
     +d11*dum1+d12*dum2+d13*dum3+d14*dum4+d15*dum5  
     +B11*(Y1+BETA*Y2)+B1+LAMBDA1*(Y1*Y1+2*D1*Y1*Y2+D2*Y2*Y2);
FRML E2 Q2=(S12*Pd1+S22*Pd2+S23*Pd3
          -THETA2*(S11*P11+S12*P12+S13*P13
                          +S22*P22+S23*P23
                                  +S33*P33))*(Y1+BETA*Y2)
     +d21*dum1+d22*dum2+d23*dum3+d24*dum4+d25*dum5     
     +B22*(Y1+BETA*Y2)+B2+LAMBDA2*(Y1*Y1+2*D1*Y1*Y2+D2*Y2*Y2);
FRML E3 Q3=(S13*Pd1+S23*Pd2+S33*Pd3
          -THETA3*(S11*P11+S12*P12+S13*P13
                          +S22*P22+S23*P23
                                  +S33*P33))*(Y1+BETA*Y2)
     +d31*dum1+d32*dum2+d33*dum3+d34*dum4+d35*dum5  
     +B33*(Y1+BETA*Y2)+B3+LAMBDA3*(Y1*Y1+2*D1*Y1*Y2+D2*Y2*Y2);
FRML E4 Q4=-((S11+S12+S13)*Pd1+(S12+S22+S23)*Pd2+(S13+S23+S33)*Pd3
           +THETA4*(S11*P11+S12*P12+S13*P13
                           +S22*P22+S23*P23
                                   +S33*P33))*(Y1+BETA*Y2)
     +d41*dum1+d42*dum2+d43*dum3+d44*dum4+d45*dum5   
     +B44*(Y1+BETA*Y2)+B4+LAMBDA4*(Y1*Y1+2*D1*Y1*Y2+D2*Y2*Y2);
     
? Restrictions yielding the single-output model: beta=d1=d2=0
frml c1 d1;
frml c2 d2;
frml c3 beta;
######
底下例子取自analyzr.tsp
?Y = a0 + b1*K + b2*L + b3*M + b4*Q + b5*F + b6*B + b7*Q2 + b8*F2 +
?b9*B2 + b10*QF + b11*QB + b12*FB

olsq(silent) Y C K L M Q F B Q2 F2 B2 QF QB FB ;

? write restrictions in terms of variable names, using unnormalized FRMLs
?R1: b4 + b5 + b6 = 1
frml R1 Q + F + B - 1;  ? expression equals zero
?R2: b7 + b10 + b11 = 0
frml R2 Q2 + QF + QB;
?R3: b8 + b10 + b12 = 0
frml R3 F2 + QF + FB;
?R4: b9 + b11 + b12 = 0
frml R4 B2 + QB + FB;

analyz(silent) r1-r4;  ? test restrictions and compute restricted coefs

其它问题再逐步回答
二维码

扫码加我 拉你入群

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

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

2011-12-24 12:51:23
epoh 发表于 2011-12-24 10:29
FRML 有两种表示法
tsp50rm.pdf page 173/466
  normalized form
谢谢epoh老师,这个问题经过您解释明白了
的确这是为后文做wald检验做铺垫的


但是我仔细看了程序,都没有看到哪里可以通过设置1%,5%,10%水平而获得临界值。


二维码

扫码加我 拉你入群

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

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

2011-12-25 15:25:29
epoh 发表于 2011-12-24 10:29
FRML 有两种表示法
tsp50rm.pdf page 173/466
  normalized form
epoh老师,如果是在常见的分布中,要计算某一分位水平下的临界值有现成的函数,比如:
在卡方分布表中,若自由度为10, alpha=0.975,求临界值Lambda。
解:因为表中给出的值满足p(x2>lambda),而逆累积分布函数icdf求满足p(x2<lambda)的临界值 。所以,这里的alpha取为0.025,即
Lambda=icdf('chi2',0.025,10)

而我现在的wald检验是经过自己构造出来的,现在有的只是cdfwald,对应抽样中每个观察值都对应一个值。比如抽样100次,得出100个cdfwald。
这时候如何求解临界值呢?


二维码

扫码加我 拉你入群

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

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

2011-12-25 21:09:54
ywh19860616 发表于 2011-12-25 15:25
epoh老师,如果是在常见的分布中,要计算某一分位水平下的临界值有现成的函数,比如:
在卡方分布表中, ...
Bootstrap panel Granger-causality between
government spending and revenue in the EU.pdf
   
老兄看完这篇文献的Table1~ Table2b
大概就清楚了.
程序好像是跟文献配合的一样

二维码

扫码加我 拉你入群

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

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

2011-12-25 22:12:38
epoh 发表于 2011-12-25 21:09
Bootstrap panel Granger-causality between
government spending and revenue in the EU.pdf
epoh老师,谢谢您,方法是和这篇文章一样的
但是我还是没有明白如何获得table1a-table2b
二维码

扫码加我 拉你入群

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

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

2011-12-25 22:57:40
ywh19860616 发表于 2011-12-25 22:12
epoh老师,谢谢您,方法是和这篇文章一样的
但是我还是没有明白如何获得table1a-table2b
Table 1
                                                      Bootstrap critical values
             coefficient  Test Statistic  1%        5%        10%
Austria   0.1351        1.2361       26.5043 15.5115 11.5606

......
应该如同 ,我设 ntrial =1000
                                                           Bootstrap critical values
                   coefficient  Test Statistic      1%         5%        10%
CONSTR1    .315681   0.401457291     53.3929 31.4970  22.3084
......
详细可以参照
ywhwald_result
   
ywhwald_result.txt
大小:(729.85 KB)

 马上下载


补充说明过程:

假设 ntrial=10000,

bootstrap(bs)10000次之后的wald stat

存放在waldstatm  (waldstatm(bs,i)=@WALD;)

然后

UNMAKE waldstatm ws1-ws29;

SORT ws.;

建立cdfwald

   DO pc=1 TO 100;

     SELECT ttrial=ceil(pc*ntrial/100);

     SET cdfwald(pc,i)=ws.;

   ENDDO;

  所以排序后(由小大大)

  第100个是 1%

  第9500个是 95%

  第9900个是 99%


二维码

扫码加我 拉你入群

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

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

2011-12-27 07:25:32
epoh 发表于 2011-12-25 22:57
Table 1
                                                      Bootstrap critical values
      ...
epoh老师,谢谢您
昨天一天都在忙别的事情,没有及时回复,呵呵
我先自己试加了一下后面的几句命令,就是sort以后取出第950个值,不懂再请教您
二维码

扫码加我 拉你入群

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

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

2011-12-27 22:06:26
ywh19860616 发表于 2011-12-27 07:25
epoh老师,谢谢您
昨天一天都在忙别的事情,没有及时回复,呵呵
我先自己试加了一下后面的几句命令,就 ...
建议你存成文件,可能比较方便.
write(file ='cdfwald.xls') cdfwald;

二维码

扫码加我 拉你入群

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

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

2011-12-28 12:09:58
epoh 发表于 2011-12-27 22:06
建议你存成文件,可能比较方便.
write(file ='cdfwald.xls') cdfwald;
好的,谢谢epoh老师
二维码

扫码加我 拉你入群

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

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

2011-12-31 09:07:21
epoh 发表于 2011-12-27 22:06
建议你存成文件,可能比较方便.
write(file ='cdfwald.xls') cdfwald;
epoh老师,我还想问下程序中的@NOB表示什么意思
按参考手册上解释是Number of observations,而我的运行结果不会这样
我有N=12,t=20,即20个个体,12年数据,为什么运行出@NOB=19.0000呢?我哪里理解错了?
还有,利用FRML命令时,若我想在方程中加入时间趋势变量t,应该如何存在?

程序中有这样一句(133 line):select ttrial=ceil(pc*ntrial/100),此处100是不是指抽样的间隔?
二维码

扫码加我 拉你入群

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

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

2011-12-31 10:27:31
epoh老师,还有一个问题,您有空帮忙看看,谢谢了
原程序只是滞后一阶,而我现在想修改滞后阶数,出现了很多错误,找不出原因
比如原来程序: FRML eq. ys.=a.+b1.*ys(-1)+c1.*xs.(-1)
我想取滞后2阶: FRML eq. ys.=a.+b1.*ys(-1)+b2.*ys(-2)+c1.*xs.(-1)+c2.*xs(-2)   
后面有些也对应修改
比如在  FRML constr. c1.;  
后加了一句:
FRML constr1. c2.;  
epoh老师,这两个个约束应该不是同时c1=c2=0吧?
如果要2个同时为0,应该如何?

但是出现很多错误,您看可以解决吗?



附件是您上传的文章一段,想要实现two-way granger causality



附件列表
未命名.jpg

原图尺寸 33.46 KB

未命名.jpg

二维码

扫码加我 拉你入群

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

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

2011-12-31 11:56:26
ywh19860616 发表于 2011-12-31 10:27
epoh老师,还有一个问题,您有空帮忙看看,谢谢了
原程序只是滞后一阶,而我现在想修改滞后阶数,出现了很 ...
select ttrial=ceil(pc*ntrial/100),
此处100跟抽样无关
假设ntrial=100,而pc=1 TO 100;
当pc=95,pc*ntrial/100=95,第95个就是wald95
当pc=99,pc*ntrial/100=99,第99个就是wald99
同理ntrial=10000,而pc=1 TO 100;
当pc=95,pc*ntrial/100=9500,第9500个就是wald95
当pc=99,pc*ntrial/100=9900,第9900个就是wald99

假设你不想把cdfwald存成文件,
也可以改成底下:
MAT tcdfwald=cdfwald';
print tcdfwald;
unmake tcdfwald wald1-wald100;
print wald90;
print wald95;
print wald99;
#######
滞后阶数修改,这两天我试试.
#######
freq(panel,n=29,t=38,id=country,time=year) n;
set nt = 29*38;
smpl 1,nt;
read(file='c:\wald1.xls');
print @NOB;
  @NOB =  1102.00000

######
TREND generates a series with a linear growth trend.
      The trend may be repeated under the control of several options.

ywhtsp.tsp
FREQ A;
SMPL fyear lyear;
TREND t;

                     T
1960           1.00000
1961           2.00000
1962           3.00000
1963           4.00000
1964           5.00000
1965           6.00000
1966           7.00000
1967           8.00000
1968           9.00000
1969          10.00000
1970          11.00000
1971          12.00000
1972          13.00000
1973          14.00000
1974          15.00000
1975          16.00000
1976          17.00000
1977          18.00000
1978          19.00000
1979          20.00000
1980          21.00000
1981          22.00000
1982          23.00000
1983          24.00000
1984          25.00000
1985          26.00000
1986          27.00000
1987          28.00000
1988          29.00000
1989          30.00000
1990          31.00000
1991          32.00000
1992          33.00000
1993          34.00000
1994          35.00000
1995          36.00000
1996          37.00000
1997          38.00000
二维码

扫码加我 拉你入群

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

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

2011-12-31 12:18:32
epoh 发表于 2011-12-31 11:56
select ttrial=ceil(pc*ntrial/100),
此处100跟抽样无关
假设ntrial=100,而pc=1 TO 100;
好的,谢谢epoh老师。真是打扰您了.现在的主要问题应该是如何实现c1=c2=0的设置以及最后面的PROC BSTR的ystarm(u,v)表达式

祝epoh老师元旦快乐。
二维码

扫码加我 拉你入群

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

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

2011-12-31 20:25:49
ywh19860616 发表于 2011-12-31 12:18
好的,谢谢epoh老师。真是打扰您了.现在的主要问题应该是如何实现c1=c2=0的设置以及最后面的PROC BSTR的 ...
老兄这个也应该可以由
R package “plm"来做
因为只要有coefficient covariance matrix
就可带入Wald statistic formula
详细可看section 3.4. Inference in the panel model
from the plm vignette. Page 7/48
你只要把plm式子列出即可
顺祝 老兄新年快乐!
二维码

扫码加我 拉你入群

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

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

2011-12-31 21:44:42
epoh 发表于 2011-12-31 20:25
老兄这个也应该可以由
R package “plm"来做
因为只要有coefficient covariance matrix
恩,好的,谢谢epoh老师,我再仔细研究下TSP程序
因为TSP和R交叉使用,可能比较困难
二维码

扫码加我 拉你入群

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

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

2011-12-31 22:03:04
ywh19860616 发表于 2011-12-31 21:44
恩,好的,谢谢epoh老师,我再仔细研究下TSP程序
因为TSP和R交叉使用,可能比较困难
我的意思是直接在R 作

二维码

扫码加我 拉你入群

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

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

2011-12-31 23:54:16
epoh 发表于 2011-12-31 22:03
我的意思是直接在R 作
epoh老师,嗯,我明白您的想法了
经过您修改,1阶滞后的已经可以运行了,我想在原程序基础上修改到2阶滞后应该也是可行的,比如您发的那篇文章就用到了
如果完全用R实现,那所有程序可能都得重新编写过,也许更困难些
二维码

扫码加我 拉你入群

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

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

2012-1-1 20:41:54
epoh 发表于 2011-12-31 22:03
我的意思是直接在R 作
epoh老师,我想问下这两段语句有何区别呢?
FRML constr1. c1.;
FRML constr2 c2.;
ANALYZ  constr1.;
ANALYZ constr2.;

FRML constr1. c1.;
FRML constr2 c2.;
ANALYZ  constr1. constr2.;


二维码

扫码加我 拉你入群

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

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

2012-1-2 11:08:40
ywh19860616 发表于 2012-1-1 20:41
epoh老师,我想问下这两段语句有何区别呢?
FRML constr1. c1.;
FRML constr2 c2.;
TSP最麻烦的地方在
产生错误时没说出原因
很不方便除错
LAG2修改后,经检查SAMPLE,EQ都对
但是SUR说找不到YS,XS,困扰思考中.
#######
Current sample:  1962 to 1997
                                   
EQ = EQ1 EQ2 EQ3 EQ4 EQ6 EQ7 EQ8 EQ10 EQ12 EQ13 EQ14 EQ15 EQ16
    EQ17 EQ18 EQ19 EQ20 EQ21 EQ23 EQ24 EQ25 EQ26 EQ28 EQ29

EQUATION: EQ1

     FRML EQ1 YS1 = A1 + B11*YS(-1) + B21*YS(-2) + C11*XS1(-1) + C21*XS(-2)

EQUATION: EQ2

     FRML EQ2 YS2 = A2 + B12*YS(-1) + B22*YS(-2) + C12*XS2(-1) + C22*XS(-2)

C:\Program Files\TSP 5.0\ywhtsp_lag2.tsp (113):  ERROR
*** ERROR in command 81 Procedure SUR: Illegal variable type or variable
    missing ====>  YS

C:\Program Files\TSP 5.0\ywhtsp_lag2.tsp (113):  ERROR
*** ERROR in command 81 Procedure SUR: Illegal variable type or variable
    missing ====>  YS

C:\Program Files\TSP 5.0\ywhtsp_lag2.tsp (113):  ERROR
*** ERROR in command 81 Procedure SUR: Illegal variable type or variable
    missing ====>  XS
二维码

扫码加我 拉你入群

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

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

2012-1-2 12:02:41
epoh 发表于 2012-1-2 11:08
TSP最麻烦的地方在
产生错误时没说出原因
很不方便除错
epoh老师,您上面写的XS1 是表示什么?
可能我没有想到那么复杂,我的想法是现在要同时检验XS(-1)和XS(-2)前的系数是否同时为0,如c1==c2==0,如果在winrats,都有直接test函数可用,而TSP中不知道如何使用。因此,我才问
ANALYZ constr1 constr2
与分开两个命令写
ANALYZ constr1
ANALYZ constr2

有何区别?
二维码

扫码加我 拉你入群

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

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

2012-1-2 14:51:32
ywh19860616 发表于 2012-1-2 12:02
epoh老师,您上面写的XS1 是表示什么?
可能我没有想到那么复杂,我的想法是现在要同时检验XS(-1)和XS(- ...
哈哈!老兄原来你才是罪魁祸首.
42楼你写错了
FRML eq. ys.=a.+b1.*ys(-1)+b2.*ys(-2)+c1.*xs.(-1)+c2.*xs(-2)
应该是  
FRML eq. ys.=a.+b1.*ys.(-1)+b2.*ys.(-2)+c1.*xs.(-1)+c2.*xs.(-2);
这才折腾了半天.呵呵
现在程序可以运行了,你试试.
请注意短信息

不过如果要做bootstrap
PROC RESIDNUL beta eps;   
还需要修改
因为原来是
FRML eq0. ys.=a0.+b01.*ys.(-1);
PARAM a0. b01.;

*******************
FRML constr1. c1.;
FRML constr2 c2.;
ANALYZ  constr1.;
ANALYZ constr2.;

FRML constr1. c1.;
FRML constr2 c2.;
ANALYZ  constr1. constr2.;
是不一样的
底下例子(MSGM_modified)可以说明:
MSGM_modified.tsp,HedPrice1.dat,Team.dat
  
MSGM_modified.rar
大小:(12.44 KB)

 马上下载


frml c1 d1;
frml c2 d2;
frml c3 beta;

? Benchmark estimation with Wald test of separability
LSQ(MAXIT=500,TOL=0.0001) E1 E2 E3 E4;
set loglu=@logl;
analyz c1 c2 c3;

Results of Parameter Analysis
=============================

                         Standard
Parameter  Estimate        Error       t-statistic   P-value
C1         13.4943       18.7142       .721072       [.471]
C2         292.864       214.771       1.36361       [.173]
C3         18.7071       5.50715       3.39687       [.001]

Wald Test for the Hypothesis that the given set of Parameters are jointly zero:

CHISQ(3) =   19.941252     ; P-value = 0.00017

*******************

LSQ(MAXIT=500,TOL=0.0001) E1 E2 E3 E4;
set loglu=@logl;
analyz c1 ;
analyz c2 ;
analyz c3 ;

Results of Parameter Analysis
=============================


                         Standard
Parameter  Estimate        Error       t-statistic   P-value
C1         13.4940       18.7140       .721067       [.471]

Wald Test for the Hypothesis that the given set of Parameters are jointly zero:

CHISQ(1) =  0.51993751     ; P-value = 0.47087

Results of Parameter Analysi
=============================

                         Standard
Parameter  Estimate        Error       t-statistic   P-value
C2         292.867       214.769       1.36363       [.173]

Wald Test for the Hypothesis that the given set of Parameters are jointly zero:

CHISQ(1) =   1.8594975     ; P-value = 0.17268

Results of Parameter Analysis
=============================

                         Standard
Parameter  Estimate        Error       t-statistic   P-value
C3         18.7071       5.50717       3.39687       [.001]

Wald Test for the Hypothesis that the given set of Parameters are jointly zero:

CHISQ(1) =   11.538697     ; P-value = 0.00068
二维码

扫码加我 拉你入群

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

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

2012-1-2 15:48:33
epoh 发表于 2012-1-2 14:51
哈哈!老兄原来你才是罪魁祸首.
42楼你写错了
FRML eq. ys.=a.+b1.*ys(-1)+b2.*ys(-2)+c1.*xs.(-1)+c2 ...
哈哈,epoh老师,是我自己写错了,罪过罪过了
如果
FRML c1;
FRML c2;
FRML c3;
然后
ANALYZ  c1 c2 c3;
这句命令是否可以认为是test c1==c2==c3==0?
二维码

扫码加我 拉你入群

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

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

2012-1-2 16:16:45
ywh19860616 发表于 2012-1-2 15:48
哈哈,epoh老师,是我自己写错了,罪过罪过了
如果
FRML c1;
ANALYZ  c1 c2 c3;
这句命令是否可以认为是test c1==c2==c3==0?
对!就是.
另短信息提到的
我再帮你修改.
二维码

扫码加我 拉你入群

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

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

2012-1-2 16:23:45
epoh 发表于 2012-1-2 16:16
ANALYZ  c1 c2 c3;
这句命令是否可以认为是test c1==c2==c3==0?
对!就是.
好的,谢谢epoh老师
您帮忙看看最后面的这句,如果lag length 改变,是否需要变化?
SET ystarm(u,v)=beta(v,1)+beta(v,2)*ystarm(u1,v)+epsstarm(u,v);
二维码

扫码加我 拉你入群

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

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

2012-1-2 19:50:59
ywh19860616 发表于 2012-1-2 16:23
好的,谢谢epoh老师
您帮忙看看最后面的这句,如果lag length 改变,是否需要变化?
SET ystarm(u,v)=b ...
SET u1=u-1;
SET u2=u-2;
    DO v=1 TO ncountry;
SET ystarm(u,v)=beta(v,1)+beta(v,2)*ystarm(u1,v)+beta(v,3)*ystarm(u2,v)+epsstarm(u,v);
    ENDDO;

#####Results:
Current sample:  1962 to 1997
                  @NOB        @NCOEF             P
Value         36.00000     120.00000       5.00000
                  UAIC         USBIC           AIC          SBIC
Value        -19.62815     -14.34976      -0.10309       0.11685
                   C11           C21
Value          1.84125      -1.50339
                   C12           C22
Value         0.094831       0.26454
                   C13           C23
Value         -0.66898      -0.36604
.....
.....
                  C128          C228
Value          0.55224      0.048596
                  C129          C229
Value         -0.23415       0.83328
*************
Results of Parameter Analysis
=============================
                                     Standard
Parameter    Estimate        Error          t-statistic      P-value
CONSTR11   1.84125       .837120        2.19950      [.028]
CONSTR21   -1.50339      .798314       -1.88321      [.060]

Wald Test for the Hypothesis that the given set of Parameters are jointly zero:

CHISQ(2) =   4.9251172     ; P-value = 0.08522
.....
.....
                                       Standard
Parameter       Estimate        Error         t-statistic   P-value
CONSTR129  -.234145      .433760      -.539804      [.589]
CONSTR229  .833276       .435543       1.91319       [.056]
二维码

扫码加我 拉你入群

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

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

2012-1-2 20:25:08
epoh 发表于 2012-1-2 19:50
SET u1=u-1;
SET u2=u-2;
    DO v=1 TO ncountry;
epoh老师,我运行您修改的程序会出现错误,为何?
错误提示:
183      ENDDO;
*** ERROR in command 183: Operator required before item ====>   ENDDO
184  ENDPROC;
*** WARNING in command 184: ENDP statement encountered but there are
    unfinished DO loop(s) ====>    1
185  ;END;
*** WARNING in command 185: END  statement encountered but there are
    unfinished DO loop(s) ====>    1

您修改程序有加上
SET u1=u-1;
SET u2=u-2;
    DO v=1 TO ncountry;
SET ystarm(u,v)=beta(v,1)+beta(v,2)*ystarm(u1,v)+beta(v,3)*ystarm(u2,v)+epsstarm(u,v);
    ENDDO;
这段再运行吗?
二维码

扫码加我 拉你入群

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

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

2012-1-2 20:57:22
ywh19860616 发表于 2012-1-2 20:25
epoh老师,我运行您修改的程序会出现错误,为何?
错误提示:
183      ENDDO;
可能贴在短信息被截断了
我发在这里好了
ywhtsp_lag2.tsp (need pw)
  
ywhtsp_lag2.rar
大小:(1.66 KB)

 马上下载



二维码

扫码加我 拉你入群

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

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

2012-1-2 21:07:39
恩,epoh老师,我试试看
好像还是有错误,也许是我运行操作有问题
我是直接input 程序名

二维码

扫码加我 拉你入群

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

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

2012-1-3 08:30:00
ywh19860616 发表于 2012-1-2 21:07
恩,epoh老师,我试试看
好像还是有错误,也许是我运行操作有问题
我是直接input 程序名
DOT(VALUE=i) 1-29;
  IF choice(i)=1; THEN;
    DO;
    .....
    .....
      IF i=1; THEN;
         LIST eq eq.;
      ELSE;
         LIST eq eq eq.;
    ENDDO;
ENDDOT;

漏掉回答这个问题
i  1 ~ 29
when i = 1 ,
执行这一句: LIST eq eq.;
        eq eq1
否则执行 LIST eq eq eq.;
i=2 , eq eq1 eq2
i=3 , eq eq1 eq2 eq3
...
i=5 , 都不执行
二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群