我正在处理一些RNA-seq计数数据,我有大约60,000个包含基因名称的列和24个包含样本名称的行.当我做了一些基因名称转换时,我留下了一堆被命名的列NA.我知道R的处理NA方式与典型的列名不同,我的问题是如何删除这些列.这是我的数据的一个例子.
"Gene1" "Gene2" "Gene3" NA "Gene4"
1 10 11 12 10 15
2 13 12 50 40 30
3 34 23 23 21 22
我想它最终会像
"Gene1" "Gene2" "Gene3" "Gene4"
1 10 11 12 15
2 13 12 50 30
3 34 23 23 22
我确实找到了一些适用于其他人但不适合我的R代码
df<-df[, grep("^(NA)", names(df), value = TRUE, invert = TRUE)]