24小时热门版块排行榜    

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

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

英鹤鸣

铜虫 (初入文坛)

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 16:25 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 16:00 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:57 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见