我使用genefilter作为步骤1对R进行差异表达分析。
我正在使用该pOverA功能以及在使用之前设置过滤器功能genefilter()（我已经查看了有关 Bioconductor 的教程）。
我的 miRNA Expression 数据没有NAs 或NaNs。它由规范化的日志数据组成，并0在其间散布有值。
我创建了一个表达式数据的 eset，然后：
ff1<-pOverA(0.2,-3.2,na.rm=TRUE)
ff2<-function(x) IQR(x) > 0.2
ff<-filterfun(ff1,ff2)
我现在运行genefilter( exprs(miREset), ff)并收到一条错误消息：
quantile.default(as.numeric(x), c(0.25, 0.75), na.rm = na.rm, 错误：如果 'na.rm' 为 FALSE，则不允许缺失值和 NaN 另外：警告消息：在quantile(as.numeric(x), c(0.25, 0.75), na.rm = na.rm, names = FALSE, : 强制引入的 NA
我花了很多时间试图理解为什么会这样。我对R相当陌生。