#! /垃圾桶/桶#这是一个shell存档。去掉这行之前的任何东西,然后打开包装#将其保存到文件中并键入“sh文件”。覆盖现有#文件,键入“sh file-c”。您也可以通过以下方式将其作为标准输入#取消shar,或键入“sh'ch5eg1/detail.bak'<<'END_OF_FILE'X在t=142时未观测到。X在t=170时未观测到。X在t=172时观察不到。X在t=193时的观测缺失。XWarning,状态::to,var未初始化,tmp.dat不完整。X(X)X(X)X(X)**********************************************************************X轴**X*nlmdl 2.0版*X轴**X(X)**********************************************************************X(X)X(X)X(X)**********************************************************************X轴**X*参数设置*X轴**X(X)**********************************************************************X(X)X Gallant,“非线性统计模型”,第5章,图3,第324-326页X SUR什么估算方法?代码SUR、TSLS或GMM。X 224观测次数,t=1。。。,n.(名词)。X2方程的数量,即e的维数。X 0仪器数量,即Z的尺寸。X 8参数的数量,即θ的维数。X 50高斯-纽顿迭代的上限。X 1 Number var iterates,ivar=0表示无。X同型码同型或异型码。X 0用于var估计的移动平均项MA的数量。X none代码none或Parzen,当MA>0时为none是不明智的。X 1.000000e-13收敛公差,tol=1.0e-8是合理的。X 1.000000e-10反转公差,eps=1.0e-13是合理的X最小输出?无代码、最小代码或完整代码。X(X)X(X)X(X)**********************************************************************X轴**X*起始θ*X轴**X(X)**********************************************************************X(X)X(X)X第1列X(X)X行1-2.90000X行2-1.30000X行3 0.820000X行4 0.360000X行5-1.50000X行6-1.00000X行7-0.0300000X行8-0.470000X(X)X(X)X(X)**********************************************************************X轴**Xθ*X轴**X(X)**********************************************************************X(X)X(X)X第1列X(X)X行1-2.92458X行2-1.28675X行3 0.818570X行4 0.361158X行5-1.53759X第6排-1.004896X第7排0.0300867X第8排-0.467420X(X)X(X)X(X)X(X)**********************************************************************X轴**X*变量*X*(无自由度修正)*X轴**X(X)**********************************************************************X(X)X(X)X第1列第2列X(X)X行1 0.164925 0.0920057X行2 0.0920057 0.0896426X(X)X(X)X(X)X(X)**********************************************************************X轴**X*V(X*V)*X*(等级=8)*X轴**X(X)**********************************************************************X(X)X(X)X第1栏第2栏第3栏第4栏第5栏第6栏X(X)X行1 0.0760485-0.0394416-0.00813632 0.00225515 0.0151611-0.00086387X行2-0.0394416 0.0506098 0.00219396 0.00071347-0.00244522 0.0122632X行3-0.00813632 0.00219396 0.00644155 0.00078138-0.00603532-0.00238599X行4 0.00225515 0.00071347 0.00078138 0.00090344-0.00048087 0.00101291X行5 0.0151611-0.00244522-0.00603532-0.00048087 0.00832137 0.00200847X行6-0.00086387 0.0122632-0.00238599 0.00101291 0.00200847 0.00689444X行7-0.00293607 0.00321764-0.00075345-0.00079074 0.00110681 0.00015596X行8-0.00218205 0.00172663 0.00080830 0.00027549-0.00115374 0.00049917X(X)X(X)X第7列第8列X(X)X行1-0.00293607-0.00218205X行2 0.00321764 0.00172663X行3-0.00075345 0.00080830X行4-0.00079074 0.00027549X行5 0.00110681-0.00115374X行6 0.00015596 0.00049917X行7 0.00128616-0.00024276X行8-0.00024276 0.00036532X(X)文件结束如果测试5719-ne `wc-c<'ch5eg1/detail.bak'';然后echo-shar:\“'ch5eg1/detail.bak'\”解压缩时大小错误!fi(菲涅耳)#“ch5eg1/detail.bak”结束fi(菲涅耳)如果测试-f'ch5eg1/end3l.bak'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch5eg1/end3l.bak'\”其他的echo shar:提取“'ch5eg1/end3l.bak'”(5775个字符)sed“s/^X//”>'ch5eg1/end3l.bak'XGallant,“非线性统计模型”,第5章,图3,第324-326页XSUR什么估算方法?代码SUR、TSLS或GMM。X224观察次数,t=1。。。,n.(名词)。X2方程式数量,即e的尺寸。X0仪器数量,即Z的尺寸。X8参数的数量,即θ的维数。X50高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。异方差代码同方差或异方差。X3用于var估计的移动平均项MA的数量。XParzen Code none或Parzen,当MA>0时为none是不明智的。X1.000000e-13收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xfull输出多少?无代码、最小代码或完整代码。X8 V的计算秩。Xuncorrected var的除数,修正或未修正。X空行。X空行。X空行。X-2.9245814501616803e+00-2.92458145θ(1)X-1.2867460434063589e+00-1.28674604θ(2)X 8.1856990357452608e-01 0.81856990θ(3)X 3.6115785200404138e-01 0.36115785θ(4)X-1.5375885808695555e+00-1.53758858θ(5)X-1.0489590970263374e+00-1.04895910θ(6)X 3.0086700241162653e-02 0.03008670θ(7)X-4.6742011797582117e-01-0.46742012θ(8)X 1.6492462994112111e-01 0.16492463变量(1,1)X 9.2005729031321579e-02 0.09200573变量(2,1)X 9.2005729031321579e-02 0.09200573变量(1,2)X 8.964242316875812e-02 0.08964264变量(2,2)X 6.4151871786190293e-02 0.06415187伏(1,1)X-2.6808435093990508e-02-0.02680844伏(2.1)X-1.0481302083617307e-02-0.01048130伏(3.1)X 2.4179792658493493e-03 0.00241798伏(4.1)X 1.5271744257123141e-02 0.01527174伏(5.1)X 3.6217216495842341e-03 0.00362172伏(6.1)X-2.2745052058113957e-03-0.00227451伏(7,1)X-2.2282353784641093e-03-0.00222824伏(8,1)X-2.6808435093990942e-02-0.02680844伏(1,2)X 3.9204445770723149e-02 0.03920445伏(2,2)X 3.5239211150354388e-03 0.00352392伏(3.2)X 1.9044116642894958e-03 0.00190441伏(4.2)X-2.5090660931227499e-03-0.00250907伏(5.2)X 9.9860387081928022e-03 0.00998604伏(6.2)X 1.0763522857786497e-03 0.00107635伏(7.2)X 1.4608657346846976e-03 0.00146087伏(8.2)X-1.0481302083617045e-02-0.01048130伏(1,3)X 3.523921150353313e-03 0.00352392伏(2,3)X 5.9499339084749959e-03 0.00594993伏(3,3)X 3.8379282539160748e-04 0.00038379伏(4.3)X-5.5488697133901426e-03-0.00554887伏(5.3)X-2.3181283618853221e-03-0.00231813伏(6.3)X-3.0995472128008387e-04-0.00030995伏(7,3)X 7.9619528555448433e-04 0.00079620伏(8.3)X 2.4179792658493345e-03 0.00241798伏(1.4)X 1.9044116642895203e-03 0.00190441伏(2,4)X 3.8379282539158981e-04 0.00038379伏(3,4)X 1.1185997609060533e-03 0.00111860伏(4,4)X-8.2948235409772357e-05-0.00008295伏(5.4)X 1.8128096548391446e-03 0.00181281伏(6.4)X-8.6907934492451726e-04-0.00086908伏(7.4)X 2.3704685064746336e-04 0.00023705伏(8.4)X 1.5271744257122990e-02 0.01527174伏(1,5)X-2.5090660931226623e-03-0.00250907伏(2,5)X-5.5488697133901886e-03-0.00554887伏(3,5)X-8.2948235409779458e-05-0.00008295伏(4,5)X 7.2303542854823596e-03 0.00723035伏(5,5)X 2.4292056516279819e-03 0.00242921伏(6.5)X 6.6272949286796750e-04 0.00066273伏(7.5)X-1.0628753376875112e-03-0.00106288伏(8.5)X 3.6217216495840446e-03 0.00362172伏(1.6)X 9.9860387081929288e-03 0.00998604伏(2,6)X-2.3181283618853659e-03-0.00231813伏(3,6)X 1.8128096548391255e-03 0.00181281伏(4.6)X 2.4292056516280035e-03 0.00242921伏(5.6)X 6.9725212588786700e-03 0.00697252伏(6,6)X-8.0090995756335992e-04-0.00080091伏(7.6)X 2.8957023926462434e-04 0.00028957伏(8,6)X-2.274505205811465e-03-0.00227451伏(1,7)X 1.07635228577877e-03 0.00107635伏(2,7)X-3.0995472128007016e-04-0.00030995伏(3,7)X-8.6907934492452008e-04-0.00086908伏(4,7)X 6.6272949286796273e-04 0.00066273伏(5,7)X-8.0090995756337260e-04-0.00080091伏(6,7)X 1.1126181202496821e-03 0.00111262伏(7,7)X-2.579931664200757e-04-0.00025799伏(8,7)X-2.2282353784641279e-03-0.00222824伏(1.8)X 1.4608657346847011e-03 0.00146087伏(2,8)X 7.961952855549908e-04 0.00079620伏(3,8)X 2.37046850646418e-04 0.00023705伏(4,8)X-1.0628753376875236e-03-0.00106288伏(5,8)X 2.89570239264623e-04 0.00028957伏(6,8)X-2.579931664200887e-04-0.00025799伏(7,8)X 3.6558708385757817e-04 0.00036559伏(8,8)X 7.1212458944209244e-09 0.00000001 D(1)X-9.5291010557382897e-09-0.00000001 D(2)X-1.5778368693623276e-09-0.00000000 D(3)X-6.4433263552306315e-10-0.00000000 D(4)X 1.5722264246975941e-09 0.00000000 D(5)X-2.4312439036197547e-09-0.00000000 D(6)X-1.5845329475213192e-10-0.00000000 D(7)X-8.7707819833213637e-10-0.00000000 D(8)X 4.4685695288184939电子+02 446.85695288物镜文件结束如果测试5775-ne“wc-c<”ch5eg1/end3l.bak“”;然后echo shar:\“'ch5eg1/end3l.bak'”以错误的大小解包!fi(菲涅耳)#“ch5eg1/end3l.bak”的末尾fi(菲涅耳)如果测试-f'ch5eg1/ending.bak'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch5eg1/ending.bak'\”其他的echo shar:提取“'ch5eg1/ending.bak'”(5775个字符)sed“s/^X//”>'ch5eg1/ending.bak'<<'END_OF_FILE'XGallant,“非线性统计模型”,第5章,图3,第324-326页XSUR什么估算方法?代码SUR、TSLS或GMM。X224观察次数,t=1。。。,n.(名词)。X2方程式数量,即e的尺寸。X0仪器数量,即Z的尺寸。X8参数的数量,即θ的维数。X50高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。Xnone代码none或Parzen,当MA>0时为none是不明智的。X1.000000e-13收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xminimal输出多少?无代码、最小代码或完整代码。X8 V的计算秩。Xuncorrected var的除数,修正或未修正。X空行。X空行。X空行。X-2.9245814501616803e+00-2.92458145θ(1)X-1.2867460434063589e+00-1.28674604θ(2)X 8.1856990357452608e-01 0.81856990θ(3)X 3.6115785200404138e-01 0.36115785θ(4)X-1.5375885808695555e+00-1.53758858θ(5)X-1.0489590970263374e+00-1.0048995910θ(6)X 3.0086700241162653e-02 0.03008670θ(7)X-4.6742011797582117e-01-0.46742012θ(8)X 1.649246299411211e-01 0.16492463变量(1,1)X 9.2005729031321579e-02 0.09200573变量(2,1)X 9.2005729031321579e-02 0.09200573变量(1,2)X 8.964242316875812e-02 0.08964264变量(2,2)X 7.6048540298630313e-02 0.07604854伏(1,1)X-3.9441553460601381e-02-0.03944155伏(2.1)X-8.1363244658375945e-03-0.00813632伏(3.1)X 2.2551454088996347e-03 0.00225515伏(4.1)X 1.5161071040894071e-02 0.01516107伏(5.1)X-8.6386978244942293e-04-0.00086387伏(6.1)X-2.9360670383682128e-03-0.00293607伏(7.1)X-2.1820501983515709e-03-0.00218205伏(8,1)X-3.9441553460601381e-02-0.03944155伏(1,2)X 5.0609776285336930e-02 0.05060978伏(2,2)X 2.1939598380520749e-03 0.00219396伏(3.2)X 7.1346736728890531e-04 0.00071347伏(4.2)X-2.4452170498858623e-03-0.00244522伏(5.2)X 1.2263248643774908e-02 0.01226325伏(6.2)X 3.2176432662395357e-03 0.00321764伏(7.2)X 1.7266330377801262e-03 0.00172663伏(8.2)X-8.1363244658375945e-03-0.00813632伏(1,3)X 2.1939598380520749e-03 0.00219396伏(2,3)X 6.4415535695106418e-03 0.00644155伏(3,3)X 7.8138028913347761e-04 0.00078138伏(4.3)X-6.0353186597764865e-03-0.00603532伏(5.3)X-2.3859940684536310e-03-0.00238599伏(6.3)X-7.5345196703219372e-04-0.00075345伏(7,3)X 8.0829924747999706e-04 0.00080830伏(8.3)X 2.2551454088996347e-03 0.00225515伏(1.4)X 7.1346736728890531e-04 0.00071347伏(2,4)X 7.8138028913347761e-04 0.00078138伏(3,4)X 9.0343959846035443e-04 0.00090344伏(4,4)X-4.8087021200371650e-04-0.00048087伏(5.4)X 1.0129070210761102e-03 0.00101291伏(6.4)X-7.9073548924408663e-04-0.00079074伏(7.4)X 2.7549163836087631e-04 0.00027549伏(8.4伏)X 1.5161071040894071e-02 0.01516107伏(1,5)X-2.4452170498858623e-03-0.00244522伏(2,5)X-6.0353186597764865e-03-0.00603532伏(3,5)X-4.8087021200371650e-04-0.00048087伏(4,5)X 8.3213669176229327e-03 0.00832137伏(5,5)X 2.0084691585193232e-03 0.00200847伏(6.5)X 1.1068054871987483e-03 0.00110681伏(7.5)X-1.1537431897854821e-03-0.00115374伏(8.5)X-8.6386978244942293e-04-0.00086387伏(1.6)X 1.2263248643774908e-02 0.01226325伏(2,6)X-2.3859940684536310e-03-0.00238599伏(3,6)X 1.0129070210761102e-03 0.00101291伏(4,6)X 2.0084691585193232e-03 0.00200847伏(5.6)X 6.8944369024945473e-03 0.00689444伏(6,6)X 1.5596012364267730e-04 0.00015596伏(7.6)X 4.9917336963915457e-04 0.00049917伏(8.6)X-2.9360670383682128e-03-0.00293607伏(1,7)X 3.2176432662395357e-03 0.00321764伏(2,7)X-7.5345196703219372e-04-0.00075345伏(3,7)X-7.9073548924408663e-04-0.00079074伏(4,7)X 1.1068054871987483e-03 0.00110681伏(5,7)X 1.5596012364267730e-04 0.00015596伏(6,7)X 1.2861607343363046e-03 0.00128616伏(7,7)X-2.4275726609718795e-04-0.00024276伏(8,7)X-2.1820501983515709e-03-0.00218205伏(1,8)X 1.7266330377801262e-03 0.00172663伏(2,8)X 8.0829924747999706e-04 0.00080830伏(3,8)X 2.7549163836087631e-04 0.00027549伏(4,8)X-1.1537431897854821e-03-0.00115374伏(5,8)X 4.9917336963915457e-04 0.00049917伏(6,8)X-2.4275726609718795e-04-0.00024276伏(7,8)X 3.6531994032499421e-04 0.00036532伏(8,8)X 7.1212458944209244e-09 0.00000001 D(1)X-9.5291010557382897e-09-0.00000001 D(2)X-1.5778368693623276e-09-0.00000000 D(3)X-6.4433263552306315e-10-0.00000000 D(4)X 1.5722264246975941e-09 0.00000000 D(5)X-2.4312439036197547e-09-0.00000000 D(6)X-1.5845329475213192e-10-0.00000000 D(7)X-8.7707819833213637e-10-0.00000000 D(8)X 4.4685695288184939e+02 446.85695288对象文件结束如果测试5775-ne `wc-c<'ch5eg1/ending.bak'';然后echo shar:\“'ch5eg1/ending.bak'”以错误的大小解包!fi(菲涅耳)#“ch5eg1/ending.bak”的末尾fi(菲涅耳)如果测试-f'ch6eg1/ending.bak'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg1/ending.bak'\”其他的echo shar:提取“'ch6eg1/ending.bak'”\(5772个字符\)sed“s/^X//”>'ch6eg1/ending.bak'<<'END_OF_FILE'XGallant,“非线性统计模型”,第6章,图1,第436页。XTSLS什么估算方法?代码SUR、TSLS或GMM。X220观察次数,t=1。。。,n.(名词)。X2方程式数量,即e的尺寸。X17仪器数量,即Z的尺寸。X8参数的数量,即θ的维数。X20高斯-纽顿迭代的上限。X1数值var迭代,ivar=0表示无。同型或异型编码。X0用于var估计的移动平均项MA的数量。Xnone代码none或Parzen,当MA>0时为none是不明智的。X1.000000e-08收敛公差,tol=1.0e-8是合理的。X1.000000e-10反转公差,eps=1.0e-13是合理的Xfull输出多少?无代码、最小代码或完整代码。X8 V的计算秩。Xuncorrected var的除数,修正或未修正。X空行。X空行。X空行。X-2.1378778275860273e+00-2.13787783θ(1)X-1.9893911556433312e+00-1.98939116θ(2)X 7.0939153916745556e-01 0.70939154θ(3)X 3.36633999323270310e-01 0.33663399θ(4)X-1.40019968694481941e+00-1.40199687θ(5)X-1.1388969627624823e+00-1.13889696θ(6)X 2.9131488559477250e-02 0.02913149θ(7)X-5.0050229673416280e-01-0.50050230θ(8)X 1.6847307354905444e-01 0.16847307变量(1,1)X 9.4991006557633997e-02 0.09499101变量(2,1)X 9.4991006557633997e-02 0.09499101变量(1,2)X 9.3711564655296287e-02 0.09371156变量(2,2)X 3.4123357521137065e-01 0.34123358伏(1,1)X-4.0515026192173231e-01-0.40515026伏(2.1)X-7.4443670981270527e-02-0.07444367伏(3.1)X-1.7075137202523048e-02-0.01707514伏(4,1)X 7.3616408974018430e-02 0.07361641伏(5.1)X-8.0520944995638929e-02-0.08052094伏(6.1)X-1.3988025621517654e-02-0.01398803伏(7.1)X-2.2077927314748533e-02-0.02207793伏(8,1)X-4.0515026192173231e-01-0.40515026伏(1,2)X 5.6591782097119825e-01 0.56591782伏(2,2)X 9.3113487046268872e-02 0.09311349伏(3.2)X 2.8410019934882371e-02 0.02841002伏(4.2)X-8.0617530766411080e-02-0.08061753伏(5.2)X 1.2644653136927941e-01 0.12644653伏(6.2)X 1.9325013502098160e-02 0.01932501伏(7.2)X 2.9896132263088925e-02 0.02989613伏(8,2)X-7.4443670981270527e-02-0.07444367伏(1,3)X 9.3113487046268872e-02 0.09311349伏(2,3)X 2.4067454001684946e-02 0.02406745伏(3.3)X 5.7835699332500186e-03 0.00578357伏(4.3)X-2.1581815564116969e-02-0.02158182伏(5.3)X 1.7403224408106693e-02 0.01740322伏(6.3)X 1.9150094956931823e-03 0.00191501伏(7,3)X 5.9076727860237991e-03 0.00590767伏(8.3)X-1.7075137202523048e-02-0.01707514伏(1,4)X 2.8410019934882371e-02 0.02841002伏(2,4)X 5.7835699332500186e-03 0.00578357伏(3,4)X 2.5487432759845568e-03 0.00254874伏(4,4)X-4.7950486640943122e-03-0.00479505伏(5.4伏)X 7.3439708547267122e-03 0.00734397伏(6,4)X-5.1094810414861448e-05-0.00005109伏(7,4)X 1.8674253117946520e-03 0.00186743伏(8.4)X 7.3616408974018430e-02 0.07361641伏(1.5)X-8.0617530766411080e-02-0.08061753伏(2,5)X-2.1581815564116969e-02-0.02158182伏(3,5)X-4.7950486640943122e-03-0.00479505伏(4,5)X 2.2762057613803403e-02 0.02276206伏(5,5)X-1.5195066143983841e-02-0.01519507伏(6.5)X-9.9415605596577862e-04-0.00099416伏(7,5)X-5.6251801394123823e-03-0.00562518伏(8.5)X-8.0520944995638929e-02-0.08052094伏(1,6)X 1.2644653136927941e-01 0.12644653伏(2,6)X 1.7403224408106693e-02 0.01740322伏(3.6)X 7.3439708547267122e-03 0.00734397伏(4,6)X-1.5195066143983841e-02-0.01519507伏(5.6伏)X 3.3346276344180567e-02 0.03334628伏(6,6)X 3.7331467729901690e-03 0.00373315伏(7.6)X 6.9269333322694495e-03 0.00692693伏(8,6)X-1.3988025621517654e-02-0.01398803伏(1.7)X 1.9325013502098160e-02 0.01932501伏(2,7)X 1.9150094956931823e-03 0.00191501伏(3,7)X-5.1094810414861448e-05-0.00005109伏(4.7伏)X-9.9415605596577862e-04-0.00099416伏(5,7)X 3.7331467729901690e-03 0.00373315伏(6,7)X 2.0417788436814010e-03 0.00204178伏(7,7)X 6.0905830809060847e-04 0.00060906伏(8,7)X-2.2077927314748533e-02-0.02207793伏(1.8)X 2.9896132263088925e-02 0.02989613伏(2,8)X 5.9076727860237991e-03 0.00590767伏(3,8)X 1.8674253117946520e-03 0.00186743伏(4,8)X-5.6251801394123823e-03-0.00562518伏(5,8)X 6.9269333322694495e-03 0.00692693伏(6,8)X 6.0905830809060847e-04 0.00060906伏(7,8)X 2.0036106220501352e-03 0.00200361伏(8,8)X 1.6787206308357473e-07 0.00000017 D(1)X-2.4625529712808257e-07-0.00000025 D(2)X-4.6504414197150224e-08-0.00000005 D(3)X-1.4610224422905096e-08-0.00000001 D(4)X 4.0400010629413197e-08 0.00000004 D(5)X-5.8821932254927893e-08-0.00000006 D(6)X-8.2632941559150507e-09-0.00000001 D(7)X-1.6640001385280539e-08-0.00000002 D(8)X 3.5611514196184793e+01 35.61151420对象文件结束如果测试5772-ne `wc-c<'ch6eg1/ending.bak'`;然后echo shar:\“'ch6eg1/ending.bak'”以错误的大小解包!fi(菲涅耳)#“ch6eg1/ending.bak”的末尾fi(菲涅耳)如果测试-f'ch6eg2/detail3l.bak'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“'ch6eg2/detail3l.bak'\”其他的echo shar:提取“'ch6eg2/detail3l.bak'”\(7043个字符\)sed“s/^X//”>'ch6eg2/detail3l.bak'X(X)X(X)X(X)**********************************************************************X轴**X*nlmdl 2.0版*X轴**X(X)**********************************************************************X(X)X(X)X(X)**********************************************************************X轴**X*参数设置*X轴**X(X)**********************************************************************X(X)X Gallant,“非线性统计模型”,第6章,图2,第449-450页。X GMM什么估算方法?代码SUR、TSLS或GMM。X 238观察次数,t=1。。。,n.(名词)。X 1方程式数量,即e的尺寸。X 3仪器数量,即Z的尺寸。X2参数的数量,即θ的维数。X 20高斯-纽顿迭代的上限。X 1 Number var iterates,ivar=0表示无。X异方差代码同方差或异方差。X 3 var估计的移动平均项MA的数量。当MA>0时,X Parzen Code none或Parzen,none是不明智的。X 1.00000e-05收敛公差,tol=1.0e-8是合理的。X 1.000000e-13反转公差,eps=1.0e-13是合理的X满输出多少?无代码、最小代码或完整代码。X(X)X(X)X(X)**********************************************************************X轴**X*起始θ*X轴**X(X)**********************************************************************X(X)X(X)X第1列X(X)X行1-0.400000X行2 0.900000X(X)Xvar_循环0X 2.38000000000000e+02 23800000000 var(1,1)X 2.3848513503508150e+02 238.48513504变量(2,1)X 2.3866807009536583e+02 238.66807010 var(3,1)X 2.3848513503508150e+02 238.48513504变量(1,2)X 2.3897620940526255e+02 238.97620941 var(2,2)X 2.3916481122626234e+02 239.16481123变量(3,2)X 2.3866807009536583e+02 238.66807010 var(1,3)X 2.3916481122626234e+02 239.16481123变量(2,3)X 2.3975964092305995e+02 239.75964092 var(3,3)X(X)Xθ_回路0X-4.000000000000002e-01-0.40000000θ(1)X 9.000000000000002e-01 0.90000000θ(2)X 2.2977593855825589e+00 2.2977593目标X-4.9910043429446205e-01-0.49910043 D(1)X 9.8933364963875864e-02 0.09893336 D(2)X步长=1X(X)Xθ_回路1X-8.9910043429446207e-01-0.89910043θ(1)X 9.98933336496387586e-01 0.99893336θ(2)X 2.7751760221041621e-03 0.00277518对象X 5.0366043831498272e-02 0.05036604 D(1)X-4.4414610501568677e-06-0.00000444 D(2)X步长=1X(X)Xθ_回路2X-8.4873439046296384e-01-0.84873439θ(1)X 9.9892892350282569e-01 0.99892892θ(2)X 2.7717868638761308e-03 0.00277179对象X-1.1734797022229638e-04-0.00011735 D(1)X 2.1096605255709332e-07 0.00000021 D(2)X步长=1X(X)Xθ_回路3X-8.4885173843318618e-01-0.84885174θ(1)X 9.9892913446887821e-01 0.99892913θ(2)X 2.7717868580610532e-03 0.00277179对象X 2.6798358348564710e-07 0.00000027 D(1)X-5.5350439549602871e-10-00000000 D(2)XTolence检查通过。X(X)Xvar_回路1X 4.397820854904504e-01 0.43978221变量(1,1)X 4.4050039895600307e-01 0.44050040变量(2,1)X 4.3420956263910121e-01 0.43420956 var(3.1)X 4.4050039895600307e-01 0.44050040变量(1,2)X 4.4122700259798669e-01 0.44122700 var(2,2)X 4.349363977747488e-01 0.43493664 var(3,2)X 4.3420956263910121e-01 0.43420956 var(1,3)X 4.349363977747488e-01 0.43493664变量(2,3)X 4.2967081964723475e-01 0.42967082 var(3,3)X(X)Xθ_回路0X-8.4885173843318618e-01-0.84885174θ(1)X 9.9892913446887821e-01 0.99892913θ(2)X 1.2014636301992196e+00 1.20146633对象X-6.4020486065461402e-03-0.00640205 D(1)X-8.8353145570584189e-04-0.00088353 D(2)X步长=1X(X)Xθ_回路1X-8.5525378703973232e-01-0.85525379θ(1)X 9.980456031317236e-01 0.99804560θ(2)X 1.0861341198092420e+00 1.08613412对象X 2.6747016119385693e-06 0.00000267 D(1)X-1.7073275785211286e-08-0.00000002 D(2)XTolence检查通过。X(X)X(X)X**********************************************************************X**X*θ*X轴**X(X)**********************************************************************X(X)X(X)X第1列X(X)X行1-0.855254X行2 0.998046X(X)X(X)X(X)X(X)**********************************************************************X轴**X*变量*X*(无自由度修正)*X轴**X(X)**********************************************************************X(X)X(X)X第1列第2列第3列X(X)X行1 0.439782 0.440500 0.434210X行2 0.440500 0.441227 0.434937X行3 0.434210 0.434937 0.429671X(X)X(X)X(X)X(X)**********************************************************************X轴**X*V(X*V)*X*(秩=2)*X轴**X(X)**********************************************************************X(X)X(X)X第1列第2列X(X)X行1 3.60663-0.00715768X行2-0.00715768 2.117e-05X(X)文件结束如果测试7043-ne“wc-c<‘ch6eg2/detail3l.bak’”;然后echo shar:\“'ch6eg2/detail3l.bak'”打开时大小错误!fi(菲涅耳)#“ch6eg2/detail3l.bak”结束fi(菲涅耳)如果测试-f“display.cc”-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“display.cc”其他的echo shar:提取“display.cc”(4567个字符)sed“s/^X//”>'显示.cc'X/*----------------------------------------------------------------------------X(X)Xnlmdl:显示.ccX(X)Xnlmdl是a.Ronald中统计方法的C++实现XGallant,“非线性统计模型”,纽约:John Wiley and Sons,X1987,ISBN 0-471-80260-3,使用分布式的矩阵类realmat头文件nlmdl.h和realmat.h描述了Xprogram和矩阵类。X(X)X版权所有(C)1990。X(X)XA公司。罗纳德·格兰特XP。邮政信箱5513XRaleigh北卡罗来纳州27650-5513XUSA公司X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且X版权声明和本许可声明出现在支持X文档。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X----------------------------------------------------------------------------*/X/*display.cc是一段被分割的状态.cc*/X/*输出以减少Turbo C++内存需求*/X(X)X#包含“status.h”X(X)Xint状态::display(display_mode模式)X(X){X(X)X整数i,j,len;X字符温度[MAX_STATUS_LINE];X char*pad=“”;X(X)X(X)X开关(模式){X(X)X箱启动:X(X)X cout<<starbox(“/Parameter settings//_”)<<“\n”;X(X)X cout公司<0是不明智的。\n“;X(X)X sprintf(温度,“%-19.6e”,tol);X cout公司<0列(var)>0){X if(strcmp(method,“SUR”)==0||strcmp{X代表(j=1;j<=M;j++){X表示(i=1;i<=M;i++){X sprintf(温度,“%25.16e%20.8f”,变量检查2(i,j),变量长度(i,j));X高级定制<0&&列(θ)>0){X代表(i=1;i<=p;i++){X冲刺(温度,“%25.16e%20.8f”,θ[i],θ[i]);X cout公司<0列(D)>0){X代表(i=1;i<=p;i++){X冲刺(温度,“%25.16e%20.8f”,D[i],D[i);X cout公司<“nlmdl.h”<<“END_OF_FILE”X/*---------------------------------------------------------------------------X(X)Xnlmdl:nlmdl.hX(X)Xnlmdl是a.Ronald中统计方法的C++实现XGallant,“非线性统计模型”,纽约:John Wiley and Sons,X1987,ISBN 0-471-80260-3,使用分布式的矩阵类realmat头文件nlmdl.h和realmat.h描述了Xprogram和矩阵类。X(X)X版权所有(C)1990。X(X)XA公司。罗纳德·格兰特XP。邮政信箱5513XRaleigh北卡罗来纳州27650-5513XUSA公司X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且X版权声明和本许可声明出现在支持X文档。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X-----------------------------------------------------------------------------X(X)XT此标题描述用于估算模型θ的程序nlmdl.ccXe=q(t,θ)。用户提供一个名为starting.dat的文件,其内容为模板描述的X:X(X)X开关此行作为字符串传递给类模型。X方法什么估算方法?代码为SUR、TSLS或GMM。X n观测次数,t=1。。。,n.(名词)。X M系统中方程的数量,即e的维数。X K仪器数量,即Z尺寸,代码K=0,带SUR。X p参数的数量,即θ的维数。X itheta Gauss-Newton迭代的上限。X ivar var迭代次数,ivar=0表示无。X vartype代码同型或异型。X MA var估计的移动平均项数,如果e iid,代码MA=0。X权重“Code none”或“Parzen”,当MA>0时,“Code none”是自找麻烦。X tol收敛公差,tol=1.0e-8是合理的。X eps反转容差,eps=1.0e-13是合理的(如果REAL是双精度的)。X细节输出多少?无代码、最小代码或完整代码。X空行。X空行。X空行。X空行。X空行。Xθ(1)θ的起始值。必须提供这些。X(X)X至X(X)Xθ(p)X var(1)方差估计的起始值。这些可以省略。X矩阵按列存储。如果方法=SUR或TSLS,则X到l=M,var对应于C(e,e')。如果方法=GMM,则l=M*KX和var对应于e Kronecker t上总和的方差X var(l*l)乘积Z。X(X)X(X)XSee类状态以了解更多详细信息,尤其是状态::from。文件Xstarting.dat被读取两次。第一次阅读时,它必须至少包含X第一条线,交换机。接下来描述的类模型可以提供条目X表示二读。X(X)X用户还将类模型作为文件model.h和model.cc提供Xclass模型是在model.h中声明的,它应该与以下模板匹配:X(X)X#包含“status.h”X外部状态;X(X)X类模型X{X受保护:X//。。。X(X)X公众:X模型();X~模型();X realmat e_t;//e_t是s.M乘以1X realmat dele_t;//dele_t是s.M乘以s.pX realmat Z_t;//Z_t是s.K乘以1,如果s.K=0,则为空;X int initialize();X int terminate();X无效e(整数t);//e填写e_tX void dele(整数t);//dele填充delet=del e_t wrt thetaX无效Z(整数t);//Z填充Z_tX}(X)X(X)XF对于SUR,按如下方式对函数Z进行编码:X(X)X无效模型::Z(整数t){}X(X)XI在编写类模型时,有关程序流的相关事实Xnlmdl.cc如下:X(X)X 1。首先,s暂定为starting.dat,使用X(X)X状态;X s.from(s.starting);X(X)X此读取的主要目的是从X starting.dat,以便可以将其建模为字符串X s.交换机。尽管如此,starting.dat中的所有内容都被读取并X输入s。X(X)X 2。接下来,调用类模型的initialize()。函数initialize()X可以读取starting.dat,可以读取其他数据等X状态的数据是公开的,任何数据都可以由填写或更改X initialize()。例如,可以填写所有类状态X数据并使用写入X(X)X s.to(s.starting);X(X)X One还可以使用切换到不同的启动状态文件X(X)X s.starting=“文件名”;X(X)X或其他结束状态文件使用X(X)X s.ending=“文件名”;X(X)如果initialize()返回0,则执行停止。X(X)X 3。然后,类状态的数据由调用最终填充X(X)X s.from(s.starting);X(X)X和m.e_t、m.dele_t和m.Z_t的尺寸为X(X)X m.e_t.调整大小(s.m,1);X m.dele_t.调整大小(s.m,s.p);X,如果(s.K>0)m.Z_t.调整大小(s.K,1);X(X)X使用这些设置和维度进行估算。X(X)X 4。最后,调用类模型的terminate();terminate()可以读取数据,X写入数据等。如果terminate()返回0,则执行停止,否则再次调用X initialize(),循环重复(从步骤2开始)。一个X可以使用此功能在起始值网格上循环。X(X)X(X)X参考:X(X)XGallant,A.Ronald(1987),“非线性统计模型”,纽约:JohnXWiley父子公司。国际标准图书编号0-471-80260-3。X(X)X----------------------------------------------------------------------------*/X(X)X(X)X#ifndef __FILE_NLMDL_H_SEEN__X#杂注一次X#定义__FILE_NLMDL_H_SEEN__ 1X(X)X#include“normal.h”//此标头也处于状态.hX#include“realmat.h”//此标头也处于状态.hX#include“tools.h”//此标题也处于状态.hX(X)X#包含“status.h”X#包含“model.h”X(X)X#define版本“2.1”X(X)X#结尾X(X)文件结束如果测试5850-ne“wc-c<'nlmdl.h'”;然后echo shar:“'nlmdl.h'”打开时大小错误!fi(菲涅耳)#“nlmdl.h”的结尾fi(菲涅耳)如果测试-f'realmat.h'-a“${1}”!=“-c”;然后echo shar:不会破坏现有文件“realmat.h”其他的echo shar:提取“realmat.h”(5870个字符)sed“s/^X//”>'真实.h'<<'END_OF_FILE'X/*----------------------------------------------------------------------------X(X)Xrealmat:realmat.hX(X)Xrealmat是一个C++矩阵类。头文件realmat.h描述了它的用法。X(X)X版权所有(C)1990。X(X)XA公司。罗纳德·格兰特XP。邮政信箱5513XRaleigh北卡罗来纳州27650-5513XUSA公司X(X)X使用、复制、修改和分发本软件及其特此免费授予用于任何目的的Xdocumentation,前提是X上述版权声明出现在所有副本中,并且X版权声明和本许可声明出现在支持X文档。X(X)XT本软件按“原样”提供,无任何明示或暗示担保。X(X)X--------------------------------------------------------------------------------X(X)XT此标头定义了一个矩阵类。具有r行和c列的矩阵X是使用实数X(r,c)构造。X按列存储,没有浪费空间也就是说vec(X)是存储的。元素索引开始X为1,而不是零。realmat的受保护成员包括行、列、len和x。X(X)XA是一个示例,假设X是使用realmat X(r,c)构造的。然后二十、 rows=r,X.cols=c,X.len=r*c,指针X.X指向Xvec(X),X[i]返回vec(X)的第i个元素,i=1,。。。,X.len,X.elem(i,j)X返回X的第(i,j)个元素,i=1,。。。,X行,j=1,。。。,。X.系列。一个元素X的X为REAL类型,X.row、X.cols和X.len为INTEGER类型。这个REAL和INTEGER的Xtypedef是in-usual.h。要在运行时分配存储,请使用Xrealmat X;X尺寸(r,c);。X(X)X运算符+,-,*,=,+=,-=,++,--,==,!=,<,<=,>,>=和<<是X定义。T(X)转置X.X.check1(i)和X.check 2(i,j)访问元素vec(X)的X和X分别进行范围检查。invpsd(X)或invpsd(X,eps)X返回半正定对称矩阵X的g-逆X最小错误检查;eps是等级确定中使用的公差X默认为1.e-13。rows(X)返回X.row,cols(X)则返回X.cols。这个Xoperator++的行为类似于++X,而不是X++;类似于--。运算符==,如果比较的矩阵大小相同,则X<、<=、>和>=的计算结果为真X并且每个元素的关系为真;(A!=B)=!(A==B)。X(X)X-----------------------------------------------------------------------------*/X(X)X#ifndef __FILE_REALMAT_H_SEEN(如果文件名为__FILE_REALMAT_H_SEEN)__X#杂注一次X#定义__FILE_REALMAT_H_SEEN__1X(X)X#包括“usual.h”X#包括“tools.h”X(X)Xclass realmatX(X){X保护:X(X)X整数行;X整数列;X整数透镜;X真实*X;X(X)X实数(整数r、整数c、实数*a);X void resize(整数r,整数c,实*a)//应为。x==0X(X)X(X)Xpublic公司:X realmat();X(X)X realmat(整数r,整数c);X实数(INTEGER r,INTEGER c,REAL fill_value);X(X)X realmat(realmat&);X(X)X~实数();X(X)Xvoid大小调整(整数r,整数c);Xvoid大小调整(INTEGER r,INTEGER c,REAL fill_value);X(X)Xrealmat&operator=(realmat/a);Xrealmat&operator+=(realmat/a);Xrealmat&operator-=(realmat/a);Xrealmat&operator++();Xrealmat和操作员--();X(X)XREAL&operator[](整数i);XREAL&elem(整数i,整数j);XREAL和check1(整数i);XREAL&check2(整数i,整数j);X(X)Xvoid错误(const char*msg);X(X)Xfriend int运算符==(realmat&a,realmat/b);Xfriend int运算符=(真实材料&a、真实材料&b);X(X)Xfriend int运算符<(realmat&a,realmat/b);Xfriend int运算符>(realmat&a,realmat&b);X(X)Xfriend int运算符<=(realmat&a,realmat&b);Xfriend int运算符>=(realmat&a,realmat&b);X(X)Xfriend INTEGER行(realmat&a);Xfriend INTEGER cols(realmat&a);X(X)Xfriend ostream&operator<<(ostream&stream,realmat&a);X(X)Xfriend realmat operator+(realmat&a,realmat/b);Xfriend realmat运营商+(realmat&a);X(X)Xfriend realmat运营商-(realmat&a,realmat/b);Xfriend realmat运营商-(realmat&a);X(X)Xfriend realmat运营商*(realmat&a、realmat&b);Xfriend realmat运营商*(REAL&a,realmat&b);Xfriend realmat运营商*(INTEGER&a,realmat&b);X(X)Xfriend realmat T(realmat&a);X(X)Xfriend realmat invpsd(realmat&a,REAL eps=1.0e-13)//a是psd,对称X(X)十} ;X(X)X(X)X//受保护:X(X)Xinline realmat::realmat(INTEGER r,INTEGER c,REAL*a)X{行=r;列=c;长度=r*c;X=a;}X(X)Xinline void realmat::resize(INTEGER r,INTEGER c,REAL*a)X{rows=r;cols=c;len=r*c;X=a;}//注意:没有删除,应该是this.X==0。X(X)X(X)X//公开:X(X)Xinline int运算符=(真实材料&a,真实材料&b)X{返回!(a==b);}X(X)Xinline REAL和realmat::operator[](INTEGER i)X{返回X[i-1];}X(X)Xinline REAL和realmat::elem(INTEGER i,INTEGER j)X{return X[i+行*j-行-1];}//返回X[行*(j-1)+i-1]X(X)Xinline INTEGER行(realmat&a)X{返回a.rows;}X(X)Xinline INTEGER cols(realmat&a)X{返回a.cols;}X(X)Xinline realmat::~realmat()X{删除X;}X(X)Xinline realmat::realmat()X{行=0;列=0;长度=0;X=0;}X(X)Xinline realmat::realmat(整数r,整数c)X{X if(r<=0)error(“error,realmat,realmat,Number of rows not positive”(错误,实数,实数,非正数行数));X if(c<=0)error(“error,realmat,realmath,Number of columns not positive”);X行=r;X列=c;X len=r*c;X X=新REAL[len];X,如果(X==0)错误(“错误,realmat,realmet,操作员新建失败”);X}(X)X(X)Xinline void realmat::resize(INTEGER r,INTEGER c)X(X){X如果(r<=0)错误(“错误,realmat,调整大小,行数非正数”);X if(c<=0)error(“error,realmat,resize,Number of columns not positive”);X删除X;X行=r;X列=c;X len=r*c;X X=新REAL[len];X,如果(X==0)错误(“错误,realmat,调整大小,操作员新建失败”);X}(X)X(X)Xextern void default_realmat_error_handler(const char*msg);Xextern ONE_ARG_ERROR_HANDLER_T realmat_ERROR_HANDLER;Xextern ONE_ARG_ERROR_HANDLER_T公司X set_realmat_error_handler(ONE_ARG_error_handler_T f);X(X)X#结尾文件结束如果测试5870-ne‘wc-c<‘realmat.h’’;然后echo shar:“realmat.h”的大小不正确!fi(菲涅耳)#“realmat.h”结束fi(菲涅耳)echo shar:存档2结束(共6个)。cp/dev/null ark2已完成缺少=“”对于1 2 3 4 5 6中的I;如果测试-方舟${一} 完成了 ; 然后MISSING=“${MISSING}${I}”fi(菲涅耳)完成如果测试“${MISSING}”=“”;然后echo您已打开所有6个档案。rm-f方格[1-9]完成其他的echo您仍然需要解压缩以下存档文件:echo“”${MISSING}fi(菲涅耳)##shell存档结束。退出0