具有 2 个或更多有效数字的地板和天花板

Floor and ceiling with 2 or more significant digits

可以使用 signif:

将结果四舍五入为两位有效数字
> signif(12500,2)
[1] 12000
> signif(12501,2)
[1] 13000

但是是否有同样方便的函数,例如 signif.floorsignif.ceiling 下方的 fictitious 函数,以便我可以获得两位或更多位有效数字有地板还是天花板?

> signif.ceiling(12500,2)
[1] 13000
> signif.floor(12501,2)
[1] 12000

编辑:

现有的 signif 函数适用于负数和小数。

因此,可能的解决方案最好也适用于负数:

> signif(-125,2)
[1] -120
> signif.floor(-125,2)
[1] -130

和小数:

> signif(1.23,2)
[1] 1.2
> signif.ceiling(1.23,2)
[1] 1.3

作为特例,0也应该return 0:

> signif.floor(0,2)
[1] 0

编辑 远不及@storaged 的​​回答那么好,但我已经开始了,所以我不妨完成:

基本贯穿每一个case(正负小数与否)

signif.floor=function(x,n){
  if(x==0)(out=0)
  if(x%%round(x)==0 & sign(x)==1){out=as.numeric(paste0(el(strsplit(as.character(x),''))[1:n],collapse=''))*10^(nchar(x)-n)}
  if(x%%round(x) >0 & sign(x)==1){out=as.numeric(paste0(el(strsplit(as.character(x),''))[1:(n+1)],collapse=''))}
  if(x%%round(x)==0 & sign(x)==-1){out=(as.numeric(paste0(el(strsplit(as.character(x),''))[1:(n+1)],collapse=''))-1)*10^(nchar(x)-n-1)}
  if(x%%round(x) <0 & sign(x)==-1){out=as.numeric(paste0(el(strsplit(as.character(x),''))[1:(n+2)],collapse=''))-+10^(-n+1)}
  return(out)
}

signif.ceiling=function(x,n){
  if(x==0)(out=0)
  if(x%%round(x)==0 & sign(x)==1){out=(as.numeric(paste0(el(strsplit(as.character(x),''))[1:n],collapse=''))+1)*10^(nchar(x)-n)}
  if(x%%round(x) >0 & sign(x)==1){out=as.numeric(paste0(el(strsplit(as.character(x),''))[1:(n+1)],collapse=''))+10^(-n+1)}
  if(x%%round(x)==0 & sign(x)==-1){out=(as.numeric(paste0(el(strsplit(as.character(x),''))[1:(n+1)],collapse='')))*10^(nchar(x)-n-1)}
  if(x%%round(x) < 0 & sign(x)==-1){out=as.numeric(paste0(el(strsplit(as.character(x),''))[1:(n+2)],collapse=''))}
  return(out)
}

我认为这种方法适用于所有类型的数字(即整数、负数、小数)。

地板函数

signif.floor <- function(x, n){
  pow <- floor( log10( abs(x) ) ) + 1 - n
  y <- floor(x / 10 ^ pow) * 10^pow
  # handle the x = 0 case
  y[x==0] <- 0
  y
}

上限函数

signif.ceiling <- function(x, n){
  pow <- floor( log10( abs(x) ) ) + 1 - n
  y <- ceiling(x / 10 ^ pow) * 10^pow
  # handle the x = 0 case
  y[x==0] <- 0
  y
}

他们都做同样的事情。首先计算位数,然后使用标准的 floor/ceiling 函数。检查它是否适合您。

编辑 1 按照 Heikki 的评论中的建议添加了 x = 0 案例的处理程序。

编辑 2 再次跟随 Heikki 我添加一些示例:

测试 x

的不同值
# for negative values
> values <- -0.12151 * 10^(0:4); values
# [1]    -0.12151    -1.21510   -12.15100  -121.51000 -1215.10000
> sapply(values, function(x) signif.floor(x, 2))
# [1]    -0.13    -1.30   -13.00  -130.00 -1300.00
> sapply(values, function(x) signif.ceiling(x, 2))
# [1]    -0.12    -1.20   -12.00  -120.00 -1200.00

# for positive values
> sapply(-values, function(x) signif.floor(x, 2))
# [1]    0.12    1.20   12.00  120.00 1200.00
> sapply(-values, function(x) signif.ceiling(x, 2))
# [1]    0.13    1.30   13.00  130.00 1300.00

测试 n

的不同值
> sapply(1:5, function(n) signif.floor(-121.51,n))
# [1] -200.00 -130.00 -122.00 -121.60 -121.51
> sapply(1:5, function(n) signif.ceiling(-121.51,n))
# [1] -100.00 -120.00 -121.00 -121.50 -121.51