similar to: Revisiting multiple plots

I create a zooreg object that runs from Jan-1-2002 0:00 to Jun-1-2005 0:00... regts.start = ISOdatetime(2002, 1, 1, hour=0, min=0, sec=0, tz="") regts.end = ISOdatetime(2005, 6, 1, hour=0, min=0, sec=0, tz="") regts.zoo <- zooreg( NA, regts.start, regts.end, deltat=3600 ) Upon inspection: > regts.zoo[1:3] 2002-01-01 00:00:00 2002-01-01 01:00:00 2002-01-01 02:00:00

Dear R users, I have a quick question regarding creating an index for a zoo object. I am using R 2.12 on winxp. I have read through the R archives searching for information on dates and time series analysis and nothing seems to cover my question. I am extracting data from a time series of remote sensing data. The frequency of the data is 16 days which means I cant use ts (ts likes months single

Readers, I am trying to use the zoo package with an array of data: file1: hh:mm:ss 1 hh:mm:ss 2 hh:mm:ss 3 hh:mm:ss 4 file2: hh:mm:ss 11 55 hh:mm:ss 22 66 hh:mm:ss 33 77 hh:mm:ss 44 88 I wanted to merge these data set so I tried the following commands: library(chron) library(zoo) z1<-read.zoo("path/to/file1.csv",header=TRUE,sep=",",FUN=times)

Hi, Having a matrix F(189,6575) I want to do this: z1<-subset(F[,1], F[,1] >= 5 & F[,1] <= 10) . . . z189<-subset(F[,189], F[,189] >= 5 & F[,189] <= 10) I would prefer to have an empty matrix, say 'z' in order to fill its columns with the output of subsetting F. But each of the subsets can differ in length. It's to say: length(z1) 1189 length(z2) 1238

Hi Jakob, dat1<-read.table(text=" TIME, Value1, Value2 01.08.2011 02:30:00, 4.4, 4.7 01.09.2011 03:00:00, 4.2, 4.3 01.11.2011 01:00:00, 3.5, 4.3 01.12.2011 01:40:00, 3.4, 4.5 01.01.2012 02:00:00, 4.8, 5.3 01.02.2012 02:30:00, 4.9, 5.2 01.08.2012 02:30:00, 4.1, 4.7 01.12.2012 03:00:00, 4.7, 4.3 01.01.2013 01:00:00, 3, 4.3 01.01.2013 01:30:00, 3.8, 4.1 01.01.2013 02:00:00, 3.8,

Hi, I have attached the data files to this note. I use this code: library(zoo) z1 <- read.zoo("baltimorefludata.txt", format = "%m/%d/%Y", header = TRUE) z2 <- read.zoo("baltimorew.txt", format = "%Y%m%d", header = TRUE) z3<-merge(z1,z2) write.table(z3, "fluweatherdata_baltimore2.txt", sep="\t") R is writing the other data

Hi everyone and Achim, Achim, I appreciate your help about the function "NCOL". When I use "NCOL" instead of "ncol", I can find out the number of columns (number of time series) in the presence of only one time series (one variable, one column). Now I want to know how I can find out the column names of the zoo objects? In case of more than one time series, the

# chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) } z1 <- read.zoo(SC3.csv, sep = ",", header = TRUE, FUN = fmt.chron) z2 <- read.zoo(SC2.csv, sep = ",", header = TRUE, FUN = fmt.chron) z3<-c(z1, z2) write.table(z3, sep="," , "SC.csv") How do you

I noticed that under R 2.6.0 there is a warning about closing the connection in the code from this post: https://stat.ethz.ch/pipermail/r-help/2007-September/140601.html which is evidently related to the following from the NEWS file: o Connections will be closed if there is no R object referring to them. A warning is issued if this is done, either at garbage collection or if all the

Hi, I am interested in decomposing a time series and getting the trend, seasonal and?irregular variations, as one can get with the "stl" command. My time series is fairly regular, but it has some breaks. From the zoo manual, I gather that it should be possible to convert it to a regular time series and then fill the NA entries by interpolation. I am not able to proceed beyond a certain

Thanks for everyone's help with zoo -- I think I've got my data set ready. (The data consists of surface weather temperatures, from 2002 to 2005, one observation per hour. Some values are missing... i.e. NA) I have three goals: GOAL #1:Get the data in proper time series form, preserving frequency information: > w4.ts <- as.ts( w3.zoo, frequency=(1/3600) ) I hope that 1/3600

This is the same set of data that I have been working with for those in the know. it is a matrix of ~174 columns and ~70,000 rows. I have it as a zoo object, but I could read it in as just a matrix as long as the date time stamp won't be corrupted. here is an example of what a column would look like: 1/1/06 12:00, 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 1, 2, 3, 4, 5,6 ,7, NA #read in with the

DateTime RM61 11/30/2006 12:31 NA 11/30/2006 12:46 NA 11/30/2006 13:01 2646784125 11/30/2006 13:16 NA 11/30/2006 13:31 NA 11/30/2006 13:46 NA 11/30/2006 14:01 2666435177 11/30/2006 14:16 NA 11/30/2006 14:31 NA 11/30/2006 14:46 NA 11/30/2006 15:01 2653041914 11/30/2006 15:16 NA 11/30/2006 15:31 NA 11/30/2006 15:46 NA 11/30/2006 16:01 2693126189 11/30/2006 16:16 NA 11/30/2006 16:31 NA 11/30/2006

Hi, I have a problem with the appearance of legend in my plots. If I specify the legend positions in characters like "topright"..etc, it appears, if i specify it in terms of coordinates like "-1, 1" .. etc, it does not appear. Can anyone help me with this? script - x.date <- as.Date(paste(year, month, day, sep="-")) ts1.n.e3 <- ts(data.nemr.e3[,3]) z1.n.e3

Hello, I have been trying to plot correctly a graph for 2 month now, with no success.I want to put a grid, adjusted on the X axis tickers. Here is the way I build my X-Axis and my grid: grid(11, NULL, col="grey40") axis(2) # build the tickers on the beginning of each monthticks.at <- seq(ISOdatetime(lastYear, 01, 01, hour=0, min=0, sec=0, tz="GMT"),

Hi R Users! I would like to create an high frequency series but I am experiencing some difficulties. My series should start at 09.30 a.m. each day and end at 16.00 for, let's say, 2 years. I don't care on how many observations are for each day. It's ok also one observation each minute. In this case, I would have 390 observations each day. I have tried the following: start <-

Hi, when I run > plot.default(z1, z2, xlab = "x", ylab = "y", main = "xxxx", pch = "+") > abline(lm(z1 ~ z2)) then the plot is plotted perfectly (scatterplot), however, the lm() function doesnt appear on the plot. What could be wrong? (Yesterday it worked perfectly, with the lm() line.) Running R 2 on OS X. Mathias Hunsk??r Furevik Norway

Hi all, a problem I encounter again and again: I plot a zoo object using "plot" and then want to add lines or points to this plot. I usually circumvent this problem by adding artificial coloumns to the zoo object before plotting, but I am sure there's a better solution. To be specific: Assume I did x <-

Another newbie question I got the 1 minute spine interpolation and 15 mean aggregation working with many thanks to Gabor Grothendieck using Zoo functions. I got a tip from Hasan Diwan to look at xts but it seemed I would make better progress using code from Gabor. Now I'm having trouble plotting this zoo object. I'm thinking I want a function to "split" the zoo object back to

Hi all, I am using the zoo package to plot time series. I have a problem with formatting the axes. my zoo object (z) looks like the following. c1 1992-01-10 21 1992-01-17 34 1992-01-24 33 1992-01-31 41 1992-02-07 39 1992-02-14 38 1992-02-21 37 1992-02-28 28 1992-03-06 33 1992-03-13 40