24小时热门版块排行榜    

查看: 5209  |  回复: 10

littletingzi

新虫 (小有名气)

[求助] 用spss拟合生长曲线(非线性回归)其中参数初始值怎么确定?

我的实验是测定一种真菌的菌丝生长曲线,以下是时间(天数)和菌丝重量的数据。
time   weight
2        .0003
3        .0021
4        .0077
5        .0138
6        .0129
7        .0125
8        .0098
9        .0090
10        .0091
11        .0097
12        .0092
13        .0108

应用spss软件,采用曲线模型 Y=A/(1+B*EXP(-K*time))也可写成Y=A/(1+Be-kt)(此处-kt是e的指数)拟合,根据有人提到的三点法算出的参数初始值A=0.0109,B=0.4970,K=2.1320进行非线性拟合,但是结果却提示“EXP函数的参数太大”不知道是什么原因,向各位大虾求助,另外能不能教教我到底怎么给出参数的初始值,如何才能正确的得到拟合曲线?谢谢!!!
回复此楼

» 收录本帖的淘帖专辑推荐

科研工具、方法、技巧

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

匿名

用户注销 (小有名气)

本帖仅楼主可见
5楼2012-12-11 12:55:36
已阅   申请农林EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)

本帖仅楼主可见
7楼2012-12-11 19:08:12
已阅   申请农林EPI   回复此楼   编辑   查看我的主页
普通回帖

匿名

用户注销 (小有名气)

★ ★
感谢参与,应助指数 +1
littletingzi: 金币+2, 有帮助 2012-12-12 16:17:05
本帖仅楼主可见
2楼2012-12-09 22:02:14
已阅   申请农林EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)

本帖仅楼主可见
3楼2012-12-09 22:07:38
已阅   申请农林EPI   回复此楼   编辑   查看我的主页

littletingzi

新虫 (小有名气)

引用回帖:
3楼: Originally posted by guilaiguoke at 2012-12-09 22:07:38
补充一个,要算出参数的初始值应该先做出散点图,这样有助于计算,我明天晚上电脑能拿回来的话就给你试试。

不好意思啊,我忘记说了,Y代表的就是weight.ABK就像你说的那样都是待定参数。
我也看了一些用SPSS 做这个的但是都没有很清楚的说明怎么确定初始值,首次接触这个曲线拟合,所以不是很清楚原理,谢谢你的回复。那就麻烦你帮我试试,静等回复!
4楼2012-12-10 09:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

littletingzi

新虫 (小有名气)

引用回帖:
5楼: Originally posted by guilaiguoke at 2012-12-11 12:55:36
我用spss,得出B的系数特别大,我不知道对不对;

首先,从散点图上可以看出,大致是个S曲线形,
A代表了最大值,初始值:0.0138
B,取t趋近于0时,weight有最小值,取t=0,得到B=45
随便带一组数字进去,我选 ...

方差分析我跟你得到的是一样的,我当时也是觉得B的值特别大,我也搞不清楚是怎么回事。
因为我想得拟合得到的是真菌的生长曲线,我看文献里面采用的都是这个模型,所以我采用这个模型试试,但是发现结果很奇怪。你说B的建议性不好,主要是什么原因造成的呢,我应该怎么避免呢?
真的很感谢你这么详细的解答,我想知道如果出现这种情况下面应该怎么做呢,感谢你给我点建议!
6楼2012-12-11 14:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

littletingzi

新虫 (小有名气)

引用回帖:
7楼: Originally posted by guilaiguoke at 2012-12-11 19:08:12
这里就一个自变量一个因变量,因变量的weight肯定受好多个因素的影响,具体的话我也不是你们专业的,你可以看下别人文献里面的散点图,看看和你那个是不是一个样子的,看看人家的模型方程是怎么样的,然后从实验过 ...

真的很谢谢你!
我自己再琢磨琢磨……
8楼2012-12-11 19:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我就是我n109

银虫 (初入文坛)

您好!不知道楼主现在会分析了吗?我不知道这个初始值怎么设置的,而且我的藻的logistic方程是:μ=μm (1-X/Xm ),μ=1/X  dx/dt,怎么跟A,B,K对应关系呢?
9楼2013-11-27 14:17:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

18520195269

新虫 (初入文坛)

我也遇到了类似的问题,求教如何算出复杂的参数。
拟合曲线已知:y=a*exp[-b*exp(-kt)]
t表示测定的时间
y表示相应的时间t测出来的实验组体重
A为极限体重参数,k为生长速度参数,
b为达到最大生长速度时的日龄
t         y
1        0.258
1        0.251
1        0.255
2        0.564
2        0.495
2        0.503
3        0.788
3        0.781
3        0.818
4        1.422
4        1.51
4        1.551
5        2.312
5        2.233
5        2.34
10楼2015-01-12 14:22:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 littletingzi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 01:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
信息提示
请填处理意见