将多元回归方程式添加到分段图中

时间:2017-11-14 18:02:21

标签: r linear-regression trendline glmulti

我发现这个例子对RDocumentation很有帮助,我想进一步探索它。以下代码改编自https://www.rdocumentation.org/packages/segmented/versions/0.5-2.1/topics/plot.segmented

GDD<- c(221.2765957446810000,
    309.2198581560280000,
    431.2056737588650000,
    483.6879432624110000,
    553.1914893617020000)

biom<-c(0.0000000001000000,
        0.8510638297872340,
        5.9574468085106400,
        15.3191489361702000,
        22.1276595744681000)
o<- glm(biom~GDD, family = gaussian); o
o.seg<-segmented(o, ~GDD) #single segmented covariate and one breakpoint:'psi' can be omitted
par(mfrow=c(2,1))
plot(o.seg, conf.level=0.95, shade=TRUE)
points(o.seg, link=FALSE, col=2)

我想找到创建黑色分段线并找到关节点(红色)的等式。请看我能做到这一点。非常感谢任何帮助。enter image description here

1 个答案:

答案 0 :(得分:1)

o.seg告诉你你想知道的一切。 “Estimated Break-Point(s):”告诉您断点的红外线(GDD值)(红点)。它是psi1.GDD = 391.3。

o.seg$coefficients告诉您线方程,大致为-2.14 + 0.0097*GDD*(GDD < psi1.GDD) + 0.12*GDD*(GDD >= psi1.GDD)