[BioC] Sequence data and filtering

[hrishikesh deshmukh]

Dear Bioconductorians,

My sincere apologies if this email is duplicated.I
have searched the mailing list and could not find a
way to solve my problem: Find out the probe intensity
(PM adn MM) for HG_U95Av2 which fall in (200 to 20000
units) range and extract their respective PM
sequences. I have been able to filter the probes but
how do i extract the sequence data?! I am including
the small script which i have written for filtering,
now how do i extract sequence information.

Thanks in advance.

Hrishi


library(affy)
Data <- ReadAffy()
pmi <- pm(Data)
mmi <- mm(Data)
pbn <- probeNames(Data)
 
rng.pmi <- apply(pmi,1,range)
rng.mmi <- apply(mmi,1,range)


out.of.boundspm <- ((rng.pmi[1,] >=200) & (rng.pmi[1,]
<=20000))
out.of.boundsmm <- ((rng.mmi[1,] >=200) & (rng.mmi[1,]
<=20000))

in.bounds <- (out.of.boundspm & out.of.boundsmm) 

length(pmi[,1])
ac1 <- 1:201800
ac2 <- ac1[in.bounds]


h.pbn <- pbn[ac2]
h.pmi <- pmi[ac2,]
h.mmi <- mmi[ac2,]

normal.data <- cbind(h.pbn,h.pmi,h.mmi)
names(normal.data) <-
c('probesetnames','PMintensity','MMintensity')
write.table(normal.data, file="normaldata.txt",quote
=F,row.names=F,col.names=T,sep = " ")


		
__________________________________ 

Send holiday email and support a worthy cause. Do good.