[R] wapply from gplots -- How do I get a local estimate of variance?

[Fredrik Karlsson]

Dear list,

I am trying to plot the local variance in a moving window on a
dataset. The function that I am trying to use for this is wapply from
gtools.

However, from the lattice panel function code:

<snip>
	cat(x)
	cat(y)
  	wapply(x,y,method="range",width=1/10,fun=sd,na.rm=TRUE) -> outvar
</snip>

I get:

<snip>
109 109 109 109 109 109 116 116 116 116 119 119 123 123 123 123 127
127 133 133 133 133 133 138 138 138 138 138 138 138 138138 138 138 142
142 142 142 142 142 142 147 147 147 147 147 158 158 158 158 158 158
158 158 158 158 158 158 158
20.234 28.431 10.762 45.613 15.028 10.775 18.516 13.371 18.516 13.371
20.61 31.08 9.338 210.77 11.927 39.438 -35.079 15.872 12.272
9.1122.665 7.355 22.057 -49.289 11.212 16.236 20.654 16.236 20.654
15.181 57.271 37.513 57.271 37.513 8.518 -108.517 -14.444 10.702 4.482
16.422 23.003 42.451 12.998 42.451 12.998 14.292 5.945 7.115 16.079
-13.172 10.62 16.079 -13.172 10.62 16.079 -13.172 10.62 20.477 22.656
Error in var(x, na.rm = na.rm) : 'x' is empty
</snip>

What am I doing wrong with this? The variables that I think is
supplied to the var function does not seem to be empty like the error
message suggests.

The gplots package is from the 2.0.7 version of the gregmisc bundle.


/Fredrik Karlsson