[BioC] removing probes from affybatch error
James MacDonald
jmacdon at med.umich.edu
Wed Mar 31 19:31:34 CEST 2010
The AnnotationDbi package contains a function makeProbePackage() that
you can use to make probe packages.
Best,
Jim
Jenny Drnevich wrote:
> Hi George,
>
> I've never made a probe package either? Did you try searching the
> BioC archives for "make probe package"?
>
> http://dir.gmane.org/gmane.science.biology.informatics.conductor
>
> Also, you forgot to post to the BioC list as well!
>
> Cheers,
> Jenny
>
> At 08:40 AM 3/31/2010, George Roberts III wrote:
>> Hi Jenny,
>>
>> I installed the .CDF package for my chips but still don't know how
>> to make a cdfprobe package and couldn't find any further info.
>>
>> Thanks for your help,
>> George
>>
>> On Sat, Jan 23, 2010 at 11:58 PM, George Roberts
>> <<mailto:groberts0429 at gmail.com>groberts0429 at gmail.com> wrote:
>> Hi Jenny,
>>
>> Sorry for not posting in the BioC list, I will in the future; I
>> found the URL where you posted the previous email but don't know how
>> to reply to that thread so I replied to this email instead.
>>
>> I spoke with the designers of the custom chip to find some insight
>> into the "NA" issue. All of the "NA" values were traced to a series
>> of probes that consist of only a single PM probe for each gene. They
>> were analyzed with GCOS and none gave "NA"; all gave a number. The
>> individual who performed the first analysis had no further insight
>> into whey R might have behaved this way. I want to exclude these
>> probes from my analysis in any event.
>>
>> Thanks for pointing out that this error had to do with not
>> installing the .CDF package. I had loaded the .CDF environment in
>> the belief that it would be enough. Installing the .CDF package and
>> .CDFprobe package will be my next step.
>>
>> Thank you very much for your help,
>> George
>>
>>
>> On Tue, Jan 19, 2010 at 9:38 AM, Jenny Drnevich
>> <<mailto:drnevich at illinois.edu>drnevich at illinois.edu> wrote:
>> Hi George,
>>
>> Sorry I haven't been able to respond - I've been extremely busy.
>> It's actually better to post questions like this to the BioC list
>> instead of to one poster directly - you might get a more timely
>> response, and it's good to have the response in the archives so
>> other people can access it. Your error message is pretty simple -
>> you don't have package mpauta520274cdfprobe. I'm guessing this is a
>> custom Affy array, so in addition to the mpaut1a520274fcdf package,
>> you'll need a mpaut1a520274fprobe package to be able to use the
>> RemoveProbes function "as is" and other functions like gcrma(). I
>> forget off the top of my head what you need to do to make the probe
>> package, so perhaps others will comment. If you can't make the probe
>> package and only want to remove entire probe sets, then the code of
>> RemoveProbes might be able to be modified to work without the probe
>> package. Unfortunately, I won't have time for a couple of weeks to
>> look at it. Feel free to hack it all you want if you can! Sorry I
>> can't be of more help right now...
>>
>> Jenny
>>
>> P.S. What's going on with your scans that you are getting "NA"
>> values for individual probes? That's very odd!
>>
>>
>> At 05:48 PM 1/12/2010, George Roberts wrote:
>>> Dr. Drnevich,
>>>
>>> I would like to use your code
>>> <https://stat.ethz.ch/pipermail/bioconductor/2008-September/024296.html>https://stat.ethz.ch/pipermail/bioconductor/2008-September/024296.html
>>> to remove some probesets that are causing "NA" values resulting in
>>> errors. My CDF is "MPAUT1a520274F.CDF". I get the following error messages:
>>>
>>>> RemoveProbes(listOutProbes=NULL, listOutProbeSets="
>>> + 1415670_at
>>> + 1415671_at
>>> + 1415672_at
>>> + ", cleancdf)
>>> Loading required package: mpaut1a520274fcdfcdf
>>> Loading required package: mpaut1a520274fcdfprobe
>>> Error in get(probepackagename) :
>>> object 'mpaut1a520274fcdfprobe' not found
>>> In addition: Warning messages:
>>> 1: In library(package, lib.loc = lib.loc, character.only = TRUE,
>>> logical.return = TRUE, :
>>> there is no package called 'mpaut1a520274fcdfcdf'
>>> 2: In library(package, lib.loc = lib.loc, character.only = TRUE,
>>> logical.return = TRUE, :
>>> there is no package called 'mpaut1a520274fcdfprobe'
>>>
>>> Thanks,
>>> George Roberts
>>> Wayne State University
>> 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: <mailto:drnevich at illinois.edu>drnevich at illinois.edu
>>
>>
>
> 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
> [[alternative HTML version deleted]]
>
> _______________________________________________
> 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
--
James W. MacDonald, M.S.
Biostatistician
Hildebrandt Lab
8220D MSRB III
1150 W. Medical Center Drive
Ann Arbor MI 48109-5646
734-936-8662
**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues
More information about the Bioconductor
mailing list