

Dear All,

I'm trying to use dynlm to fit a time series.

I have 3 seasonal terms. Here is an example of the problem.


This is my time variable, hourly data:
timeSeries <- seq(as.POSIXct("2011-01-01 00:00:00"), as.POSIXct("2011-12-31 23:00:00"), by="hour") 


My response 
is:
y <- rnorm(length(t), 1000, 500) 

There are 3 seasonal factors:
t.h <- as.POSIXlt(t)$hour      # hours of the day

t.d <- as.POSIXlt(t)$wday      # days of the week
t.m <- as.POSIXlt(t)$mon            # months of the year

I can create a time series object for the hour, by doing

t.hh <- ts(t.h, freq=24)

but how to I create time series objects for the day (t.d) and the month (t.m)?

I tried 

t.dd <- ts(t.d, freq=t.d)

but this doesn't work.

I would appreciate any help.

Regards,
Dave

	[[alternative HTML version deleted]]

