[R] How to extract "specific"/"last" intercept value from segmented package.

ruipbarradas at sapo.pt ruipbarradas at sapo.pt
Tue Jul 5 21:28:24 CEST 2016


Hello,

Try

dimnames(inter$x)[[1]]

You could have seen this by inspecting 'inter':

str(inter)

Hope this helps,

Rui Barradas
 

Citando Narendra Modi <bjpmodi2016 at gmail.com>:

> I am able to perform regression on a dataset as below:
>
> plot(x,y)
> lin.mod <- lm(y~x)
> m <- mean(x)
> m
>
> segmented.mod <- segmented(lin.mod, seg.Z = ~x, psi= m)
>
> plot(segmented.mod, add=T)
> sl <- slope(segmented.mod)
> inter <- intercept(segmented.mod)
>
> summary(segmented.mod)    # Show Summary
> sl                        # show all the slopes
> inter                     # show all the intercepts
>
> In my dataset, the above method correctly identifies the breakpoints and
> hence I get two intercepts.
>> inter
>
> $x
>              Est.
> intercept1  -3.269
> intercept2 -19.980
>
> What I am interested is the "intercept2" value. How can I obtain this?
>
> The method needs to be dynamic as in if the next dataset has 3 intercepts,
> I would like to get "intercept3 value.
>
> PD
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide  
> http://www.R-project.org/posting-guide.htmland provide commented,  
> minimal, self-contained, reproducible code.

 

	[[alternative HTML version deleted]]



More information about the R-help mailing list