[Bio] / Sprout / Sprout.pm Repository:
ViewVC logotype

Diff of /Sprout/Sprout.pm

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.42, Tue Oct 18 06:58:09 2005 UTC revision 1.55, Mon Apr 3 18:33:58 2006 UTC
# Line 92  Line 92 
92                                                          # database type                                                          # database type
93                         dataDir      => $FIG_Config::sproutData,                         dataDir      => $FIG_Config::sproutData,
94                                                          # data file directory                                                          # data file directory
95                         xmlFileName  => "$FIG_Config::sproutData/SproutDBD.xml",                         xmlFileName  => "$FIG_Config::fig/SproutDBD.xml",
96                                                          # database definition file name                                                          # database definition file name
97                         userData     => "$FIG_Config::dbuser/$FIG_Config::dbpass",                         userData     => "$FIG_Config::dbuser/$FIG_Config::dbpass",
98                                                          # user name and password                                                          # user name and password
# Line 389  Line 389 
389    
390  =head3 LoadUpdate  =head3 LoadUpdate
391    
392  C<< my %stats = $sprout->LoadUpdate($truncateFlag, \@tableList); >>  C<< my $stats = $sprout->LoadUpdate($truncateFlag, \@tableList); >>
393    
394  Load updates to one or more database tables. This method enables the client to make changes to one  Load updates to one or more database tables. This method enables the client to make changes to one
395  or two tables without reloading the whole database. For each table, there must be a corresponding  or two tables without reloading the whole database. For each table, there must be a corresponding
# Line 740  Line 740 
740  should be of the form returned by L</featureLocation> when in a list context. In other words,  should be of the form returned by L</featureLocation> when in a list context. In other words,
741  each location is of the form I<contigID>C<_>I<begin>I<dir>I<end>.  each location is of the form I<contigID>C<_>I<begin>I<dir>I<end>.
742    
743    For example, the following would return the DNA sequence for contig C<83333.1:NC_000913>
744    between positions 1401 and 1532, inclusive.
745    
746        my $sequence = $sprout->DNASeq('83333.1:NC_000913_1401_1532');
747    
748  =over 4  =over 4
749    
750  =item locationList  =item locationList
751    
752  List of location specifiers, each in the form I<contigID>C<_>I<begin>I<dir>I<end> (see  List of location specifiers, each in the form I<contigID>C<_>I<begin>I<dir>I<len> or
753  L</FeatureLocation> for more about this format).  I<contigID>C<_>I<begin>C<_>I<end> (see L</FeatureLocation> for more about this format).
754    
755  =item RETURN  =item RETURN
756    
# Line 1153  Line 1158 
1158    
1159  =item RETURN  =item RETURN
1160    
1161  Returns a hash mapping the functional assignment IDs to user IDs.  Returns a hash mapping the user IDs to functional assignment IDs.
1162    
1163  =back  =back
1164    
# Line 1163  Line 1168 
1168      # Get the parameters.      # Get the parameters.
1169      my ($self, $featureID) = @_;      my ($self, $featureID) = @_;
1170      # Get all of the feature's annotations.      # Get all of the feature's annotations.
1171      my @query = $self->GetAll(['IsTargetOfAnnotation', 'Annotation'],      my @query = $self->GetAll(['IsTargetOfAnnotation', 'Annotation', 'MadeAnnotation'],
1172                              "IsTargetOfAnnotation(from-link) = ?",                              "IsTargetOfAnnotation(from-link) = ?",
1173                              [$featureID], ['Annotation(time)', 'Annotation(annotation)']);                              [$featureID], ['Annotation(time)', 'Annotation(annotation)',
1174                                               'MadeAnnotation(from-link)']);
1175      # Declare the return hash.      # Declare the return hash.
1176      my %retVal;      my %retVal;
     # Declare a hash for insuring we only make one assignment per user.  
     my %timeHash = ();  
1177      # Now we sort the assignments by timestamp in reverse.      # Now we sort the assignments by timestamp in reverse.
1178      my @sortedQuery = sort { -($a->[0] <=> $b->[0]) } @query;      my @sortedQuery = sort { -($a->[0] <=> $b->[0]) } @query;
1179      # Loop until we run out of annotations.      # Loop until we run out of annotations.
1180      for my $annotation (@sortedQuery) {      for my $annotation (@sortedQuery) {
1181          # Get the annotation fields.          # Get the annotation fields.
1182          my ($timeStamp, $text) = @{$annotation};          my ($timeStamp, $text, $user) = @{$annotation};
1183          # Check to see if this is a functional assignment.          # Check to see if this is a functional assignment.
1184          my ($user, $function) = _ParseAssignment($text);          my ($actualUser, $function) = _ParseAssignment($user, $text);
1185          if ($user && ! exists $timeHash{$user}) {          if ($actualUser && ! exists $retVal{$actualUser}) {
1186              # Here it is a functional assignment and there has been no              # Here it is a functional assignment and there has been no
1187              # previous assignment for this user, so we stuff it in the              # previous assignment for this user, so we stuff it in the
1188              # return hash.              # return hash.
1189              $retVal{$function} = $user;              $retVal{$actualUser} = $function;
             # Insure we don't assign to this user again.  
             $timeHash{$user} = 1;  
1190          }          }
1191      }      }
1192      # Return the hash of assignments found.      # Return the hash of assignments found.
# Line 1200  Line 1202 
1202  The functional assignment is handled differently depending on the type of feature. If  The functional assignment is handled differently depending on the type of feature. If
1203  the feature is identified by a FIG ID (begins with the string C<fig|>), then a functional  the feature is identified by a FIG ID (begins with the string C<fig|>), then a functional
1204  assignment is a type of annotation. The format of an assignment is described in  assignment is a type of annotation. The format of an assignment is described in
1205  L</ParseLocation>. Its worth noting that we cannot filter on the content of the  L</ParseAssignment>. Its worth noting that we cannot filter on the content of the
1206  annotation itself because it's a text field; however, this is not a big problem because  annotation itself because it's a text field; however, this is not a big problem because
1207  most features only have a small number of annotations.  most features only have a small number of annotations.
1208    
# Line 1262  Line 1264 
1264              }              }
1265          }          }
1266          # Build a query for all of the feature's annotations, sorted by date.          # Build a query for all of the feature's annotations, sorted by date.
1267          my $query = $self->Get(['IsTargetOfAnnotation', 'Annotation'],          my $query = $self->Get(['IsTargetOfAnnotation', 'Annotation', 'MadeAnnotation'],
1268                                 "IsTargetOfAnnotation(from-link) = ? ORDER BY Annotation(time) DESC",                                 "IsTargetOfAnnotation(from-link) = ? ORDER BY Annotation(time) DESC",
1269                                 [$featureID]);                                 [$featureID]);
1270          my $timeSelected = 0;          my $timeSelected = 0;
1271          # Loop until we run out of annotations.          # Loop until we run out of annotations.
1272          while (my $annotation = $query->Fetch()) {          while (my $annotation = $query->Fetch()) {
1273              # Get the annotation text.              # Get the annotation text.
1274              my ($text, $time) = $annotation->Values(['Annotation(annotation)','Annotation(time)']);              my ($text, $time, $user) = $annotation->Values(['Annotation(annotation)',
1275                                                         'Annotation(time)', 'MadeAnnotation(from-link)']);
1276              # Check to see if this is a functional assignment for a trusted user.              # Check to see if this is a functional assignment for a trusted user.
1277              my ($user, $function) = _ParseAssignment($text);              my ($actualUser, $function) = _ParseAssignment($user, $text);
1278              if ($user) {              Trace("Assignment user is $actualUser, text is $function.") if T(4);
1279                if ($actualUser) {
1280                  # Here it is a functional assignment. Check the time and the user                  # Here it is a functional assignment. Check the time and the user
1281                  # name. The time must be recent and the user must be trusted.                  # name. The time must be recent and the user must be trusted.
1282                  if ((exists $trusteeTable{$user}) && ($time > $timeSelected)) {                  if ((exists $trusteeTable{$actualUser}) && ($time > $timeSelected)) {
1283                      $retVal = $function;                      $retVal = $function;
1284                      $timeSelected = $time;                      $timeSelected = $time;
1285                  }                  }
# Line 1291  Line 1295 
1295      return $retVal;      return $retVal;
1296  }  }
1297    
1298    =head3 FunctionsOf
1299    
1300    C<< my @functionList = $sprout->FunctionOf($featureID, $userID); >>
1301    
1302    Return the functional assignments of a particular feature.
1303    
1304    The functional assignment is handled differently depending on the type of feature. If
1305    the feature is identified by a FIG ID (begins with the string C<fig|>), then a functional
1306    assignment is a type of annotation. The format of an assignment is described in
1307    L</ParseAssignment>. Its worth noting that we cannot filter on the content of the
1308    annotation itself because it's a text field; however, this is not a big problem because
1309    most features only have a small number of annotations.
1310    
1311    If the feature is B<not> identified by a FIG ID, then the functional assignment
1312    information is taken from the B<ExternalAliasFunc> table. If the table does
1313    not contain an entry for the feature, an empty list is returned.
1314    
1315    =over 4
1316    
1317    =item featureID
1318    
1319    ID of the feature whose functional assignments are desired.
1320    
1321    =item RETURN
1322    
1323    Returns a list of 2-tuples, each consisting of a user ID and the text of an assignment by
1324    that user.
1325    
1326    =back
1327    
1328    =cut
1329    #: Return Type @@;
1330    sub FunctionsOf {
1331        # Get the parameters.
1332        my ($self, $featureID) = @_;
1333        # Declare the return value.
1334        my @retVal = ();
1335        # Determine the ID type.
1336        if ($featureID =~ m/^fig\|/) {
1337            # Here we have a FIG feature ID. We must build the list of trusted
1338            # users.
1339            my %trusteeTable = ();
1340            # Build a query for all of the feature's annotations, sorted by date.
1341            my $query = $self->Get(['IsTargetOfAnnotation', 'Annotation', 'MadeAnnotation'],
1342                                   "IsTargetOfAnnotation(from-link) = ? ORDER BY Annotation(time) DESC",
1343                                   [$featureID]);
1344            my $timeSelected = 0;
1345            # Loop until we run out of annotations.
1346            while (my $annotation = $query->Fetch()) {
1347                # Get the annotation text.
1348                my ($text, $time, $user) = $annotation->Values(['Annotation(annotation)',
1349                                                                'Annotation(time)',
1350                                                                'MadeAnnotation(user)']);
1351                # Check to see if this is a functional assignment for a trusted user.
1352                my ($actualUser, $function) = _ParseAssignment($user, $text);
1353                if ($actualUser) {
1354                    # Here it is a functional assignment.
1355                    push @retVal, [$actualUser, $function];
1356                }
1357            }
1358        } else {
1359            # Here we have a non-FIG feature ID. In this case the user ID does not
1360            # matter. We simply get the information from the External Alias Function
1361            # table.
1362            my @assignments = $self->GetEntityValues('ExternalAliasFunc', $featureID,
1363                                                     ['ExternalAliasFunc(func)']);
1364            push @retVal, map { ['master', $_] } @assignments;
1365        }
1366        # Return the assignments found.
1367        return @retVal;
1368    }
1369    
1370  =head3 BBHList  =head3 BBHList
1371    
1372  C<< my $bbhHash = $sprout->BBHList($genomeID, \@featureList); >>  C<< my $bbhHash = $sprout->BBHList($genomeID, \@featureList); >>
# Line 1411  Line 1487 
1487      my $genomeData = $self->GetEntity('Genome', $genomeID);      my $genomeData = $self->GetEntity('Genome', $genomeID);
1488      if ($genomeData) {      if ($genomeData) {
1489          # The genome exists, so get the completeness flag.          # The genome exists, so get the completeness flag.
1490          ($retVal) = $genomeData->Value('complete');          ($retVal) = $genomeData->Value('Genome(complete)');
1491      }      }
1492      # Return the result.      # Return the result.
1493      return $retVal;      return $retVal;
# Line 3150  Line 3226 
3226      return $retVal;      return $retVal;
3227  }  }
3228    
3229    =head3 DeleteGenome
3230    
3231    C<< my $stats = $sprout->DeleteGenome($genomeID, $testFlag); >>
3232    
3233    Delete a genome from the database.
3234    
3235    =over 4
3236    
3237    =item genomeID
3238    
3239    ID of the genome to delete
3240    
3241    =item testFlag
3242    
3243    If TRUE, then the DELETE statements will be traced, but no deletions will occur.
3244    
3245    =item RETURN
3246    
3247    Returns a statistics object describing the rows deleted.
3248    
3249    =back
3250    
3251    =cut
3252    #: Return Type $%;
3253    sub DeleteGenome {
3254        # Get the parameters.
3255        my ($self, $genomeID, $testFlag) = @_;
3256        # Perform the delete for the genome's features.
3257        my $retVal = $self->{_erdb}->Delete('Feature', "fig|$genomeID.%", $testFlag);
3258        # Perform the delete for the primary genome data.
3259        my $stats = $self->{_erdb}->Delete('Genome', $genomeID, $testFlag);
3260        $retVal->Accumulate($stats);
3261        # Return the result.
3262        return $retVal;
3263    }
3264    
3265  =head2 Internal Utility Methods  =head2 Internal Utility Methods
3266    
3267  =head3 ParseAssignment  =head3 ParseAssignment
# Line 3160  Line 3272 
3272    
3273  A functional assignment is always of the form  A functional assignment is always of the form
3274    
3275      I<XXXX>C<\nset >I<YYYY>C< function to\n>I<ZZZZZ>      C<set >I<YYYY>C< function to\n>I<ZZZZZ>
3276    
3277  where I<XXXX> is the B<assigning user>, I<YYYY> is the B<user>, and I<ZZZZ> is the  where I<YYYY> is the B<user>, and I<ZZZZ> is the actual functional role. In most cases,
3278  actual functional role. In most cases, the user and the assigning user will be the  the user and the assigning user (from MadeAnnotation) will be the same, but that is
3279  same, but that is not always the case.  not always the case.
3280    
3281    In addition, the functional role may contain extra data that is stripped, such as
3282    terminating spaces or a comment separated from the rest of the text by a tab.
3283    
3284  This is a static method.  This is a static method.
3285    
3286  =over 4  =over 4
3287    
3288    =item user
3289    
3290    Name of the assigning user.
3291    
3292  =item text  =item text
3293    
3294  Text of the annotation.  Text of the annotation.
# Line 3185  Line 3304 
3304    
3305  sub _ParseAssignment {  sub _ParseAssignment {
3306      # Get the parameters.      # Get the parameters.
3307      my ($text) = @_;      my ($user, $text) = @_;
3308      # Declare the return value.      # Declare the return value.
3309      my @retVal = ();      my @retVal = ();
3310      # Check to see if this is a functional assignment.      # Check to see if this is a functional assignment.
3311      my ($user, $type, $function) = split(/\n/, $text);      my ($type, $function) = split(/\n/, $text);
3312      if ($type =~ m/^set ([^ ]+) function to$/i) {      if ($type =~ m/^set function to$/i) {
3313          # Here it is, so we return the user name (which is in $1), the functional role text,          # Here we have an assignment without a user, so we use the incoming user ID.
3314          # and the assigning user.          @retVal = ($user, $function);
3315          @retVal = ($1, $function, $user);      } elsif ($type =~ m/^set (\S+) function to$/i) {
3316            # Here we have an assignment with a user that is passed back to the caller.
3317            @retVal = ($1, $function);
3318        }
3319        # If we have an assignment, we need to clean the function text. There may be
3320        # extra junk at the end added as a note from the user.
3321        if (@retVal) {
3322            $retVal[1] =~ s/(\t\S)?\s*$//;
3323      }      }
3324      # Return the result list.      # Return the result list.
3325      return @retVal;      return @retVal;
# Line 3283  Line 3409 
3409  }  }
3410    
3411    
   
3412  1;  1;

Legend:
Removed from v.1.42  
changed lines
  Added in v.1.55

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3