医学统计:随机变量的概率分布、概率计算的SAS实现

概率分布 (Probability Distribution),用于描述随机变量取某个特定的值或取某一区间范围内值的概率。

可以这样理解:概率1以一定的规律,分布在各个可能值或可能范围上

概率分布可以分为离散型概率分布和连续型概率分布。对于离散型随机变量,概率分布描述变量取某一可能值的概率;对于连续型随机变量,概率分布描述变量取某一范围内值的概率。

1. 相关概念

1.1. 概率计算

概率计算是指计算随机变量在某个区间或某个具体取值上的概率。不同类型的随机变量使用的概率计算方法不同:

  • 对于离散型随机变量,可以直接求取其概率质量函数 (PMF)。
  • 对于连续型随机变量,则需要求取其概率密度函数 (PDF),然后通过积分计算概率。

概率质量函数 (Probability Mass Function) 是离散型随机变量的概率分布函数。它描述了离散型随机变量取每一个特定值的概率。

概率密度函数 (Probability Density Function) 是连续型随机变量的概率分布函数。它描述了随机变量在某一特定值附近的概率密度,而非直接描述随机变量取某一特定值的概率。

1.2. 累积分布函数 (Cumulative Distribution Function)

累积分布函数,也称分布函数,是用于描述随机变量在某个值以下的概率的函数。对于随机变量 X,其累积分布函数 F(x)定义为:

F(x) = P(X \leq x)

这表示随机变量X小于或等于某个值x的概率。

对于离散型随机变量,其累积分布函数是一个阶梯状的函数,常见的离散型分布有二项分布、泊松分布等。

对于连续型随机变量,其累积分布函数是一个连续的、光滑的函数。常见的连续型分布有正态分布、指数分布等。

1.3. 分位数/分位点 (Quantile)

分位数,是将数据分布划分成若干等份的点,表示在统计学中,将一组数据按大小排序后,某一百分比例位置上的数值。分位数在数据分析中用于描述数据的分布情况,衡量数据的集中趋势和离散程度。

设随机变量X的累积分布函数为F(x),对于给定的概率p(0到1之间),分位数p_q是满足以下条件的数值:

F(q_p) = P(X \leq q_p) = p

这意味着随机变量X 取值小于或等于p_q的概率为p

下面介绍几个常见的随机变量,以及用SAS实现相关的概率计算。

2. 离散型随机变量

离散型随机变量,是指全部可能取到的值是有限个或可列无限多个

常见的的离散型随机变量有:

  1. (0-1)分布 (两点分布)
  2. 二项分布 (伯努利试验)
  3. 泊松分布

2.1. (0-1)分布

对于(0-1)分布,其概率分布规律如下,即,随机变量取0或1这两个值的概率。

P\{X=k\} = p^k (1-p)^{1-k}, \; k=0,1 \; (0 < p < 1)

也可以写成:

X 0 1
p_k 1-p p

(0-1)分布的概率分布很明确,无需SAS处理。

2.2. 二项分布 (Binomial Distribution)

二项分布,描述在n次独立试验中,某事件发生k次的概率。每次试验只有两种可能结果:成功或失败,成功的概率为p,失败的概率为1-p

二项分布的概率质量函数为:

P(X = k) = \binom{n}{k} p^k (1 - p)^{n - k}, \quad 0 < p < 1, \quad k = 0, \, 1, \, \ldots, \, n.

SAS函数probbnml返回二项分布的概率,在成功概率为0.5的试验中, 重复10次试验,4次成功的概率的SAS程序,计算如下:

data binomial;
    p = 0.5; 
    n = 10;    
    k = 4;    
    probability = probbnml(p, n, k);  
run;

2.3. 泊松分布 (Poisson Distribution)

泊松分布,通常用于描述在固定时间间隔或空间区域内,某事件发生的次数。泊松分布假设事件发生的概率是恒定的,并且事件是独立发生的。

设随机变量X的取值为0, 1, 2, \ldots, n, \ldots, 相应的概率质量函数为

P(X = k) = \frac{\lambda^k e^{-\lambda}}{k!}, \quad \lambda > 0, \quad k = 0, \, 1, \, 2, \, \ldots, \, n, \, \ldots

称随机变量X服从参数为\lambda的泊松分布,记为 X \sim P(\lambda)

其中,X是在固定时间间隔或空间区域内发生的事件次数;k是事件发生的具体次数(0, 1, 2, ...);\lambda是单位时间间隔或空间区域内事件的平均发生率;e是自然常数,约等于2.71828。

SAS函数poisson返回泊松分布的概率,在平均发生率为3的情况下,事件发生2次的概率的SAS程序如下:

data poisson;
    lambda = 3;  
    k = 2;       
    probability = poisson(lambda, k);  
run;

3. 连续型随机变量

概率密度函数

如果对于随机变量X的分布函数F(x), 存在非负函数f(x), 使对于任意实数x

F(x) = \int_{-\infty}^{x} f(t) \, dt

则称X为连续型随机变量,其中函数f(x)X的概率密度函数,简称概率密度。即,通过概率密度函数的定积分来求得分布函数值

于是有,

P\{x_1 < X \leq x_2\} = F(x_2) - F(x_1) = \int_{x_1}^{x_2} f(x) \, dx

由定义可知,f(x)个别点的函数值不影响分布函数F(x)的取值,因为点的定积分值为0。

下面介绍下几个常见连续型随机变量分布的概率密度。

3.1. 正态分布 (Normal Distribution)

正态分布,又称为高斯分布 (Gaussian Distribution),是一种连续型概率分布,在许多自然现象和统计模型中广泛应用。正态分布的概率密度函数呈钟形曲线,对称且集中在均值附近。

正态分布的概率密度为,

f(x) = \frac{1}{\sqrt{2\pi\sigma}} e^{-\frac{(x-\mu)^2}{2\sigma^2}}, \; -\infty < x < \infty

一般对于正态分布的随机变量,可以通过线性变换,将其转化为标准正态分布:

分布函数是对概率密度的定积分,SAS中可以通过函数probnorm计算标准正态分布的分布函数值。例如,对于标准正态分布随机变量Z,获取 P\{Z \leq -1.96\}P\{Z > 1.96\}P\{-1.96 < Z \leq 1.96\}程序如下:

data normal;
  a = probnorm(-1.96);
  b = 1 - probnorm(1.96);
  c = probnorm(1.96) - probnorm(-1.96);
run;

其含义是,对于服从标准正态分布的随机变量,在一次抽样中,其值≤-1.96的概率为0.025,>1.96的概率为0.025,在两者之间的概率为0.095。

-1.96又称为标准正态分布的0.025分位数。

一般地,当X \sim N(0, 1)时,满足概率表达式P(X \leq u_\alpha) = \alphau_\alpha称为标准正态分布的\alpha 分位数。

标准正态分布常见的分位数有:

\alpha 0.9 0.95 0.975 0.99 0.995 0.999
u_\alpha 1.282 1.645 1.96 2.326 2.576 3.090

3.2. 卡方分布 (Chi-Square Distribution)

卡方分布,用于描述从标准正态分布中独立抽取的平方和的分布。

假设Z_1, Z_2, \ldots, Z_kk个相互独立的标准正态分布随机变量,即Z_i \sim N(0, 1),则这些随机变量的平方和:X = \sum_{i=1}^{k} Z_i^2服从自由度为 k 的卡方分布,记作 X \sim \chi^2(k)

卡方分布对应的概率密度如下,n为卡方分布的自由度,

f(y) = \begin{cases} \frac{1}{2^{n/2} \Gamma(n/2)} y^{n/2-1} e^{-y/2}, & \text{if } y > 0, \\ 0, & \text{otherwise}. \end{cases}

卡方分布的分布函数值,可以通过SAS函数probchi获取,获取自由度为10的卡方分布分布函数F(2) = P(X \leq 2),程序如下:

data chi;
    chi = 2; 
    df = 10;    
    probability = probchi(chi, df);  
run;

3.3. t分布 (Student's t-Distribution)

t分布,也称为学生t分布,是一种在样本量较小或总体标准差未知的情况下,用于进行假设检验和构建置信区间的概率分布。t分布由威廉·戈塞特 (William Sealy Gosset) 以“学生” (Student) 的笔名发表,因此得名。

设随机变量XY相互独立,且X \sim N(0, 1)Y \sim \chi^2(n),则称t = \frac{X}{\sqrt{Y/n}}服从自由度为nt分布,记为t \sim t(n)

t(n)分布的概率密度函数如下,Γ是伽玛函数,n是自由度

f(x) = \frac{\Gamma\left(\frac{n+1}{2}\right)}{\sqrt{n\pi} \, \Gamma\left(\frac{n}{2}\right)} \left(1 + \frac{x^2}{n}\right)^{-\frac{n+1}{2}}, \quad x \in \mathbb{R}.

t分布的分布函数值,可以通过SAS函数probt获取,获取自由度为10的t分布分布函数F(2) = P(X \leq 2),程序如下:

data t;
    t = 2; 
    df = 10;    
    probability = probt(t, df);  
run;

3.4. F分布 (F-Distribution)

F分布,用于检验两个样本的方差是否相等,在方差分析 (ANOVA)、回归分析和假设检验中广泛应用。

设随机变量XY相互独立,X \sim \chi^2(m)Y \sim \chi^2(n),则称 F = \frac{X/m}{Y/n} 服从自由度为 (m, n)F分布,记为 F \sim F(m, n),其中m称为第一自由度,n称为第二自由度。

F(m, n)分布的概率密度函数为

f(y) = \begin{cases} \frac{\Gamma\left(\frac{m+n}{2}\right)}{\Gamma\left(\frac{m}{2}\right) \Gamma\left(\frac{n}{2}\right)} \left(\frac{m}{n}\right)^{m/2} y^{m/2-1} \left(1 + \frac{my}{n}\right)^{-(m+n)/2}, & y > 0, \\ 0, & \text{其他}. \end{cases}

F分布的分布函数值,可以通过SAS函数probf获取,获取自由度为(10, 10)的分布函数,程序如下:

data f;
    f = 2; 
    df1 = 10;    
    df2 = 10;    
    probability = probf(f, df1,df2);  
run;

总结

文章介绍了随机变量的概率分布、累积分布函数以及分位数的概念,列举了常见的离散型与连续型随机变量概率分布,以及对应概率计算的SAS程序。

对于离散型随机变量,概率分布描述变量取某一可能值的概率;对于连续型随机变量,概率分布描述变量取某一范围内值的概率。

感谢阅读, 欢迎关注:SAS茶谈!
若有疑问,欢迎评论交流!

梳理不易,转载请注明出处 (by Jihai / SAS茶谈)

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,417评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,921评论 3 387
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,850评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,945评论 1 285
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,069评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,188评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,239评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,994评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,409评论 1 304
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,735评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,898评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,578评论 4 336
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,205评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,916评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,156评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,722评论 2 363
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,781评论 2 351

推荐阅读更多精彩内容