Gene dataset turned into NA after I tried to match probe sets to genes
0
0
Entering edit mode
7.9 years ago

I am using the GitHub tutorial https://github.com/icnn/Microarray-Tutorials/wiki/Affymetrix#4-batch-correction

Every code was perfect and I got exact expected results as in the tutorial. However, on the step 8. Annotating Probes, I ran:

to_keep <- is.na(geneDat$ensembl_gene_id) == FALSE)

geneDat <- geneDat[to_keep,]

datAll <- datAll[to_keep,]

my whole data "geneDat" turned into NA https://postimg.org/image/8bm09qizr/ instead of the example provided in the tutorial of what it should be.

I don't have enough experience to know whether the tutorial missed something or something is wrong on my end. I followed every code exactly, so can't understand what's the problem.

Thanks for help!

RNA-Seq affymetrix microarray • 1.3k views
ADD COMMENT
1
Entering edit mode

Most likely it's a problem with mapping the probes to EnsEMBL genes. If all values are NA, I would suspect that you didn't have access to EnsEMBL data. Check the biomaRt user guide especially if you're behind a proxy.

ADD REPLY
0
Entering edit mode

Thank you, I fixed the problem. I believe that was the problem you described.

ADD REPLY

Login before adding your answer.

Traffic: 2015 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6