24小时热门版块排行榜    

查看: 1482  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Will-Lee

金虫 (小有名气)

[求助] 求前辈帮忙用1stOpt你个一下数据,写论文急用,先谢过

x=0.125, 0.375, 1, 3, 7, 14, 21
y1=0.121, 0.07, 0.036, 0.197, 0.157, 0.190, 0.247
y2=0.136, 0.067, 0.03, 0.175, 0.134, 0.164, 0.219;

y1=a1*exp(-p*x)+(-k1-a1)*exp(-q*x)+k1;
y2=a2*exp(-p*x)+b*exp(-q*x)+(-k2-a2-b)*exp(-c*x)+k2。

要求这两个方程共享参数p和q(即,这两个解有公共的参数p,q),同时对这两个方程进行拟合。论文急用,后辈不胜感激。

[ Last edited by Will-Lee on 2013-4-17 at 09:47 ]
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Will-Lee

金虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-04-17 17:21:34
数据1:
均方差(RMSE): 0.0366860912198703
残差平方和(SSE): 0.018842170045897
相关系数(R): 0.970917816614817
相关系数之平方(R^2): 0.942681406620082
决定系数(DC): 0.942681404186766
F统计(F-Statist ...

非常感谢前辈的帮助啊!我还有一个疑问想请教您!这样同时拟合y1,y2得到了相关系数之平方(R^2)以及各个参数值,那我现在把参数分别代入y1,y2会得到它们各自解析解的表达式,此时它们的表达式与它们的实际数据单独的相关性是不是不用考虑了啊?就是说,共享拟合后,y1的表达式和y1的数据;以及y2的表达式和y2的数据之间,各自的这个相关性就无关紧要了是吗?只要考虑这个软件输出的两个同时拟合得到的这个(R^2)是吗?请前辈指教啊!
另外,我在论文里需要写使用的软件版本,麻烦您可以告诉我您运行的1st0pt的版本吗?万分感谢!!!
5楼2013-04-18 09:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
csgt0: 金币+2, 谢谢 2013-04-17 14:08:29
公式应该可以化简下吧,7组数据,8个参数,结果不唯一:
1:
均方差(RMSE): 0.0366816872715916
残差平方和(SSE): 0.0188376465352719
相关系数(R): 0.97092490177028
相关系数之平方(R^2): 0.942695164877628
决定系数(DC): 0.942695164877628
F统计(F-Statistic): -0.407553856793672

参数                  最佳估算
--------------------        -------------
a1        -0.0715119190800517
p        19188.0541782803
k1        0.2168405771542
q        0.228631851923395
a2        -0.251784621263395
b        -0.167788449604747
k2        0.196740312405319
c        6.19744351252297

2:
均方差(RMSE): 0.0366816872715916
残差平方和(SSE): 0.0188376465352719
相关系数(R): 0.97092490177028
相关系数之平方(R^2): 0.942695164877628
决定系数(DC): 0.942695164877628
F统计(F-Statistic): -0.407553856660192

参数                  最佳估算
--------------------        -------------
a1        -0.0715119197388559
p        8709.26578127736
k1        0.216840576740686
q        0.228631852581812
a2        -0.251784628311086
b        -0.16778845043493
k2        0.196740312651761
c        6.19744356226596

3:
均方差(RMSE): 0.0366816872715916
残差平方和(SSE): 0.0188376465352719
相关系数(R): 0.97092490177028
相关系数之平方(R^2): 0.942695164877628
决定系数(DC): 0.942695164877628
F统计(F-Statistic): -0.407553855552041

参数                  最佳估算
--------------------        -------------
a1        -0.0715119191198544
p        1209.82711843774
k1        0.216840576170168
q        0.22863185804796
a2        -0.251784622436768
b        -0.167788450839218
k2        0.196740311949434
c        6.1974433956109
2楼2013-04-17 10:19:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Will-Lee

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2013-04-17 10:19:38
公式应该可以化简下吧,7组数据,8个参数,结果不唯一:
1:
均方差(RMSE): 0.0366816872715916
残差平方和(SSE): 0.0188376465352719
相关系数(R): 0.97092490177028
相关系数之平方(R^2): 0.942695164877628 ...

麻烦前辈了啊,上午比较匆忙,发的那个公式有一个参数没有简化,现在重新计算了一下,把参数c消去了,
c=k1^2*p*q/(k1+k2)*[a1*(q-p)+q*k1]
前辈帮我重新拟合一下哈,太感谢了。我还有另一套数据,公式都是用这个进行共享拟合:
y1=a1*exp(-p*x)+(-k1-a1)*exp(-q*x)+k1;
y2=a2*exp(-p*x)+b*exp(-q*x)+(-k2-a2-b)*exp(-c*x)+k2。
数据1:
x=0.125, 0.375, 1, 3, 7, 14, 21
y1=0.121, 0.07, 0.036, 0.197, 0.157, 0.190, 0.247
y2=0.136, 0.067, 0.03, 0.175, 0.134, 0.164, 0.219
数据2:
x=0.125, 0.375, 1, 3, 7, 14, 21
y1=0.042,0.081,0.003,0.014,0.007,0.01,0.007
y2=0.034,0.069,0.004,0.017,0.011,0.013,0.01
非常感谢前辈啊,另外想问一下,前辈方便告诉我,您用的1stopt是什么版本的吗?感激涕零啊!
3楼2013-04-17 17:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Will-Lee: 金币+40, ★★★★★最佳答案 2013-04-17 22:00:49
csgt0: 金币+2, 谢谢 2013-04-18 14:28:08
数据1:
均方差(RMSE): 0.0366860912198703
残差平方和(SSE): 0.018842170045897
相关系数(R): 0.970917816614817
相关系数之平方(R^2): 0.942681406620082
决定系数(DC): 0.942681404186766
F统计(F-Statistic): 0

参数                  最佳估算
--------------------        -------------
a1        -0.0716981074946548
p        58.3099186235832
k1        0.217229689427639
q        0.225630583468914
a2        -0.254590627940963
b        -0.167186991575341
k2        0.197032981018483

====== 输出结果 =====

文件:数据文件-1
No        目标y1        计算y1        目标y2        计算y2
1        0.121        0.0756963295350097        0.136        0.136388627375288
2        0.07        0.083505197585773        0.067        0.064461306895123
3        0.036        0.101093617494315        0.03        0.0640223025541952
4        0.197        0.143271393143855        0.175        0.112069533911774
5        0.157        0.187236239060414        0.134        0.162576441295238
6        0.19        0.211048164751123        0.164        0.189931632292051
7        0.247        0.215955703045328        0.219        0.195569422746543

数据2:
均方差(RMSE): 0.0102388830426604
残差平方和(SSE): 0.0014676861634579
相关系数(R): 0.95095898439747
相关系数之平方(R^2): 0.904322990006267
决定系数(DC): 0.904322935889315
F统计(F-Statistic): 0

参数                  最佳估算
--------------------        -------------
a1        -140.208201956267
p        3.82503293783713
k1        0.00685029735039128
q        3.820230451081
a2        -114.246613332228
b        114.238791246974
k2        0.0246065387831164

====== 输出结果 =====

文件:数据文件-1
No        目标y1        计算y1        目标y2        计算y2
1        0.042        0.0547962273963608        0.034        0.0455324496547414
2        0.081        0.0654316902880739        0.069        0.0551172638680705
3        0.003        0.0214263697654643        0.004        0.0199040703785599
4        0.014        0.00687135622292324        0.017        0.00858973741779155
5        0.007        0.00685029736165498        0.011        0.00952134780710019
6        0.01        0.00685029735039128        0.013        0.0110485765490737
7        0.007        0.00685029735039128        0.01        0.0124211882548926
4楼2013-04-17 17:21:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +7 爱看书的可乐 2026-08-26 8/400 2026-08-26 10:50 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 系统进不去 +3 yanglien 2026-08-26 4/200 2026-08-26 10:40 by bycb2005
[基金申请] 怎么看青基中了没有啊 +3 叶九微 2026-08-26 3/150 2026-08-26 10:38 by LemmonTr
[基金申请] 科研孤儿太难了 +19 我4大白菜 2026-08-20 20/1000 2026-08-26 09:49 by zzuzxg
[基金申请] 出来了 +8 trojank 2026-08-26 8/400 2026-08-26 09:00 by gxc
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 范进中举一文的中心思想 +8 炎黄贵胄 2026-08-22 9/450 2026-08-26 07:47 by WASM
[基金申请] 2026国自然函评费到账 +21 羊腰板 2026-08-21 24/1200 2026-08-26 07:20 by finnigan
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
请填处理意见