Normalization of Illumina Microarray Data
1
0
Entering edit mode
@paulelliot1994-12566
Last seen 6.3 years ago

If it's not too much Trouble, can someone tell me as to how do I, read in the txt file in R studio? The Raw data was downloaded from GEO. The data accession number is GSE48624.

Moreover, how do i Differentiate the control sample from the normal sample since there are only txt files available in the raw data compressed file!

Thank you very much!

 

illumina microarray normalization • 2.7k views
ADD COMMENT
2
Entering edit mode
Mark Dunning ★ 1.1k
@mark-dunning-3319
Last seen 20 months ago
Sheffield, Uk

The "raw data" for this dataset seem to be the annotation files for the arrays, and not the raw intensities at all. So not helpful for analysis.

However, you can use the GEOquery package to import the normalised intensities

library(GEOquery)
remotefile <- "ftp://ftp.ncbi.nlm.nih.gov/geo/series/GSE48nnn/GSE48624/matrix/GSE48624_series_matrix.txt.gz"

destfile <- "GSE48624_series_matrix.txt.gz"
if(!file.exists(destfile)) download.file(remotefile, destfile)
mydata<- getGEO(filename=destfile)

Phenotypical data can be obtained using the pData function

pd <- pData(mydata)

pd

The vector pd$source_name_ch1 will tell you whether each sample was in the before or after group

pd$source_name_ch1

The function exprs will extract the normalised values

exprs(mydata)

Hope this helps,

ADD COMMENT
0
Entering edit mode

Thank You so much Mark, For helping me out. God bless!

ADD REPLY
0
Entering edit mode

If it's not too much trouble, can you tell me as to how to perform differential expression analysis on the same and get the list of upregulated and downregulated genes using the same Dataset(GSE48624)Thank You !

ADD REPLY

Login before adding your answer.

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