How to normalize in R Illumina HumanHT-12 V3.0 expression beadchip data
4
0
Entering edit mode
kanacska ▴ 10
@kanacska-7375
Last seen 8.6 years ago
Hungary

Hi,

Im doing annotation on this article: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE38376

Ive read in the txt file in R studio, but i dont know how to normalize it.

Thanks,

 

Anna

illumina normalization microarray • 4.0k views
ADD COMMENT
0
Entering edit mode

If it's not too much trouble can you tell me as to how did you read in the txt file in R studio? Thank you!

ADD REPLY
0
Entering edit mode
Mark Dunning ★ 1.1k
@mark-dunning-3319
Last seen 13 months ago
Sheffield, Uk

Not sure which text file you are referring to, but if you just want to have a look at the data and start exploring it, the easiest thing would be to load it using GEOquery.

library(GEOquery)

data <- getGEO("GSE38376")[[1]]

ADD COMMENT
0
Entering edit mode
kanacska ▴ 10
@kanacska-7375
Last seen 8.6 years ago
Hungary

Thank you. This one on the page: GSE38376_non-normalized.txt

I want to do a QC and to get a gene list of diferencialy expressed genes

ADD COMMENT
0
Entering edit mode
kanacska ▴ 10
@kanacska-7375
Last seen 8.6 years ago
Hungary

but ive never used illumina just affymetrix

ADD COMMENT
0
Entering edit mode
Mark Dunning ★ 1.1k
@mark-dunning-3319
Last seen 13 months ago
Sheffield, Uk

I'm pretty sure that file isn't in the standard format exported by Illumina's software. Usually you expect to see columns with the prefix AVG_Signal. The following is a bit of a hack, but will give you an object to play around with;

 

library(limma)

idata <- read.ilmn("GSE38376_non-normalized.txt",probeid = "PROBE_ID",expr="SKBR")

 

The limma user guide has examples of processing Illumina data and doing differential expression with these kinds of objects.

ADD COMMENT
0
Entering edit mode

Dear Mark,

Sorry to bother you but if you can take a look at my post and can give me any suggestion, should be helpful. https://support.bioconductor.org/p/125225/

Thank you,

ADD REPLY

Login before adding your answer.

Traffic: 769 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