R中的重要数字四舍五入

Question

所以基本上我正在做一个物理实验，在我的表中我希望我的数据四舍五入到与错误相同的精度，这四舍五入为1 sig fig。

例如，如果我有以下内容：

angle <- c(4, 4.1, 4.2)
angle.error <- c(0.024, 0.3, 0.113)
data <- data.frame(angle,angle.error)

想要最终得到的是将角度误差四舍五入为1 sig fig，并将角度四舍五入到相应的小数位数作为角度误差，给出

angle <- c(4.00, 4.1, 4.2)
angle.error <- c(0.02,0.3,0.1)
data <- data.frame(angle, angle.error)

希望这有道理！ 这是我们被教导提供数据的标准方式，因此我很惊讶我找到一种正确的方法来实现这一目标是多么困难。 任何贡献都会很棒!!

编辑：

如果我想将其转换为UDR，其中数据在第5列，而数据框第6列中的错误，我写道：

Conv2 <- function(x) {
x[6] <- signif(x[6],1)
exp <- floor(log10(x[6]))
man <- x[6]/10^exp         
x[5] <- round(x[5], -exp)   # Error associated with this line I think
sapply(seq_along(x[5]), function(i) format(x[i,5], nsmall=-exp[i]))

return(x)
}

当我实现它时，我收到错误'FUN中的错误（X [[i]]，...）：数学函数的非数字参数'

Answer 1

这可以用signif()来完成。 与round()到多个小数位的round()相反， signif() round()到特定数量的重要位置。

x <- c(1,1.0001,0.00001)
round(x,1) # rounds to the specified number of decimal places
signif(x,1) # rounds to the specified number of significant digits

让我们将此应用于您的案例。 请注意，signif在两个参数上都是矢量化的。

> signif(c(1.111,2.222,3.333),c(1,2,3))
[1] 1.00 2.20 3.33

这意味着您的情况如下：

angle.error <- c(0.024, 0.3, 0.113)
cor <- signif(angle.error,1)

圆角越圆越困难。 我们需要将角度舍入到与​​angle.error_corr相同的小数位数。 因此，我们将提取该数字的尾数和指数。

exp <- floor(log10(cor)) 
man <- cor/10^exp

现在可以使用指数来圆角。

angle <- c(4, 4.1, 4.2)
angle_cor <- round(angle, -exp)

但是，这并没有给出预期的结果，因为R表示最精确的双精度数字和丢失尾随零的数字。 我们可以用格式的nsmall参数解决这个问题，遗憾的是它没有矢量化。

sapply(seq_along(angle), function(i) format(angle[i], nsmall = -exp[i]))

哪个应该给你你想要的答案。