Re: [R] Error with Segmented package

2015-10-10 Thread andrew haywood
Thanks Dave, when I run traceback() the following output Error in if (psi == Inf) psi <- median(XREGseg) : missing value where TRUE/FALSE needed > traceback() 2: segmented.lm(lm(y ~ x, data = data), seg.Z = ~x, psi = NA, control = seg.control(K = 1)) 1: segmented(lm(y ~ x, data = data), seg.Z =

Re: [R] Error with Segmented package

2015-10-10 Thread David Winsemius
On Oct 10, 2015, at 6:58 AM, andrew haywood wrote: > Thanks Dave, > > when I run traceback() the following output > > Error in if (psi == Inf) psi <- median(XREGseg) : > missing value where TRUE/FALSE needed > > traceback() > 2: segmented.lm(lm(y ~ x, data = data), seg.Z = ~x, psi = NA, con

Re: [R] Error with Segmented package

2015-10-07 Thread andrew haywood
Thanks David, I apologies for not posting the versions. I am running R version 3.2.2 (2015-08-14) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1 locale: [1] LC_COLLATE=English_United Kingdom.1252 [2] LC_CTYPE=English_United Kingdom.1252 [3] LC

Re: [R] Error with Segmented package

2015-10-07 Thread David Winsemius
On Oct 7, 2015, at 6:10 PM, andrew haywood wrote: > Thanks David, > > I apologies for not posting the versions. > > I am running > > R version 3.2.2 (2015-08-14) > Platform: x86_64-w64-mingw32/x64 (64-bit) > Running under: Windows 7 x64 (build 7601) Service Pack 1 > > locale: > [1] LC_COLLAT

Re: [R] Error with Segmented package

2015-10-07 Thread David Winsemius
On Oct 7, 2015, at 6:50 AM, andrew haywood wrote: > Dear List, > > I am trying to run a simple pieewise regression using the segmented package. > > When running the following code > > library(segmented) > data = data.frame(x=c(50,60,70,80,90,100,110) , y= > c(703.786,705.857,708.153,711.056,70

[R] Error with Segmented package

2015-10-07 Thread andrew haywood
Dear List, I am trying to run a simple pieewise regression using the segmented package. When running the following code library(segmented) data = data.frame(x=c(50,60,70,80,90,100,110) , y= c(703.786,705.857,708.153,711.056,709.257, 707.4, 705.6)) model.lm = segmented(lm(y~x,data = data),seg.Z