Egambia: Gene expression in TB patients and Healthy controls

Description Details Examples

Description

Gene expression in TB patients and Healthy controls

Details

This data set has been constructed from the gene expression data set accessible in the Gene Expression Omnibus (GEO) at the accession number GSE28623. Ten healthy donors (NID, non-infected donors) and 10 tubercolosis patients (TB) have been randomly selected from the full data set, and top 25 genes with the highest IQR have been selected for further analysis. Genes without an Entrez gene (EG) identifier have likewise been omitted.

The Egambia object is a data frame. The first three columns are the gene symbol, gene name and Entrez gene (EG) identifier. The remaining columns correspond to the gene expression data.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
## Not run: 
# The data set has been generated as follows:
# get the data set from GEO
library( GEOquery )
gambia <- getGEO( "GSE28623" )[[1]]

# Convert to limma and normalize
library( limma )
e <- new( "EListRaw", list( E= exprs( gambia ), genes= fData( gambia ), targets= pData( gambia ) ) )
e.bg <- backgroundCorrect( e, method= "normexp" )
en <- normalizeBetweenArrays( e.bg, method= "q" )
en <- avereps( en, ID= en$genes$NAME )
en <- en[ en$genes$CONTROL_TYPE == "FALSE", ]
en$targets$group <- factor( gsub( " whole blood RNA *", "", en$targets$description ) )

# Fill in Entrez Gene IDs
library( org.Hs.eg.db )
en$genes$EG <- ""
sel <- en$genes$REFSEQ %in% ls( org.Hs.egREFSEQ2EG )
en$genes$EG[sel] <- mget( as.character( en$genes$REFSEQ[sel] ), org.Hs.egREFSEQ2EG )

# Filter by IQR and missing EG's
iqrs <- apply( en$E, 1, IQR )
en2 <- en[ iqrs > quantile( iqrs, 0.75 ) & en$genes$EG != "", ]

# Select 10 random samples from NID and TB groups
en2 <- en2[ , c( sample( which( en2$targets$group == "NID" ), 10 ), 
                 sample( which( en2$targets$group == "TB" ), 10 ) ) ]
colnames( en2$E ) <- en2$targets$group
Egambia <- cbind( en2$genes[ , c( "GENE_SYMBOL", "GENE_NAME", "EG" ) ], en2$E )

## End(Not run)

tmod documentation built on March 18, 2018, 1:59 p.m.