[R] cor.test observations limit

[Prof Brian Ripley]

Yes.  About 1290 for the Spearman test.

You need (n^3 - n) representable as an integer.

On Tue, 1 Apr 2003, Ming-Chung Li wrote:

> Is there a limit on the number of observations for using cor.test. For
> example,
> 
> > library(ctest)

That's not needed, and has not been for about 3 years.

> > cor.test(rnorm(3000), rnorm(3000), method="spearman")
> Error in if (q > (n^3 - n)/6) pspearman(q - 1, n, lower.tail = FALSE) else
> pspearman(q,  : 
>         missing value where logical needed
> In addition: Warning message: 
> NAs introduced by coercion 

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595