[R] function running in package gregmisc

[Don MacQueen]

At 9:23 AM +0200 7/23/02, Xavier.Abulker at fimat.com wrote:
>Hello,
>I've got a problem using the function "running" in package gregmisc
>
>For example:
>
>  test<-c(1,2,3,4,5)
>running(test,fun=var,width=3) gives
>1:1 1:2 1:3 2:4 3:5
>  NA  NA   2   3   4
>
>which is wrong because
>
>  var(test[1:3])
>[1] 1
>
>Where am I wrong?

Look at the definition of running:

>  running
function (X, fun = mean, width = min(length(X), 20), allow.fewer = FALSE,
     ...)
running2(X = X, fun = mean, width = width, allow.fewer = allow.fewer,
     ...)

You can see that running() calls running2(), but does not pass the 
value of its argument `fun' to running2(). So try running2() instead:

>  running2(1:5,fun=var,width=3)
1:1 1:2 1:3 2:4 3:5
  NA  NA   1   1   1


>Thanks
>Xavier
>

Trying some other functions besides var() can help understand what's going on:

>  running(1:5,fun=max,width=3)
1:1 1:2 1:3 2:4 3:5
  NA  NA   2   3   4
>  running(1:5,fun=function(x) 0,width=3)
1:1 1:2 1:3 2:4 3:5
  NA  NA   2   3   4
>  args(running)
function (X, fun = mean, width = min(length(X), 20), allow.fewer = FALSE,
     ...)
NULL

At this point the results suggest, even without inspecting the 
definition of running(), that running() ignores `fun' and always 
returns the mean.

Looks like the author needs to fix running(), but fortunately 
running2() is available.

-Don

-- 
--------------------------------------
Don MacQueen
Environmental Protection Department
Lawrence Livermore National Laboratory
Livermore, CA, USA
--------------------------------------
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._