[BioC] questions about agilent 44k array

xuejun qin xqin at chg.duhs.duke.edu
Wed Feb 23 20:10:17 CET 2005


Hi, dear all,


I need you guys help, I got array data in Agilent 44k format, it has  number
for columns  and rows but no number for grid,   there is a lot of missing
data, that means different numbers of  spot for each row, but all array have
same number of spots. I try to clean the data by package of marray. as
following

NTD.layout <- read.marrayLayout(fname = file.path(datadir,"agi.txt"),
ngr = 1, ngc = 1, nsr =43931, nsc = 1, skip = 1) ## ngr=1, ngc=1, which can
import all data, as there is different number  of spots for each row, I set
nsc=1

NTD.gnames <- read.marrayInfo(fname=file.path(datadir,"agi.txt"),
info.id=4:5,labels=4)
NTD.samples <- read.marrayInfo(fname=file.path(datadir,"tempgets1.txt"),
info.id=1:5,labels=3)
nameTab <- read.table("tempnames.txt",header=TRUE,sep="\t",quote="\""
fnames <- as.character(nameTab$FileName)
NTD.raw <- read.marrayRaw(fnames,path=datadir,skip=NULL,
name.Gf="gMeanSignal",name.Gb="gBGMedianSignal",name.Rf="rMeanSignal",name.R
b="rBGMedianSignal",layout=NTD.layout,
gnames=NTD.gnames,targets=NTD.samples)
Rcol <- maPalette(low="white",high="red")
Gcol <- maPalette(low="white",high="green")

image(NTD.raw[,1],yvar="maRb",col=Rcol,main=maLabels(maTargets(NTD.raw[1,j])
),contours=TRUE,bar=TRUE)

Error in image.default(1:nc, 1:nr, z, axes = FALSE, col = col, ...) :
        `z' must be a matrix

I know if I set number for ngr, ngc other than 1, it works, but that means I
lost some data.

I try to get the graph  to check  quality of data, such as Rb, Gb, Rf, Gf,
also I want to check those Rb>Rf , Gb>Gf, but marrayTwo doesn't exist
anymore .

Any suggestion about  Agilent 44k data will be welcome.

Thank You!!



More information about the Bioconductor mailing list