24小时热门版块排行榜    

查看: 3849  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 8/400 2026-08-26 11:03 by zjlp
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +24 木水思豆 2026-08-25 27/1350 2026-08-26 08:26 by yanchen918
[基金申请] 范进中举一文的中心思想 +8 炎黄贵胄 2026-08-22 9/450 2026-08-26 07:47 by WASM
[基金申请] 2026国自然函评费到账 +21 羊腰板 2026-08-21 24/1200 2026-08-26 07:20 by finnigan
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见