File Coverage

blib/lib/XMS/MotifSet.pm

Criterion	Covered	Total	%
statement	13	15	86.6
branch			n/a
condition			n/a
subroutine	5	5	100.0
pod			n/a
total	18	20	90.0

package XMS::MotifSet;

use 5.008008;
use strict;
use warnings;
use Carp;
require Exporter;

our @ISA = qw(Exporter);


# This allows declaration	use XMS ':all';


our %EXPORT_TAGS = ( 'all' => [ qw(

) ] );

our @EXPORT_OK = ( @{ $EXPORT_TAGS{'all'} } );

our @EXPORT = qw(

);

our $VERSION = '0.01';

use XML::Writer;
use XML::Writer::String;
use IO::File;
use XMS::WeightMatrix;
use XMS::Motif;
use XML::DOM;



sub new {

    my $class = shift;
    my $self = {};
    my(@motifs) = @_;
    my ($temp) = @motifs;

    if(ref($temp) eq "XMS::Motif"){
	@{$self->{motifs}} = @motifs;

    }else{

	my $parser = new XML::DOM::Parser;
	my $doc = $parser->parsefile($temp);
	my $root = $doc->getDocumentElement();

	my @motifnodes = $root->getElementsByTagName("motif");

	if (scalar @motifnodes == 0){
	    die "Corrupt input XMS file";
	}

	my $m=0;
	my @motifarray=();

	foreach my $motif (@motifnodes){


	    my $motifname = $motif->getElementsByTagName("name")->item(0)->getFirstChild()->getData;

	    my $threshold = $motif->getElementsByTagName("threshold")->item(0)->getFirstChild->getData;

	    my @props = $motif->getElementsByTagName("prop");

	    ####### Begin reading annotation key value pairs ##########

	    my %annotations;
	    foreach my $prop (@props){
		my ($keynode) = $prop->getElementsByTagName("key");
		my $key="";
		if ( $keynode->getFirstChild() ) {
		    $key = $keynode->getFirstChild()->getData;
		}
		my ($valuenode) = $prop->getElementsByTagName("value");
	 	my $value="";
		if ( $valuenode->getFirstChild() ) {
		    $value = $valuenode->getFirstChild()->getData;
		}
		if ( $key ne "" ){
		    $annotations{$key} = $value;
		}
	    }

	    ####### End reading annotation key value pairs ##########

	    my @wmnodes = $motif->getElementsByTagName("weightmatrix");
	    my @columnsarray=();
	    foreach my $wmnode (@wmnodes) {

		my @columns=$motif->getElementsByTagName("column");

		my %columnhash;
		foreach my $column (@columns){

		   my @weights=$column->getElementsByTagName("weight");

		    foreach my $weight (@weights){

			my $weightsymbol = $weight->getAttributeNode("symbol");
			my $symbolvalue = $weightsymbol->getValue;

			my $weightvalue=$weight->getFirstChild->getData;

			$columnhash{$symbolvalue} = $weightvalue;
		    }
		   my $count=0;
		   my @wmvalues=();
		   foreach my $key (sort keys %columnhash){
		       $wmvalues[$count]=$columnhash{$key};
		       $count++;
		   }
		   push(@columnsarray,[@wmvalues]);
		}
	    }

	    my $wmobj = XMS::WeightMatrix->new(@columnsarray);
	    my $motifobj = XMS::Motif->new($wmobj,$motifname,$threshold,%annotations);
	    $motifarray[$m] = $motifobj;
	    $m++;
	}
	@{$self->{motifs}} = @motifarray;
    }

    $self->{output} =  XML::Writer::String->new();
    $self->{writer} = new XML::Writer(OUTPUT => $self->{output}, DATA_MODE => 'TRUE', DATA_INDENT=>3);

    bless($self,$class);
    return $self;

}


sub toXML {

    my $self = shift;
    my $writer = $self->{writer};
    my $output = $self->{output};

    my @motifs = @{$self->{motifs}};


    $writer->startTag("motifset");
    for(my $m=0;$m<@motifs;$m++){
	$self->{motifs}[$m]->toXML($writer);
    }

    $writer->endTag("motifset");

    $writer->end();
    return $self->{output}->value();
}


sub toString {

    my $self = shift;
    my @motifs = @{$self->{motifs}};

    my $rawstring = "";

    for(my $m=0;$m<@motifs;$m++){
        $rawstring = $rawstring.$self->{motifs}[$m]->toString();
        if ($m<(@motifs-1)){
	    $rawstring = $rawstring."\n";
	}
    }
    return $rawstring;
}


1;
__END__