Question: Normalization agilent microarray data?
0
gravatar for mathavanbioinfo
8 months ago by
India
mathavanbioinfo50 wrote:

Hello All I am analyzing Agilent microarray data to study the infection condition, I used networkanalyst.ca tool for analysis. The normalization was done using Variance Stabilizing Normalization. Herewith attached image produced after normalization. Please comment is this correctly normalized.

Than you enter image description here

microarray normalization • 466 views
ADD COMMENTlink modified 8 months ago • written 8 months ago by mathavanbioinfo50

We are [mostly] analysts here. Your description of your data processing steps are quite ambiguous. Please explain exactly how you processed your data, providing code used, where applicable.

You appear to have used networkanalyst.ca - is that correct?

ADD REPLYlink modified 8 months ago • written 8 months ago by Kevin Blighe60k

Yes, I have used networkanalyst.ca tool, I have used microarray Series matrix file, that contain the gene expression value of the sample. enter image description here

ADD REPLYlink written 8 months ago by mathavanbioinfo50

This image represents before normalization

ADD REPLYlink written 8 months ago by mathavanbioinfo50

Okay, it is just not common to use networkanalyst.ca to process microarray data.

You can also obtain normalised data by simply doing (in R):

library(Biobase)
library(GEOquery)

# load series and platform data from GEO
gset <- getGEO("GSE93861", GSEMatrix =TRUE, getGPL=FALSE)
if (length(gset) > 1) idx <- grep("GPL6480", attr(gset, "names")) else idx <- 1
gset <- gset[[idx]]

# set parameters and draw the plot
dev.new(width=4+dim(gset)[[2]]/5, height=6)
par(mar=c(2+round(max(nchar(sampleNames(gset)))/2),4,2,1))
title <- paste ("GSE93861", '/', annotation(gset), " selected samples", sep ='')
boxplot(exprs(gset), boxwex=0.7, notch=T, main=title, outline=FALSE, las=2)

hh

ADD REPLYlink written 8 months ago by Kevin Blighe60k

I have done it, can please share the code for DEG analysis

ADD REPLYlink written 8 months ago by mathavanbioinfo50
1

Are you asking Kevin for exact working code that will perform differential expression analysis on your data? Please do not expect volunteers here to do your work for you.

ADD REPLYlink written 8 months ago by RamRS27k
1

Thank you, I will do it

ADD REPLYlink written 8 months ago by mathavanbioinfo50

How to solve this

sigmatrix <- project.NormData$E[probeset.list,] Error in project.NormData$E[probeset.list, ] : invalid subscript type 'list'

ADD REPLYlink modified 8 months ago • written 8 months ago by mathavanbioinfo50

probeset.list should be a character or integer vector, not a list object. If character, it should contain rownames of your object being subset; if integer, these should correspond to row indices of your object being subset.

Please take some time to review subsetting in R.

ADD REPLYlink modified 8 months ago • written 8 months ago by Kevin Blighe60k

Sorry Kevin

I perused your answer here and I finished with something names gset with class of Expressionset but I don't know how to get normalized matrix for `` GEO accession through this code

Actually the ultimate goal is a survival Kaplan meier plot for this list of genes for this data sets

MAGEA4
SPP1
ALB
CCNE1
CXCL8
IL1R2
IL1B
IL11
OSM
GAL
CCL21
CCL4
LAMC2
CCL20
AQP9
ISG15
IL6
FOS
CXCL10
BIRC3
NR4A3
IFIT2
CLCF1
CYP1B1
GZMB
OR10J3
C5AR1
CCL2
ATG12
EGR1
TNFAIP3
CTSL
CD40
SLC7A7
ETV7
LTB
IL2RB
CEBPB
ITGB2
ZNF502
POLD4
DDX10
ADD REPLYlink modified 7 months ago • written 8 months ago by A3.8k

Hey, you can obtain the data like this:

library(Biobase)
library(GEOquery)

gset <- getGEO("GSE19417", GSEMatrix =TRUE, getGPL=FALSE)
if (length(gset) > 1) idx <- grep("GPL4372", attr(gset, "names")) else idx <- 1
gset <- gset[[idx]]
ADD REPLYlink written 8 months ago by Kevin Blighe60k

Thank you your code gave me

>

Is this matrix normalized now? I am not sure which accession are in rows :( and how I could gene symbol instead

ADD REPLYlink modified 7 months ago • written 8 months ago by A3.8k
1

Yes, it seems to be normalised.

They are Rosetta IDs - you can obtain an ID mapping here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL4372 (download the 'Annotation SOFT table...')

ADD REPLYlink modified 8 months ago • written 8 months ago by Kevin Blighe60k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 908 users visited in the last hour