一般p维正态随机向量的二次型

ops/2024/12/19 15:50:13/

内容来源
应用多元统计分析 北京大学出版社 高惠璇编著


1

X ∼ N p ( μ , Σ ) , Σ > 0 X\sim N_p(\mu,\Sigma),\Sigma>0 XNp(μ,Σ),Σ>0

X ′ Σ − 1 X ∼ χ 2 ( p , δ ) X'\Sigma^{-1}X\sim\chi^2(p,\delta) XΣ1Xχ2(p,δ),其中 δ = μ ′ Σ − 1 μ \delta=\mu'\Sigma^{-1}\mu δ=μΣ1μ

证明

因为 Σ > 0 \Sigma>0 Σ>0,所以存在可逆矩阵 C C C

Σ = C C ′ \Sigma=CC' Σ=CC

Y = C − 1 X Y=C^{-1}X Y=C1X,则

Y ∼ N p ( C − 1 μ , C − 1 Σ ( C − 1 ) ′ ) = N p ( C − 1 μ , I ) Y\sim N_p(C^{-1}\mu,C^{-1}\Sigma(C^{-1})')= N_p(C^{-1}\mu,I) YNp(C1μ,C1Σ(C1))=Np(C1μ,I)

X ′ Σ − 1 X = Y ′ C ′ Σ − 1 C Y = Y ′ Y ∼ χ 2 ( p , δ ) X'\Sigma^{-1}X=Y'C'\Sigma^{-1}CY=Y'Y\sim\chi^2(p,\delta) XΣ1X=YCΣ1CY=YYχ2(p,δ)

其中 δ = ( C − 1 μ ) ′ ( C − 1 μ ) = μ ′ Σ − 1 μ \delta=(C^{-1}\mu)'(C^{-1}\mu)=\mu'\Sigma^{-1}\mu δ=(C1μ)(C1μ)=μΣ1μ

2

X ∼ N p ( μ , Σ ) , Σ > 0 X\sim N_p(\mu,\Sigma),\Sigma>0 XNp(μ,Σ),Σ>0

A A A 为对称矩阵, r a n k ( A ) = r rank(A)=r rank(A)=r

( X − μ ) ′ A ( X − μ ) ∼ χ 2 ( r ) ⟺ Σ A Σ A Σ = Σ A Σ (X-\mu)'A(X-\mu)\sim\chi^2(r)\Longleftrightarrow \Sigma A\Sigma A\Sigma=\Sigma A\Sigma (Xμ)A(Xμ)χ2(r)ΣAΣAΣ=ΣAΣ

证明

因为 Σ > 0 \Sigma>0 Σ>0

所以它的平方根矩阵 Σ 1 / 2 > 0 \Sigma^{1/2}>0 Σ1/2>0

所以存在正交矩阵 Γ \Gamma Γ 使得

Σ 1 / 2 = Γ d i a g ( λ 1 , ⋯ , λ r ) Γ ′ \Sigma^{1/2}= \Gamma diag\Big(\sqrt{\lambda_1},\cdots,\sqrt{\lambda_r}\Big)\Gamma' Σ1/2=Γdiag(λ1 ,,λr )Γ

其中 λ i \lambda_i λi Σ \Sigma Σ 的特征值

Y = Σ − 1 / 2 ( X − μ ) ∼ N p ( 0 , I ) Y=\Sigma^{-1/2}(X-\mu)\sim N_p(0,I) Y=Σ1/2(Xμ)Np(0,I),则

( X − μ ) ′ A ( X − μ ) = Y ′ Σ 1 / 2 A Σ 1 / 2 Y (X-\mu)'A(X-\mu)=Y'\Sigma^{1/2} A\Sigma^{1/2}Y (Xμ)A(Xμ)=YΣ1/2AΣ1/2Y

C = Σ 1 / 2 A Σ 1 / 2 C=\Sigma^{1/2} A\Sigma^{1/2} C=Σ1/2AΣ1/2,结合上篇笔记-分量独立的n维正态随机向量的二次型的结论有

Y ′ C Y ∼ χ 2 ( p ) ⟺ C 2 = C Y'CY\sim\chi^2(p)\Longleftrightarrow C^2=C YCYχ2(p)C2=C

Σ 1 / 2 A Σ 1 / 2 ⋅ Σ 1 / 2 A Σ 1 / 2 = Σ 1 / 2 A Σ 1 / 2 \Sigma^{1/2} A\Sigma^{1/2}\cdot\Sigma^{1/2} A\Sigma^{1/2}=\Sigma^{1/2} A\Sigma^{1/2} Σ1/2AΣ1/2Σ1/2AΣ1/2=Σ1/2AΣ1/2

3

设 X ∼ N p ( μ , Σ ) , Σ > 0 设 X\sim N_p(\mu,\Sigma),\Sigma>0 XNp(μ,Σ),Σ>0

A , B A,B A,B p p p 阶非对称矩阵,则

( X − μ ) ′ A ( X − μ ) 与 ( X − μ ) ′ B ( X − μ ) 独立 ⟺ Σ A Σ B Σ = 0 (X-\mu)'A(X-\mu)与(X-\mu)'B(X-\mu)独立\\\Longleftrightarrow \Sigma A\Sigma B\Sigma=0 (Xμ)A(Xμ)(Xμ)B(Xμ)独立ΣAΣBΣ=0


http://www.ppmy.cn/ops/143216.html

相关文章

【AutoDL】通过【SSH远程连接】【vscode】

小帅碎碎念 0. 起因1. SSH信息获取2. 给你的vscode安装支持SSH远程连接的插件3. SSH远程连接入口4. 输入密码登陆5. 总结 0. 起因 之前使用AutoDL和Jupyter进行代码编辑和执行确实很方便,尤其是对于交互式数据分析项目。然而,也存在一些限制和不便之处&…

RabbitMQ安装延迟消息插件(mq报错)

之前启动一个springboot的单体项目,一直mq的错误,即便我更新了最新版本的mq,还是报错。 后来才发现,项目使用了延时队列,是需要单独下载延时插件的。 1如果判断mq有没有延时队列插件【没有x-delayed-message】 2下载…

生活小妙招之UE CaptureRT改

需求,四个不同的相机拍摄结果同屏分屏显示 一般的想法是四个Capture拍四张RT,然后最后在面片/UI上组合。这样的开销是创建4张RT,材质中采样4次RT。 以更省的角度,想要对以上流程做优化,4个相机拍摄是必须的&#xff…

Spring Mvc面试题(常见)

1 Spring MVC的执行流程 用户发起请求,请求先被Servlet拦截以后,转发给SpringMVC框架SpringMVC 里面的DispatcherServlet(核心控制器) 接收到请求,并转发给HandlerMappingHandlerMapping负责解析请求,根据请求信息和配置信息找到匹配的Controller类(当这里有配置拦截器,会…

Vue Web开发(十)

1. 用户管理新增,搜索,编辑,删除 本节课完成用户列表表单设计,使用table组件,同样模块化组件,CommonTable.vue组件,并且在User页面中引入,mock实现数据模拟,最终完成用户…

golang操作sqlite3加速本地结构化数据查询

目录 摘要Sqlite3SQLite 命令SQLite 语法SQLite 数据类型列亲和类型——优先选择机制 SQLite 创建数据库SQLite 附加数据库SQLite 分离数据库 SQLite 创建表SQLite 删除表 SQLite Insert 语句SQLite Select 语句SQLite 运算符SQLite 算术运算符SQLite 比较运算符SQLite 逻辑运算…

基于单片机的蔬菜大棚环境监测系统

本设计以单片机作为主控,通过温湿度检测模块监测、二氧化碳浓度采集模块、光照强度监测大棚中的环境数据,可以通过按键模块设置阈值,当超过阈值后将通过声光报警模块提示人们此时环境并不适合作物生长,并将以上监测数据通过液晶显…

C# WinForm移除非法字符的输入框

C# WinForm移除非法字符的输入框 文章目录 namespace System.Windows.Forms {using System.ComponentModel;/// <summary>/// 支持移除 非法字符 的输入框。/// </summary>public class RemoveInvalidCharTextBox : TextBox{/// <summary>/// 测试代码&#…