我有一个p值向量,我想对-log10
转换进行转换。我了解log2
和log10
等功能,但似乎无法找到执行负日志转换的任何内容。我发现的不起作用。这是我目前的代码:
x <- -log(dataTable$Pvalues, base = 10)
执行时,这会给我一个包含数值和无限(Inf
)值的向量。
如何在不获取无限值的情况下使用-log10
正确执行转换?
答案 0 :(得分:1)
不知道您的数据,但您可以这样做:
library(dplyr)
dataTable.1 <- dataTable %>%
filter(Pvalues != 0)
x <- -log(dataTable.1$Pvalues, base = 10)