python显著性检验

fanpy

相关系数(回归系数)检验

其实在(2)中已经提到了相关系数和回归系数，在计算过程中，直接返回了对应的p-value，因此可以直接使用p-value。

1.相关系数from scipy.stats import pearsonrr,p = pearsonr(x, y)2.线性回归系数from sklearn.linear_model import LinearRegressionslope, intercept, r_value, p_value, std_err = stats.linregress(x,y)
均值检验t-test

计算两个独立样本得分均值的T检验。
这是对两个独立样本具有相同平均值（预期值）的零假设的双边检验。此测试假设默认情况下总体具有相同的方差。在合成分析中通常用到t-test。

from scipy.stats.mstats import ttest_indt,p = ttest_ind(a, b, 0, equal_var=True, nan_policy='propagate')当a,b为变量场时，即[time,lat,lon]时，a,b两个数组的经纬度需相同。
nan_policy 可选{‘propagate’, ‘raise’, ‘omit’}
“propagate”：返回nan
“raise”：报错
“omit”：执行忽略nan值的计算
绘图计算得到的P值用于绘图，当p<0.01时，通过99%显著性检验，p<0.05，通过95%显著性检验，以此类推。
图形绘制只需在原有填色图上叠加打点图层,实际上打点也是特殊的图色，只不过将颜色换成了点，实际上用到的还是contourf函数。

proj = ccrs.PlateCarree(central_longitude=95)leftlon, rightlon, lowerlat, upperlat = (0,160,30,90)img_extent = [leftlon, rightlon, lowerlat, upperlat]fig3 = plt.figure(figsize=(12,8),dpi=1000)f3_ax1 = fig3.add_axes([1.24, 0.8, 0.5, 0.5],projection = proj)contour_map(f3_ax1,img_extent,20)c1 = f3_ax1.contourf(lon,lat, z500_cs_0_0, zorder=0,levels =range(-80,90,10) , extend = 'both', transform=ccrs.PlateCarree(), cmap=plt.cm.RdBu_r)c1b = f3_ax1.contourf(lon,lat, p_0_500_0,[np.min(p_0_500_0),0.05,np.max(p_0_500_0)], zorder=1,hatches=['.', None],colors="none", transform=ccrs.PlateCarree())












作者：摸鱼咯

来源：简书
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。