24小时热门版块排行榜    

查看: 3703  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 英鹤鸣 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] E0414, 我的本子有没有希望? +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[硕博家园] 一篇论文同时出现在两个期刊,一模一样,这算不算学术不端,请各位老师斧正。 +11 mahl005 2026-06-14 12/600 2026-06-18 17:01 by 770904134
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
请填处理意见