帮忙算一组数据

handb

我的软件是试用版的，数据较多，不能粘贴，请帮帮忙，谢谢！
Parameters a,b,c,d,e,f,g,h,i,j;
Variable x,y,z,C;
Function C=a+b*x+c*y+d*z+e*x*y+f*x*z+g*y*z+h*x^2+i*y^2+j*z^2;
Data;
  x   y   z     C
0.5        10        20        0.010902
0.5        15        20        0.0033786
0.5        20        20        0.0025
0.5        25        20        0.0012
0.5        30        20        0.00050257
1        10        20        0.0087
1        15        20        0.0036
1        20        20        0.0017
1        25        20        0.00068731
1        30        20        0.00032019
1.5        10        20        0.0071
1.5        15        20        0.003
1.5        20        20        0.0012
1.5        25        20        0.00082019
1.5        30        20        0.0003993
2        10        20        0.0085
2        15        20        0.0029
2        20        20        0.0013
2        25        20        0.00092489
2        30        20        0.00058607
2.5        10        20        0.0061
2.5        15        20        0.0026
2.5        20        20        0.0013
2.5        25        20        0.00058824
2.5        30        20        0.00036031
3        10        20        0.0052
3        15        20        0.0029
3        20        20        0.0018
3        25        20        0.00090369
3        30        20        0.00051407
0.5        10        30        0.0085
0.5        15        30        0.0042
0.5        20        30        0.0024
0.5        25        30        0.0014
0.5        30        30        0.001
1        10        30        0.0077
1        15        30        0.0031
1        20        30        0.0017
1        25        30        0.0013
1        30        30        0.00022884
1.5        10        30        0.0076
1.5        15        30        0.0043
1.5        20        30        0.0016
1.5        25        30        0.00069833
1.5        30        30        0.00032979
2        10        30        0.0073
2        15        30        0.0029
2        20        30        0.0016
2        25        30        0.00094093
2        30        30        0.00069881
2.5        10        30        0.0071
2.5        15        30        0.0026
2.5        20        30        0.0014
2.5        25        30        0.0011
2.5        30        30        0.00053801
3        10        30        0.0071
3        15        30        0.0032
3        20        30        0.0014
3        25        30        0.00086919
3        30        30        0.00045052
0.5        10        40        0.0096
0.5        15        40        0.0039
0.5        20        40        0.0018
0.5        25        40        0.0016
0.5        30        40        0.0013
1        10        40        0.0068
1        15        40        0.0032
1        20        40        0.0016
1        25        40        0.00067646
1        30        40        0.00071598
1.5        10        40        0.0078
1.5        15        40        0.003
1.5        20        40        0.0018
1.5        25        40        0.00070634
1.5        30        40        0.00052355
2        10        40        0.0061
2        15        40        0.0027
2        20        40        0.0015
2        25        40        0.00097435
2        30        40        0.00044496
2.5        10        40        0.0066
2.5        15        40        0.0027
2.5        20        40        0.0013
2.5        25        40        0.0007444
2.5        30        40        0.00033344
3        10        40        0.0054
3        15        40        0.003
3        20        40        0.0015
3        25        40        0.00068942
3        30        40        0.00026619
0.5        10        50        0.0101
0.5        15        50        0.0035
0.5        20        50        0.0018
0.5        25        50        0.00097358
0.5        30        50        0.00041595
1        10        50        0.0079
1        15        50        0.0037
1        20        50        0.0019
1        25        50        0.0013
1        30        50        0.00074984
1.5        10        50        0.007
1.5        15        50        0.0028
1.5        20        50        0.0015
1.5        25        50        0.00078604
1.5        30        50        0.00038018
2        10        50        0.0056
2        15        50        0.0027
2        20        50        0.0014
2        25        50        0.00087912
2        30        50        0.00034373
2.5        10        50        0.0049
2.5        15        50        0.0025
2.5        20        50        0.0016
2.5        25        50        0.0012
2.5        30        50        0.00075411
3        10        50        0.0062
3        15        50        0.0031
3        20        50        0.0013
3        25        50        0.00075392
3        30        50        0.00053634
0.5        10        60        0.0091
0.5        15        60        0.0034
0.5        20        60        0.0019
0.5        25        60        0.0015
0.5        30        60        0.00080194
1        10        60        0.0074
1        15        60        0.0033
1        20        60        0.0017
1        25        60        0.0012
1        30        60        0.00073223
1.5        10        60        0.006
1.5        15        60        0.0031
1.5        20        60        0.0015
1.5        25        60        0.00077907
1.5        30        60        0.00056884
2        10        60        0.0059
2        15        60        0.0026
2        20        60        0.0014
2        25        60        0.00093124
2        30        60        0.00061023
2.5        10        60        0.0051
2.5        15        60        0.0024
2.5        20        60        0.0015
2.5        25        60        0.00089189
2.5        30        60        0.00050778
3        10        60        0.0058
3        15        60        0.0026
3        20        60        0.0014
3        25        60        0.001
3        30        60        0.00040474

piaomiao2080

参数        最佳数据
a                  0.0216147637656143
b                  -0.00248619392665948
d                  -2.68298197946595E-5
e                  5.52289691682725E-5
f                  6.6289204408169E-8
g                  1.51821767280956E-6
h                  0.000265944941437366
i                  2.32708400258303E-5
j                  -1.27300706633786E-7

计算用时00:01:21:953
(RMSE): 0.000576041854122351
(SSE): 4.97736326551074E-5
(R): 0.973428872952514
(R^2): 0.947563770697601
(DC): 0.947563770697602
(Chi-Square): 0.0114337012336212
(F-Statistic): 2674.47602409562

handb

多谢！结果我分析分析，再问一下，系数c是0吗？

handb

这个结果为什么反不回去呢？例如，当在z=20时，得到的C值与我给出的值差很多（100多倍），怎么回事，请指教，等待中...

handb

是不是我的问题说的不清楚，为什么用1stopt拟合后，输出的数据与试验值相差很小，(SSE): 4.97736326551074E-5，但把x=0.5,y=10,z=20...数据代人拟合的表达式中，得到的C值与试验的结果明显不符呢，想知道这是什么原因？

handb

怎么没人关注呢，难道没有人遇到这样的问题吗？用1stopt拟合后，得到最佳值，把原来数据代到模型中，结果怎么就不吻合了呢？？？

shamohu

1stOpt是不区分大小写的，所以你的拟合公式：
Function C=a+b*x+c*y+d*z+e*x*y+f*x*z+g*y*z+h*x^2+i*y^2+j*z^2;
应改为：
Function C=a+b*x+c1*y+d*z+e*x*y+f*x*z+g*y*z+h*x^2+i*y^2+j*z^2;
可得结果：
均方差(RMSE): 0.000576041854122353
残差平方和(SSE): 4.97736326551077E-5
相关系数(R): 0.973428872952513
相关系数之平方(R^2): 0.9475637706976
决定系数(DC): 0.947563770697602
卡方系数(Chi-Square): 0.0114337011772157
F统计(F-Statistic): 296.656639169212

参数        最佳估算
----------        -------------
a                  0.0216147637644857
b                  -0.00248619394776179
c1                 -0.00139583098929429
d                  -2.68298195785315E-5
e                  5.52289691494157E-5
f                  6.62896124767356E-8
g                  1.51821766705512E-6
h                  0.000265944942934823
i                  2.32708399983663E-5
j                  -1.27300715656312E-7

再把你的数据代回去验证，就不会出错了。

handb

谢谢shamohu版主，还是版主厉害，这个问题折磨我很长时间了，真的很感谢，谢谢......