[BioC] How to filter genes according to the intensities of both channels for cDNA array?

Xiao Shi bioconductor.cn at gmail.com
Fri May 6 08:46:32 CEST 2005


Hi everybody,
I am working on 10 GPR files from a time course cDNA array
experiment.Aftersome normalization steps,i got the M and A for 10
arrays.And then i want to apply the quality control procedures to flag 
logratio(M) for those spots to Na whose inteinsities form either channel is 
lower than 200.what can i do?

The procedures i did are as follows:
library(marray)
file=c('2000009478.gpr','2000009471.gpr','2000009472.gpr','2000009473.gpr','
2000009479.gpr','2000009466.gpr','2000009467.gpr','2000009468.gpr','
2000009480.gpr','2000009481.gpr')
target<-read.marrayInfo(fname="e:/gpr/zpz/zpz-target.txt")
##Read probe related information
galinfo<-read.Galfile("GAL-13K.txt",path="e:/gpr/zpz")
##Read expression data
zraw=read.GenePix
(fnames=file,path="e:/gpr/zpz/results",layout=galinfo$layout,gnames=galinfo$gnemes,target=target)

##normalization on print tips
zraw.norm=maNormMain(zraw)

##Between slides normalization
zraw.normg=maNormScale(zraw.norm,norm="g")

##And then the genefilter procedures,
## I want to flag the logratio(M) to NA according to the intensity of both 
channels of a spot,and then select the genes if the 
### sbsolute logratio(|M|) is grater than one.What can i do,what's the code 
doing that?
Thanks in advance.
Shi Jiantao

	[[alternative HTML version deleted]]



More information about the Bioconductor mailing list