You can get a slightly better X axis by converting your Date variable to chron:
Replace this: plot(mydata[,1], ...) with this: library(chron) plot(chron(unclass(mydata[,1])), ...) and ignore the warning. On Nov 13, 2007 9:08 AM, John Kane <[EMAIL PROTECTED]> wrote: > I clearly spoke too soon. > > With the actual data I am not getting sensible x-axis > units. The program with the actual data below. Graph > output is here: > http://ca.geocities.com/jrkrideau/R/hd.png . > > I seem to be getting only a single entry for the > x-axis of "2007". However dates range from > First Date Last Date > "2006-09-26" "2007-11-10" > > I must be missing something blindingly obvious but I > don't see it. > > Thanks for any suggestions. > > > Actual data and test program > ================================================ > mydata <- > read.table("http://ca.geocities.com/jrkrideau/R/heartdata.txt";, > sep="\t", > header=FALSE) > mydata[,1] <- as.Date(mydata[,1],"%m/%d/%y") > names(mydata) <- Cs(dates, sy,dys,pulse, weight) > minmax <- c(max(mydata[,2]),min(mydata[,2]), > max(mydata[,3]),min(mydata[,3])) > names(minmax) <- c("maxsys", "minsys", "maxdsy", > "mindys") ; minmax > > min.max.dates <- c(min(mydata[,1]), max(mydata[,1])) > names(min.max.dates) <- c("First Date", "Last Date"); > min.max.dates > > ss <- lm(mydata[,2]~mydata[,1]) > plot(mydata[,1],mydata[,2], xlab="Dates", ylab="Blood > pressure", > ylim=c(minmax[4], minmax[1]), col= "red", > type="l") > abline (ss, col ="yellow") > dd <- lm(mydata[,3]~mydata[,1]) > points(mydata[,1], mydata[,3], type="l", col="blue" ) > abline(dd, col= "yellow", lwd=2) > > > > > > --- John Kane <[EMAIL PROTECTED]> wrote: > > > Thanks to Gabor and Jim. I am not sure if the first > > entry year = 2009 is all the problem I'm getting but > > it is certainly seems like the worst of it. > > > > My stupidity: Someone sent me the data set in Excel > > > > and I didn't do the basic data checks on. I _KNEW_ > > the > > data went from 2006 to 2007. > > > > --- Gabor Grothendieck <[EMAIL PROTECTED]> > > wrote: > > > > > In your examples the first line of your data > > refers > > > to the > > > year 2009 and Oct 1st is repeated. Is that really > > > what > > > you meant? > > > > > > I can't tell what your problem is from your > > > description > > > other than the data problems cited but there are > > > lots of > > > examples of plotting with zoo in the following > > which > > > may > > > help you: > > > vignette("zoo") > > > vignette("zoo-quickref") > > > ?plot.zoo > > > ?xyplot.zoo > > > > > > Note that zoo series must be time series, i.e. > > they > > > must > > > have unique times. > > > > > > On Nov 12, 2007 1:47 PM, John Kane > > > <[EMAIL PROTECTED]> wrote: > > > > I am completely misunderstanding how to handle > > > dates. > > > > I want to plot a couple of data series against > > > some > > > > dates. Simple example 1 below works fine. > > > > Unfortunately I have multiple observations per > > day > > > (no > > > > time breakdowns) and observations across years. > > > > (example 2 very simplistic version ) > > > > > > > > Can anyone suggest a quick fix or point me to > > > > something to read? I thought that zoo might do > > it > > > but > > > > I seem to be missing something there too. > > > > > > > > Any suggestions gratefully recieved. > > > > > > > > > > > > Example 1 consecutive dates same year. > > > > > > ================================================= > > > > x <- "days > > > > 9/26/09 > > > > 9/27/06 > > > > 9/28/06 > > > > 9/29/06 > > > > 9/29/06 > > > > 9/29/06 > > > > 10/1/06 > > > > 10/1/06 > > > > 10/2/06 > > > > 10/3/06" > > > > > > > > mydata <- read.table(textConnection(x), > > > header=TRUE, > > > > as.is=TRUE); mydata > > > > > > > > mydates <- as.Date(mydata[,1], "%m/%d/%y"); > > > mydates > > > > mynums <- rnorm(10) > > > > plot(mydates, mynums) > > > > ================================================ > > > > Example 2 (things go blooy!) > > > > non-consecutive dates different years. > > > > > > > > > > ================================================= > > > > x <- "days > > > > 9/26/09 > > > > 9/27/06 > > > > 9/28/06 > > > > 9/29/06 > > > > 9/29/06 > > > > 9/29/06 > > > > 10/1/07 # <- year changes > > > > 10/1/07 > > > > 10/2/07 > > > > 10/3/07" > > > > > > > > mydata <- read.table(textConnection(x), > > > header=TRUE, > > > > as.is=TRUE); mydata > > > > > > > > mydates <- as.Date(mydata[,1], "%m/%d/%y"); > > > mydates > > > > mynums <- rnorm(10) > > > > plot(mydates, mynums) > > > > > > > > ______________________________________________ > > > > R-help@r-project.org mailing list > > > > 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 > > 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. > > > > > > Be smarter than spam. See how smart SpamGuard is at giving junk email > the boot with the All-new Yahoo! Mail. Click on Options in Mail and switch > to New Mail today or register for free at http://mail.yahoo.ca > ______________________________________________ R-help@r-project.org mailing list 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.
