24小时热门版块排行榜    

查看: 3922  |  回复: 6

英鹤鸣

铜虫 (初入文坛)

[求助] 求助SPSS进行光合光响应曲线的拟合的问题 已有1人参与

数据如下:
I          Pn(因变量)
1200  9.25
1000  8.75
800         7.55
600        6.2
400        4.3
200        1.5
100        -0.45
50        -1.6
0        -2.9
利用SPSS进行非线性拟合,公式如下
(Q * I + Pm - SQRT((Q * I + Pm) * (Q * I + Pm) - 4 * Q * I * K * Pm)) / (2 * K) - Rd
参数设置:Pm 30
                 K  0.5
                 Q 0.05
                 Rd  2
约束 ,定义参数约束 Pm《=50
                               0《=K《=1
                              0《Q《=0.125

进行曲线拟合,得到结果
怎么输出有三个结果呢?一个R平方为1的 Pm K Q Rd值(按道理不能为1吧?拟合就这么准??)
另外两个为R平方为0.997的Pm Pm K Q Rd值
我知道这些差异是迭代次数不同造成的,应该用那些值呢?
还有一些其他数据(跟上面的数据类似的)
,软件拟合后都不一样啊?有些输出只有R平方为1的Pm K Q Rd值 ,有些输出R平方为0.999、0.998,该用那些值呢???(这些值都不同……)
大虾门 求助。。。

[ Last edited by 英鹤鸣 on 2013-11-21 at 10:43 ]
回复此楼
见到你,我就高兴了!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

英鹤鸣

铜虫 (初入文坛)

R平方为1啊
求助SPSS进行光合光响应曲线的拟合的问题
捕获.PNG

见到你,我就高兴了!
2楼2013-11-21 10:49:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
英鹤鸣: 金币+5, ★★★很有帮助, 谢谢。很有帮助! 2013-11-26 22:47:26
用1stOpt试试:
CODE:
Parameter Pm<=50, 0<=K<=1, 0<=Q<=0.125;
Variable I, Pn;
Function Pn=(Q * I + Pm - SQRT((Q * I + Pm) * (Q * I + Pm) - 4 * Q * I * K * Pm)) / (2 * K) - Rd;
Data;
//I          Pn(因变量)
1200  9.25
1000  8.75
800         7.55
600        6.2
400        4.3
200        1.5
100        -0.45
50        -1.6
0        -2.9

均方差(RMSE):0.0665301868052333
残差平方和(SSE):0.0398363918070531
相关系数(R): 0.99988378138897
相关系数之平方(R^2): 0.999767576284706
确定系数(DC): 0.999767576284706
卡方系数(Chi-Square): 0.00227340160046171
F统计(F-Statistic): 7170.81156955291

参数        最佳估算
----------        -------------
pm        17.6529635884422
k        0.262942188623271
q        0.0271813330822702
rd        2.88963088638502

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

No        实测值pn        计算值pn
1        9.25        9.34760810868276
2        8.75        8.58486926486495
3        7.55        7.58601822020897
4        6.2        6.23165235277238
5        4.3        4.31763667363598
6        1.5        1.48458562613415
7        -0.45        -0.457850617675411
8        -1.6        -1.60488873434242
9        -2.9        -2.88963088638502
3楼2013-11-22 11:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

英鹤鸣

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by dingd at 2013-11-22 11:42:39
用1stOpt试试:

Parameter Pm<=50, 0<=K<=1, 0<=Q<=0.125;
Variable I, Pn;
Function Pn=(Q * I + Pm - SQRT((Q * I + Pm) * (Q * I + Pm) - 4 * Q * I * K * Pm)) / (2 * K) - Rd;
Data;
// ...

此软件的分析跟R平方为1的 Pm ,  K  ,Q , Rd值一样的,不明白为啥R平方为1?
见到你,我就高兴了!
4楼2013-11-22 18:59:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

应该是精度设定问题。
5楼2013-11-22 20:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kkzhang

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 英鹤鸣 at 2013-11-22 18:59:14
此软件的分析跟R平方为1的 Pm ,  K  ,Q , Rd值一样的,不明白为啥R平方为1?...

请问初始值怎样确定。谢谢
6楼2015-03-13 15:16:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kkzhang

新虫 (初入文坛)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 英鹤鸣 at 2013-11-22 18:59:14
此软件的分析跟R平方为1的 Pm ,  K  ,Q , Rd值一样的,不明白为啥R平方为1?...

你好,请问小麦旗叶的光响应曲线在spss中参数初始值怎样确定、、、急用!!!!!
7楼2015-03-13 17:10:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 英鹤鸣 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 14:13 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:53 by VOMexwnDzOFZ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 10:53 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见