[R] All sub-summands of a vector
Gabor Grothendieck
ggrothendieck at gmail.com
Fri Apr 2 20:43:50 CEST 2010
There is also rollmean in the zoo package which might be slightly
faster since its optimized for that operation.
k * rollmean(x, k)
e.g.
> 2 * rollmean(1:4, 2)
[1] 3 5 7
will give a rolling sum. runmean in the caTools package is even faster.
On Fri, Apr 2, 2010 at 2:31 PM, Jorge Ivan Velez
<jorgeivanvelez at gmail.com> wrote:
> Hi Andy,
>
> Take a look at the rollapply function in the zoo package.
>
>> require(zoo)
> Loading required package: zoo
>> x <- 1:4
>> rollapply(zoo(x), 1, sum)
> 1 2 3 4
> 1 2 3 4
>> rollapply(zoo(x), 2, sum)
> 1 2 3
> 3 5 7
>> rollapply(zoo(x), 3, sum)
> 2 3
> 6 9
>> rollapply(zoo(x), 4, sum)
> 2
> 10
>
> # all at once
> sapply(1:4, function(r) rollapply(zoo(x), r, sum))
>
>
> HTH,
> Jorge
>
>
> On Fri, Apr 2, 2010 at 2:24 PM, Andy Rominger <> wrote:
>
>> Hello,
>>
>> I'd like to take all possible sub-summands of a vector in the quickest and
>> most efficient way possible. By "sub-summands" I mean for each sub-vector,
>> take its sum. Which is to say: if I had the vector
>>
>> x<-1:4
>>
>> I'd want the "sum" of x[1], x[2], etc. And then the sum of x[1:2], x[2:3],
>> etc. And then...so on.
>>
>> The result would be:
>> 1 2 3 4
>> 2 5 7
>> 6 9
>> 10
>>
>> I can do this with for loops (code below) but for long vectors (10^6
>> elements) looping takes more time than I'd like. Any suggestions?
>>
>> Thanks very much in advance--
>> Andy
>>
>>
>> # calculate sums of all sub-vectors...
>> x <- 1:4
>>
>> sub.vect <- vector("list",4)
>>
>> for(t in 1:4) {
>> maxi <- 4 - t + 1
>> this.sub <- numeric(maxi)
>> for(i in 1:maxi) {
>> this.sub[i] <- sum(x[i:(i+t-1)])
>> }
>> sub.vect[[t]] <- this.sub
>> }
>>
>> sub.vect
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> 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.
>>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.
>
More information about the R-help
mailing list