[R] Help with segmented package

Filoche pmassicotte at hotmail.com
Thu Mar 1 17:57:39 CET 2012


Hi everyone.

I'm using segmented package to find break point in a bi-linear relationship.
In a particular case, I find 1 pointcut (so 2 slopes).

I would like to know if it is possible to retrieve information in the
segmented object that could let me to plot 1 particular segment with a
different color.

For example, in that folowing example, I would like to plot the second
segment in red.

library(segmented)
set.seed(1234)
z<-runif(100)
y<-rpois(100,exp(2+1.8*pmax(z-.6,0)))
o<-glm(y~z,family=poisson)
o.seg<-segmented(o,seg.Z=~z,psi=list(z=.5))
plot(o.seg)

Thank in advance,
Phil



--
View this message in context: http://r.789695.n4.nabble.com/Help-with-segmented-package-tp4435550p4435550.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list