24小时热门版块排行榜    

CyRhmU.jpeg
查看: 943  |  回复: 4

070166017

木虫 (正式写手)

[求助] 求助!求助!小女子快崩溃了,关于1STOPT的拟合问题,谢谢各位大神啦,求解答啊!已有1人参与

各位大神,求助!急求!急求!急求!急求!最近要进行数据拟合,得到其方程,我最近看了可以使用1STOPT,就自己尝试了一下,完全是新手啊,我的数据还都是离散的三维数据,压根没有方程直接可以拿过来用然后求系数,也没有可以参考的方程。别人说可以用1STOPT的自动搜索功能,就用了,但问题有三个:
1,得出的方程很复杂,参数有10个左右,我就两个自变量,一个时间,一个温度,因变量是在不同温度和时间的组合下测得的指标含量(举个例子:蛋白质含量等),是不是自动搜索功能得到的方程都是这么复杂?能不能精简方程呢?
2,尝试做了一下,结果如下,但是不同算法得出的结果不一样,我又不知道哪种算法好,怎么选择?
3,用其中一种算法得到的R2为0.852048,但看三维图时候有太多的点都不在曲面上,怎么回事?为什么R2还会这么大?
小女子求助大神了,已经查了好几天资料了,没有头绪,快崩溃了,谢谢啦!!!
数据如下:第一列x是时间,第二列y是温度,第三列z是指标
10        75        91.78
20        75        90.28
30        75        92.94
40        75        85.88
50        75        85.43
60        75        84.39
90        75        78.60
120        75        77.04
10        80        87.39
20        80        92.60
30        80        79.82
40        80        82.68
50        80        78.95
60        80        77.04
90        80        74.89
120        80        75.85
10        85        86.13
20        85        86.61
30        85        80.97
40        85        67.10
50        85        79.18
60        85        72.05
90        85        73.52
120        85        71.17
10        90        89.16
20        90        80.05
30        90        78.31
40        90        73.13
50        90        80.74
60        90        74.80
90        90        76.87
120        90        68.13
10        95        88.93
20        95        81.60
30        95        75.31
40        95        71.52
50        95        58.31
60        95        65.60
90        95        66.79
120        95        66.52
10        100        77.47
20        100        74.34
30        100        65.94
40        100        64.97
50        100        62.98
60        100        61.58
90        100        59.09
做出的一张图如下:附件中
问题是这么多的点都不在曲面上,为什么R2还这么高?

求助!求助!小女子快崩溃了,关于1STOPT的拟合问题,谢谢各位大神啦,求解答啊!
44.jpg


求助!求助!小女子快崩溃了,关于1STOPT的拟合问题,谢谢各位大神啦,求解答啊!-1
44.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

留住的是幸福,丢失的是遗憾!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
070166017: 金币+10, ★★★很有帮助, 做出来的效果很好,但是希望解答小女子的疑惑啊,毕竟后面还有大量数据要处理呀,谢谢啦,!!! 2015-10-17 10:13:32
原数据就比较离散,下面供参考:
z = p1+p2*ln(x)+p3*y+p4*y^2+p5*y^3+p6*(x+p8*y)^p7

均方差(RMSE):3.63072801752595
残差平方和(SSE):619.562739050652
相关系数(R): 0.915134781140587
相关系数之平方(R^2): 0.83747166765323
确定系数(DC): 0.83747166765323
卡方系数(Chi-Square): 4.16038253332308
F统计(F-Statistic): 34.2797375139579

参数        最佳估算
----------        -------------
p1        2063.99025358502
p2        -12.7732050903734
p3        -65.8668598387318
p4        0.756750204657691
p5        -0.00292102719737957
p6        -15636.9896782746
p8        1.08916742249694
p7        -1.4379572915664
求助!求助!小女子快崩溃了,关于1STOPT的拟合问题,谢谢各位大神啦,求解答啊!-2
c1.jpg


求助!求助!小女子快崩溃了,关于1STOPT的拟合问题,谢谢各位大神啦,求解答啊!-3
c2.jpg

2楼2015-10-16 10:03:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

070166017

木虫 (正式写手)

你好,昨天有事情没有及时回复,为什么咱们做出来的方程不一样呢?我的是Function: z = (p1+p2*Ln(x)+p3*(Ln(x))^2+p4*Ln(y)+p5*(Ln(y))^2+p6*(Ln(y))^3)/(1+p7*Ln(x)+p8*(Ln(x))^2+p9*Ln(y)),方程很复杂,我也不知道做出的这个方程有没有意义,该如何确定方程是否有意义?然后三维图就是上面上传的那个,请问你的方程是怎么做出的,不是自动搜索做出来的吗?如果是自动搜索,那咱们应该搜索结果一样啊,从三维图上看比我的好很多,大部分数据点都在曲面上,我的是什么问题,查资料也没有找到解答,求解答,麻烦了,谢谢!!!!!!!!!!!!
留住的是幸福,丢失的是遗憾!
3楼2015-10-17 10:05:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

070166017

木虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2015-10-16 10:03:19
原数据就比较离散,下面供参考:
z = p1+p2*ln(x)+p3*y+p4*y^2+p5*y^3+p6*(x+p8*y)^p7

均方差(RMSE):3.63072801752595
残差平方和(SSE):619.562739050652
相关系数(R): 0.915134781140587
相关系数之平方(R^ ...

你好,昨天有事情没有及时回复,为什么咱们做出来的方程不一样呢?我的是Function: z = (p1+p2*Ln(x)+p3*(Ln(x))^2+p4*Ln(y)+p5*(Ln(y))^2+p6*(Ln(y))^3)/(1+p7*Ln(x)+p8*(Ln(x))^2+p9*Ln(y)),方程很复杂,我也不知道做出的这个方程有没有意义,该如何确定方程是否有意义?然后三维图就是上面上传的那个,请问你的方程是怎么做出的,不是自动搜索做出来的吗?如果是自动搜索,那咱们应该搜索结果一样啊,从三维图上看比我的好很多,大部分数据点都在曲面上,我的是什么问题,查资料也没有找到解答,求解答,麻烦了,谢谢!!!!!!!!!!!!
留住的是幸福,丢失的是遗憾!
4楼2015-10-17 10:10:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浮云几度

银虫 (小有名气)

你好,请问你之前发的关于用1stopt拟合参数的问题解决了么?我也想用1stopt拟合参数,参数有9个,但是只能找到5.0破解版的,这个版本最多只能有4个参数。所以你最后用什么进行拟合的?matlab么?
5楼2015-11-18 16:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 070166017 的主题更新
信息提示
请填处理意见