search for: ndeltat

Hi, I am getting unexpected behaviour from aggregate.ts(). The 'ndeltat' argument is effectively being reduced by 1 in some cases, even when it is an integer, with the result that the blocks to be aggregated are not of the expected size, and also that the end() of the aggregated series is much later than the end() of the original series. rawts <- ts(rep(1:10, e...

I've stumbled across an issue with aggregate.ts that either is due to a misuse of %/% or something deeper relating to numerical precision on Windows. The test code is x <- rep(6:10, 1:5) as.vector(aggregate(as.ts(x), FUN = mean, ndeltat = 5)) On Linux and Mac I get the correct answer > x <- rep(6:10, 1:5) > as.vector(aggregate(as.ts(x), FUN = mean, ndeltat = 5) [1] 7.2 8.8 10.0 > sessionInfo() R version 2.11.0 RC (2010-04-18 r51771) i386-apple-darwin9.8.0 locale: [1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_...

...d like to proceed with evaluating the function; otherwise, I would like NA returned. I'm not wedded to aggregate.ts and don't need ts class if there is another fast, vectorized way to handle this. Here is what I am trying to do, with the invalid na.rm thrown in: as.vector(aggregate.ts(x, ndeltat=24, FUN=min, na.rm=F)) Thanks, Scott Waichler Pacific Northwest National Laboratory scott.waichler at pnl.gov

...on circular data, but I can't figure out how to get the circular data components to pass through: library(circular) x <- circular(c(20, 30, 355, 5, 345, 25), units = "degrees") % I want to get the mean angle for each successive pair; answer should be c(25, 0, 5): aggregate.ts(x, ndeltat=2, FUN="mean.circular") Error in as.circular(x) : (converted from warning) an object is coerced to the class 'circular' using default value for the following components: type: 'angles' units: 'radians' template: 'none' modulo: 'asis' zer...

...the codes right to do what I want. Here are examples of unsuccessful attempts: try = aggregate(test.irts$value, by= list(time=test.ts$time), FUN= mean) ERROR MESS: Error in FUN(X[[1L]], ...) : arguments must have same length try = aggregate.ts(test.irts$value,frequency(test.ts$time),nfrequency=1, ndeltat=1) PROBLEM: gives me only one vector for my values (no time vector) but they are NOT aggregated by hour. try = list(x=test.ts$time, y=aggregate.ts(test.irts$value, nfrequency=1)) PROBLEM: Gives 2 separate vectors :hourly time step and values - but again, values are not aggregated by hour. I woul...

.../usr/local/lib), [DY]LIB_LIBRARY_PATH is now set without a trailing colon. (PR#13637) o More efficient utf8ToInt() on long multi-byte strings with many multi-byte characters. (PR#14262) o aggregate.ts() gave platform-depedent results due to rounding error for ndeltat != 1. o package.skeleton() sometimes failed to fix filenames for .R or .Rd files to start with an alphanumeric. (PR#14253) It also failed when only an S4 class without any methods was defined. (PR#14280) o splinefun(*, method = "monoH.FC") was not...

.../usr/local/lib), [DY]LIB_LIBRARY_PATH is now set without a trailing colon. (PR#13637) o More efficient utf8ToInt() on long multi-byte strings with many multi-byte characters. (PR#14262) o aggregate.ts() gave platform-depedent results due to rounding error for ndeltat != 1. o package.skeleton() sometimes failed to fix filenames for .R or .Rd files to start with an alphanumeric. (PR#14253) It also failed when only an S4 class without any methods was defined. (PR#14280) o splinefun(*, method = "monoH.FC") was not...

I'm having trouble with column-subsetting a multivariate time series. The underlying problem is with start(). R> x <- ts(matrix(rnorm(100), ncol=2)) R> y <- aggregate(x, ndeltat=2) R> y[,1] Error in ts(y, start = start(x), freq = frequency(x)) : invalid start R> start(y) [1] 2 1 R> frequency(y) [1] 0.5 The point where the error occurs in ts() is: if(length(start) > 1) { if(start[2] > frequency) stop("invalid start") .....