基于python的显著性检验

需要用到numpy库

import numpy as np
import scipy.stats as stats
import scipy.optimize as opt

首先我们来创造两个数组作为测试数据

n = 200
norm_dist = stats.norm(loc=0.5, scale=10)  #构造一个正态分布，均值为0.5，标准差为10  “标准差”也称“均方差”，是“方差”开根号
dat = norm_dist.rvs(size=n)         #随机取200个点
print ("mean of data is: " + str(np.mean(dat)))
print ("median of data is: " + str(np.median(dat)))
print ("standard deviation of data is: " + str(np.std(dat)))  #因为这200个点是随机取得，所以跟原先的正态分布可能有一些不同

norm_dist2 = stats.norm(loc=0.2, scale=1)
dat2 = norm_dist2.rvs(size=n/2)#随机取100个点
print ("mean of data is: " + str(np.mean(dat2)))
print ("median of data is: " + str(np.median(dat2)))
print ("standard deviation of data is: " + str(np.std(dat2)))

对这两个数组分析差异---双样本的t检验

stat_val, p_val = stats.ttest_ind(dat, dat2, equal_var=False)
#看看两个分布在均值上有没有显著差异
#注意，这里我们生成的第二组数据样本大小、方差和第一组均不相等，在运用t检验时需要使用Welch's t-test
#即指定ttest_ind中的equal_var=False。
print ('Two-sample t-statistic D = %6.3f, p-value = %6.4f' % (stat_val, p_val))

计算两个序列的相关性，并做显著性检验

import scipy.stats as stats  
x = [76,81,78,76,76,78,76,78,98,88,76,66,44,67,65,59,87,77,79,85,68,76,77,98,99,98,87,67,78]    
y = [43,33,23,34,31,51,56,43,44,45,32,33,28,39,31,38,21,27,43,46,41,41,48,56,55,45,68,54,33] 
r, p=stats.pearsonr(x,y) 
[out]:(0.39341862097439129, 0.034735931329532836) 

相关系数为0.39，说明这两个序列存在一定的相关性
p-value为0.035，说明结果是统计显著的