27

假设我们有一个数字 1.000633，我想计算小数点后的零个数，直到分数中的第一个非零数字，答案应该是 3。对于 0.002，答案应该是 2。

R 中没有这样的功能可以提供帮助。我在 DescTools 包中的 Ndec 函数中进行了探索，但它没有完成这项工作。

4

7 回答 7

24

使用regexpr及其match.length参数

attr(regexpr("(?<=\\.)0+", x, perl = TRUE), "match.length")
于 2016-02-22T12:23:07.933 回答
15

这是另一种可能性：

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
于 2016-02-22T12:17:10.137 回答
10

我们可以用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))
于 2016-02-22T12:15:52.340 回答
7

使用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
于 2016-02-22T12:14:34.330 回答
7

另一种使用str_countfrom stringrpackage 的方式，

 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
于 2016-02-22T12:21:29.990 回答
5 
floor( -log10( eps + abs(x) - floor( abs( x ) ) ) )
于 2016-02-22T16:55:07.950 回答
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))
于 2018-03-05T07:34:56.700 回答