24小时热门版块排行榜    

查看: 1504  |  回复: 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的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

共享模式下的相关系数是他们的平均相关系数(总体残差平方和最小时),他们也有各自的相关系数。用的5.0版。

» 本帖已获得的红花(最新10朵)

6楼2013-04-18 09:29:02
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见