The R Project SVN R

Rev

Rev 71991 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 71991 Rev 72224
Line 3... Line 3...
3 
options(digits = 5)
 3 
options(digits = 5)
4 
 
 4 
 
5 
## UCV, BCV and SJ all gave wrong answers in example(bw.ucv)
 5 
## UCV, BCV and SJ all gave wrong answers in example(bw.ucv)
6 
c(nrd0 = bw.nrd0(precip), nrd = bw.nrd(precip),
 6 
c(nrd0 = bw.nrd0(precip), nrd = bw.nrd(precip),
7 
  ucv = bw.ucv(precip), bcv = bw.bcv(precip),
 7 
  ucv = bw.ucv(precip), bcv = bw.bcv(precip),
8 
  "SJ-dpi" = bw.SJ(precip, method="dpi"),
 8 
  "SJ-dpi" = bw.SJ(precip, method = "dpi"),
9 
  "SJ-ste" = bw.SJ(precip, method="ste"))
 9 
  "SJ-ste" = bw.SJ(precip, method = "ste"))
10 
 
 10 
 
11 
## wrong answers/errors in R < 3.3.3 for largish datasets
 11 
## wrong answers/errors in R < 3.3.3 for largish datasets
12 
set.seed(1); bw.bcv(rnorm(6000))
 12 
set.seed(1); bw.bcv(rnorm(6000))
13 
set.seed(1); x <- rnorm(47000)
 13 
set.seed(1); x <- rnorm(47000)
14 
bw.ucv(x); bw.SJ(x,  method="dpi"); bw.SJ(x,  method="ste")
 14 
bw.ucv(x); bw.SJ(x,  method = "dpi"); bw.SJ(x,  method = "ste")
15 
 
 15 
 
16 
## An extremely unbalanced example where counts exceed INT_MAX
 16 
## An extremely unbalanced example where counts exceed INT_MAX
- 
 
 17 
## Prior to R 3.4.0 was slow as O(n^2)
17 
set.seed(1); bw.SJ(c(runif(65537), 1e7))
 18 
set.seed(1); bw.SJ(c(runif(65537), 1e7))
18 
 
 19