1. 程式人生 > >分位數、上側分位數及python實現

分位數、上側分位數及python實現

分位數

定義:設隨機變數X的分佈函式為F(x),對任意給定的實數 p ,(0<p<1),若存在xpx_p 使得 P{XX\leq xpx_p } = F(xpx_p) = p成立,那麼稱xpx_p為此概率分佈的p分位數。 例: 在這裡插入圖片描述 如上圖所示,影象為標準正態分佈的概率密度函式,xpx_p 為p分位數。

上側分位數

定義:對隨機變數X和給定的 α\alpha (0<α\alpha<1),若存在 xαx_\alpha ,使得 P{X \geq xαx_\alpha } = α\alpha,那麼稱 xαx_\alpha

為X的上側分位數。 例:在這裡插入圖片描述 如上圖所示,xαx_\alpha 為標準正態分佈的上側 α\alpha 分位數。 在數理統計教程中,又把標準正態分佈的上側分位數記為:uαu_\alpha

Python實現

首先介紹一下幾個常用分佈的函式 import scipy.stats as st st.norm() #正態分佈 st.t() #t分佈 st.f() #f分佈 st.chi2() #χ2\chi^2 分佈 然後介紹一下基本的操作: pdf 概率密度函式

cdf 分佈函式 ppf 分佈函式的逆 sf 殘損函式(1-cdf) isf 逆殘損函式

下面是程式碼實現

在這裡插入圖片描述 看到結果,你可能會問,怎麼跟數理統計書附錄上的t分位數表不一樣,是因為書上求的是上側分位數表。

最後再畫個圖吧

在這裡插入圖片描述