Reference: https://stat.ethz.ch/pipermail/bioconductor/2006-September/014242.html Hi, To calculate gcrma values of HG_U133_Plus_2 arrays, I am using Jenny Drnevich's code (see the link above) (Thanks to Jim who provided that link to me). This code is to mask either probesets or individual probes while calculating gcrma values. I have successfully run the code to mask some of the probe sets. Now, I want to mask individual probes. For example, I want to calculate gcrma value of probe set 117_at without a probe 60:815 and the probe set 121_at excluding three of its probes 801:783, 228:969, 268:273 (the numbers are the probes' x:y positions). According to the code (in the above link), to call function RemoveProbes(), I need to provide list of probe to be masked but it is not mentioned in which format one should provide probe list. I used the following code that is not working: probe_list=c("117_at:60:815","121_at:801:783", "121_at:228:969", "121_at:268:273") cleancdf=cleancdfname("HG_U133_Plus_2", addcdf=FALSE) ResetEnvir(cleancdf) RemoveProbes(listOutProbes=probe_list, listOutProbeSets=NULL, cleancdf) I am getting error while masking probe. If anybody used this code to mask individual probes, please let me know the format of the probe list. Thanks for your help in advance. Regards, Farhat

Hi Memon, RemoveProbes() requires the probe name, in the format of 117_at1 or 117_at9, not the x,y coordinates. I'm not sure there's an easy way to get from x,y coordinates to probe names, but some combination of indexProbes() and indices2xy() should work. You'll need an AffyBatch object of the raw data to test it, called "rawdata" in my code below. > x <- unlist(indexProbes(rawdata,genenames="117_at")) > x # this will give you the indices of the pm probes for probe set 117_at along with their names > y <- indices2xy(x, abatch=rawdata) > y # this will give you the x,y coordinates along with the correct probe names HTH, Jenny At 08:56 AM 1/18/2010, Memon, Farhat N wrote: >Reference: >https://stat.ethz.ch/pipermail/bioconductor/2006-September/014242.htm l > >Hi, > >To calculate gcrma values of HG_U133_Plus_2 arrays, I am using Jenny >Drnevich's code (see the link above) (Thanks to Jim who provided >that link to me). This code is to mask either probesets or >individual probes while calculating gcrma values. I have >successfully run the code to mask some of the probe sets. Now, I >want to mask individual probes. > >For example, I want to calculate gcrma value of probe set 117_at >without a probe 60:815 and the probe set 121_at excluding three of >its probes 801:783, 228:969, 268:273 (the numbers are the >probes' x:y positions). > >According to the code (in the above link), to call function >RemoveProbes(), I need to provide list of probe to be masked but it >is not mentioned in which format one should provide probe list. > >I used the following code that is not working: > >probe_list=c("117_at:60:815","121_at:801:783", "121_at:228:969", >"121_at:268:273") > >cleancdf=cleancdfname("HG_U133_Plus_2", addcdf=FALSE) >ResetEnvir(cleancdf) >RemoveProbes(listOutProbes=probe_list, listOutProbeSets=NULL, cleancdf) > >I am getting error while masking probe. If anybody used this code to >mask individual probes, please let me know the format of the probe list. > >Thanks for your help in advance. > >Regards, >Farhat > >_______________________________________________ >Bioconductor mailing list >Bioconductor at stat.math.ethz.ch >https://stat.ethz.ch/mailman/listinfo/bioconductor >Search the archives: >http://news.gmane.org/gmane.science.biology.informatics.conductor Jenny Drnevich, Ph.D. Functional Genomics Bioinformatics Specialist W.M. Keck Center for Comparative and Functional Genomics Roy J. Carver Biotechnology Center University of Illinois, Urbana-Champaign 330 ERML 1201 W. Gregory Dr. Urbana, IL 61801 USA ph: 217-244-7355 fax: 217-265-5066 e-mail: drnevich at illinois.edu