充分统计量和因子分解定理
充分统计量
-
定义: 设样本 X X X的服从分布 f ( X ∣ θ ) f(X|\theta) f(X∣θ), θ ∈ Θ \theta\in\Theta θ∈Θ,设 T = T ( X ) T=T(X) T=T(X)为一统计量,若在已知 T T T的条件下,样本 X X X的条件分布与参数 θ \theta θ无关,则称 T = T ( X ) T=T(X) T=T(X)为 θ \theta θ的充分统计量
-
Example:
设 X = ( x 1 , x 2 , . . , x n ) X=(x_1,x_2,..,x_n) X=(x1,x2,..,xn)是从泊松分布 P ( λ ) P(\lambda) P(λ)中抽取的随机样本,下面将从定义出发证明 T ( X ) = ∑ i = 1 n x i T(X)=\sum_{i=1}^nx_i T(X)=∑i=1nxi是 θ \theta θ的充分统计量∵ x i ∼ P ( λ ) , ∴ ∑ i = 1 n x i ∼ P ( n λ ) \because x_i \sim P(\lambda),\therefore\sum_{i=1}^nx_i \sim P(n\lambda) ∵xi∼P(λ),∴∑i=1nxi∼P(nλ),我们将其记为 T ∼ P ( θ ) , θ = n λ T\sim P(\theta),\theta=n\lambda T∼P(θ),θ=nλ
由已知可得,样本的条件分布为 f ( X ∣ λ ) = ∏ i = 1 n e − λ λ x i x i ! = e − n λ λ ∑ i = 1 n x i ∏ i = 1 n x i ! = e − θ λ T ∏ i = 1 n x i ! f(X|\lambda)=\prod_{i=1}^n\frac{e^{-\lambda}\lambda^{x_i}}{x_i!}=\frac{e^{-n\lambda}\lambda^{\sum_{i=1}^nx_i}}{\prod_{i=1}^nx_i!}=\frac{e^{-\theta}\lambda^{T}}{\prod_{i=1}^nx_i!} f(X∣λ)=i=1∏nxi!e−λλxi=∏i=1nxi!e−nλλ∑i=1nxi=∏i=1nxi!e−θλT
此时样本 X X X的条件分布 f ( X ∣ λ ) f(X|\lambda) f(X∣λ)与参数 λ \lambda λ无关,因此 T ( X ) = ∑ i = 1 n x i T(X)=\sum_{i=1}^nx_i T(X)=∑i=1nxi是 θ \theta θ的充分统计量
因子分解定理
-
从定义出发证明充分统计量显得有些繁琐,因此我们引入因子分解定理
-
定义: 设样本 X = ( x 1 , x 2 , . . , x n ) X=(x_1,x_2,..,x_n) X=(x1,x2,..,xn)的条件分布为 f ( X ∣ θ ) f(X|\theta) f(X∣θ), θ ∈ Θ \theta\in\Theta θ∈Θ, T = T ( X ) T=T(X) T=T(X)为一统计量,则 T = T ( X ) T=T(X) T=T(X)是充分统计量的充分必要条件为条件分布为 f ( X ∣ θ ) f(X|\theta) f(X∣θ)可被分解为如下形式: f ( X ∣ θ ) = g ( T ( X ) , θ ) ⋅ h ( X ) f(X|\theta)=g(T(X),\theta)·h(X) f(X∣θ)=g(T(X),θ)⋅h(X)也就是可被分解为两部分,一部分仅与 T ( X ) T(X) T(X)和 θ \theta θ有关,另一部分为一个常数或仅与样本 X X X有关。
-
重要推论: 若 T = T ( X ) T=T(X) T=T(X)是充分统计量, S = g ( T ) S=g(T) S=g(T)是 T T T一一对应的变换,则 S S S也是 θ \theta θ的充分统计量
-
Example:
证明以下命题:设 X = ( x 1 , x 2 , . . , x n ) X=(x_1,x_2,..,x_n) X=(x1,x2,..,xn)为从正态总体 N ( a , σ 2 ) N(a,\sigma^2) N(a,σ2)中抽取的随机样本,令 θ = ( a , σ 2 ) \theta=(a,\sigma^2) θ=(a,σ2),则 T ( X ) = ( ∑ x i , ∑ x i 2 ) T(X)=(\sum{x_i},\sum{x_{i}^2}) T(X)=(∑xi,∑xi2)为充分统计量,且 ( X ‾ , S 2 ) (\overline{X},S^2) (X,S2)也是充分统计量,此处 X ‾ = 1 n ∑ x i , S 2 = 1 n − 1 ∑ ( x i − X ‾ ) 2 \overline{X}=\frac{1}{n}\sum{x_i},S^2=\frac{1}{n-1}\sum{(x_i-\overline{X})^2} X=n1∑xi,S2=n−11∑(xi−X)2由已知得,样本的条件分布为
f ( x ) = ( 1 2 π σ ) n exp ( − 1 2 σ 2 ∑ ( x i − a ) 2 ) = ( 1 2 π σ ) n exp ( − 1 2 σ 2 ( ∑ x i 2 − 2 a ∑ x i + n a 2 ) ) = g ( T ( X ) , θ ) ⋅ h ( X ) \begin{aligned} f(x) &= (\frac{1}{\sqrt{2\pi}\sigma})^n\exp(-\frac{1}{2\sigma^2}\sum{(x_i-a)^2}) \\ &=(\frac{1}{\sqrt{2\pi}\sigma})^n\exp(-\frac{1}{2\sigma^2}(\sum{x_i^2}-2a\sum{x_i}+na^2)) \\ &= g(T(X),\theta)·h(X) \end{aligned} f(x)=(2πσ1)nexp(−2σ21∑(xi−a)2)=(2πσ1)nexp(−2σ21(∑xi2−2a∑xi+na2))=g(T(X),θ)⋅h(X)
此处的 h ( X ) ≡ 1 h(X)\equiv1 h(X)≡1,至此, T ( X ) = ( ∑ x i , ∑ x i 2 ) T(X)=(\sum{x_i},\sum{x_{i}^2}) T(X)=(∑xi,∑xi2)为充分统计量得证,又因为 ( X ‾ , S 2 ) (\overline{X},S^2) (X,S2)为 T ( X ) = ( ∑ x i , ∑ x i 2 ) T(X)=(\sum{x_i},\sum{x_{i}^2}) T(X)=(∑xi,∑xi2)一一对应的变换,由推论可得, ( X ‾ , S 2 ) (\overline{X},S^2) (X,S2)也是充分统计量
理解:
- 充分统计量对于简化计算是有显著的帮助的
- 一一对应的变换可理解为一个函数
- 样本的条件分布其实就是样本似然
- 无论是从定义出发证明充分统计量,还是通过因子分解定理,都需要先求出样本的条件分布,然后再选择一种方法
- 从定义出发证明需要想方设法消除式子中原来的参数
相关文章:
充分统计量和因子分解定理
充分统计量 定义: 设样本 X X X的服从分布 f ( X ∣ θ ) f(X|\theta) f(X∣θ), θ ∈ Θ \theta\in\Theta θ∈Θ,设 T T ( X ) TT(X) TT(X)为一统计量,若在已知 T T T的条件下,样本 X X X的条件分布与参数 θ \the…...
M1 PD安装arm ubuntu及Docker
M1 PD安装arm ubuntu 下载 Ubuntu 22.04.2 LTS https://cn.ubuntu.com/download/server/arm 参考视频安装 https://www.bilibili.com/video/BV1Mu4y1f74v/?spm_id_from333.999.0.0&vd_source9056c6d3c91a117baaceb663957daa08 PD Ubuntu安装docker 删除现有的docker安装…...
TCP协议的RST标志
下文中的内容多数来自【参考】中的文章,这边进行一个整理和总结,后续会慢慢增加出现各个 RST 包的测试代码,便于理解。 TCP的 “断开连接” 标志 RST 标志 Reset,复位标志,用于非正常地关闭连接。它是 TCP 协议首部里…...
【软件质量与软件测试 白盒测试与黑盒测试】
第十章 黑盒测试 10.1 等价类划分: 10.1.1 划分等价类 等价类是指所有数据中的一组,它们具有相同的测试结果或相同的响应。等价类划分是将输入数据分为多个等价类的过程。 10.1.2 划分等价类的方法 划分等价类方法主要包括以下几种: 特…...
JavaScript教程(高级)
面向对象编程介绍 两大编程思想 (1)、 面向过程编程: (缩写 POP)( Process-oriented programming)面向过程就是分析出解决问题所需要的步骤,然后用函数把这些步骤一步一步实现&am…...
C++进阶 —— 范围for(C++11新特性)
目录 一,范围for介绍 二,范围for注意事项 一,范围for介绍 范围for(range-based for loop)是C11新引入的特性,可遍历各种序列结构的容器(如数组、vector、list等);每次循…...
ELK +Filebeat日志分析系统
一、 ELK日志分析系统概述 1、ELK简介 ELK是三个开源软件的缩写,分别表示:Elasticsearch , Logstash, Kibana , 它们都是开源软件。新增了一个FileBeat,它是一个轻量级的日志收集处理工具(Agent),Filebeat占用资源少,…...
万字解析PELT算法!
Linux是一个通用操作系统的内核,她的目标是星辰大海,上到网络服务器,下至嵌入式设备都能运行良好。做一款好的linux进程调度器是一项非常具有挑战性的任务,因为设计约束太多了: 它必须是公平的快速响应系统的throughp…...
腾讯云服务器端口怎么全开?教程来了
腾讯云服务器端口怎么全开?云服务器CVM在安全组中设置开通,轻量应用服务器在防火墙中设置,腾讯云百科来详细说下腾讯云服务器端口全开放教程: 目录 腾讯云服务器端口全部开通教程 云服务器CVM端口全开放教程 轻量应用服务器开…...
深入理解Java虚拟机:JVM高级特性与最佳实践-总结-13
深入理解Java虚拟机:JVM高级特性与最佳实践-总结-13 Java内存模型与线程Java内存模型原子性、可见性与有序性先行发生原则 Java内存模型与线程 Java内存模型 原子性、可见性与有序性 Java内存模型是围绕着在并发过程中如何处理原子性、可见性和有序性这三个特征来…...
租售keysight E8257D 50G模拟信号发生器 销售/回收
是德(Keysight) E8257D 模拟信号发生器 Keysight E8257D (Agilent) PSG 模拟信号发生器提供业界领先的输出功率、电平精度和高达 67 GHz 的相位噪声性能(工作频率可达 70 GHz)。Agilent PSG 模拟信号发生器的高输出功率和卓越的电…...
【C++】什么是函数模板/类模板?
文章目录 一、函数模板1.什么是函数模板?2.函数模板格式3.函数模板原理4.函数模板实例化(1)隐式实例化(2)显示实例化 二.类模板1.类模板定义格式2.类模板的实例化 总结 一、函数模板 1.什么是函数模板? 函…...
为什么是ChatGPT引发了AI浪潮?
目录 BERT和GPT简介 BERT和GPT核心差异 GPT的优势 GPT的劣势 总结 随着近期ChatGPT的火热,引发各行各业都开始讨论AI,以及AI可以如何应用到各个细分场景。为了不被时代“抛弃”,我也投入了相当的精力用于研究和探索。但在试验的过程中&…...
批处理文件(.bat)启动redis及任何软件(同理)
批处理文件 每次从文件根目录用配置文件格式来启动redis太麻烦了 可以在桌面上使用批处理文件(.bat)启动Redis,请按照以下步骤进行操作: 打开文本编辑器,如记事本。 在编辑器中输入以下内容: 将文件保存…...
深度学习求解稀疏最优控制问题的并行化算法
稀疏最优控制问题 问题改编自论文An FE-Inexact Heterogeneous ADMM for Elliptic Optimal Control Problems with L1-Control Cost { min y ( μ ) , u ( μ )...
牛客网项目—开发社区首页
视频连接:开发社区首页_哔哩哔哩_bilibili 代码地址:Community: msf begin 仿牛客论坛项目 (gitee.com) 本文是对仿牛客论坛项目的学习,学习本文之前需要了解Java开发的常用框架,例如SpringBoot、Mybatis等等。如果你也在学习牛…...
uniapp水文【uniapp】
文章目录 1、前言2、历史3、发展4、功能5、优缺点6、总结7、附录7.1、高频使用7.2、使用注意 1、前言 Uniapp是一种跨平台的移动应用开发框架,它允许开发者使用一套代码库,同时生成iOS、Android等多个平台的应用程序。这种技术方案可以大大降低开发成本…...
Java函数式接口
3 函数式接口 3.1 函数式接口概述 函数式接口:有且仅有一个抽象方法的接口 Java中的函数式编程体现就是Lambda表达式,所以函数式接口就是可以适用于Lambda使用的接口只有确保接口中有且仅有一个抽象方法, Java中的Lambda才能顺利地进行推导…...
安装libevent库
安装libevent库 yum install libevent libevent-devel 自动安装Memcached yum install memcached 源码安装 下载1.6.19版本 wget https://www.memcached.org/files/memcached-1.6.19.tar.gz (若证书过期yum install -y ca-certificates) 解压源码 tar -zxvf…...
vue 截取字符串的方法
vue中的字符串方法,我目前使用最多的是下面两种方法,因为 vue的字符串方法支持断言操作。 1、 vue中截取字符串的方法如下: 2、 vue中截取字符串的方法,这个方法也是需要依赖于 vue库提供的支持。 3、 vue中截取字符串的方法&…...
FPGA新手入门:用Verilog手搓一个交通灯控制器(附完整代码与仿真)
FPGA实战:从零构建智能交通灯控制系统的Verilog全流程指南 引言 第一次接触FPGA开发时,我被硬件描述语言的独特思维方式所吸引。与软件编程不同,Verilog让我们能够直接描述硬件电路的行为。交通灯控制系统作为数字电路设计的经典案例…...
KityMinder:可视化思维的协作引擎 | 高效工作者必备工具
KityMinder:可视化思维的协作引擎 | 高效工作者必备工具 【免费下载链接】kityminder 百度脑图 项目地址: https://gitcode.com/gh_mirrors/ki/kityminder 在信息爆炸的时代,如何将零散的想法系统化、复杂的项目结构化?作为一款开源免…...
NCM格式突破全攻略:从解密到跨平台播放的自由解锁方案
NCM格式突破全攻略:从解密到跨平台播放的自由解锁方案 【免费下载链接】ncmdump 项目地址: https://gitcode.com/gh_mirrors/ncmd/ncmdump 音乐作为数字生活的重要组成部分,却常常受到格式限制的困扰。网易云音乐的NCM加密格式就是其中典型代表&…...
Z-Image-Turbo_Sugar脸部Lora问题排查:常见错误403 Forbidden与连接问题解决
Z-Image-Turbo_Sugar脸部Lora问题排查:常见错误403 Forbidden与连接问题解决 部署和调用AI模型服务,就像组装一台新电脑,硬件都插好了,但开机时屏幕就是不亮,或者提示你密码错误。最近不少朋友在折腾Z-Image-Turbo_Su…...
OpenSSL实战:从零构建私有CA体系及多级证书签发指南
1. 为什么需要私有CA体系? 在日常开发中,我们经常遇到需要HTTPS加密通信的场景。比如微服务之间的API调用、内部系统的数据传输、物联网设备的安全连接等。虽然可以使用公共CA机构颁发的证书,但在以下场景中,自建CA体系会更加灵活…...
vue路由跳转打开新窗口并携带参数(vue2/vue3)
概要 在一些需求中经常遇到跳转页面,但是产品让跳转页面的同时打开一个新窗口方便用户进行对比数据,接下来就是跳转页面打开新窗口的方法 vue2的写法 const routeUrl this.$router.resolve({path: "/页面路由",query: { id: xx参数 },});wi…...
DAMOYOLO-S惊艳效果案例集:多领域高难度场景检测展示
DAMOYOLO-S惊艳效果案例集:多领域高难度场景检测展示 今天咱们不聊枯燥的理论和复杂的部署,直接来看点“硬货”。如果你正在寻找一个能在各种刁钻场景下都表现稳定的目标检测模型,那么DAMOYOLO-S绝对值得你花几分钟了解一下。它不是什么新概…...
千问3.5-2B集成IDEA开发环境:Java大模型应用快速构建指南
千问3.5-2B集成IDEA开发环境:Java大模型应用快速构建指南 1. 为什么要在IDEA中集成大模型? 作为Java开发者,我们经常需要在项目中处理各种文本处理任务。传统方式要么需要调用外部API(有网络延迟和费用问题)…...
基于FPGA的伺服驱动系统:电流环控制与多环路反馈、SVPWM及编码器协议实现的研究
伺服驱动FPGA电流环,包含坐标变换,电流环,速度环,位置环,电机反馈接口,SVPWM,编码器协议,电流环和编码器协议都是FPGA里实现的伺服驱动系统里玩FPGA可不是闹着玩的,尤其是…...
Python打包神器大PK:Nuitka vs PyInstaller,谁才是你的菜?(附实测数据)
Python打包工具深度评测:Nuitka与PyInstaller的终极对决 当开发者需要将Python项目分发给没有Python环境的用户时,打包工具的选择往往成为关键决策。本文将深入分析两大主流工具Nuitka和PyInstaller在多个维度的表现,帮助开发者根据项目需求做…...
