在miRNA表达数据上使用genefilter

Question

我使用genefilter作为步骤1对R进行差异表达分析。

我正在使用pOverA函数以及在使用genefilter()之前设置过滤函数（我已经查看了Bioconductor上的教程）。

我的miRNA表达数据没有NA s或NaN s。它由标准化的日志数据组成，其中包含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 In，则不允许缺少值和NaN   另外：警告信息：在分位数（as.numeric（x），c（0.25,0.75），   na.rm = na.rm，names = FALSE，：由强制引入的NAs

我花了很多时间试图理解为什么会这样。我是R的新手。