[R] help with time series
Steve_Friedman at nps.gov
Steve_Friedman at nps.gov
Fri Jan 30 21:17:24 CET 2009
Hello everyone
I'm working with R 2.8.1 on a windows machine
I have a question regarding time series analysis
The first question is how does R expect the input file to be structured?
I'm working with a *.txt file similar to the abbreviated one here:
Date,stage
4/2/1953,7.56
4/3/1953,7.56
4/4/1953,7.54
4/5/1953,7.53
4/6/1953,7.5
4/7/1953,7.47
4/8/1953,7.44
4/9/1953,7.41
4/10/1953,7.37
4/11/1953,7.33
4/12/1953,7.3
4/13/1953,7.26
4/14/1953,7.28
4/15/1953,7.28
4/16/1953,7.23
4/17/1953,7.47
4/18/1953,7.59
4/19/1953,7.58
4/20/1953,7.57
4/21/1953,7.56
4/22/1953,7.55
4/23/1953,7.53
4/24/1953,7.51
4/25/1953,7.48
4/26/1953,7.46
4/27/1953,7.5
4/28/1953,7.56
The data record is substantially longer - 50 years worth of daily
hydrologic water stage data (column 2).
R seems to get confused by the format of this input not knowing what to do
with the date field, and also deciding to treat everything as a level.
I'm reading the data as follows:
mystage <- read.table("C:\\Documents and
Settings\\skfriedman\\Desktop\\R-scripts\\stage.txt", header = TRUE)
looking at the data I get the following:
mystage[1:4,]
[1] 4/2/1953,7.56 4/3/1953,7.56 4/4/1953,7.54 4/5/1953,7.53
20195 Levels: 1/1/1954,8.72 1/1/1955,8.48 1/1/1956,7.94 1/1/1957,7.88
1/1/1958,8.5 ... 9/9/2007,8.84
>
What I'd like is a time series with a starting data of April 21, 1953.
ending December 30, 2008. data are daily records, so the frequency should
be 365 (?) not counting leap year nuisances.
So the first question is how should I build the input file to correctly
import it to a time series with an odd beginning date?
The analysis I'm really trying to get to will involve calculating the mean
monthly stage, the mean seasonal (aggregated over several months) stage,
the annual maximum period with a continuous stage greater than 0.
Thanks in advance I will summary solutions.
Much appreciated
Steve
Steve Friedman Ph. D.
Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
Steve_Friedman at nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
More information about the R-help
mailing list