Re: [R] Gregmisc 'running' question

[Don MacQueen]

It says "not enough x observations". That's pretty clear.
So ask yourself, what is the minimum number of observations one needs 
to do a t test?

Try this to see how many observations t.test was given each time.

   running(dat,width=50,fun=function(x) length(x),allow=T)

or this
  running(dat,width=50,fun=function(x) cat('x:',x,'\n\n'),allow=T)
-Don

At 9:31 AM -0500 3/26/04, Sean Davis wrote:
Just a quick and probably simple question:

 dat <- rnorm(500,sd=1+(1:500)/500)
 fun <- function(x) t.test(x)$p.value
 running(dat,width=50,fun=fun,allow=T)
Error in t.test.default(x) : not enough x observations
 running(dat,width=50,fun=fun,allow=F)
Error in t.test.default(x) : not enough x observations
 fun2 <- function(x) mean(x)
 running(dat,width=50,fun=fun2,allow=T)
         1:1          1:2          1:3          1:4          1:5
1:6
-0.334134613 -0.626595581 -0.368967457 -0.113737178 -0.057448771
0.228643936
         1:7          1:8          1:9         1:10         1:11
1:12
-0.058807689  0.021762463 -0.063805657  0.031931121  0.080465708
0.087062800
....
However, this works fine.

 t.test(dat[1:50])$p.value
[1] 0.1661845

Why doesn't t.test work with running here?

Thanks,
Sean
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


--
--------------------------------------
Don MacQueen
Environmental Protection Department
Lawrence Livermore National Laboratory
Livermore, CA, USA
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html