Boxplot from GEO2R vs my own script

Hi all!

I have downloaded the .cel files to GSE79326. The boxplot I draw is different from that drawn using GEO2R. I wondered why. My code is:

library (oligo)
afbatch=read.celfiles(list.celfiles()) 
dataset=exprs (rma (afbatch, normalize=F)) 
boxplot (dataset)

Thank you in advance.

r boxplot

• 2,173 views • link

updated 8.4 years ago by Sean Davis

• written 8.4 years ago by english.server

Don't you think if you show us two boxplots would be helpful to say what's happening there?

• 0 views • link

written 8.4 years ago by venu

1 answer

The code for a GEO2R analysis is available from the GEO2R site. You might take a look at that code to see what GEO2R is doing that might be different.

One difference is that you are starting with .CEL files. GEO2R uses data as supplied by the submitter.

• 0 views • link

written 8.4 years ago by Sean Davis

More posts like this

problem with reading cel files from [HuGene-2_0-st] Affymetrix Human Gene 2.0 ST Array [HuGene20stv…
written 1.7 years ago by drwaqasakbar •

Hi, I am trying to analyze a microarray dataset with platform “[HuGene-2_0-st] Affymetrix Human Gene 2.0 ST Array [HuGene20stv1_Hs_ENTREZG_17.0.0]”. I have read almost all papers …
using RMA normalized microarray in Limma
written 2.7 years ago by Fluke •

Hi, i downloaded raw.tar file containing .CEL microarray data from GEO database. and then i use affy() package in R to extracted .CEL file and …
affy expression array data
written 2.9 years ago by liuyibin2005 •

Hi, everyone I use the oligo package to read cel files, and use the rma function to transformed them into expression matrix, but what I …
How to find upregulated and downregulated genes for my microarray data with limma
written 6.7 years ago by kai_bio •

I have annotated the .cel files and got the log2 transformation values of 12 gastric effected patients and 15 normal (control) samples patients data and …
How to do Clariom D Human analysis for 3 groups in R?
written 7.2 years ago by jiamacro •

Hi! I am a freshman who wanna analyze my clariom d human data. I have 3 groups and total 14 .CEL files and I wanna …
Affymetrix miRNA 4.0 Normalization
written 7.4 years ago by gracie •

I am analyzing a GeneChip miRNA 4.0 Array dataset. here is my code; library(oligo) library(pd.mirna.4.0) celFiles <- list.celfiles( full.names=TRUE) rawData <- read.celfiles(celFiles, pkgname="pd.mirna.4.0") eset <- …
Error in match.arg in boxplot - oligo package
written 9.3 years ago by salamandra

When running the following code: library(oligo) library(GEOquery) getGEOSuppFiles("GSE38792") list.files("GSE38792") untar("GSE38792/GSE38792_RAW.tar", exdir = "GSE38792/CEL") list.files("GSE38792/CEL") celfiles <- list.files("GSE38792/CEL", full=T) celfiles # read cel files rawData <- …
missing probeID in affymetrix using oligo package
written 11.2 years ago by benhrif.oussama •

Hello, I am working on an affymetrix project. My data are cel files, and I am using oligo package. My script is ```r celFiles<-list.celfiles(dir,full.names=TRUE) rawData<-read.celfiles(celFiles) …
missing probeID in affymetrix using oligo package
written 11.2 years ago by benhrif.oussama •

Hello, I am working on an affymetrix project. My data are cel files, and I am using oligo package. My script is ```r celFiles<-list.celfiles(dir,full.names=TRUE) rawData<-read.celfiles(celFiles) …
Getting Probe Ids For Expression Matrix Before Summarisation (In R If Possible)
written 14.7 years ago by Jimbo

Hi guys Anyone know how to get the probe set ids for the expression matrix row names BEFORE summarisation (so like 140292811, 140292812 etc.). When …