File Coverage

blib/lib/Statistics/DependantTTest.pm

Criterion	Covered	Total	%
statement	52	53	98.1
branch	1	2	50.0
condition	1	3	33.3
subroutine	8	8	100.0
pod	0	3	0.0
total	62	69	89.8

package Statistics::DependantTTest;
use strict;
use Carp;
use vars qw($VERSION);
use Statistics::PointEstimation;
no strict 'refs';
$VERSION='0.03';

##############
sub new
{
    my $proto = shift;
    my $class = ref($proto) || $proto;
    my $self= {};

    $self->{sample_data}  = undef;
    $self->{s}            = undef;



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

##############
sub load_data {
    my $self = shift;
    my $sample_name = shift;

    my (@sample_data)=@_;

    my $s= new Statistics::PointEstimation;

    $s->add_data(\@sample_data);

    $self->{sample_data}->{$sample_name}=\@sample_data;
    $self->{s}->{$sample_name}=$s;

    return $self;
} # end sub load_data
##############

##############
sub perform_t_test {
    my $self=shift;
    my $first_sample_name = shift;
    my $second_sample_name = shift;
    my $s1=$self->{s}->{$first_sample_name};
    my $s2=$self->{s}->{$second_sample_name};


    if(@{$self->{sample_data}->{$first_sample_name}} != @{$self->{sample_data}->{$second_sample_name}})
    {
	croak "The two results sets are of different length.\n For the paired T-test, the two results must be from two identical sets of subjects tested in a reference and test condition.\n";
    }

    my $s = new Statistics::PointEstimation;

    my @sample_difference;

    my $count=0;
    while ($count < @{$self->{sample_data}->{$first_sample_name}})
    {
	$sample_difference[$count] = ${$self->{sample_data}->{$first_sample_name}}[$count] - ${$self->{sample_data}->{$second_sample_name}}[$count];
        $count++;
    } # end while

    $s->add_data(\@sample_difference);

    my $t_value = $s->t_statistic();

    my $deg_freedom = scalar @{$self->{sample_data}->{$first_sample_name}} - 1;

    return ($t_value, $deg_freedom );
} # end sub perform_t_test
##############


1;
__END__