How to convert dataframe into ExpressionSetIllumina / summarized data with beadarray package?
1
0
Entering edit mode
salamandra ▴ 20
@salamandra-12825
Last seen 2.2 years ago
Portugal

In documentation (beadsummary.pdf) of beadarray package it explains how to get a ExpressionSetIllumina from data deposited on GEO. Something like:

library(GEOquery)
url <- "ftp://ftp.ncbi.nih.gov/pub/geo/DATA/SeriesMatrix/GSE33126/"
filenm <- "GSE33126_series_matrix.txt.gz"
if(!file.exists("GSE33126_series_matrix.txt.gz")) download.file(paste(url, filenm, sep=""), paste0(getwd(),'/',filenm))
gse <- getGEO(filename=filenm)
# Convert into ExpresssionSetIllumina:
summaryData <- as(gse, "ExpressionSetIllumina")
summaryData

However, this is already normalized data. It also explains how to get "ExpressionSetIllumina" from IDAT files, but the raw data I am trying to get is not idat file, as you see in suppl. files of following studies: GSE22247 GSE26677

I can get those supplementary txt files with:

dataRAW <- read.delim(paste0(txtpath), sep='\t', skip = 4, header = T)

How can I convert dataframe (dataRAW) into an ExpressionSetIllumina object ?

p.s: I need raw data to test different normalisations and latter do dif. gene expression analysis.

beadarray • 838 views
ADD COMMENT
3
Entering edit mode
Mark Dunning ★ 1.1k
@mark-dunning-3319
Last seen 12 months ago
Sheffield, Uk

I would probably start by reading these data using limma

data <- limma::read.ilmn("GSE22247_non-normalized_data.txt",expr = "SAMPLE")

and then following steps from the limma user guide for normalisation

ADD COMMENT

Login before adding your answer.

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