2017-11-14 9 views
0

この例は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

o.segには、知りたいことがすべて表示されます。 「推定ブレークポイント:」は、ブレークポイントのx値(GDD値)(赤い点)を示します。それはpsi1.GDD = 391.3です。

o.seg$coefficientsは線の式を示しています。おおよそ-2.14 + 0.0097*GDD*(GDD < psi1.GDD) + 0.12*GDD*(GDD >= psi1.GDD)です。

関連する問題