程序运行问题

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是不是指抽样的间隔？

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

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

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老师元旦快乐。

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式子列出即可
顺祝 老兄新年快乐!

epoh 发表于 2011-12-31 20:25
老兄这个也应该可以由
R package “plm＂来做
因为只要有coefficient covariance matrix

恩，好的，谢谢epoh老师，我再仔细研究下TSP程序
因为TSP和R交叉使用，可能比较困难

ywh19860616 发表于 2011-12-31 21:44
恩，好的，谢谢epoh老师，我再仔细研究下TSP程序
因为TSP和R交叉使用，可能比较困难

我的意思是直接在R 作

epoh 发表于 2011-12-31 22:03
我的意思是直接在R 作

epoh老师，嗯，我明白您的想法了
经过您修改，1阶滞后的已经可以运行了，我想在原程序基础上修改到2阶滞后应该也是可行的，比如您发的那篇文章就用到了
如果完全用R实现，那所有程序可能都得重新编写过，也许更困难些

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

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