方法和詳細的操作步驟如下:
創新互聯建站專注為客戶提供全方位的互聯網綜合服務,包含不限于做網站、成都做網站、榕江網絡推廣、微信小程序定制開發、榕江網絡營銷、榕江企業策劃、榕江品牌公關、搜索引擎seo、人物專訪、企業宣傳片、企業代運營等,從售前售中售后,我們都將竭誠為您服務,您的肯定,是我們最大的嘉獎;創新互聯建站為所有大學生創業者提供榕江建站搭建服務,24小時服務熱線:028-86922220,官方網址:www.yijiale78.com
1、第一步,新建Excel文檔,見下圖,轉到下面的步驟。
2、第二步,執行完上面的操作之后,輸入x軸值(計算分布度),例如區間[-1,1],間隔為0.1,見下圖,轉到下面的步驟。
3、第三步,執行完上面的操作之后,由AVERAGE函數計算的平均值為0,見下圖,轉到下面的步驟。
4、第四步,執行完上面的操作之后,選擇函數STDEV并計算標準偏差,見下圖,轉到下面的步驟。
5、第五步,執行完上面的操作之后,選擇正態分布函數NORMDIST并計算返回概率密度分布值,見下圖,轉到下面的步驟。
6、第六步,執行完上面的操作之后,選擇“圖表”--“折線圖”選項,然后完成分布圖,見下圖。這樣,就解決了這個問題了。 ?
Python正態分布概率計算方法,喜歡算法的伙伴們可以參考學習下。需要用到math模塊。先了解一下這個模塊方法,再來寫代碼會更好上手。
def st_norm(u):
'''標準正態分布'''
import math
x=abs(u)/math.sqrt(2)
T=(0.0705230784,0.0422820123,0.0092705272,
0.0001520143,0.0002765672,0.0000430638)
E=1-pow((1+sum([a*pow(x,(i+1))
for i,a in enumerate(T)])),-16)
p=0.5-0.5*E if u0 else 0.5+0.5*E
return(p)
def norm(a,sigma,x):
'''一般正態分布'''
u=(x-a)/sigma
return(st_norm(u))
while 1:
'''輸入一個數時默認為標準正態分布
輸入三個數(空格隔開)時分別為期望、方差、x
輸入 stop 停止'''
S=input('please input the parameters:\n')
if S=='stop':break
try:
L=[float(s) for s in S.split()]
except:
print('Input error!')
continue
if len(L)==1:
print('f(x)=%.5f'%st_norm(L[0]))
elif len(L)==3:
print('f(x)=%.5f'%norm(L[0],L[1],L[2]))
else:
print('Input error!')
1.直方圖
由于正態分布具有非常典型的中間高,兩邊低的圖形特征,如果樣本數據并不服從正態分布,我們可以通過直方圖很快地分辨出來。更進一步地,Python可以輔助生成基于樣本數據估計的正態曲線,這樣就容易輔助我們進行判斷。
圖形觀察雖然直觀,但是部分研究者認為單純觀察圖形過于主觀,因此我們也可以選擇使用統計檢驗的方法去研究數據是否服從正態分布。
操作步驟:
導入相關的包及數據
2 P-P圖及Q-Q圖
直方圖是最長用于觀察數據分布的常用圖形選項,尤其是帶正態曲線的直方圖,可以非常直觀地看到實際數據分布和正態曲線的對比,而P-P圖及Q-Q圖則是另一種選擇,它可以直觀給出實際數據分布和理論的差距。
值得注意的是,雖然P-P圖及Q-Q圖常用用于判斷數據樣本是否服從正態分布,但實際上它們也能判斷數據樣本是否服從其他的分布
P-P圖:反映的是數據的實際累積概率與假定所服從分布的理論累積概率的符合程度。在此處,我們所假定的分布就是正態分布,如果數據樣本是服從正態分布的話,那么實際的累積概率與理論的累積概率應該是相對一致的,放映在圖形中就是數據點應該沿著圖形的對角線分布。
Q-Q圖的原理與P-P圖幾乎一致。P-P圖考察的是實際分布與理論分布的累積概率分布差異,而Q-Q圖考察的是實際百分位數與理論百分位數的差異。同理在此處,我們所假定的分布就是正態分布,如果數據樣本是服從正態分布的話,那么實際的分布應該是相對一致的,反映在圖形中就是數據點應該沿著圖形的對角線分布。
在Python中,statsmodels包中目前主要提供的是Q-Q圖的繪制
柯爾莫戈洛夫-斯米諾夫檢驗(Kolmogorov-Smirnov test),一般又稱K-S檢驗,是一種基于累計分布函數的非參數檢驗,用以檢驗兩個經驗分布是否不同或一個經驗分布與另一個理想分布是否不同。
K-S檢驗的原假設是“樣本數據來自的分布與正態分布無顯著差異”,因此一般來說,KS檢驗最終返回兩個結果,分別是檢驗統計量及P值,檢驗結果P0.05才是我們的目標。
實際上,GraphPad不推薦使用單純的Kolmogorov-Smirnov test方法
夏皮洛-威爾克檢驗(Shapiro—Wilk test),一般又稱W檢驗。W檢驗是一種類似于利用秩進行相關性檢驗的方法。同樣需要注意的是,W檢驗與K-S檢驗一樣,原假設是“樣本數據來自的分布與正態分布無顯著差異”,因此一般來說,W檢驗最終返回兩個結果,分別是檢驗統計量及P值。,檢驗結果P0.05才是我們的目標。
當數據集中的數據無重復值時,該方法的檢驗效果比較好,但是當數據集中有些數據不是獨一無二的,即有些數據的數值是相同的,那么該方法的檢驗效果就不是很好
GraphPad官方推薦使用該方法。
首先計算 偏度和峰度以便在不對稱和形狀方面量化分布離高斯分布的距離。然后,其計算這些值中的每一個與高斯分布的預期值之間的差異,并基于這些差異的總和,計算各P值。這是一種通用和強大的正態性檢驗,推薦使用。請注意,D'Agostino開發了幾種正態性檢驗。Prism使用的其中一個是“綜合K2”檢驗。
安德森-達令檢驗樣本數據是否來自特定分布,包括分布:'norm', 'expon', 'gumbel', 'extreme1' or 'logistic'.
原假設 H0:樣本服從特定分布; 備擇假設 H1:樣本不服從特定分布
實際上,從已有的文獻表明,對于數據分布的正態性研究,首選方法是圖形觀察,即利用直方圖、P-P圖或Q-Q圖進行觀察,如果分布嚴重偏態和尖峰分布則建議進行進一步的假設檢驗。如果圖形分布結果不好判斷,則再進行正態性檢驗。
實際上,從已有的文獻表明,對于數據分布的正態性研究,首選方法是圖形觀察,即利用直方圖、P-P圖或Q-Q圖進行觀察,如果分布嚴重偏態和尖峰分布則建議進行進一步的假設檢驗。如果圖形分布結果不好判斷,則再進行正態性檢驗。
其次,對于檢驗方法來說,對于K-S檢驗及W檢驗結果來說,有文獻采用蒙特卡羅模擬方法進行多次驗證,結果表明W檢驗結果相比于大部分方法都有較大的檢驗功效,而K-S方法的檢驗結果相對不佳。并且部分學者認為,K-S檢驗的實用性遠不如圖形工具,因為在樣本量少時,該檢驗不太敏感,但是在樣本量大時,該檢驗卻過于敏感。因此正常情況下,我們更常采用W檢驗的結果。
值得注意的是,雖然說K-S檢驗結果相對不佳,但是不同檢驗方法對于樣本量的敏感度是不一樣的。在樣本量較小的情況下(小于50個樣本的情況下),請優先選擇W檢驗;在樣本量50-5000的情況下,可以酌情使用W檢驗及K—S檢驗;在樣本量大于5000的情況下,請使用K-S檢驗結果,尤其是在SPSS中,當樣本量大于5000的情況下,將只顯示K-S檢驗結果,而不顯示W檢驗結果。
import scipy.stats as sta
import math
def option_call(s,x,r,sigma,t):
d1=(math.log(s/x)+(r+sigma**2/2)*t)/(math.sqrt(t)*sigma)
d2=d1-sigma*math.sqrt(t)
c=s*sta.norm.cdf(d1,0,1)-x*sta.norm.cdf(d2,0,1)*math.exp(-r*t)
return c
網站名稱:累計正態函數python 正態累積函數
分享路徑:http://www.yijiale78.com/article24/hhscce.html
成都網站建設公司_創新互聯,為您提供做網站、定制網站、微信公眾號、服務器托管、營銷型網站建設、自適應網站
聲明:本網站發布的內容(圖片、視頻和文字)以用戶投稿、用戶轉載內容為主,如果涉及侵權請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網站立場,如需處理請聯系客服。電話:028-86922220;郵箱:631063699@qq.com。內容未經允許不得轉載,或轉載時需注明來源: 創新互聯