packrat/lib-R/stats/tests/bandwidth.R

In UBC-MDS/Karl:

### Sanity checks for some of the Jan 2017 breakage.

options(digits = 5)

## UCV, BCV and SJ all gave wrong answers in example(bw.ucv)
c(nrd0 = bw.nrd0(precip), nrd = bw.nrd(precip),
  ucv = bw.ucv(precip), bcv = bw.bcv(precip),
  "SJ-dpi" = bw.SJ(precip, method = "dpi"),
  "SJ-ste" = bw.SJ(precip, method = "ste"))

## wrong answers/errors in R < 3.3.3 for largish datasets
set.seed(1); bw.bcv(rnorm(6000))
set.seed(1); x <- rnorm(47000)
bw.ucv(x); bw.SJ(x,  method = "dpi"); bw.SJ(x,  method = "ste")

## An extremely unbalanced example where counts exceed INT_MAX
## Prior to R 3.4.0 was slow as O(n^2)
set.seed(1); bw.SJ(c(runif(65537), 1e7))