GSP
Quick Navigator

Search Site

Unix VPS
A - Starter
B - Basic
C - Preferred
D - Commercial
MPS - Dedicated
Previous VPSs
* Sign Up! *

Support
Contact Us
Online Help
Handbooks
Domain Status
Man Pages

FAQ
Virtual Servers
Pricing
Billing
Technical

Network
Facilities
Connectivity
Topology Map

Miscellaneous
Server Agreement
Year 2038
Credits
 

USA Flag

 

 

Man Pages


Manual Reference Pages  -  BIO::DB::SOAPEUTILITIES (3)

.ds Aq ’

NAME

Bio::DB::SoapEUtilities - Interface to the NCBI Entrez web service *BETA*

CONTENTS

SYNOPSIS



 use Bio::DB::SoapEUtilities;

 # factory construction

 my $fac = Bio::DB::SoapEUtilities->new()

 # executing a utility call

 #get an iteratable adaptor
 my $links = $fac->elink(
               -dbfrom => protein,
               -db => taxonomy,
               -id => \@protein_ids )->run(-auto_adapt => 1);

 # get a Bio::DB::SoapEUtilities::Result object
 my $result = $fac->esearch(
               -db => gene,
               -term => sonic and human)->run;

 # get the raw XML message
 my $xml = $fac->efetch(
             -db => gene,
             -id => \@gids )->run( -raw_xml => 1 );

 # change parameters
 my $new_result = $fac->efetch(
                   -db => gene,
                   -id => \@more_gids)->run;
 # reset parameters
 $fac->efetch->reset_parameters( -db => nucleotide,
                                 -id => $nucid );
 $result = $fac->efetch->run;
               
 # parsing and iterating the results

 $count = $result->count;
 @ids = $result->ids;

 while ( my $linkset = $links->next_link ) {
    $submitted = $linkset->submitted_id;
 }

 ($taxid) = $links->id_map($submitted_prot_id);
 $species_io = $fac->efetch( -db => taxonomy,
                             -id => $taxid )->run( -auto_adapt => 1);
 $species = $species_io->next_species;
 $linnaeus = $species->binomial;



DESCRIPTION

This module allows the user to query the NCBI Entrez database via its SOAP (Simple Object Access Protocol) web service (described at <http://eutils.ncbi.nlm.nih.gov/entrez/eutils/soap/v2.0/DOC/esoap_help.html>). The basic tools (einfo, esearch, elink, efetch, espell, epost) are available as methods off a SoapEUtilities factory object. Parameters for each tool can be queried, set and reset for each method through the Bio::ParameterBaseI standard calls (available_parameters(), set_parameters(), get_parameters(), reset_parameters()). Returned data can be retrieved, accessed and parsed in several ways, according to user preference. Adaptors and object iterators are available for efetch, egquery, elink, and esummary results.

USAGE

The SoapEU system has been designed to be as easy (few includes, available parameter facilities, reasonable defaults, intuitive aliases, built-in pipelines) or as complex (accessors for underlying low-level objects, all parameters accessible, custom hooks for builder objects, facilities for providing local copies of WSDLs) as the user requires or desires. (To the extent that it does not succeed in either direction, it is up to the user to report to the mailing list (FEEDBACK)!)

    Factory

To begin, make a factory:



 my $fac = Bio::DB::SoapEUtilities->new();



From the factory, utilities are called, parameters are set, and results or adaptors are retrieved.

If you have your own copy of the wsdl, use



 my $fac = Bio::Db::SoapEUtilities->new( -wsdl_file => $my_wsdl );



otherwise, the correct one will be obtained over the network (by Bio::DB::ESoap and friends).

    Utilities and parameters

To run any of the standard NCBI EUtilities (einfo, esearch, esummary, elink, egquery, epost, espell), call the desired utility from the factory. To use a utility, you must set its parameters and run it to get a result. TMTOWTDI:



 # verbose
 my $fetch = $fac->efetch();
 $fetch->set_parameters( -db => gene, -id => [828392, 790]);
 my $result = $fetch->run;

 # compact
 my $result = $fac->efetch(-db =>gene,-id => [828392,790])->run;

 # change ids
 $fac->efetch->set_parameters( -id => 470338 );
 $result = $fac->run;

 # another util
 $result = $fac->esearch(-db => protein, -term => BRCA and human)->run;

 # the utilities are kept separate
 %search_params = $fac->esearch->get_parameters;
 %fetch_params = $fac->efetch->get_parameters;
 $search_param{db}; # is protein
 $fetch_params{db}; # is gene



The factory is Bio::ParameterBaseI compliant: that means you can find out what you can set with



 @available_search = $fac->esearch->available_parameters;
 @available_egquery = $fac->egquery->available_parameters;



For more information on parameters, see <http://www.ncbi.nlm.nih.gov/entrez/query/static/eutils_help.html>.

    Results

The intermediate object for SoapEU query results is the Bio::DB::SoapEUtilities::Result. This is a BioPerly parsing of the SOAP message sent by NCBI when a query is run(). This can be very useful on it’s own, but most users will likely want to proceed directly to Adaptors, which take a Result and turn it into more intuitive/familiar BioPerl objects. Go there if the following details are too gory.

Results can be highly- or lowly-parsed, depending on the parameters passed to the factory run() method. To get the raw XML message with no parsing, do



 my $xml = $fac->$util->run(-raw_xml => 1); # $xml is a scalar string



To retrieve a Bio::DB::SoapEUtilities::Result object with limited parsing, but with accessors to the SOAP::SOM message (provided by SOAP::Lite), do



 my $result = $fac->$util->run(-no_parse => 1);
 my $som = $result->som;
 my $method_hash = $som->method; # etc...



To retrieve a Result object with message elements parsed into accessors, including count() and ids(), run without arguments:



 my $result = $fac->esearch->run()
 my $count = $result->count;
 my @Count = $result->Count; # counts for each member of
                             # the translation stack
 my @ids = $result->IdList_Id; # from automatic message parsing
 @ids = $result->ids; # a convenient alias



See Bio::DB::SoapEUtilities::Result for more, even gorier details.

    Adaptors

Adaptors convert EUtility Results into convenient objects, via a handle that usually provides an iterator, in the spirit of Bio::SeqIO. These are probably more useful than the Result to the typical user, and so you can retrieve them automatically by setting the run() parameter -auto_adapt = 1>.

In general, retrieve an adaptor like so:



 $adp = $fac->$util->run( -auto_adapt => 1 );
 # iterate...
 while ( my $obj = $adp->next_obj ) {
    # do stuff with $obj
 }



The adaptor itself occasionally possesses useful methods besides the iterator. The method next_obj always works, but a natural alias is also always available:



 $seqio = $fac->esearch->run( -auto_adapt => 1 );
 while ( my $seq = $seqio->next_seq ) {
    # do stuff with $seq
 }



In the above example, -auto_adapt = 1> also instructs the factory to perform an efetch based on the ids returned by the esearch (if any), so that the adaptor returned iterates over Bio::SeqI objects.

Here is a rundown of the different adaptor flavors:
o efetch, Fetch Adaptors, and BioPerl object iterators

The FetchAdaptor creates bona fide BioPerl objects. Currently, there are FetchAdaptor subclasses for sequence data (both Genbank and FASTA rettypes) and taxonomy data. The choice of FetchAdaptor is based on information in the result message, and should be transparent to the user.



 $seqio = $fac->efetch( -db =>nucleotide,
                        -id => \@ids,
                        -rettype => gb )->run( -auto_adapt => 1 );
 while (my $seq = $seqio->next_seq) {
    my $taxio = $fac->efetch(
        -db => taxonomy,
        -id => $seq->species->ncbi_taxid )->run(-auto_adapt => 1);
    my $tax = $taxio->next_species;
    unless ( $tax->TaxId == $seq->species->ncbi_taxid ) {
      print "more work for MAJ"
    }
 }



See the pod for the FetchAdaptor subclasses (e.g., Bio::DB::SoapEUtilities::FetchAdaptor::seq) for more detail.

o elink, the Link adaptor, and the linkset iterator

The LinkAdaptor manages LinkSets. In SoapEU, an elink call <B>alwaysB> preserves the correspondence between submitted and retrieved ids. The mapping between these can be accessed from the adaptor object directly as id_map()



 my $links = $fac->elink( -db => protein,
                          -dbfrom => nucleotide,
                          -id => \@nucids )->run( -auto_adapt => 1 );

 # maybe more than one associated id...
 my @prot_0 = $links->id_map( $nucids[0] );



Or iterate over the linksets:



 while ( my $ls = $links->next_linkset ) {
    @ids = $ls->ids;
    @submitted_ids = $ls->submitted_ids;
    # etc.
 }



o esummary, the DocSum adaptor, and the docsum iterator

The DocSumAdaptor manages docsums, the esummary return type. The objects returned by iterating with a DocSumAdaptor have accessors that let you obtain field information directly. Docsums contain lots of easy-to-forget fields; use item_names() to remind yourself.



 my $docs = $fac->esummary( -db => taxonomy,
                            -id => 527031 )->run(-auto_adapt=>1);
 # iterate over docsums
 while (my $d = $docs->next_docsum) {
    @available_items = $docsum->item_names;
    # any available item can be called as an accessor
    # from the docsum object...watch your case...
    $sci_name = $d->ScientificName;
    $taxid = $d->TaxId;
 }



o egquery, the GQuery adaptor, and the query iterator

The GQueryAdaptor manages global query items returned by calls to egquery, which identifies all NCBI databases containing hits for your query term. The databases actually containing hits can be retrieved directly from the adaptor with found_in_dbs:



 my $queries = $fac->egquery(
     -term => BRCA and human
    )->run(-auto_adapt=>1);
 my @dbs = $queries->found_in_dbs;



Retrieve the global query info returned for <B>anyB> database with query_by_db:



 my $prot_q = $queries->query_by_db(protein);
 if ($prot_q->count) {
    #do something
 }



Or iterate as usual:



 while ( my $q = $queries->next_query ) {
    if ($q->status eq Ok) {
      # do sth
    }
 }



    Web environments and query keys

To make large or complex requests for data, or to share queries, it may be helpful to use the NCBI WebEnv system to manage your queries. Each EUtility accepts the following parameters:



 -usehistory
 -WebEnv
 -QueryKey



for this purpose. These store the details of your queries serverside.

SoapEU attempts to make using these relatively straightforward. Use Result objects to obtain the correct parameters, and don’t forget -usehistory:



 my $result1 = $fac->esearch(
     -term => BRCA and human,
     -db => nucleotide,
     -usehistory => 1 )->run( -no_parse=>1 );

 my $result = $fac->esearch(
     -term => AND early onset,
     -QueryKey => $result1->query_key,
     -WebEnv => $result1->webenv )->run( -no_parse => 1 );

 my $result = $fac->esearch(
    -db => protein,
    -term => sonic,
    -usehistory => 1 )->run( -no_parse => 1 );

 # later (but not more than 8 hours later) that day...

 $result = $fac->esearch(
    -WebEnv => $result->webenv,
    -QueryKey => $result->query_key,
    -RetMax => 800 # get em all
    )->run; # note were parsing the result...
 @all_ids = $result->ids;



    Error checking

Two kinds of errors can ensue on an Entrez SOAP run. One is a SOAP fault, and the other is an error sent in non-faulted SOAP message from the server. The distinction is probably systematic, and I would welcome an explanation of it. To check for result errors, try something like:



 unless ( $result = $fac->$util->run ) {
    die $fac->errstr; # this will catch a SOAP fault
 }
 # a valid result object was returned, but it may carry an error
 if ($result->count == 0) {
    warn "No hits returned";
    if ($result->ERROR) {
      warn "Entrez error : ".$result->ERROR;
    }
 }



Error handling will be improved in the package eventually.

SEE ALSO

Bio::DB::EUtilities, Bio::DB::SoapEUtilities::Result, Bio::DB::ESoap.

FEEDBACK

    Mailing Lists

User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the Bioperl mailing list. Your participation is much appreciated.



  bioperl-l@bioperl.org                  - General discussion
http://bioperl.org/wiki/Mailing_lists  - About the mailing lists



    Support

Please direct usage questions or support issues to the mailing list:

bioperl-l@bioperl.org

rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible.

    Reporting Bugs

Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the web:



  http://redmine.open-bio.org/projects/bioperl/



AUTHOR - Mark A. Jensen

Email maj -at- fortinbras -dot- us

APPENDIX

The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _

    new



 Title   : new
 Usage   : my $eutil = new Bio::DB::SoapEUtilities();
 Function: Builds a new Bio::DB::SoapEUtilities object
 Returns : an instance of Bio::DB::SoapEUtilities
 Args    :



run()



 Title   : run
 Usage   : $fac->$eutility->run(@args)
 Function: Execute the EUtility
 Returns : true on success, false on fault or error
           (reason in errstr(), for more detail check the SOAP message
            in last_result() )
 Args    : named params appropriate to utility
           -auto_adapt => boolean ( return an iterator over results as
                                    appropriate to util if true)
           -raw_xml => boolean ( return raw xml result; no processing )
           Bio::DB::SoapEUtilities::Result constructor parms



    Useful Accessors

response_message()



 Title   : response_message
 Aliases : last_response, last_result
 Usage   : $som = $fac->response_message
 Function: get the last response message
 Returns : a SOAP::SOM object
 Args    : none



webenv()



 Title   : webenv
 Usage   :
 Function: contains WebEnv key referencing the session
           (set after run() )
 Returns : scalar
 Args    : none



errstr()



 Title   : errstr
 Usage   : $fac->errstr
 Function: get the last error, if any
 Example :
 Returns : value of errstr (a scalar)
 Args    : none



    Bio::ParameterBaseI compliance

available_parameters()



 Title   : available_parameters
 Usage   :
 Function: get available request parameters for calling
           utility
 Returns :
 Args    : -util => $desired_utility [optional, default is
           caller utility]



set_parameters()



 Title   : set_parameters
 Usage   :
 Function:
 Returns : none
 Args    : -util => $desired_utility [optional, default is
            caller utility],
           named utility arguments



get_parameters()



 Title   : get_parameters
 Usage   :
 Function:
 Returns : array of named parameters
 Args    : utility (scalar string) [optional]
           (default is caller utility)



reset_parameters()



 Title   : reset_parameters
 Usage   :
 Function:
 Returns : none
 Args    : -util => $desired_utility [optional, default is
            caller utility],
           named utility arguments



parameters_changed()



 Title   : parameters_changed
 Usage   :
 Function:
 Returns : boolean
 Args    : utility (scalar string) [optional]
           (default is caller utility)



_soap_facs()



 Title   : _soap_facs
 Usage   : $self->_soap_facs($util, $fac)
 Function: caches Bio::DB::ESoap factories for the
           eutils in use by this instance
 Example :
 Returns : Bio::DB::ESoap object
 Args    : $eutility, [optional on set] $esoap_factory_object



_caller_util()



 Title   : _caller_util
 Usage   : $self->_caller_util($newval)
 Function: the utility requested off the main SoapEUtilities
           object
 Example :
 Returns : value of _caller_util (a scalar string, a valid eutility)
 Args    : on set, new value (a scalar string [optional])



Search for    or go to Top of page |  Section 3 |  Main Index


perl v5.20.3 BIO::DB::SOAPEUTILITIES (3) 2016-04-05

Powered by GSP Visit the GSP FreeBSD Man Page Interface.
Output converted with manServer 1.07.