Предположим, если у нас есть число 1,000633, я хочу подсчитать количество нулей после запятой до первой ненулевой цифры в дроби, ответ должен быть 3. Для 0,002 ответ должен быть 2.
В R нет такой функции, которая могла бы помочь. Я изучил функцию Ndec в пакете DescTools, но она не работает.
Использование regexpr и его аргумента match.length
attr(regexpr("(?<=\\.)0+", x, perl = TRUE), "match.length")
regexpr вместо gregexpr, чтобы избежать sapply. Теперь он полностью векторизован и намного быстрее.
- person David Arenburg; 22.02.2016
x <- 10.2 это возвращает -1 вместо 0. Мне пришлось вставить оператор ifelse в свое решение, чтобы зафиксировать случай, когда без него не получится. Это может быть причиной того, что вы считаете мою реализацию сложной. С другой стороны, возможно, вы могли бы также рассмотреть возможность захвата таких случаев, чтобы ваше решение также работало для любого числа.
- person RHertel; 22.02.2016
-1 при отсутствии совпадения. Это обозначение regexpr для несоответствия. Мое решение работает для любого числа.
- person David Arenburg; 22.02.2016
(?<=\\.)0+|$ как регулярное выражение, если вы хотите получить 0 вместо -1.
- person maaartinus; 23.02.2016
Вот еще одна возможность:
zeros_after_period <- function(x) {
if (isTRUE(all.equal(round(x),x))) return (0) # y would be -Inf for integer values
y <- log10(abs(x)-floor(abs(x)))
ifelse(isTRUE(all.equal(round(y),y)), -y-1, -ceiling(y))} # corrects case ending with ..01
Пример:
x <- c(1.000633, 0.002, -10.01, 7.00010001, 62.01)
sapply(x,zeros_after_period)
#[1] 3 2 1 3 1
ifelse(round(y) == y, -y-1, -ceiling(y)) ?
- person David Arenburg; 22.02.2016
x <- c(0.1, 1.0, 1.001)
- person David Arenburg; 22.02.2016
y = log10(abs(x) %% 1) тоже работает. Думаю, чтобы сделать его векторизованным, y = -log10(abs(x) %% 1); ceiling(y) - ( (y %% 1) < 10^-options()$digits ) или использовать какой-то другой порог. Вероятно, все еще есть крайний случай или два
- person Frank; 22.02.2016
Мы можем использовать sub
ifelse(grepl("\\.0", str1),
nchar(sub("[^\\.]+\\.(0+)[^0]+.*", "\\1", str1)), NA)
#[1] 3 2 3 3 2
Или с помощью stringi
library(stringi)
r1 <- stri_extract(str1, regex="(?<=\\.)0+")
ifelse(is.na(r1), NA, nchar(r1))
#[1] 3 2 3 3 2
Просто чтобы проверить, работает ли он со странными случаями
str2 <- "0.00A-Z"
nchar(sub("[^\\.]+\\.(0+)[^0]+.*", "\\1", str2))
#[1] 2
str1 <- as.character(c(1.000633, 0.002, 0.000633,
10.000633, 3.0069006))
+
- person akrun; 22.02.2016
"[^\\.]+\\.(0+)[^0]{1}.*", и он будет найден (хотя я все же предпочитаю подход numeric RHertel). Дело в точном решении, а не в голосовании
- person Cath; 22.02.2016
you are online. Я сказал что-то плохое?
- person akrun; 22.02.2016
stringi, так как оно также считает NA за два символа... Извините за назойливость :)
- person David Arenburg; 22.02.2016
-1
- person akrun; 22.02.2016
Использование функции rle:
#test values
x <- c(0.000633,0.003,0.1,0.001,0.00633044,10.25,111.00012,-0.02)
#result
sapply(x, function(i){
myNum <- unlist(strsplit(as.character(i), ".", fixed = TRUE))[2]
myNumRle <- rle(unlist(strsplit(myNum, "")))
if(myNumRle$values[1] == 0) myNumRle$lengths[1] else 0
})
#output
# [1] 3 2 0 2 2 0 3 1
Другой способ, используя str_count из пакета stringr,
x <- as.character(1.000633)
str_count(gsub(".*[.]","",x), "0")
#[1] 3
РЕДАКТИРОВАТЬ: это подсчитывает все нули после десятичного числа и до первого ненулевого значения.
y <- c(1.00215, 1.010001, 50.000809058, 0.1)
str_count(gsub(".*[.]","",gsub("(?:(0+))[1-9].*","\\1",as.character(y))),"0")
#[1] 2 1 3 0
Вы можете использовать sub, так как нам не нужно прыгать. Таким образом, нет необходимости в gsub
nchar(sub(".*\\.(0*).*","\\1",str1))
[1] 3 2 3 3 2
куда
str1 <- as.character(c(1.000633, 0.002, 0.000633,
10.000633, 3.0069006))
count0 <- function(x, tol = .Machine$double.eps ^ 0.5) { x <- abs(x); y <- -log10(x - floor(x)); floor(y) - (y %% 1 < tol) }
- person Roland; 23.02.2016
