[R] Running sum

[Philippe Grosjean]

Ooops! Sorry. I did not look carefully enough on the code and focused on the
term! My proposition is for a "running sum with a window width of 2
observations" (the title of the question is "Running sum", isn't it?), while
the true question was indeed about a *cumulative sum*, which is something
totally different.
Best,

Philippe Grosjean 

> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch 
> [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of 
> Philippe Grosjean
> Sent: Friday, November 19, 2004 9:11 PM
> To: MSchwartz at MedAnalytics.com; 'Sean Davis'
> Cc: 'R-Help'
> Subject: RE: [R] Running sum
> 
> ?cumsum is not exactly the answer (as I understand it), but a 
> part of it.
> I propose:
> 
> runSum2 <- function(x)
> 	cumsum(x)[-1] - c(0, cumsum(x[1:(length(x) - 2)]))
> 
> # Example
> a <- round(runif(10, 0, 10))
> a
> runSum2(a)
> max(runSum2(a)) # To get only the max
> 
> Best,
> 
> Philippe
> 
> ..............................................<°}))><........
>  ) ) ) ) )
> ( ( ( ( (    Prof. Philippe Grosjean
>  ) ) ) ) )
> ( ( ( ( (    Numerical Ecology of Aquatic Systems
>  ) ) ) ) )   Mons-Hainaut University, Pentagone
> ( ( ( ( (    Academie Universitaire Wallonie-Bruxelles
>  ) ) ) ) )   6, av du Champ de Mars, 7000 Mons, Belgium  
> ( ( ( ( (       
>  ) ) ) ) )   phone: + 32.65.37.34.97, fax: + 32.65.37.33.12
> ( ( ( ( (    email: Philippe.Grosjean at umh.ac.be
>  ) ) ) ) )      
> ( ( ( ( (    web:   http://www.umh.ac.be/~econum
>  ) ) ) ) )
> ..............................................................
> 
>  
> 
> > -----Original Message-----
> > From: r-help-bounces at stat.math.ethz.ch 
> > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Marc Schwartz
> > Sent: Friday, November 19, 2004 7:57 PM
> > To: Sean Davis
> > Cc: R-Help
> > Subject: Re: [R] Running sum
> > 
> > On Fri, 2004-11-19 at 13:08 -0500, Sean Davis wrote:
> > > I have vector X of length N that I want to have a running sum for 
> > > (called Y).  I just need max(Y).  I do this with a "for"
> > loop like so:
> > > 
> > >      Y <- vector(length=N)
> > >      Y[1] <- X[1]
> > >      for (i in 2:N) {
> > >        Y[i] <- Y[i-1]+X[i]
> > >      }
> > >      return(max(Y))
> > > 
> > > Is there a faster way to do this?
> > > 
> > > Thanks,
> > > Sean
> > 
> > 
> > Something like:
> > 
> > > cumsum(1:10)
> >  [1]  1  3  6 10 15 21 28 36 45 55
> > 
> > > max(cumsum(1:10))
> > [1] 55
> > 
> > Does that help?
> > 
> > Marc Schwartz
> > 
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide! 
> > http://www.R-project.org/posting-guide.html
> > 
> >
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! 
> http://www.R-project.org/posting-guide.html
> 
>