使用edgeR包进行差异表达分析。出现这样的问题：
data=data[rowMeans(data)>1,]nrow(data)[1] 3226group=c(rep("normal",50),rep("tumor",374))design <- model.matrix(~group）y <- DGEList(counts=data,group=group)
错误: NA counts not allowed
看了别人的回答：
These are array data, edgeR is for RNA-seq, just saying... For this error, well the error is clear, NAs are not allowed and your data have NAs. If this was RNA-seq this would not even be possible since counts cannot be absent but only zero. Check the limma package for microarray analysis.
First try to determine why the NAs are in your data. It would be mal-practice to simply replace them with a value without understanding why they exist in the first place.
Also, as ATpoint says, these are microarray samples, so, you should be using Limma.
我读取文件命令是rt=read.table("mRNA.symbol.txt", sep = "\t",header = T, check.names = F，fill=NA)
但是还是一头雾水，想请教详细的解决方法是什么，具体怎么做