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}
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
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
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:'
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]
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:'