1stopt非线性拟合时多次运行结果相差太大的疑问
本人采用1stopt进行参数拟合,代码如下Parameters E1[1,50000],eta1[1,1000],E2[1,10000],eta2[1,2000];
Variable x,y(6:12);
//t 6 7 8 9 10 11 12
function y6=1/E1+x/eta1+(1-exp(-E2/eta2*x))/E2 ;
Data;
0.047 1.53E-04 1.53E-04 1.47E-04 1.53E-04 1.54E-04 1.49E-04 1.52E-04
0.097 2.66E-04 2.66E-04 2.54E-04 2.66E-04 2.67E-04 2.59E-04 2.66E-04
0.149 3.72E-04 3.71E-04 3.58E-04 3.74E-04 3.75E-04 3.62E-04 3.72E-04
0.203 4.77E-04 4.74E-04 4.59E-04 4.79E-04 4.79E-04 4.63E-04 4.77E-04
0.26 5.83E-04 5.77E-04 5.60E-04 5.86E-04 5.83E-04 5.65E-04 5.83E-04
0.32 6.90E-04 6.81E-04 6.62E-04 6.94E-04 6.87E-04 6.69E-04 6.90E-04
0.382 8.00E-04 7.86E-04 7.69E-04 8.05E-04 7.94E-04 7.75E-04 7.99E-04
0.447 9.12E-04 8.93E-04 8.78E-04 9.18E-04 9.02E-04 8.84E-04 9.11E-04
0.516 1.03E-03 1.00E-03 9.91E-04 1.03E-03 1.01E-03 9.97E-04 1.03E-03
0.587 1.14E-03 1.12E-03 1.11E-03 1.27E-03 1.13E-03 1.11E-03 1.14E-03
0.662 1.26E-03 1.23E-03 1.23E-03 1.39E-03 1.24E-03 1.24E-03 1.26E-03
0.741 1.39E-03 1.35E-03 1.36E-03 1.52E-03 1.37E-03 1.36E-03 1.39E-03
0.823 1.52E-03 1.48E-03 1.49E-03 1.64E-03 1.49E-03 1.49E-03 1.52E-03
0.91 1.65E-03 1.61E-03 1.63E-03 1.78E-03 1.62E-03 1.63E-03 1.65E-03
1 1.79E-03 1.74E-03 1.77E-03 1.91E-03 1.76E-03 1.77E-03 1.78E-03
在采用同组数据进行拟和时,多次运行得到的结果有时差别比较大,但是,其相关系数、残差平方和等都表明回归具有很高的质量,这种情况下怎么判断哪一个结果更合理
另外采用不同的回归方法得到的参数结果相差就更大,如以上代码采用简面体爬山法与麦夸特法得到的回归结果差别更大
想请教如何选择合适的回归方法及对回归结果的合理性进行识别?
哪位有关于1stopt中几种回归方法介绍的资料啊,想看一看,希望热心人提供共享 本人采用同一种方法回归曾得到一下两个结果
E1 21817.99942
eta1 657.4972261
E2 4359.326634
eta2 1166.419372
残差平方和: 6.56921262993701E-11
相关系数(R^2): 0.99998241999163
E1 12750.91385
eta1 190.709309
E2 954.2788789
eta2 388.8222842
残差平方和: 7.25872833499016E-11
相关系数(R^2): 0.999981338901178
但是可以看到二者差别很大,这时候怎么确定哪一个更合理,很令人头疼 用差分或爬山法可得较稳定得结果:
优化算法: 标准差分进化算法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 2.09248811116761E-6
残差平方和(SSE): 6.56775974306669E-11
相关系数(R): 0.999991211713161
相关系数之平方(R^2): 0.999982423503557
决定系数(DC): 0.999982423503557
卡方系数(Chi-Square): 5.15594463345163E-8
F统计(F-Statistic): 208611.91981562
参数 最佳估算
---------- -------------
e1 21216.5578502987
eta1 657.299336720303
e2 4386.64987777283
eta2 1182.8148458846 本人经过对方程y6进行多次回归,请楼上的看看这个结果:
均方差(RMSE): 2.13976694540969E-6
残差平方和(SSE): 6.86790387100187E-11
相关系数((R): 0.999991177676713
相关系数之平方(R^2): 0.999982355431258
决定系数(DC): 0.999981620264339
卡方系数(Chi-Square): 5.89339541530299E-8
F统计(F-Statistic): 736757.62842478
参数 最佳估算
---------- -------------
E1 21188.6705988168
eta1 719.437129927333
E2 1322.65092706657
eta2 1355.87514879997
这个结果与楼上的结果差别就比较大。
楼上的有没有经历过回归的效果基本相同的情况下,得到的回归参数差别较大,
我关注的是这种情况下怎么判断哪个更合理 我得到得参数结果至少整数部分都是一样得,小数部分有微小差别。 [quote]原帖由 [i]shamohu[/i] 于 2007-1-10 17:01 发表
我得到得参数结果至少整数部分都是一样得,小数部分有微小差别。 [/quote]
请问shamohu你用的是哪个版本啊,我用的是1.0的
页:
[1]
