[Bioperl-l] [Gmod-gbrowse] bioperl-live & annotator plugin
Lincoln Stein
lincoln.stein at gmail.com
Tue Dec 2 10:51:06 EST 2008
I take it back. You should add the matrix in add_type(), assuming that
$matrix is some sort of descriptive term.
Lincoln
On Tue, Dec 2, 2008 at 10:49 AM, Lincoln Stein <lincoln.stein at gmail.com>wrote:
> Hi Sean,
>
> Did your plugin break when you upgraded to GBrowse 1.7, or when you
> upgraded bioperl? Knowing this will help me figure out where to look for the
> problem. It frankly looks like a database problem in the sense that the
> xyplot glyph is trying to render 100 features rather than a feature with 100
> parts.
>
> Are you sure this step is correct?
>
> $feature_list->add_type($matrix=>{glyph => 'xyplot',
>> graph_type => 'boxes',
>> scale => 'auto',
>> height => 40,
>> bump => 0,
>> min_score=>0,
>> max_score=>$maxscore,
>> key => "$matrix affinities ($maxscore)",
>> fgcolor => $COLORS[$i % @COLORS],
>> bgcolor => $COLORS[$i % @COLORS],
>> });
>
>
> It seems to me you intend to add the $feature, not the $matrix variable.
>
> Lincoln
>
> On Tue, Dec 2, 2008 at 10:32 AM, Sean O'Keeffe <okeeffe at molgen.mpg.de>wrote:
>
>> Hi all,
>> I'm experiencing a little issue when upgrading to the latest bioperl-live
>> build.
>> We needed the new glyph libraries as older versions can't display wiggle
>> density plots and also experience versioning issues for Gbrowse 1.7.
>> On upgrading to the latest nightly build, all seemed fine until we tried
>> our TF affinity annotator plugin.
>> The plugin basically generates GFF output of a selected matrix binding to
>> the sequence of the current window, and displays the score values as an
>> xyplot.
>> We've never had any real issues 'til the new libraries were loaded. I
>> tried switching the old xyplot.pm with the new one but to no avail.
>> I've attached 2 images old and new of the generated output, plus some code
>> and feature file output.
>> Any insight would be great,
>>
>> Sean.
>>
>> +++++++++++++++++++++++++++++++++++++++++++++++++++
>>
>> the code of interest in the plugin
>>
>> ......
>> my $i = 0;
>> my (@features,$feature,$start,$end,$score,$maxscore,$strand);
>> $feature_list = Bio::Graphics::FeatureFile->new();
>>
>> foreach $matrix (@matrices){
>> $maxscore=0;
>> # Usually the track is created first but here we want the maxscore
>> passed to
>> # the track so we process the features first and create the track after.
>> foreach $_ (@result){
>> next if (/\#/);
>> next unless (/$matrix/);
>> @features = split("\t", $_);
>>
>> $start=$features[3];$end=$features[4];$score=$features[5];$strand=$features[6];
>> if ($score>$maxscore) {
>> $maxscore=$score;
>> }
>> # Fill the feature list with features
>> $feature = Bio::Graphics::Feature->new(-start=>$abs_start + $start,
>> -stop=>$abs_start + $end,
>> -score=>$score,
>> -strand=>$strand,
>> -ref=>$ref,
>> -name=>$matrix,
>> -type=>'region',
>> -source=>'ANNOTATE');
>> $feature_list->add_feature($feature,$matrix);
>> }
>> # Create the track and add features
>> $feature_list->add_type($matrix=>{glyph => 'xyplot',
>> graph_type => 'boxes',
>> scale => 'auto',
>> height => 40,
>> bump => 0,
>> min_score=>0,
>> max_score=>$maxscore,
>> key => "$matrix affinities ($maxscore)",
>> fgcolor => $COLORS[$i % @COLORS],
>> bgcolor => $COLORS[$i % @COLORS],
>> });
>> $i++;
>> }
>> return $feature_list;
>> .........
>>
>> +++++++++++++++++++++++++++++++++++++++++++
>>
>> the feature file which is passed to the plugin
>>
>> #command_line= /project/gbrowse/Pipeline/Annotation/ANNOTATE_v2.07 -s
>> /project/general/tmp/TIwySJV5kx --pscm /project/general/tmp/FkztCDpBvb -name
>> SRF_C -a -g 0.5 -t 0 -R 200 -S 100 -p 1.0
>> segment ANNOTATE region 1 200 0.000103017 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 101 300 0.000101877 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 201 400 6.76318e-05 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 301 500 7.0194e-05 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 401 600 3.36309e-05 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 501 700 7.15977e-05 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 601 800 0.000639066 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 701 900 0.000589631 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 801 1000 0.000580096 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> segment ANNOTATE region 901 1100 0.00139282 + .
>> matrix_id=M00215;matrix_name='V$SRF_C';seq_id=segment;def_L=200
>> ...
>>
>>
>>
>>
>> -------------------------------------------------------------------------
>> This SF.Net email is sponsored by the Moblin Your Move Developer's
>> challenge
>> Build the coolest Linux based applications with Moblin SDK & win great
>> prizes
>> Grand prize is a trip for two to an Open Source event anywhere in the
>> world
>> http://moblin-contest.org/redirect.php?banner_id=100&url=/
>> _______________________________________________
>> Gmod-gbrowse mailing list
>> Gmod-gbrowse at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/gmod-gbrowse
>>
>>
>
>
> --
> Lincoln D. Stein
>
> Ontario Institute for Cancer Research
> 101 College St., Suite 800
> Toronto, ON, Canada M5G0A3
> 416 673-8514
> Assistant: Stacey Quinn <Stacey.Quinn at oicr.on.ca>
>
> Cold Spring Harbor Laboratory
> 1 Bungtown Road
> Cold Spring Harbor, NY 11724 USA
> (516) 367-8380
> Assistant: Sandra Michelsen <michelse at cshl.edu>
>
--
Lincoln D. Stein
Ontario Institute for Cancer Research
101 College St., Suite 800
Toronto, ON, Canada M5G0A3
416 673-8514
Assistant: Stacey Quinn <Stacey.Quinn at oicr.on.ca>
Cold Spring Harbor Laboratory
1 Bungtown Road
Cold Spring Harbor, NY 11724 USA
(516) 367-8380
Assistant: Sandra Michelsen <michelse at cshl.edu>
More information about the Bioperl-l
mailing list