This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Help! Select Differentially Expressed Genes Using Bioconductor

I am trying to use filtered_R to select differential expressed genes from Filtall (a non-specific filtering subset ) by BH. An error occurred. How should I do this correctly? Thanks a lot.

# Filtall
ExpressionSet (storageMode: lockedEnvironment)
assayData: 379 features, 111 samples 

rawp <- esApply(Filtall, MARGIN = 1, FUN = function(x) {
t.test(x[BCR], x[!NEG], alternative="two.sided",mu=0, paired = FALSE, var.equal = FALSE)$p.value
})

sum(rawp<=0.05)
[1] 312

fwer <- mt.rawp2adjp(rawp, c("BH", "BY"))
mt.reject(fwer$adjp, alpha = 0.05)$r

 rawp  BH  BY 
 0.05  312 310 286

theta = seq(from=0, to=0.5, by=0.1)

rejBH = filtered_R(alpha=0.05, filter= Filtall, test=rawp, theta=theta, method=BH)

Error in (1 - h) * qs[i] : non-numeric argument to binary operator
In addition: Warning messages:
1: In is.na(x) : is.na() applied to non-(list or vector) of type 'S4'
2: In is.na(x) : is.na() applied to non-(list or vector) of type 'S4'

There is not enough information in this question for us to help you. What is Filtall? Where is this filtered_R function coming from? Is this from a package or your own function?

f1<-pOverA(p=0.25,A=100, na.rm=TRUE) f2<-cv(a = 0.7, b = 10) ff<-filterfun(f1,f2) filtall <- genefilter(2^exprs(all), ff) sum(filtall) [1] 379 Filtall<-all[filtall,] Filtall ExpressionSet (storageMode: lockedEnvironment) assayData: 379 features, 111 samples element names: exprs protocolData: none phenoData sampleNames: 01005 01010 ... LAL4 (111 total) varLabels: cod diagnosis ... date last seen (21 total) varMetadata: labelDescription featureData: none experimentData: use 'experimentData(object)' pubMedIds: 14684422 16243790 Annotation: hgu95av2

So I plan to use filtered_R (from genefilter package in Bioconductor ) to select differential genes with BH<=0.5. filtered_R(alpha, filter, test, theta, data, method = "none") #http://127.0.0.1:19971/library/genefilter/html/filtered_p.html However, it doesn't work. please see the error information in the original post.

Please provide a reproducible example and the output of sessionInfo().

0 answers

No answers yet.

Log in to answer this question.