原文来师兄的博客:http://blog.csdn.net/wjj5881005/article/details/53320403

  • 多元正态分布
  • 多元正态分布的条件密度

多元正态分布

多元正态分布的密度函数如下 :

fx(x1,...xn)=1(2π)k√|Σ|1/2exp(−12(x−μ)TΣ−1(x−μ))f_{x}(x_{1},...x_{n})=\frac{1}{\sqrt{(2\pi)^{k}}|\Sigma|^{1/2}}exp(-\frac{1}{2}(x-\mu)^{T}\Sigma^{-1}(x-\mu)) (1)
其对应的矩母函数(也有称动差函数)为 exp(μTt+12tTΣt)exp(\mu^{T}t+\frac{1}{2}t^{T}\Sigma t)。事实上,如果随机向量 [X1,...Xn][X_{1},...X_{n}]满足上面的动差函数,那么我们就称随机向量 [X1,...Xn][X_{1},...X_{n}]服从多元高斯分布。具体地证明可以看 这里。

多元正态分布的条件密度

令随机向量[X1,...Xn][X_{1},...X_{n}]服从多元高斯分布。我们可以推导XnX_{n}在给定X1,...Xn−1X_{1},...X_{n-1}的情况下的条件密度分布:

f(xn|x1,...,xn−1)=f(x1,...,xn−1,xn)f(x1,...,xn−1)f(x_{n}|x_{1},...,x_{n-1})=\frac{f(x_{1},...,x_{n-1},x_{n})}{f(x_{1},...,x_{n-1})} (2),
其中 f(x1,...,xn)=(2π)−n/2(|Σ|−1/2)exp[−12∑ni,j=1yiqijyj]f(x_{1},...,x_{n})=(2\pi)^{-n/2}(|\Sigma|^{-1/2})exp[-\frac{1}{2}\sum_{i,j=1}^{n}y_{i}q_{ij}y_{j}] (3)
其中 Q=Σ−1=[qij],yi=xi−μiQ=\Sigma^{-1}=[q_{ij}],y_{i}=x_{i}-\mu_{i}。同样地,
f(x1,...,xn−1)=∫∞∞f(x1,...,xn−1,xn)dxn=B(y1,...,yn−1)f(x_{1},...,x_{n-1})=\int_\infty^\infty {f(x_{1},...,x_{n-1},x_{n})}\,dx_{n}=B(y_{1},...,y_{n-1}) (4).
现在我们将公式(3)中的求和项进行分解,有:
∑ni,j=1yiqijyj=∑n−1i,j=1yiqijyj+yn∑n−1j=1qnjyj+yn∑n−1i=1qinyj+qnny2n\sum_{i,j=1}^{n}y_{i}q_{ij}y_{j}=\sum_{i,j=1}^{n-1}y_{i}q_{ij}y_{j}+y_{n}\sum_{j=1}^{n-1}q_{nj}y_{j}+y_{n}\sum_{i=1}^{n-1}q_{in}y_{j}+q_{nn}y_{n}^{2}(5)
因此,最终地条件分布具有如下的形式:
A(y1,...,yn−1)B(y1,...,yn−1)exp[−(Cy2n+D(y1,...,yn−1)yn)]\frac{A(y_{1},...,y_{n-1})}{B(y_{1},...,y_{n-1})}exp[-(Cy_{n}^{2}+D(y_{1},...,y_{n-1})y_{n})] (6)
其中 C=(1/2)qnnC=(1/2)q_{nn},因为 Q=Σ−1Q=\Sigma^{-1}是对称矩阵,所以 D=∑n−1j=1qnjyj=∑n−1i=1qinyiD=\sum_{j=1}^{n-1}q_{nj}y_{j}=\sum_{i=1}^{n-1}q_{in}y_{i}.(6)式又可以进一步表示称如下的式子:
[ABexp(DD24C)]exp[−(yn+D2C)2]1C[\frac{A}{B}exp(\frac{DD^{2}}{4C})]exp[-\frac{(y_{n}+\frac{D}{2C})^{2}]}{\frac{1}{C}} (7)
从公式(7)很容看出 xnx_{n}的条件密度函数是服从正态分布的。
所以条件分布的方差为: 2Var(Xn|X1,...,Xn−1)=1/C2Var(X_{n}|X_{1},...,X_{n-1})=1/C,进一步有: Var(Xn|X1,...,Xn−1)=12C=1qnnVar(X_{n}|X_{1},...,X_{n-1})=\frac{1}{2C}=\frac{1}{q_{nn}}
均值为:
E(Xn|X1,...,Xn−1)=μn−D2C=μn−1qnn∑n−1j=1qnj(Xj−μj)E(X_{n}|X_{1},...,X_{n-1})=\mu_{n}-\frac{D}{2C}=\mu_{n}-\frac{1}{q_{nn}}\sum_{j=1}^{n-1}q_{nj}(X_{j}-\mu_{j})
这就说明了再抽样多元正态分布时,如果已知了其它维度的随机变量值, 剩下的那个维度的随机变量也是服从正态分布

关于多元正态分布的条件概率密度相关推荐

  1. 多元正太分布条件密度

    原文来自师兄的博文:http://blog.csdn.net/wjj5881005/article/details/53320403 多元正态分布 多元正态分布的条件密度 多元正态分布 多元正态分布的 ...

  2. UA MATH564 概率论VI 数理统计基础2 多元正态分布

    UA MATH564 概率论VI 数理统计基础2 多元正态分布 矩母函数 概率密度 多元正态分布的矩 条件分布 独立性 抽样分布简单地说就是统计量服从的分布,正态分布时最常用的总体分布,因此研究正态总 ...

  3. 【ML学习笔记】17:多元正态分布下极大似然估计最小错误率贝叶斯决策

    简述多元正态分布下的最小错误率贝叶斯 如果特征的值向量服从d元正态分布,即其概率密度函数为: 即其分布可以由均值向量和对称的协方差矩阵 唯一确定. 如果认为样本的特征向量在类内服从多元正态分布: 即对 ...

  4. 【应用多元统计分析】CH3 多元正态分布

    目录 一.多元正态分布的定义 1.定义 2.二元正态分布 二.多元正态分布的性质 [property1*] [property2] [property3] [property4] [property5 ...

  5. 透彻理解多元正态分布

    文章目录 前言 多元标准正态分布 多元一般正态分布 缩放与位移相同尺度 缩放与位移不同尺度 旋转变换 多元正态分布的概率密度函数 多元正态分布的性质 本篇内容主要是对于基本书籍教材多元正态分布相关章节 ...

  6. 两个多元正态分布的KL散度、巴氏距离和W距离

    ©PaperWeekly 原创 · 作者 | 苏剑林 单位 | 追一科技 研究方向 | NLP.神经网络 正态分布是最常见的连续型概率分布之一.它是给定均值和协方差后的最大熵分布(参考<&quo ...

  7. 【多元统计分析】02.多元正态分布

    文章目录 二.多元正态分布 1.多元正态分布的定义 2.多元正态分布分量独立性 3.多元正态分布的条件分布 4.最佳预测 回顾总结 二.多元正态分布 1.多元正态分布的定义 由大数定律,自然界中许多随 ...

  8. R语言学习——一元与多元正态分布检验(也可以用于其他分布的检验)

    文章目录 1 一元正态的评估 1.1 图像法 1.1.1 直方图 1.1.2 Q-Q图 1.2 峰度和偏度 1.3 统计检验 1.3.1 Shapiro-Wilks检验 1.3.2 Kolmogoro ...

  9. 多元正态分布(Multivariate normal distribution)

    多元正态分布(Multivariate normal distribution) 前言 我们通常讨论正态分布都是在一元(univariate)的情况下,相信下面的定义大家都很熟悉了:假设随机变量XXX ...

  10. 协方差矩阵与多元正态分布

    文章目录 协方差矩阵 协方差 协方差矩阵 多元正态分布 协方差矩阵的特征值分解 协方差矩阵 协方差 在统计学中,方差用来度量单个随机变量的离散程度,而协方差用来刻画两个随机变量的相似程度,方差的计算公 ...

最新文章

  1. 如何从文件的完整路径获取目录?
  2. [攻防世界 pwn]——monkey
  3. iOS 在UILabel显示不同的字体和颜色(转)
  4. oracle 怎么调整emctl,emca和emctl的日常用法及实例(经验保存)
  5. 什么是ie浏览器_?IE 浏览器为什么不招人待见?
  6. (2/2)Canvas的交互存为图片-爬坑篇
  7. 小型机oracle命令,IBM 小型机基本知识及AIX操作系统基本命令.doc
  8. MYSQL 更新字段,向字段追加字符串
  9. (十五)洞悉linux下的Netfilteriptables:开发自己的hook函数【实战】
  10. 计算机接口中的shl,西安交通大学19年5月补考《微机原理与接口技术(高起专)》作业考核试题【参考答案】...
  11. MATLAB疲劳检测系统
  12. mysql怎么设主键和外键_MySQL基础之 主键外键设置
  13. 遭遇nat.exe,socks.exe,USP10.dll,BOSC.dll,kb080387.CNT,~ctwxw.txt等1
  14. Web Api Token验证
  15. AT32 XMC驱动PC卡/CF卡
  16. 方维团购V3.07版本短信接口开发
  17. 开上新能源车之后,如何摆脱“充电焦虑”?
  18. vue.js的快速入门使用
  19. 最详细的Android图片压缩攻略
  20. 表达式转换(中缀转后缀)

热门文章

  1. C#UDP广域网,局域网通信-原理分析
  2. 学习linux方向,学习linux方向
  3. 汉诺塔_-Chaz-_新浪博客
  4. word 参考文献插入整理
  5. 医学病理图像:细胞间质与间质细胞的区别
  6. 官方太空射击游戏总结
  7. 企业联合体的形式_企业组织形式的概念与类型怎么分?
  8. Android仿人人客户端(v5.7.1)——项目框架新做的调整描述(项目中基类java源码)
  9. 交换机及IP地址(概念)
  10. 什么是磁盘?磁盘的组成?接口和分区?