[BioC] Error in extracting genes from Goterms using GOStats package
James W. MacDonald
jmacdon at uw.edu
Fri Oct 25 16:16:42 CEST 2013
Hi Sandy,
You cannot do a GO test for unannotated reporters. In other words, if
you have a probeset that measures some unknown thing, then you cannot
use that probeset because there won't be any GO terms associated with
it.
What happens if you do
list1 <- list1[!is.na(list1)]
first?
Best,
Jim
On Friday, October 25, 2013 8:47:21 AM, Sandy wrote:
>
>
> Hi James,
>
> probeSetSummary2(hgOver) works well as long as there are no NA's in this
> particuar step:
>
> list1<- (unlist(mget(selectedgenes,ath1121501ACCNUM,)))
>
> In the presence of NA's, then probeSetSummary2(hgOver) gives me this warning
>
> Warning message:
> The vector of geneIds used to create the GOHyperGParamsobject was not a
> named vector.
> If you want to know the probesets that contributed to this result
> either use a named vector for gene Ids, or pass a vector of probeset IDs
> via sigProbesets.
>
>
> And then most GO terms go missing in the result.
>
> What do I do in such a case?
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
--
James W. MacDonald, M.S.
Biostatistician
University of Washington
Environmental and Occupational Health Sciences
4225 Roosevelt Way NE, # 100
Seattle WA 98105-6099
More information about the Bioconductor
mailing list