Example #1
0
 def test_ranksums(self):
     "Testing ranksums"
     
     data1 = [ self.L, self.A ]
     data2 = [ self.M, self.B ]
     results = (-1.677105520481424, 0.09352184967262378)
     
     i = 0
     for d in data1:
         self.assertEqual( stats.ranksums( d, data2[i] )[i], results[i] )
         i += 1        
Example #2
0
    def test_ranksums(self):
        "Testing ranksums"

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

        i = 0
        for d in data1:
            self.assertEqual(stats.ranksums(d, data2[i])[i], results[i])
            i += 1
Example #3
0
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(stats.wilcoxont(l, m))
print(stats.wilcoxont(a, b))
print('kruskalwallish:')
print(stats.kruskalwallish(l, m, l))
print(len(l), len(m))
print(stats.kruskalwallish(a, b, a))
print('friedmanchisquare:')
print(stats.friedmanchisquare(l, m, l))
print(stats.friedmanchisquare(a, b, a))

print('\nANOVAs')
#execfile('test_anova.py')
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 stats.wilcoxont(l,m)
print stats.wilcoxont(a,b)
print 'kruskalwallish:'
print stats.kruskalwallish(l,m,l)
print len(l), len(m)
print stats.kruskalwallish(a,b,a)
print 'friedmanchisquare:'
print stats.friedmanchisquare(l,m,l)
print stats.friedmanchisquare(a,b,a)

print '\nANOVAs'
#execfile('test_anova.py')
Example #5
0
print stats.chisquare(fo, [5, 45])
print stats.achisquare(array(fo), array([5, 45], 'f'))

print '\n\nMann Whitney U'

red = [540, 480, 600, 590, 605]
black = [760, 890, 1105, 595, 940]
print '\nSHOULD BE 2.0, 0.01<p<0.05 (N=5,5) ... Basic Stats 1st ed, p.473-4'
print stats.mannwhitneyu(red, black)
print stats.amannwhitneyu(array(red), array(black))

print '\n\nRank Sums'

#(using red and black from above)
print '\nSHOULD BE -2.19, p<0.0286 (slightly) ... Basic Stats 1st ed, p.474-5'
print stats.ranksums(red, black)
print stats.aranksums(N.array(red), N.array(black))

print '\n\nWilcoxon T'

red = [540, 580, 600, 680, 430, 740, 600, 690, 605, 520]
black = [760, 710, 1105, 880, 500, 990, 1050, 640, 595, 520]
print '\nSHOULD BE +3.0, 0.01<p<0.05 (N=9) ... Basic Stats 1st ed, p.477-8'
print stats.wilcoxont(red, black)
print stats.awilcoxont(array(red), array(black))

print '\n\nKruskal-Wallis H'

short = [10, 28, 26, 39, 6]
medium = [24, 27, 35, 44, 58]
tall = [68, 71, 57, 60, 62]
 def evaluate( self, *args, **params):
     return _stats.ranksums(*args, **params)
print stats.achisquare(array(fo),array([5,45],'f'))


print '\n\nMann Whitney U'

red = [540,480,600,590,605]
black = [760,890,1105,595,940]
print '\nSHOULD BE 2.0, 0.01<p<0.05 (N=5,5) ... Basic Stats 1st ed, p.473-4'
print stats.mannwhitneyu(red,black)
print stats.amannwhitneyu(array(red),array(black))

print '\n\nRank Sums'

#(using red and black from above)
print '\nSHOULD BE -2.19, p<0.0286 (slightly) ... Basic Stats 1st ed, p.474-5'
print stats.ranksums(red,black)
print stats.aranksums(N.array(red),N.array(black))


print '\n\nWilcoxon T'

red   = [540,580, 600,680,430,740, 600,690,605,520]
black = [760,710,1105,880,500,990,1050,640,595,520]
print '\nSHOULD BE +3.0, 0.01<p<0.05 (N=9) ... Basic Stats 1st ed, p.477-8'
print stats.wilcoxont(red,black)
print stats.awilcoxont(array(red),array(black))

print '\n\nKruskal-Wallis H'

short = [10,28,26,39,6]
medium = [24,27,35,44,58]
Example #8
0
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 stats.wilcoxont(l, m)
print stats.wilcoxont(a, b)
print 'kruskalwallish:'
print stats.kruskalwallish(l, m, l)
print len(l), len(m)
print stats.kruskalwallish(a, b, a)
print 'friedmanchisquare:'
print stats.friedmanchisquare(l, m, l)
print stats.friedmanchisquare(a, b, a)

print '\nANOVAs'
#execfile('test_anova.py')
print(stats.achisquare(array(fo),array([5,45],'f')))


print('\n\nMann Whitney U')

red = [540,480,600,590,605]
black = [760,890,1105,595,940]
print('\nSHOULD BE 2.0, 0.01<p<0.05 (N=5,5) ... Basic Stats 1st ed, p.473-4')
print(stats.mannwhitneyu(red,black))
print(stats.amannwhitneyu(array(red),array(black)))

print('\n\nRank Sums')

#(using red and black from above)
print('\nSHOULD BE -2.19, p<0.0286 (slightly) ... Basic Stats 1st ed, p.474-5')
print(stats.ranksums(red,black))
print(stats.aranksums(N.array(red),N.array(black)))


print('\n\nWilcoxon T')

red   = [540,580, 600,680,430,740, 600,690,605,520]
black = [760,710,1105,880,500,990,1050,640,595,520]
print('\nSHOULD BE +3.0, 0.01<p<0.05 (N=9) ... Basic Stats 1st ed, p.477-8')
print(stats.wilcoxont(red,black))
print(stats.awilcoxont(array(red),array(black)))

print('\n\nKruskal-Wallis H')

short = [10,28,26,39,6]
medium = [24,27,35,44,58]