[Bioperl-l] [Gmod-gbrowse] bioperl-live & annotator plugin
Lincoln Stein
lincoln.stein at gmail.com
Tue Dec 2 10:49:45 EST 2008
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>
More information about the Bioperl-l
mailing list