Example #1
0
 def t_test(self,A,B):
     """
     Perform a Student's t-test, for dependent paired two-tailed case.
     """
     from statlib import stats
     (tval,p) = stats.ttest_rel(A,B,1)
     return {"tval":tval,"prob":p}
Example #2
0
    def test_ttest_rel(self):
        "Testing ttest_rel"

        data1 = [self.L, self.A]
        data2 = [self.M, self.B]
        results = (-4.0, 0.00076619233678407726)

        i = 0
        for d in data1:
            self.assertEqual(stats.ttest_rel(d, data2[i])[i], results[i])
            i += 1
Example #3
0
 def test_ttest_rel(self):
     "Testing ttest_rel"
     
     data1 = [ self.L, self.A ]
     data2 = [ self.M, self.B ]
     results = (-4.0, 0.00076619233678407726)
     
     i = 0
     for d in data1:
         self.assertEqual( stats.ttest_rel( d, data2[i] )[i], results[i] )
         i += 1
Example #4
0
print('kendalltau:')
print(stats.kendalltau(l, m))
print(stats.kendalltau(a, b))
print('linregress:')
print(stats.linregress(l, m))
print(stats.linregress(a, b))

print('\nINFERENTIAL')
print('ttest_1samp:')
print(stats.ttest_1samp(l, 12))
print(stats.ttest_1samp(a, 12))
print('ttest_ind:')
print(stats.ttest_ind(l, m))
print(stats.ttest_ind(a, b))
print('ttest_rel:')
print(stats.ttest_rel(l, m))
print(stats.ttest_rel(a, b))
print('chisquare:')
print(stats.chisquare(l))
print(stats.chisquare(a))
print('ks_2samp:')
print(stats.ks_2samp(l, m))
print(stats.ks_2samp(a, b))

print('mannwhitneyu:')
print(stats.mannwhitneyu(l, m))
print(stats.mannwhitneyu(a, b))
print('ranksums:')
print(stats.ranksums(l, m))
print(stats.ranksums(a, b))
print('wilcoxont:')
print 'kendalltau:'
print stats.kendalltau(l,m)
print stats.kendalltau(a,b)
print 'linregress:'
print stats.linregress(l,m)
print stats.linregress(a,b)

print '\nINFERENTIAL'
print 'ttest_1samp:'
print stats.ttest_1samp(l,12)
print stats.ttest_1samp(a,12)
print 'ttest_ind:'
print stats.ttest_ind(l,m)
print stats.ttest_ind(a,b)
print 'ttest_rel:'
print stats.ttest_rel(l,m)
print stats.ttest_rel(a,b)
print 'chisquare:'
print stats.chisquare(l)
print stats.chisquare(a)
print 'ks_2samp:'
print stats.ks_2samp(l,m)
print stats.ks_2samp(a,b)

print 'mannwhitneyu:'
print stats.mannwhitneyu(l,m)
print stats.mannwhitneyu(a,b)
print 'ranksums:'
print stats.ranksums(l,m)
print stats.ranksums(a,b)
print 'wilcoxont:'
Example #6
0
stats.attest_1samp(array(x), 75, 1)

print '\n\nIndependent Samples t-test'

a = [11, 16, 20, 17, 10, 12]
b = [8, 11, 15, 11, 11, 12, 11, 7]
print '\n\nSHOULD BE ??? <p< (df=) ... '
stats.ttest_ind(a, b, 1)
stats.attest_ind(array(a), array(b), 0, 1)

print '\n\nRelated Samples t-test'

before = [11, 16, 20, 17, 10]
after = [8, 11, 15, 11, 11]
print '\n\nSHOULD BE t=+2.88, 0.01<p<0.05 (df=4) ... Basic Stats 1st ed, p.359'
stats.ttest_rel(before, after, 1, 'Before', 'After')
stats.attest_rel(array(before), array(after), 1, 'Before', 'After')

print "\n\nPearson's r"

x = [0, 0, 1, 1, 1, 2, 2, 3, 3, 4]
y = [8, 7, 7, 6, 5, 4, 4, 4, 2, 0]
print 'SHOULD BE -0.94535 (N=10) ... Basic Stats 1st ed, p.190'
print stats.pearsonr(x, y)
print stats.apearsonr(array(x), array(y))

print "\n\nSpearman's r"

x = [4, 1, 9, 8, 3, 5, 6, 2, 7]
y = [3, 2, 8, 6, 5, 4, 7, 1, 9]
print '\nSHOULD BE +0.85 on the dot (N=9) ... Basic Stats 1st ed, p.193'
 def evaluate( self, *args, **params):
     return _stats.ttest_rel(*args, **params)
print '\n\nIndependent Samples t-test'

a = [11,16,20,17,10,12]
b = [8,11,15,11,11,12,11,7]
print '\n\nSHOULD BE ??? <p< (df=) ... '
stats.ttest_ind(a,b,1)
stats.attest_ind(array(a),array(b),0,1)


print '\n\nRelated Samples t-test'

before = [11,16,20,17,10]
after = [8,11,15,11,11]
print '\n\nSHOULD BE t=+2.88, 0.01<p<0.05 (df=4) ... Basic Stats 1st ed, p.359'
stats.ttest_rel(before,after,1,'Before','After')
stats.attest_rel(array(before),array(after),1,'Before','After')


print "\n\nPearson's r"

x = [0,0,1,1,1,2,2,3,3,4]
y = [8,7,7,6,5,4,4,4,2,0]
print 'SHOULD BE -0.94535 (N=10) ... Basic Stats 1st ed, p.190'
print stats.pearsonr(x,y)
print stats.apearsonr(array(x),array(y))


print "\n\nSpearman's r"

x = [4,1,9,8,3,5,6,2,7]
Example #9
0
print 'kendalltau:'
print stats.kendalltau(l, m)
print stats.kendalltau(a, b)
print 'linregress:'
print stats.linregress(l, m)
print stats.linregress(a, b)

print '\nINFERENTIAL'
print 'ttest_1samp:'
print stats.ttest_1samp(l, 12)
print stats.ttest_1samp(a, 12)
print 'ttest_ind:'
print stats.ttest_ind(l, m)
print stats.ttest_ind(a, b)
print 'ttest_rel:'
print stats.ttest_rel(l, m)
print stats.ttest_rel(a, b)
print 'chisquare:'
print stats.chisquare(l)
print stats.chisquare(a)
print 'ks_2samp:'
print stats.ks_2samp(l, m)
print stats.ks_2samp(a, b)

print 'mannwhitneyu:'
print stats.mannwhitneyu(l, m)
print stats.mannwhitneyu(a, b)
print 'ranksums:'
print stats.ranksums(l, m)
print stats.ranksums(a, b)
print 'wilcoxont:'