Re: read.GenePix() function
1
0
Entering edit mode
Hua Weng ▴ 60
@hua-weng-460
Last seen 10.2 years ago
Hi, Jean and everybody: In read.GenePix() function, is there any way I can specify Background Standard Deviation for Cy5 and Cy3( such as B647 SD and B546 SD) in GPR files? user.raw <- read.GenePix(fnames=NULL, path=".", name.Gf="F546 Median", name.Gb="B546 Median", name.Rf="F647 Median", name.Rb="B647 Median", name.W="Flags",layout=user.layout, gnames=user.gnames, targets=user.samples, skip=30) I want to specify these two columns is because I want to write a function to filter spots whose RF and GF are less than RB+2*RBSD and GB+2*GBSD. Do you know an approach I can do that? Thank you very much, Hua [[alternative HTML version deleted]]
• 1.1k views
ADD COMMENT
0
Entering edit mode
@jean-yee-hwa-yang-104
Last seen 10.2 years ago
Hi Hua, There is no ready made function for this at the moment. I will suggest you modified version of the function "read.marrayRaw", and read the extra two columns of values per slide. If you like further help with that, let me know. Cheers Jean > I want to specify these two columns is because I want to write a > function to filter spots whose RF and GF are less than RB+2*RBSD and > GB+2*GBSD. Do you know an approach I can do that?
ADD COMMENT
0
Entering edit mode
This can be done by providing a user-specified function to the wt.fun argument of read.maimages() in the limma package. Just use something like mywtfun <- function(dat) as.numeric( (dat$RF>dat$RB+2*dat$RBSD) & (dat$GF>dat$GB+2*dat$GBSD)) RG <- read.maimages(files, columns=list(Gf="GF,Rf="RF",Gb="GB",Rb="RB"),wt.fun=mywtfun) Gordon > > I want to specify these two columns is because I want to write a > > function to filter spots whose RF and GF are less than RB+2*RBSD and > > GB+2*GBSD. Do you know an approach I can do that?
ADD REPLY

Login before adding your answer.

Traffic: 645 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6