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 <jholt...@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 <vibhavasrivast...@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@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.
______________________________________________
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.
