育种数据分析之放飞自我分享 http://blog.sciencenet.cn/u/yijiaobai 关注:生物统计,数量遗传,混合线性模型,生物信息,R,Perl,Python,GWAS,GS相关方法,文章及代码

博文

如何用logistic回归曲线拟合日照和生物量积累的模型

已有 9855 次阅读 2017-12-22 17:18 |个人分类:农学统计|系统分类:科研笔记

数据

data

SPSS错误的做法

很遗憾,spss对于logistic默认的方法不是我们想要的!

错误的做法步骤:

1,打开模型,回归-→ 曲线回归 —-→ logistic

wrong model

2,查看结果

wrong result
这里y变量变成了 ln(1/生物量),
系数为:0.999,常量为:0.013,标准化的系数0.397.
很显然,怎么转换不能够和我们熟悉的logistic回归曲线联系到一起

正确的做法

用GenStat分析,GenStat有logistic回归曲线模型:
GenStat result

结果
GenStat result

可以看到各个参数都给出估计,模型可以解释97.6%的变异。

方程可以写为:
functions

或者将方程变化一下:
functions2

用spss正确的打开方式

需要设置方程组,需要设置初始值,麻烦的步骤如下:
1,打开模型
分析 —→ 回归 —→  非线性
spss model

2,设置拟合方程
这里我们需要手动设置方程组,设置初始值

参考文献:戴国俊, 王金玉, 杨建生,等. 应用统计软件SPSS拟合生长曲线方程[J]. 畜牧与兽医, 2006, 38(9):28-30.

方程组:y = C /(1 + (A EXP(-B GDD)))
初始值参考GenStat的结果:
start-values

3,结果
我设置了很多初始值,结果还是不理想,最后用文章中的参数,

结果还是不太令人满意,和GenStat相差较大

本来是想写出spss的解决方案,很遗憾,没有搞定。高手可以在留言板写出解决方案。目前看来,GenStat比较靠谱。

生物统计与数量遗传学公众号:
wechat

欢迎关注。



https://blog.sciencenet.cn/blog-2577109-1091005.html

上一篇:一年多点品种的回归系数和相关系数
下一篇:用矩阵的方法计算回归分析参数
收藏 IP: 111.202.84.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-23 17:20

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部