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'
• 668 views
•
link
0 answers
No answers yet.
Log in to answer this question.
There is not enough information in this question for us to help you. What is
Filtall? Where is thisfiltered_Rfunction coming from? Is this from a package or your own function?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.
Hi and welcome to Biostar. I agree with Steve Lianoglou that your question does not contain enough information for people to be able to help you. Please consider reading this post: How to ask Good Questions on Technical and Scientific Forums to help you improve this and future questions. Cheers
Please provide a reproducible example and the output of sessionInfo().