Hi John, You seem to have 1569 days of data, so perhaps you can get around your axis problem like this:
plot(Nile,xaxt="n",xlab="Week") ... axis(1,at=seq(0,200,50),labels=seq(0,200,50)*7) (untested) Jim On Thu, Mar 7, 2019 at 10:46 AM Sparks, John <jspa...@uic.edu> wrote: > > Hi R Helpers, > > I am doing some work at identifying change points in time series data. A > very nice example is given in the R Bloggers post > > https://www.r-bloggers.com/a-look-at-strucchange-and-segmented/ > > The data for the aswan dam in that example is yearly. My data is weekly. I > ran the code switching the data for the analysis to my data and it worked, > but the scale of the line chart is not sensible. I have 225 weekly > observations and the x-axis of the line graph shows numbers from 0 to over > 1500. The information on the ts object is > > Start=1 > End=1569 > Frequency=0.1428... > > I can't share the data because it is proprietary. > > Wanting to be a good member of the list, I attempted to put weekly increments > on the Nile data so I could reproduce the x axis of the chart with the axis > scale that I am seeing. Unfortunately, in doing so I got another error that > I don't understand. > > > > library(strucchange) > library(lubridate) > library(xts) > > # example from R-Blog runs fine > data(“Nile”) > plot(Nile) > bp.nile <- breakpoints(Nile ~ 1) > ci.nile <- confint(bp.nile, breaks = 1) > lines(ci.nile) > > #problem comes in here > dfNile<-data.frame(Nile) > dfNile$week<-seq(ymd('2012-01-01'),ymd('2013-11-30'),by='weeks') > tsNile<-as.xts(x=dfNile[,-2],order.by=dfNile$week) > > Error in xts(x.mat, order.by = order.by, frequency = frequency(x), ...) : > formal argument "order.by" matched by multiple actual arguments > > > Can somebody help me to put together the ts object with weeks so that I can > demonstrate the problem with the scale on the x-axis and then try to get some > help with that original problem? > > Much appreciated. > --John Sparks > > > > > > > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@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.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@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.html and provide commented, minimal, self-contained, reproducible code.
