非常感谢各位高手的回答,向你们致敬

6# iloveyou21
按您的方法程序有所改进，可是还是有问题，如下：错误是
red= log(rd[.,2:23]);
frd= frd(rd[.,2:23]);
error G0064 : Operand missing

library pcoint coint pgraph;
_ker_fun = &fejer;

loadm f[20,23]=c:\examples\tfp.dat;

loadm f[20,23]=c:\examples\r&d.dat;

loadm f[20,23]=c:\examples\fr&d.dat;

loadm f[20,23]=c:\examples\impor.dat;
loadm m[20,23]=c:\examples\tfp.dat;
loadm m[20,23]=c:\examples\r&d.dat;
loadm m[20,23]=c:\examples\fr&d.dat;

loadm m[20,23]=c:\examples\impor.dat;

tfp=f[.,1:9]~f[.,20]~f[.,10]~f[.,12:17]~f[.,19]~f[.,18]~f[.,21:23]~f[.,11];

impt= m[.,1:9]~m[.,20]~m[.,10]~m[.,12:17]~m[.,19]~m[.,18]~m[.,21:23]~m[.,11];

impt1={1970 0.0552 0.0958 0.1914 0.1528 0.1560 0.2145 0.2001 0.1475 0.3081 0.4387 0.3090 0.2678 0.1700 0.4212 0.5002 0.4513 0.2548 0.4535 0.3359 0.1466 0.2282 0.3907};
impt_m=impt1|impt[1:19,.];
tfp=log(tfp[.,2:23]);
red= log(rd[.,2:23]);
frd= frd(rd[.,2:23]);
y=tfp;x=rd~frd;
ht_1(tfp);
ht_1(rd);
ht_1(frd);
{ols_res,b,c,d,e}=ols_rd(y,x);
{fm_res,b,c,d,e}=fm_rd(y,x);
{dols_res,b,c,d,e}=dols_rd(y,x,1,2);
kao_t(ols_res,y,x,1);
pedroni_t(ols_res,y,x);
end;
先谢谢了！

还是那个问题，你用的变量要先定义。前面loadm 都重复了。你的m数据由哪个文件来，则loadm m=。。。
而不是loadm m[20,23]=c:\examples\tfp.dat;
loadm m[20,23]=c:\examples\r&d.dat;
loadm m[20,23]=c:\examples\fr&d.dat;
loadm m[20,23]=c:\examples\impor.dat;
只要一个，就是你存放m数据的那个。我猜测是这样：
loadm tfp[20,23]=c:\examples\tfp.dat;
loadm rd[20,23]=c:\examples\r&d.dat;
loadm f[20,23]=c:\examples\fr&d.dat;
loadm m[20,23]=c:\examples\impor.dat;

9# iloveyou21
非常感谢您的及时指点，我再试试看

9# iloveyou21
非常感谢您的指点！我已经完成了一部分内容，可是程序运行到后面出现了如下的错误：
运行{ols_res,b,c,d,e}=ols_rd(y,x);出现以下错误：（需要说明的是：我已经把该要的程序包已经coint2和npt1.3按说明装到所需要的地方，Gauss右边的模块里面显示有dols_rd,fm_rd,ols_rd)请您再帮忙仔细看下，先谢谢您了！
(0) : error G0159 : Wrong number of parameters
Currently active call: ols_rd
Stack trace:
   ols_rd



运行{fm_res,b,c,d,e}=fm_rd(y,x);出现：
(0) : error G0159 : Wrong number of parameters
Currently active call: FM_RD
Stack trace:
   FM_RD



运行{dols_res,b,c,d,e}=dols_rd(y,x,1,2);出现：
(0) : error G0159 : Wrong number of parameters
Currently active call: DOLS_RD
Stack trace:
   DOLS_RD


运行kao_t(ols_res,y,x,1);出现
(0) : error G0152 : Variable not initialized
运行pedroni_t(ols_res,y,x);出现
pedroni_t(ols_res,y,x);





$ library pcoint coint pgraph panel;
$ _ker_fun = &fejer;
$ loadm f[20,23]=c:\examples\tfp.dat;
$ loadm rd[20,23]=c:\examples\r&d.dat;
$ loadm frd[20,23]=c:\examples\fr&d.dat;
$ loadm fi[20,23]=c:\examples\impor.dat;
$ loadm m[20,23]=c:\examples\tfp.dat;
$
tfp=f[.,1:9]~f[.,20]~f[.,10]~f[.,12:17]~f[.,19]~f[.,18]~f[.,21:23]~f[.,11];
$ impt= m[.,1:9]~m[.,20]~m[.,10]~m[.,12:17]~m[.,19]~m[.,18]~m[.,21:23]~m[.,11];
$ impt= m[.,1:9]~m[.,20]~m[.,10]~m[.,12:17]~m[.,19]~m[.,18]~m[.,21:23]~m[.,11];
$
impt1={1970 0.0552 0.0958 0.1914 0.1528 0.1560 0.2145 0.2001 0.1475 0.3081 0.4387 0.3090 0.2678 0.1700 0.4212 0.5002 0.4513 0.2548 0.4535 0.3359 0.1466 0.2282 0.3907};
$
impt_m=impt1|impt[1:19,.];
$ tfp=log(tfp[.,2:23]);
$ red= log(rd[.,2:23]);
$ frd= log(frd[.,2:23]);
$ y=tfp;x=rd~frd;
$ ht_1(tfp);

----------- The Harris and Tzaralis (1999) Model 1a Unit Root Test -----------

-----------
Without Intercept and Time Trend ---------------
nomalized ht_1_stat
-8.39447 :critical Probability=
0.00000

-----------------------------------------------------------------------
$
ht_1(rd);

----------- The Harris and Tzaralis (1999) Model 1a Unit Root Test -----------

-----------
Without Intercept and Time Trend ---------------
nomalized ht_1_stat
0.03340 :critical Probability=
0.48668

-----------------------------------------------------------------------
$ ht_1(frd);

----------- The Harris and Tzaralis (1999) Model 1a Unit Root Test -----------

-----------
Without Intercept and Time Trend ---------------
nomalized ht_1_stat
-5.11868 :critical Probability=
0.00000

-----------------------------------------------------------------------

{ols_res,b,c,d,e}=ols_rd(y,x);
{fm_res,b,c,d,e}=fm_rd(y,x);

{dols_res,b,c,d,e}=dols_rd(y,x,1,2);
kao_t(ols_res,y,x,1);
pedroni_t(ols_res,y,x);
end;

loadm fi[20,23]=c:\examples\impor.dat;变为loadm m[20,23]=c:\examples\impor.dat;
loadm m[20,23]=c:\examples\tfp.dat;删去，再试试；

楼上的高手，按你的做法弄了不行，我改天把我详细的情况说明给你发Email再问，谢谢了