I ' m пытается проверить, все ли элементы вектора равны друг другу. Решения, которые я придумал, кажутся несколько окольными, оба включают проверку length ()
.
x <- c(1, 2, 3, 4, 5, 6, 1) # FALSE
y <- rep(2, times = 7) # TRUE
С помощью unique ()
:
length(unique(x)) == 1
length(unique(y)) == 1
С rle ()
]:
length(rle(x)$values) == 1
length(rle(y)$values) == 1
Решение, которое позволило бы мне включить значение допуска для оценки «равенства» между элементами, было бы идеальным, чтобы избежать проблем FAQ 7.31 .
Есть ли встроенная функция для типа тест, который я полностью пропустил? идентичные ()
и all.equal ()
сравнивают два объекта R, поэтому они здесь не работают.
Edit 1
Вот некоторые результаты тестирования. Используя код:
library(rbenchmark)
John <- function() all( abs(x - mean(x)) < .Machine$double.eps ^ 0.5 )
DWin <- function() {diff(range(x)) < .Machine$double.eps ^ 0.5}
zero_range <- function() {
if (length(x) == 1) return(TRUE)
x <- range(x) / mean(x)
isTRUE(all.equal(x[1], x[2], tolerance = .Machine$double.eps ^ 0.5))
}
x <- runif(500000);
benchmark(John(), DWin(), zero_range(),
columns=c("test", "replications", "elapsed", "relative"),
order="relative", replications = 10000)
С результатами:
test replications elapsed relative
2 DWin() 10000 109.415 1.000000
3 zero_range() 10000 126.912 1.159914
1 John() 10000 208.463 1.905251
Таким образом, похоже, что diff (range (x)) <.machine double.eps> самый быстрый.