24小时热门版块排行榜    

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

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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-28 01:23 by xxniao123
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见