similar to: Time series

Hi everyone, I'm trying to plot some data across time. I have a list of articles, ranks, date/times, authors, etc. Someone suggested using zoo and someone suggested using ts. I'm pretty new at this and have been trying a simple if() plot() statement, but it doesn't seem to work. I keep getting an error that Title does not exist. Basically its saying that the headers do not exist as

Does anyone know of a test for stationarity of a time series, or like all ordination techniques it is a qualitative assessment of a quantitative result. Books, papers, etc. suggestions welcome. thanks Stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are

Hello R-help, I'm studying an example in the R book.? The data file is available from the link below.http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/fertilizer.txt Could you explain Why the results from lme() and lmer() are different in the following case? In other examples, I can get the same results using the two functions, but not here...? Thank you.Miya library(lme4)library(nlme)#

I have a time series object with two columns. a <- ts(x, frequency=1/15) I would like to split this into two time series of the same length (they are the same number of observations lined up time for time already). -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are

day<-structure(c(7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7, 7.69, 7.68, 7.68, 7.67, 7.67, 7.67, 7.66, 7.65, 7.65, 7.65, 7.64, 7.64, 7.63, 7.63, 7.63, 7.62, 7.62, 7.62, 7.62, 7.63, 7.63, 7.63, 7.63, 7.63, 7.64, 7.64, 7.65, 7.65, 7.65, 7.66, 7.66, 7.67, 7.67, 7.67, 7.68, 7.68, 7.69, 7.69, 7.69, 7.69, 7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7,

#this is a daily series of precipitation data. I would like to condense it into weekly means. How can I do this #as a side note I would like to do this same thing to two years worth of fifteen minute interval data and make it into #a series of daily averages (there are 96 readings per day) #is aggregate the right function? or... y <- c(1.23, 0, 0, 0, 0, 0, 0, 0, 0, 0.27, 0, 0.29, 0, 0, 0,

Hi all: I'm trying to plot two time series created in Rmetrics and label the x-axis with dates. I tried the following: dates <- as.Date(seriesPositions(x.agg)) r <- as.Date(range(dates)) ts.plot(x.agg@Data[,c(1,5)],gpars=list(ylab="Volume", lty=c(1:2),xaxt="n",main="Plot of Volume")) axis.Date(1, at=seq(r[1], r[2], length.out=10),

Certain timeseries I have had outliers, which I removed by assigning NA to their positions. Now acf() refuses to go to work. What's the right way to remove outliers from ts objects, and what are teh standard ways to interpolate NAs in them? Cheers, Alexy

Hello R List, I seem to have a peculiar problem. When using time series data, I get the following error when running the acf and pacf function. Using the function acf(dtxts,plot= TRUE,xaxt = "n",col="red",na.action = na.omit) (where dtxts is a time series object created with package "xts" ) results in the error below. Error in na.omit.ts(as.ts(x)) : time series

Hi R users, I am doing cross correlation analysis on 2 time series (call them y-series and x-series) where I need the use the model developed on the x-series to prewhiten the yseries.. Can someone point me to a function/filter in R that would allow me to do that? Thanks in advance for any help! -- View this message in context:

Hi all and thanks in advance. I am regressing Time and Weight, and then predicting Weight at different Time. The format of the Time data is day/month/year. How can I get R to use time series data such as this? Keith -- M. Keith Cox, Ph.D. Alaska NOAA Fisheries, National Marine Fisheries Service Auke Bay Laboratories 17109 Pt. Lena Loop Rd. Juneau, AK 99801 Keith.Cox at noaa.gov marlinkcox at

Hi, I have several time series that I need to plot on the same plot. There are 3 problems with these series: 1. they do not start or end at the same times 2. they have different time intervals (seconds, minutes or hours) 3. they all have random missing time steps of a few step to longer periods of several days The data is imported from a database as POSIXct so that the first column holds the

Hello, I'd like to use the add=TRUE parameter to add a second data line (series) to an existing plot, but R is giving me an error (see below). This code: > rap<-plot(aspen_sort,ylim=c(1,10000),log="y") ...produces the plot to which I'd like to add the second line. But this code: > rap<-plot(pine_sort,add = TRUE) ...produces this error: Warning messages: 1: In

Readers, I am having difficulty understanding how to enter commands into r. I have data arranged as: 100, 200, 300 5.6, 6.7, 7.8 8.9, 9.0, 0.1 1.2, 2.3, 3.4 The data is saved in csv format and I use the command 'read.table' to import into r. The values 5.6...3.4 are a function of values 100,...300, i.e. 100,...300 are independent variables (x). The function is a power series:

I would have expected the function getAnywhere to have behaved differently in the following: > search() [1] ".GlobalEnv" "file:C:/R/Rdata/miya/.Rdata" [3] "package:boot" "package:methods" [5] "package:ctest" "package:mva" [7] "package:modreg"

I have data that I read in using: data<-read.table("RAVK2.obs.data",sep="\t") 'data' looks like this: V1 V2 1 2009-03-25 06:00:00 12.86 2 2009-03-25 12:00:00 12.80 3 2009-03-25 18:00:00 12.76 4 2009-03-26 00:00:00 12.68 5 2009-03-26 06:00:00 12.66 6 2009-03-26 12:00:00 12.64 7 2009-03-26 18:00:00 12.83 8 2009-03-27 00:00:00 13.33 9

Hi, All I have a data frame like this. > names(x) [1] "month" "alloc_gb" One of the columns contains the dates I want to use: > data.class(x$month) [1] "yearmon" I would like to make a time series object out of this. However the ts function seems to require me to spell out the starting date. Is there a way to tell ts to use the date column I already

hello, i cant find a solution on this (might be) easy problem: i have a time serie by carlandar weeks, so for every carlendar week i have a value. now i would like to use the functions for time series, so i change structur to a time serie with cam <- ts(number,start=c(2001,1),deltat=7/365) or cam <- ts(number,start=c(2001,1),frequency=52) the problem on it is, that 2004 had 53 calendar

Hi and thanks in advance, I am fairly new with R so I hope this problem isn't too amateur. I have a vector of count data which correspond to vectors of date (%m/%d/%Y) and time of day (%H:%M:%S). I am trying to compute various statistics (e.g. daily max) by lumping the data together by day. I have been able to utilize tapply() and group the counts together, but with the method I use I end

Hello, I have a time-series of standards measured for Refractive index. They are daily standards, however, I didn't run one everyday so some days have no data. I can plot the values, but the x-axis does not represent the correct time series (i.e. it's just an evenly spaced 1,2,3 type axis). I want to plot the points with some form of representitive date line on the x-axis. I don't