need help with a time series plotting problem
Try this; this should put the labels on the x-axis like Excel does:
flow <- read.table(text = "Date USGS700 USGS1000 USGS1500
USGS1898 USGS1975 USGS2500 USGS2700 USGS2800
10/1/2001 0.05 0.57 2.32 2.27 4.11 29.45 29.45 29.45
10/1/2002 0.04 0.54 2.12 1.70 4.05 29.17 29.17 29.17
10/1/2003 0.03 0.48 1.93 1.98 3.96 28.88 28.88 28.88
10/1/2004 0.03 0.45 1.76 1.42 3.91 28.60 28.60 28.60
10/1/2005 0.03 0.42 1.64 1.27 3.82 28.32 28.32 28.32
10/1/2006 0.03 0.42 1.53 1.13 3.74 28.26 28.26 28.26
10/1/2007 0.11 0.51 1.59 5.66 3.68 28.23 28.23 28.23
10/1/2008 0.16 0.45 1.70 3.40 3.62 27.84 27.84 27.84
10/1/2009 0.10 0.42 1.78 2.55 3.54 26.56 26.56 26.56"
, as.is = TRUE
, header = TRUE
)
# setup margins
par(mar = c(5, 5, 3, 3))
# create x-axis with yearly tick marks starting in Oct
x_range <- seq(from = as.Date("2001-10-1")
, to = as.Date("2011-10-1")
, by = '1 year'
)
# now create the format of month-year
x_label <- sprintf("%s-%02d", format(x_range, "%b"),
as.POSIXlt(x_range)$year %% 100)
# plot each column of data
for (i in names(flow)[-1]){ # ignore the "Date" column
plot(as.Date(flow$Date, "%m/%d/%Y")
, flow[[i]] # column to plot
,type="l"
,xlab="date"
,ylab=expression("daily discharge (" * m^3/s * ")")
,main=i
,yaxs="i"
, xaxs="i"
, xaxt = 'n'
, xlim = range(x_range)
)
axis(1, at = x_range, labels = x_label, las = 2)
}
On Sat, Dec 24, 2011 at 3:23 PM, jim holtman <jholtman at gmail.com> wrote:
This should plot all the columns for you:
flow <- read.table(text = "Date ? ?USGS700 USGS1000 ? ? ? ?USGS1500
? ?USGS1898 ? ? ? USGS1975 ? ? ? ?USGS2500 ? ? ? ?USGS2700
USGS2800
10/1/2000 ? ? ? 0.05 ? ?0.57 ? ?2.32 ? ?2.27 ? ?4.11 ? ?29.45 ? 29.45 ? 29.45
10/2/2000 ? ? ? 0.04 ? ?0.54 ? ?2.12 ? ?1.70 ? ?4.05 ? ?29.17 ? 29.17 ? 29.17
10/3/2000 ? ? ? 0.03 ? ?0.48 ? ?1.93 ? ?1.98 ? ?3.96 ? ?28.88 ? 28.88 ? 28.88
10/4/2000 ? ? ? 0.03 ? ?0.45 ? ?1.76 ? ?1.42 ? ?3.91 ? ?28.60 ? 28.60 ? 28.60
10/5/2000 ? ? ? 0.03 ? ?0.42 ? ?1.64 ? ?1.27 ? ?3.82 ? ?28.32 ? 28.32 ? 28.32
10/6/2000 ? ? ? 0.03 ? ?0.42 ? ?1.53 ? ?1.13 ? ?3.74 ? ?28.26 ? 28.26 ? 28.26
10/7/2000 ? ? ? 0.11 ? ?0.51 ? ?1.59 ? ?5.66 ? ?3.68 ? ?28.23 ? 28.23 ? 28.23
10/8/2000 ? ? ? 0.16 ? ?0.45 ? ?1.70 ? ?3.40 ? ?3.62 ? ?27.84 ? 27.84 ? 27.84
10/9/2000 ? ? ? 0.10 ? ?0.42 ? ?1.78 ? ?2.55 ? ?3.54 ? ?26.56 ? 26.56 ? 26.56"
? ?, as.is = TRUE
? ?, header = TRUE
? ?)
# setup margins
par(mar = c(5, 5, 3, 3))
# plot each column of data
for (i in names(flow)[-1]){ ?# ignore the "Date" column
? ?plot(as.Date(flow$Date, "%m/%d/%y")
? ? ? ?, flow[[i]] ?# column to plot
? ? ? ?,type="l"
? ? ? ?,xlab="date"
? ? ? ?,ylab=expression("daily discharge (" * m^3/s * ")")
? ? ? ?,main=i
? ? ? ?,yaxs="i"
? ? ? ?, xaxs="i"
? ? ? ?)
}
On Sat, Dec 24, 2011 at 1:20 PM, vibhava <vibhavasrivastava at gmail.com> wrote:
sorry about the change in variable name (data has column named usgs700 but code has usgs700). actually i am trying to get familiar with R as i need to make more complex time series plots in near future (stackplots, scatterplot etc.). Let me try to explain what i am intending to achieve here. below i have copied first few lines of a large dataset (daily flow records for 4017 days at 8 usgs locations). first column is date and rest of them are flow data measured at these 8 different location. now all i wish to do is to read them in and make time series plots for the entire period of record for each stations. i have attached some figures that i made in excel and i wish to do something like this using R. I know my code might have some errors and that's the reason why i am requesting for help from people who know R better than I do. i am R user for less than a day but i know what i am trying to do is really simple and all i need is to read 9 columns and make a simple time series plot. i would appreciate if anyone can correct the code that i have written below or if they have some alternate way of doing this i would be happy to learn something new regards vibhava Date ? ?USGS700 USGS1000 ? ? ? ?USGS1500 ? ? ? ?USGS 1898 ? ? ? USGS1975 ? ? ? ?USGS2500 ? ? ? ?USGS2700 ? ? ? ?USGS2800 10/1/2000 ? ? ? 0.05 ? ?0.57 ? ?2.32 ? ?2.27 ? ?4.11 ? ?29.45 ? 29.45 ? 29.45 10/2/2000 ? ? ? 0.04 ? ?0.54 ? ?2.12 ? ?1.70 ? ?4.05 ? ?29.17 ? 29.17 ? 29.17 10/3/2000 ? ? ? 0.03 ? ?0.48 ? ?1.93 ? ?1.98 ? ?3.96 ? ?28.88 ? 28.88 ? 28.88 10/4/2000 ? ? ? 0.03 ? ?0.45 ? ?1.76 ? ?1.42 ? ?3.91 ? ?28.60 ? 28.60 ? 28.60 10/5/2000 ? ? ? 0.03 ? ?0.42 ? ?1.64 ? ?1.27 ? ?3.82 ? ?28.32 ? 28.32 ? 28.32 10/6/2000 ? ? ? 0.03 ? ?0.42 ? ?1.53 ? ?1.13 ? ?3.74 ? ?28.26 ? 28.26 ? 28.26 10/7/2000 ? ? ? 0.11 ? ?0.51 ? ?1.59 ? ?5.66 ? ?3.68 ? ?28.23 ? 28.23 ? 28.23 10/8/2000 ? ? ? 0.16 ? ?0.45 ? ?1.70 ? ?3.40 ? ?3.62 ? ?27.84 ? 27.84 ? 27.84 10/9/2000 ? ? ? 0.10 ? ?0.42 ? ?1.78 ? ?2.55 ? ?3.54 ? ?26.56 ? 26.56 ? 26.56
setwd("J:/Rstuff/flow")
# defining the working directory
flow=read.delim("flow.dat",header=TRUE,sep="\t")
# opening the above tab separated data file
plot(flow$USGS1500~as.Date(flow$Date,
"%m/%d/%y"),type="l",xlab="date",ylab="daily discharge (m3/s)
",main="USGS1500",yaxs="i", xaxs="i")
#just to test my code, i am trying to make time series
#plot of 1 variable USGS1500 http://r.789695.n4.nabble.com/file/n4231737/flow.docx flow.docx -- View this message in context: http://r.789695.n4.nabble.com/need-help-with-a-time-series-plotting-problem-tp4230672p4231737.html Sent from the R help mailing list archive at Nabble.com.
______________________________________________ R-help at 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.
-- Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it.
Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it.