Question: How to deal with this problem when normalize different microarrays platform?
1
gravatar for BioMed
3.3 years ago by
BioMed40
BioMed40 wrote:

Dear all,

I tried to use below code to normalize my data before using combat to deal with batch effects:

  • library(affy)
  • norm <- ReadAffy()
  • norm.rma <-rma(norm)

Unluckily, my data include 3 data sets, 2 Affymetrix Human Genome U133 Plus 2.0 Array and 1 Affymetrix Human Gene 1.0 ST Array. I received this error message "Error in affyio::read_abatch(filenames, rm.mask, rm.outliers, rm.extra,: Cel file D:////batchdemo/GSMxxxx.CEL.gz does not seem to have the correct dimensions".

I tried to use dchip instead but again, I failed. I even couldn't input the data to dchip.

Any advice and suggestions will be greatly appreciated.

Thank you.

ADD COMMENTlink modified 3.3 years ago by andrew.j.skelton735.8k • written 3.3 years ago by BioMed40
3
gravatar for andrew.j.skelton73
3.3 years ago by
London
andrew.j.skelton735.8k wrote:

Using arrays from different platforms is not as simple as just reading them in and normalising them together. You'll need to read in each dataset individually, normalise, use only probes that are on both platforms (check out nuIDs for this), then account for the variation across platforms in Limma (use an additive model), when performing differential gene expression. This is not a trivial task.

ADD COMMENTlink written 3.3 years ago by andrew.j.skelton735.8k

Thank you for your insightful comment.

ADD REPLYlink written 3.3 years ago by BioMed40

hii i am trying to normalize microarray data from two different platform. when i normalize GSE2361 indiviually using code: library(affy) norm <- ReadAffy() norm.rma <-rma(norm) Error in exprs(object)[index, , drop = FALSE] : subscript out of bounds give me suggestion.

ADD REPLYlink written 2.5 years ago by au.rinki.bio10
0
gravatar for 1234anjalianjali1234
3.3 years ago by
India
1234anjalianjali123430 wrote:

remove that particular .cel file from your data that gives you the error.

"Cel file D:////batchdemo/GSMxxxx.CEL.gz does not seem to have the correct dimensions"

ADD COMMENTlink written 3.3 years ago by 1234anjalianjali123430

Dear, the error comes from all data which belong to "Affymetrix Human Gene 1.0 ST Array".

ADD REPLYlink written 3.3 years ago by BioMed40
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: 1837 users visited in the last hour