Question: probles with loess normalization in affy
1
gravatar for Lila M
17 months ago by
Lila M 860
UK
Lila M 860 wrote:

Hi all, I'm working in published array expression data sets (affymetrix u133b) . I've tried to do some loess normalization but I found some problems. Here the code:

library(affy)
celfiles <- list.files()
rawData <- ReadAffy(filenames =celfiles)
exprs(rawData)
boxplot(rawData)

normData <-normalize.AffyBatch.loess(rawData)
exprs(normData)
boxplot(normData)

So the code compiles properly but when I plot normData the boxplot doesn't look like normalized. Any suggestion? Thanks!

ADD COMMENTlink modified 17 months ago • written 17 months ago by Lila M 860

Can you show the plot?

ADD REPLYlink written 17 months ago by ATpoint40k

I don't know how to upload the boxplot , but this is a data subset

#raw data
GSM1525198.CEL.gz GSM1525199.CEL.gz GSM1525200.CEL.gz  GSM1525201.CEL.gz
1   403.3   211.0   417.5   310.5
2   16062.0 16321.8 17672.5 13178.8
3   393.8   216.3   372.3   323.8
4   16560.0 16944.3 17771.8 13375.0

#normalized data 
GSM1525198.CEL.gz GSM1525199.CEL.gz  GSM1525200.CEL.gz GSM1525201.CEL.gz
1   373.2181    299.3475    349.5832    282.4489
2   20177.9132  23678.5436  14341.0039  8911.0882
3   363.7140    306.7556    312.7128    294.3092
4   20843.4204  24557.0469  14425.2356  9033.1897
ADD REPLYlink modified 17 months ago • written 17 months ago by Lila M 860
2

How to add images to a Biostars post

ADD REPLYlink written 17 months ago by genomax91k
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: 1631 users visited in the last hour