我有38400个正态分布的测量值。平均值为-4.10e-11,西格玛为0.0229。毫不奇怪,概率图的斜率为0.0229。为了绘制和检测离群值,我想用1的斜率绘制它。理论分位数的范围应与观测值大致相同,即在-0.09和+0.09之间。 我该怎么做?在
(对于代码示例,我使用随机数而不是实际测量值)。在
import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as stats
mean = -4.10e-11
sigma = 0.0229
np.random.seed(7654321)
e_array = np.random.normal(mean, sigma, 38400)
minvalue = np.min(e_array)
maxvalue = np.max(e_array)
qq = stats.probplot(e_array, dist="norm", plot=plt)
qq_slope, qq_intercept, qq_r = qq[1]
print('Minimum value: ', minvalue)
print('Maximum value: ', maxvalue)
print('QQ slope:', qq_slope)
print('QQ intercept:', qq_intercept)
print('QQ r:', qq_r)
plt.show()
我想这就是你的意思。在
我根据理论cdf绘制样本cdf,结果散点图应该围绕一条斜率为1的直线聚集。在
这是一个结果图。别装腔作势。在
相关问题 更多 >
编程相关推荐