| 查看: 354 | 回复: 0 | ||
swordsmen木虫 (小有名气)
|
[求助]
请高手帮我看下这个GAUSS程序
|
|
new; FORMAT /M1 /ROS 8,4; library pgraph; _pdate=""; _pcolor = 1; _pmcolor = ZEROS(8,1)|15; /* path1 = "Z:\\Stephan.Popp\\Projekte\\perronmultiple\\gauss\\output.out"; OUTPUT FILE = ^path1 ON; */ "******************************************"; "date: " datestr(date) " time: " timestr(time); "*******************************************"; /* DATA ** the first column should contain the years or quarters (i.e. 1973.3) ** and the second the observations */ load yy[1456,2]=d:\workgauss\brentw.txt; @ change path settings and yy[here include # of observations ,2]@ /* Randomly chosen dataset yy=SEQA(1900,1,100)~RNDN(100,1); */ XY(yy[.,1],yy[.,2]); yyy= yy[.,2]; kmax = 5; @ max lag @ ttt = rows(yyy); tau = 0.2; @ trimming factor @ "Series: Log(.)"; "Sample:" yy[1,1]~yy[rows(yyy),1]; "# observations:" ttt; "maximum lag " kmax; "trimming factor " tau; "program: popp2break.prg" ; "*******************************"; tbunter = MAXC(3+kmax|CEIL(tau*ttt)); @ lower break date @ "break date floor, effective tau " tbunter~yy[tbunter,1]~tbunter/ttt; tbober = MINC(ttt-3-kmax|FLOOR((1-tau)*ttt)); @ upper break date @ "break date ceiling, effective (1-tau) " tbober~yy[tbober,1]~tbober/ttt; "*******************************************"; /************** Model 0 *****************/ ergeb = ZEROS(2,1); @ Vektor zum Speichern der Ergebnisse @ tbopt = ZEROS(2,1); ttb1 = tbunter; ttb2 = 0; do while ttb1 <= tbober; @ Schleife der durchlaufenden Bruchzeitpunkte T_B @ {rho1,trho1,ttheta1,pp1,rres1,varres1,tstat1,kk1} = mio2break0b(yyy,ttb1,ttb2,kmax); if ABS(ttheta1) > ABS(ergeb[1,1]); ergeb[1,1] = ttheta1; ergeb[2,1] = ttb1; endif; ttb1 = ttb1 + 1; endo; tbopt[1] = ergeb[2,1]; ergeb = ZEROS(2,1); @ Vektor zum Speichern der Ergebnisse @ ttb2 = tbunter; do while ttb2 <= tbober; if ABS(ttb2 - tbopt[1]) < 2; ttb2 = tbopt[1] + 2; endif; {rho1,trho1,ttheta1,pp1,rres1,varres1,tstat1,kk1} = mio2break0b(yyy,tbopt[1],ttb2,kmax); if ABS(ttheta1) > ABS(ergeb[1,1]); ergeb[1,1] = ttheta1; ergeb[2,1] = ttb2; endif; ttb2 = ttb2 + 1; endo; tbopt[2] = ergeb[2,1]; tbopt = SORTC(tbopt,1); {rho1,trho1,ttheta1,pp1,rres1,varres1,tstat1,kk1} = mio2break0b(yyy,tbopt[1],tbopt[2],kmax); "****************"; "output"; "****************"; "model type M0"; "first break " tbopt[1]~yy[tbopt[1],1]~tbopt[1]/ROWS(yy); "second break " tbopt[2]~yy[tbopt[2],1]~tbopt[2]/ROWS(yy); "phi = rho-1 " rho1; "t value " trho1; "optimal lag " kk1; "variance " varres1; "------"; "coeff tstat "; "yyverz constant du1verz du2verz dtb1 dtb2"; pp1~tstat1; 我在运行过程中，系统提示Undefined symbols: mio2break0b。 这表示变量mio2break0b既没有被定义又没有被赋值。但是我不知道该怎么改，请教高手，谢谢。 |
回复此楼» 猜你喜欢
化工专硕求调剂
已经有3人回复
-
274求调剂
已经有5人回复
-
调剂
已经有3人回复
-
08工学调剂
已经有15人回复
-
招08考数学
已经有16人回复
-
340求调剂
已经有3人回复
-
材料专硕找调剂
已经有4人回复
-
0854电子信息求调剂
已经有6人回复
-
0805 316求调剂
已经有4人回复
-
0854 考研调剂 招生了！AI 方向
已经有17人回复
找到一些相关的精华帖子，希望有用哦~
小白求助怎么更改高斯计算的收敛条件
已经有6人回复
- 求助：请问各位高手在高斯09中优化单个分子时如何把一部分的参数固定了，优化其他部分
已经有12人回复
- 有个matlab程序，但是运行不下来，求高手
已经有23人回复
- Fortran编程过程中遇到错误，求大侠帮忙看下
已经有5人回复
- GaussView计算中断后，如何操作可以接着中断前的继续计算？
已经有16人回复
- 求高手解答--高斯软件计算化学反应机理
已经有6人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助
已经有5人回复
- 卟啉（100—150个原子）高斯计算的计算机配置问题
已经有14人回复
- 请问牛顿法和高斯牛顿法有什么不同？非常感谢
已经有6人回复
- 高斯irc计算出错，请高手帮忙看一下怎么解决
已经有8人回复
- gaussian 09 做好输入文件后总是关联不上，总说找不到文件，请高手指点！
已经有16人回复
- 帮忙讲解一下 文件的读取
已经有26人回复
- 用高斯计算，想买一台小型服务器自家用，请大家帮忙参考
已经有5人回复
- 【求助】高斯计算出现Fatal Problem: The smallest alpha delta epsilon is -0.785373
已经有7人回复
- 【求助】请高手分析下高斯2070错误，
已经有7人回复
- 【求助】帮忙改写展宽程序
已经有11人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
