library(segmented)
olm<-lm(walevel~day)
#specify number and starting values for the breakpoints..
oseg<-segmented(olm, seg.Z=~day, psi=c(20,50,80))
plot(oseg,add=TRUE,col=2)
oseg$psi
#or you can avoid to specify starting values for psi
oseg<-segmented(olm, seg.Z=~day,
psi=NA,control=seg.control(stop.if.error=FALSE,K=30))
plot(oseg,add=TRUE,col=2)
oseg$psi
扫码加好友，拉您进群
全部版块
我的主页
收藏
