我有一个用于实验室分析几个参数的数据框架。许多值低于检测极限并且被报告为(<检测极限;例如<0.005)
数据帧如下
varA varB varC varD varE varF varG
0.024 0.69 . . . . <0.01
0.012 0.23 0.26
<0.005 0.54 0.28
0.017 0.62 0.30
0.016 <0.10 0.31
0.018 0.34 0.24
0.020 0.44 <0.01
<0.005 0.52 0.24
0.080 <0.10 0.25
我想使用这个公式转换所有低于检测极限的值
<detection limit ====> detection limit/sqare root (2)
我为varA做了这件事,如下
df$varA1<-as.character(df$varA)
df$varA1<-ifelse(df$varA1=="<0.005", 0.005/sqrt(2), df$varA1)
df$varA2<-as.numeric(df$varA1)
我可以对其他列和低于检测极限的其他值执行同样的操作。然而,考虑到我的数据很大,这将需要时间。
如何快速更改不同列中低于检测限的所有值,有什么建议吗?
dLimit <- function(v){
v <- as.character(v)
isLimit <- grepl("<", v)
v[isLimit] <- as.numeric(gsub("<(.*)", "\1", v[isLimit]))/sqrt(2)
as.numeric(v)
}
apply(df, 2, dLimit)