【概率论】3-7:多变量分布(Multivariate Distributions Part II)

原文地址1:https://www.face2ai.com/Math-Probability-3-7-Multivariate-Distributions-P2转载请标明出处

Abstract: 本文继续上文的内容,讲解多变量分布的条件分布,全概率公式和贝叶斯公式。提出直方图的概念。
Keywords: Conditional Distributions,Law of Total Probability,Bayes’ Theorem,Histograms

多变量分布

向往自由的人更懂得约束自己,约束自己的行为,自由自己的思想,而我们现在周围的人更多的是反过来的,约束自己的思想,告诉自己不许瞎想,但是自由自己的行为,想干嘛干嘛,毕竟改革开放的红利使得大部分人从温饱过渡到有车有房还可以随意出门旅游看看世界了,而思维思想,包括知识,还停留在改革开放之前,这就导致了很诡异的一系列行为–有钱没文化。
本文继续上文,将前面的内容扩展至多变量的联合分布,最后引出关于直方图的一些知识。

多变量条件分布 Conditional Distributions

条件概率,是3.6中讲解的内容,我们反复的强调所有概率都是条件概率,所有分布也都是条件分布,所以与其说条件概率(分布)与一般的概率(分布)行为一致,不如说所有的所有的行为都是在条件概率的基础上进行的,只是我们省略了某些必然达成的条件。
假设 n n n 个随机变量 X 1 , … , X n X_1,\dots,X_n X1,,Xn 有一个连续的联合分布,其联合分布式p.d.f. 是 f f f 并且 f 0 f_0 f0 定义了其中 k < n k < n k<n 个随机变量的边缘分布有 f 0 ( x 1 , … , x n ) > 0 f_0(x_1,\dots,x_n) > 0 f0(x1,,xn)>0 那么条件分布,当条件 X 1 = x 1 , … X n = x n X_1=x_1,\dots X_n=x_n X1=x1,Xn=xn 给定时 ( x k + 1 , … , x n ) (x_{k+1},\dots,x_n) (xk+1,,xn) 的p.d.f.是:

g k + 1 , … , x n ( x k + 1 , … , x n ∣ x 1 , … , x k ) = f ( x 1 , … , x n ) f 0 ( x 1 , … , x k ) g_{k+1,\dots,x_n}(x_{k+1},\dots,x_{n}|x_1,\dots,x_k)=\frac{f(x_1,\dots ,x_n)}{f_0(x_1,\dots,x_k)} gk+1,,xn(xk+1,,xnx1,,xk)=f0(x1,,xk)f(x1,,xn)

Definition 3.3.7 Conditional p.f. or p.d.f. Suppose that the random vector X ⃗ = ( X 1 , … , X n ) \vec{X}=(X_1,\dots,X_n) X =(X1,,Xn) is divided into two subvectors Y ⃗ \vec{Y} Y and Z ⃗ \vec{Z} Z ,where Y ⃗ \vec{Y} Y is a k-dimensional random vector comprising k k k of the n n n random variables in X ⃗ \vec{X} X and Z ⃗ \vec{Z} Z is an ( n − k ) (n-k) (nk)-dimensional random vector comprising the other n − k n-k nk random variables in X ⃗ \vec{X} X .Suppose also that the n n n-dimensional joint p.f. or p.d.f. of ( Y ⃗ , Z ⃗ ) (\vec{Y},\vec{Z}) (Y ,Z ) is f f f and that the marginal ( n − k ) (n-k) (nk)-dimensional joint p.f. ,p.d.f. or p.f./p.d.f. of Z ⃗ \vec{Z} Z is f 2 f_2 f2 .Then for every given point z ∈ R n − k z\in\mathbb{R}^{n-k} zRnk such that f 2 ( z ) > 0 f_2(z)>0 f2(z)>0 ,the conditional k k k-dimensional p.f. p.d.f.or p.f./p.d.f. g 1 g_1 g1 of Y ⃗ \vec{Y} Y given Z ⃗ = z ⃗ \vec{Z}=\vec{z} Z =z is defined as follows:
g 1 ( y ⃗ ∣ z ⃗ ) = f ( y ⃗ , z ⃗ ) f 2 ( z ⃗ )  for  y ⃗ ∈ R k g_1(\vec{y}|\vec{z})=\frac{f(\vec{y},\vec{z})}{f_2(\vec{z})} \text{ for }\vec{y}\in \mathbb{R}^k g1(y z )=f2(z )f(y ,z ) for y Rk

这就是完整的定义,抄一遍下来还真写了不少字,但是整个思路很清晰,首先就是把多变量形成向量的形式,再把向量拆成两个小的空间当然上面这个定义的公式也可以写成:

f ( y ⃗ , z ⃗ ) = g 1 ( y ⃗ ∣ z ⃗ ) f 2 ( z ⃗ ) f(\vec{y},\vec{z})=g_1(\vec{y}|\vec{z}) f_2(\vec{z}) f(y ,z )=g1(y z )f2(z )

写成乘法原理的形式,解决分母是0的尴尬局面。也可以看出,通过条件分布和边缘分布得到联合分布的方法是正确的(乘法原理的正确性)。

多变量全概率公式,贝叶斯定理 Law of Total Probability & Bayes’ Theorem

条件概率和乘法原则得到证明后接着就可以扩展到全概率公式和贝叶斯公式:

Theorem Multivariate Law of Total Probability and Bayes’ Theorem .Assume the conditions and notation given in Definition 3.3.7 .If Z ⃗ \vec{Z} Z has a continuous joint distribution ,the marginal p.d.f. of Y ⃗ \vec{Y} Y is
f 1 ( y ⃗ ) = ∫ − ∞ ∞  ⁣ ⋯ ∫ − ∞ ∞ ⏟ n − k g 1 ( y ⃗ ∣ z ⃗ ) f 2 ( z ⃗ ) d z ⃗ f_1(\vec{y})=\underbrace{\int^{\infty}_{-\infty} \dots \int^{\infty}_{-\infty}}_{n-k}g_1(\vec{y}|\vec{z})f_2(\vec{z})d\vec{z} f1(y )=nk g1(y z )f2(z )dz

and the conditional p.d.f. of Z ⃗ \vec{Z} Z given Y ⃗ = y ⃗ \vec{Y}=\vec{y} Y =y is
g 2 ( z ⃗ ∣ y ⃗ ) = g 1 ( y ⃗ ∣ z ⃗ ) f 2 ( z ⃗ ) f 1 ( y ⃗ ) g_2(\vec{z}|\vec{y})=\frac{g_1(\vec{y}|\vec{z})f_2(\vec{z})}{f_1(\vec{y})} g2(z y )=f1(y )g1(y z )f2(z )

上面的积分如果在离散概率情况要变成求和。

条件独立多随机变量 Conditionally Independent Random Variables

当我们定义某两个或某几个随机变量独立的时候,是看其概率分布的关系是否满足乘法关系,但是当我们说某两个或几个随机变量条件独立的时候,条件一定是同一个条件,比如 x ⃗ \vec{x} x 在条件 z ⃗ \vec{z} z 下和 y ⃗ \vec{y} y 相互独立。 意味着 f ( x ⃗ , y ⃗ ∣ z ⃗ ) = f 1 ( x ⃗ ∣ z ⃗ ) f 2 ( y ⃗ ∣ z ⃗ ) f(\vec{x},\vec{y}|\vec{z})=f_1(\vec{x}|\vec{z})f_2(\vec{y}|\vec{z}) f(x ,y z )=f1(x z )f2(y z ) ,不可能出现这种情况: x ⃗ \vec{x} x 在条件 z ⃗ \vec{z} z 下和 y ⃗ \vec{y} y 在条件 w ⃗ \vec{w} w 下相互独立。

Definition Conditionally Independent Random Variables.Let Z be a random vector with joint p.f.,p.d.f. or p.f./p.d.f. f_0(\vec{z}) .Several random variables X_1,\dots,X_n are conditionally independent given Z if,for all z such that f_0(z)>0 we have
g ( x ⃗ ∣ z ⃗ ) = ∏ i = 1 n g i x i ∣ z ⃗ ) g(\vec{x}|\vec{z})=\prod_{i=1}^ng_i{x_i}|\vec{z}) g(x z )=i=1ngixiz )

注意上面公式中 g ( x ⃗ ∣ z ⃗ ) g(\vec{x}|\vec{z}) g(x z ) 是多变量的条件随机分布,而后面的乘数 g i x i ∣ z ⃗ ) g_i{x_i}|\vec{z}) gixiz ) 条件是多变量的而 x i x_i xi 是单变量

进一步扩展,我们现在已经基本知道了前面所有分布定理的条件版本,但是只有一个条件,当条件是多个的时候,我们进行相同的操作,比如多条件的贝叶斯公式:
g 2 ( z ⃗ ∣ y ⃗ , w ⃗ ) = g 1 ( y ⃗ ∣ z ⃗ , w ⃗ ) f 2 ( z ⃗ ∣ w ⃗ ) f 1 ( y ⃗ ∣ w ⃗ ) g_2(\vec{z}|\vec{y},\vec{w})=\frac{g_1(\vec{y}|\vec{z},\vec{w})f_2(\vec{z}|\vec{w})}{f_1(\vec{y}|\vec{w})} g2(z y ,w )=f1(y w )g1(y z ,w )f2(z w )

还是前面说的,所有分布都是条件的,独立性也一样,独立是条件独立的特殊情况。
当我们证明某个条件随机变量或者某些条件独立同分布的随机变量满足某一结论的时候,那么其非条件形式也满足!

直方图 Histograms

直方图用处很多,Excel中,高考试卷上,图像处理里面,直方图是我们收集到的数据的直观显示,直方图对于显示条件独立同分布的随机变量是最有效的,换句话说,如果收集到的数据不是(条件)独立同分布的,那么显示出来没有意义,因为我们抽样的时候都是假设他们独立同分布的i.i.d

Definition Histogram. Let x 1 , … , x n x_1,\dots,x_n x1,,xn be a collection of numbers that all lie between two values a < b a < b a<b That is a ≤ x i ≤ b a\leq x_i \leq b axib for all i = 1 , … , n i=1,\dots,n i=1,,n Choose some integer k ≥ 1 k\geq 1 k1 and divide the interval [ a , b ] [a,b] [a,b] into k k k equal-legth subintervals of length b − a k \frac{b-a}{k} kba .For each subinterval,count how many of the numbers x 1 , … , x n x_1,\dots,x_n x1,,xn are in the subinterval,Let c i c_i ci be the count for subinterval i i i for i = 1 , … , k i=1,\dots,k i=1,,k ,Choose a number r > 0 r>0 r>0 (Typically, r + 1 r+1 r+1 ,or r = n r=n r=n or r = n ( b − a ) k ) r=\frac{n(b-a)}{k}) r=kn(ba)) Draw a two-dimensional graph with the horizonal axis running from a a a to b b b .For each subinterval i = 1 , … , k i=1,\dots,k i=1,,k draw a rectangular bar of width $\frac{b-a}{k} and height equal to c i r \frac{c_i}{r} rci over the midpoint of the i i ith interval.such a graph is called a histogram.

上面是直方图的基本定义,interval被分成多个等长的subinterval 当然我们也可以开发出不等长的subinterval 然后用每个直方图内的方块的面积表示这个subinterval中的数量。

总结

今天我们就把随机变量一般化了,这更加复合实际应用,而且操作也不是很复杂。
待续。。。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值