参数估计

通过取样本,并用样本构造函数,达成估计分布函数参数的目的

矩估计法

本质:用样本的各阶矩代替总体的各阶矩,即取:

E(X)=X‾=1n∑iXiE(X2)=1n∑iXi2E(X)=\overline{X}=\dfrac{1}{n}\sum_i X_i\\ E(X^2)=\dfrac{1}{n}\sum_iX_i^2 E(X)=X=n1iXiE(X2)=n1iXi2

极大似然估计法

本质:将使得样本 AAA 发生概率最大的参数值作为估计值

1、写出总体概率/密度函数

2、构造似然函数 L(λ)L(\lambda)L(λ)

3、两边取 ln⁡\lnln

4、对 λ\lambdaλ 求导

点估计的优良性准则

一、无偏性E(θ^)=θE(\hat{\theta})=\thetaE(θ^)=θ

定理:总体为 XXX,且 E(X)=μE(X)=\muE(X)=μD(X)=σ2D(X)=\sigma^2D(X)=σ2,样本为 (X1,⋯,Xn)(X_1,\cdots,X_n)(X1,,Xn),那么有:

1、X‾\overline{X}Xμ\muμ 的无偏估计

2、样本方差 S2S^2S2σ2\sigma^2σ2 的无偏估计

3、取 μ^=C1X1+⋯+CnXn\hat{\mu}=C_1X_1 + \cdots + C_nX_nμ^=C1X1++CnXn,若 C1+⋯+Cn=1C_1+\cdots + C_n=1C1++Cn=1,则 μ^\hat{\mu}μ^μ\muμ 的无偏估计

证明:

1 与 2:

已经在 上一份笔记 中证明过 E(X‾)=μE(\overline{X})=\muE(X)=μE(S2)=σ2E(S^2)=\sigma^2E(S2)=σ2

3:

E(μ^)=E(C1X1+⋯+CnXn)=E(C1X1)+⋯+E(CnXn)=C1E(X1)+⋯+CnE(Xn)=(C1+⋯+Cn)μ=μ\begin{align*} E(\hat{\mu}) &= E(C_1X_1 + \cdots + C_nX_n)\\ &=E(C_1X_1)+\cdots+E(C_nX_n)\\ &=C_1E(X_1)+\cdots+C_nE(X_n)\\ &=(C_1+\cdots+C_n)\mu\\ &=\mu \end{align*} E(μ^)=E(C1X1++CnXn)=E(C1X1)++E(CnXn)=C1E(X1)++CnE(Xn)=(C1++Cn)μ=μ

!注意:θ^\hat{\theta}θ^θ\thetaθ 的无偏估计,但是 g(θ^)g(\hat{\theta})g(θ^) 不一定是 g(θ)g(\theta)g(θ) 的无偏估计
例如:S2S^2S2σ2\sigma^2σ2 的无偏估计,而 SSS 不是 σ\sigmaσ 的无偏估计(性质)
该性质的证明:
D(S)=E(S2)−E(S)2=σ2−E(S)2⇒E(S)=σ2−D(S)⩽σ\begin{align*} &D(S) = E(S^2)-E(S)^2\\ &=\sigma^2-E(S)^2\\ &\Rightarrow E(S)=\sqrt{\sigma^2-D(S)}\leqslant\sigma \end{align*} D(S)=E(S2)E(S)2=σ2E(S)2E(S)=σ2D(S)σ

二、有效性D(θ1^)⩽D(θ2^)D(\hat{\theta_1})\leqslant D(\hat{\theta_2})D(θ1^)D(θ2^)

定理:总体为 XXX,且 E(X)=μE(X)=\muE(X)=μD(X)=σ2D(X)=\sigma^2D(X)=σ2,样本为 (X1,⋯,Xn)(X_1,\cdots,X_n)(X1,,Xn),若 a1+⋯+an=1a_1+\cdots+a_n=1a1++an=1,则现有两种 μ\muμ 的估计:a1X1+⋯+anXna_1X_1+\cdots+a_nX_na1X1++anXnX‾\overline{X}X,由有效性准则认为,X‾\overline{X}X 更优

证明:

由 上一份笔记 知 D(X‾)=σ2nD(\overline{X})=\dfrac{\sigma^2}{n}D(X)=nσ2

那么有:

D(θ^)=D(a1X1+⋯+anXn)=a12D(X1)+⋯+an2D(Xn)=σ2(a12+⋯+an2)⩾σ2n\begin{align*} D(\hat{\theta})&=D(a_1X_1+\cdots+a_nX_n)\\ &=a_1^2D(X_1)+\cdots+a_n^2D(X_n)\\ &=\sigma^2(a^2_1+\cdots+a^2_n)\\ &\geqslant \dfrac{\sigma^2}{n} \end{align*} D(θ^)=D(a1X1++anXn)=a12D(X1)++an2D(Xn)=σ2(a12++an2)nσ2

三、相合性(一致性)lim⁡n→+∞P(∣θ^−θ∣<ε)=1\lim\limits_{n\to +\infty}P(|\hat{\theta}-\theta|<\varepsilon)=1n+limP(θ^θ<ε)=1

置信区间

区间估计时,区间长度落在区间的概率 十分重要

P(θ1⩽θ⩽θ2)=1−αP(\theta_1 \leqslant \theta \leqslant \theta_2)=1-\alphaP(θ1θθ2)=1α,则 1−α1-\alpha1α 称为 置信度,而 [θ1,θ2][\theta_1,\theta_2][θ1,θ2] 则是估计区间

定义:

1、I=I(T,θ)I=I(T,\theta)I=I(T,θ),其中 θ\thetaθ 是未知参数,TTT 是已知的,随机变量 III 的分布 FFF 已知且其分布与 θ\thetaθ 无关,则将 III 称为 枢轴变量

2、给定 1−α1-\alpha1α,确定 FFF 的上 α2\dfrac{\alpha}{2}2α 分位数为 vα2v_{\frac{\alpha}{2}}v2α,上 1−α21-\dfrac{\alpha}{2}12αv1−α2v_{1-\frac{\alpha}{2}}v12α,则:

P(vα2⩽I(T,θ)⩽v1−α2)=1−αP\left(v_{\frac{\alpha}{2}}\leqslant I(T,\theta) \leqslant v_{1-\frac{\alpha}{2}}\right)=1-\alpha P(v2αI(T,θ)v12α)=1α

一个正态总体的均值和方差的区间估计

v=n(X‾−μ)σ∼N(0,1)v=\dfrac{\sqrt{n}(\overline{X}-\mu)}{\sigma}\sim N(0,1)v=σn(Xμ)N(0,1)

给定 1−α1-\alpha1α,令 P(v>uα2)=α2P(v>u_{\frac{\alpha}{2}})=\dfrac{\alpha}{2}P(v>u2α)=2α

1、σ2\sigma^2σ2 已知,对 μ\muμ 的区间估计:

P(−uα2⩽n(X‾−μ)σ⩽uα2)=1−αP(−σuα2n⩽X‾−μ⩽σuα2n)=1−αP(X‾−σuα2n⩽μ⩽X‾+σuα2n)=1−α\begin{align*} &P\left(-u_{\frac{\alpha}{2}}\leqslant \dfrac{\sqrt{n}(\overline{X}-\mu)}{\sigma} \leqslant u_{\frac{\alpha}{2}}\right)=1-\alpha\\ &P\left(-\dfrac{\sigma u_{\frac{\alpha}{2}}}{\sqrt{n}}\leqslant\overline{X}-\mu\leqslant\dfrac{\sigma u_{\frac{\alpha}{2}}}{\sqrt{n}}\right)=1-\alpha\\ &P\left(\overline{X}-\dfrac{\sigma u_{\frac{\alpha}{2}}}{\sqrt{n}}\leqslant\mu\leqslant\overline{X}+\dfrac{\sigma u_{\frac{\alpha}{2}}}{\sqrt{n}}\right)=1-\alpha \end{align*} P(u2ασn(Xμ)u2α)=1αP(nσu2αXμnσu2α)=1αP(Xnσu2αμX+nσu2α)=1α

也就是说,有 1−α1-\alpha1α 的把握,认为 μ\muμ 在区间 [X‾−σuα2n,X‾+σuα2n]\left[\overline{X}-\dfrac{\sigma u_{\frac{\alpha}{2}}}{\sqrt{n}}, \overline{X}+\dfrac{\sigma u_{\frac{\alpha}{2}}}{\sqrt{n}}\right][Xnσu2α,X+nσu2α]

2、σ2\sigma^2σ2 未知,对 μ\muμ 的区间估计:

构造枢轴变量:T=n(X‾−μ)S∼t(n−1)T=\dfrac{\sqrt{n}(\overline{X}-\mu)}{S}\sim t(n-1)T=Sn(Xμ)t(n1)

P(−tα2(n−1)⩽n(X‾−μS⩽tα2(n−1))=1−α⇒P(X‾−Sntα2(n−1)⩽μ⩽X‾+Sntα2(n−1))=1−α\begin{align*} &P\left(-t_{\frac{\alpha}{2}}(n-1)\leqslant\dfrac{\sqrt{n}(\overline{X}-\mu}{S}\leqslant t_{\frac{\alpha}{2}}(n-1)\right)=1-\alpha\\ &\Rightarrow P\left(\overline{X}-\dfrac{S}{\sqrt{n}}t_{\frac{\alpha}{2}}(n-1) \leqslant \mu \leqslant \overline{X}+\dfrac{S}{\sqrt{n}}t_{\frac{\alpha}{2}}(n-1)\right)=1-\alpha \end{align*} P(t2α(n1)Sn(Xμt2α(n1))=1αP(XnSt2α(n1)μX+nSt2α(n1))=1α

也就是说,有 1−α1-\alpha1α 的把握,认为 μ\muμ 在区间 [X‾−Sntα2(n−1),X‾+Sntα2(n−1)]\left[\overline{X}-\dfrac{S}{\sqrt{n}}t_{\frac{\alpha}{2}}(n-1), \overline{X}+\dfrac{S}{\sqrt{n}}t_{\frac{\alpha}{2}}(n-1)\right][XnSt2α(n1),X+nSt2α(n1)]

3、μ\muμ 已知,对 σ2\sigma^2σ2 的区间估计:

构造枢轴变量:χ2=1σ2∑i=1n(Xi−μ)2∼χ2(n)\chi^2=\dfrac{1}{\sigma^2}\sum\limits_{i=1}^n(X_i-\mu)^2\sim \chi^2(n)χ2=σ21i=1n(Xiμ)2χ2(n),给定 1−α1-\alpha1αχ1−α22(n)\chi^2_{1-\frac{\alpha}{2}}(n)χ12α2(n)χα22(n)\chi^2_{\frac{\alpha}{2}}(n)χ2α2(n)

注意:之所以不继续使用正态分布的枢轴变量,是因为 σ\sigmaσ 在开平方后不是无偏估计

P(χ1−α22(n)⩽1σ2∑i=1n(Xi−μ)2⩽χα22(n))=1−α⇒∑i=1n(Xi−μ)2χα22(n)⩽σ2⩽∑i=1n(Xi−μ)2χ1−α22(n)\begin{align*} &P\left(\chi^2_{1-\frac{\alpha}{2}}(n)\leqslant\dfrac{1}{\sigma^2}\sum_{i=1}^n(X_i-\mu)^2\leqslant\chi^2_{\frac{\alpha}{2}}(n)\right)=1-\alpha\\ &\Rightarrow \dfrac{\sum_{i=1}^n(X_i-\mu)^2}{\chi^2_{\frac{\alpha}{2}}(n)}\leqslant\sigma^2\leqslant\dfrac{\sum_{i=1}^n(X_i-\mu)^2}{\chi^2_{1-\frac{\alpha}{2}}(n)} \end{align*} P(χ12α2(n)σ21i=1n(Xiμ)2χ2α2(n))=1αχ2α2(n)i=1n(Xiμ)2σ2χ12α2(n)i=1n(Xiμ)2

也就是说,有 1−α1-\alpha1α 的把握,认为 σ2\sigma^2σ2 在区间 [∑i=1n(Xi−μ)2χα22(n),∑i=1n(Xi−μ)2χ1−α22(n)]\left[\dfrac{\sum_{i=1}^n(X_i-\mu)^2}{\chi^2_{\frac{\alpha}{2}}(n)},\dfrac{\sum_{i=1}^n(X_i-\mu)^2}{\chi^2_{1-\frac{\alpha}{2}}(n)}\right][χ2α2(n)i=1n(Xiμ)2,χ12α2(n)i=1n(Xiμ)2]

4、μ\muμ 未知,对 σ2\sigma^2σ2 的区间估计:

构造枢轴变量:χ2=(n−1)S2σ2∼χ2(n−1)\chi^2=\dfrac{(n-1)S^2}{\sigma^2}\sim \chi^2(n-1)χ2=σ2(n1)S2χ2(n1),给定 1−α1-\alpha1αχ1−α22(n−1)\chi^2_{1-\frac{\alpha}{2}}(n-1)χ12α2(n1)χα22(n−1)\chi^2_{\frac{\alpha}{2}}(n-1)χ2α2(n1)

P(χ1−α22(n−1)⩽(n−1)S2σ2⩽χα22(n−1))=1−α⇒P((n−1)S2χα22(n−1)⩽σ2⩽(n−1)S2χ1−α22(n−1))=1−α\begin{align*} &P\left(\chi^2_{1-\frac{\alpha}{2}}(n-1)\leqslant\dfrac{(n-1)S^2}{\sigma^2}\leqslant\chi^2_{\frac{\alpha}{2}}(n-1)\right)=1-\alpha\\ &\Rightarrow P\left(\dfrac{(n-1)S^2}{\chi^2_{\frac{\alpha}{2}}(n-1)}\leqslant\sigma^2\leqslant\dfrac{(n-1)S^2}{\chi^2_{1-\frac{\alpha}{2}}(n-1)}\right)=1-\alpha \end{align*} P(χ12α2(n1)σ2(n1)S2χ2α2(n1))=1αP(χ2α2(n1)(n1)S2σ2χ12α2(n1)(n1)S2)=1α

也就是说,有 1−α1-\alpha1α 的把握,认为 σ2\sigma^2σ2 在区间 [(n−1)S2χα22(n−1),(n−1)S2χ1−α22(n−1)]\left[\dfrac{(n-1)S^2}{\chi^2_{\frac{\alpha}{2}}(n-1)},\dfrac{(n-1)S^2}{\chi^2_{1-\frac{\alpha}{2}}(n-1)}\right][χ2α2(n1)(n1)S2,χ12α2(n1)(n1)S2]

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如若转载,请注明出处:http://www.pswp.cn/pingmian/89772.shtml
繁体地址,请注明出处:http://hk.pswp.cn/pingmian/89772.shtml
英文地址,请注明出处:http://en.pswp.cn/pingmian/89772.shtml

如若内容造成侵权/违法违规/事实不符,请联系英文站点网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

服务器后台崩溃的原因

当我们双十一活动零点拼命刷新却卡在支付完页面&#xff0c;游戏页面等不进去&#xff0c;公司系统瘫痪全体员工干瞪眼&#xff0c;服务器崩溃绝对是数字时代中的酷刑&#xff01;那服务器为什么会说崩就崩&#xff0c;用户对于这种情况该如何进行避雷呢&#xff1f;服务器主要…

线程池与ThreadPoolExecutor源码解析(上)

一、线程池线程池&#xff08;ThreadPool&#xff09;是一种线程复用的机制。它维护着若干个线程&#xff0c;任务来了就复用这些线程去执行&#xff0c;任务做完线程不会销毁&#xff0c;而是回到池中等待下一个任务。为什么要用线程池&#xff1f;降低资源消耗&#xff1a;避…

Linux内核IP分片重组机制剖析:高效与安全的艺术

在IP网络通信中,当数据包超过MTU限制时,路由器会将其拆分为多个分片。这些分片到达目标主机后,内核必须高效、安全地重组原始数据包。Linux内核的net/ipv4/inet_fragment.c实现了一套精妙的分片管理框架,完美平衡了性能和安全性需求。本文将深入剖析其设计哲学与关键技术。…

相机模型和对极几何

一、相机模型 1.针孔相机模型-外参矩阵 1.世界坐标系到相机坐标系 世界坐标系&#xff1a;可以定义空间中任意一个位置&#xff0c;原点位置三个坐标轴方向坐标系姿态&#xff08;X,Y,Z&#xff09;相机坐标系&#xff1a;定义在相机上&#xff0c;原点是相机中心&#xff0c;z…

Git 常用命令与操作步骤

以下是 Git 常用命令与操作步骤 的整理&#xff0c;涵盖日常开发中最核心的场景&#xff0c;适合快速查阅和上手&#xff1a;1. 初始化与克隆仓库操作命令本地初始化仓库git init克隆远程仓库git clone <仓库URL> &#xff08;如 git clone https://gitlab.com/user/repo…

Leetcode-.283移动零

class Solution:def moveZeroes(self, nums: List[int]) -> None:"""Do not return anything, modify nums in-place instead."""pos0for i in range(len(nums)):if nums[i]!0:nums[pos],nums[i]nums[i],nums[pos]pos1本题运用双指针来写&…

在React中做过哪些性能优化?

1. 使用 React.memo 进行组件优化 问题:当父组件重新渲染时,子组件也会重新渲染,即使它的 props 没有变化。 解决方案:使用 React.memo 包裹子组件,让其只在 props 变化时才重新渲染。 const MyComponent = React.memo((props) => {// 子组件代码 }); 2. 使用 useCa…

安装docker可视化工具 Portainer中文版(ubuntu上演示,所有docker通用) 支持控制各种容器,容器操作简单化 降低容器门槛

以下有免费的4090云主机提供ubuntu22.04系统的其他入门实践操作 地址&#xff1a;星宇科技 | GPU服务器 高性能云主机 云服务器-登录 相关兑换码星宇社区---4090算力卡免费体验、共享开发社区-CSDN博客 兑换码要是过期了&#xff0c;可以私信我获取最新兑换码&#xff01;&a…

ansible批量部署zabbix客户端

✅ansible编写剧本步骤 1️⃣创建roles目录结构2️⃣在group_vars/all/main.yml中定义变量列表3️⃣在tasks目录下编写tasks任务4️⃣在files目录下准备部署文件5️⃣在templates目录下创建j2模板文件6️⃣在handlers目录下编写handlers7️⃣在roles目录下编写主playbook8️⃣运…

蚂蚁数科AI数据产业基地正式投产,携手苏州推进AI产业落地

近日&#xff0c;蚂蚁数科AI数据产业基地在太仓智汇谷科技创新园正式投产。该基地作为苏州市首个AI数据产业基地&#xff0c;旨在通过跨行业人才与前沿技术&#xff0c;为长三角制造业、金融、医疗等领域的大模型落地提供场景化、高质量的训练数据支撑。数据被视为AI学习的核心…

计算机的网络体系及协议模型介绍

目录 1、网络协议介绍 1.1、定义 1.2、基本作用 1.3、协议的主要内容 2、网络协议分层 2.1、协议分层原因 2.2、网络协议分层的缺点 2.3、OSI协议和TCP/IP协议的联系 3、TCP/IP 协议族 3.1、定义介绍 3.2、组成 1、应用层 2、运输层 3、网络层 3.3、底层流程 4、…

密码管理安全防御

密码管理是信息安全的核心环节,其目标是通过规范密码的生成、存储、传输、验证和生命周期管理,防止未授权访问,保护用户账号和系统资源的安全。以下从核心原则、技术实践、常见问题及解决方案等方面详细说明: 一、密码管理的核心原则 密码管理需遵循“安全性”与“可用性…

Java异步日志系统性能优化实践指南:基于Log4j2异步Appender与Disruptor

Java异步日志系统性能优化实践指南&#xff1a;基于Log4j2异步Appender与Disruptor 一、技术背景与应用场景 在高并发的后端应用中&#xff0c;日志记录往往成为性能瓶颈之一。同步写日志会阻塞业务线程&#xff0c;导致响应延迟&#xff1b;而简单的异步队列实现又可能出现积压…

Mybatis07-缓存

一、缓存机制的原理计算机每次从mysql中执行sql语句&#xff0c;都是内存与硬盘的通信&#xff0c;对计算机来说&#xff0c;影响效率。因此使用缓存机制。1-1、MyBatis 的缓存机制&#xff1a;执行 DQL&#xff08;select 语句&#xff09;的时候&#xff0c;将查询结果放到缓…

【机器学习深度学习】LoRA 与 QLoRA:大模型高效微调的进阶指南

目录 前言 一、LoRA&#xff1a;低秩微调的经典之作 二、QLoRA&#xff1a;效率与精度的升级版 三、LoRA vs QLoRA&#xff1a;如何选择&#xff1f; 3.1 性能维度对比 3.2 根据「显卡资源」选择 3.3 根据「任务类型与目标」选择 3.4 根据「模型规模」选择 3.5 根据…

教育行业网络升级最佳实践:SD-WAN、传统方案与混合方案对比分析

随着教育行业的数字化转型不断深入&#xff0c;网络的稳定性、灵活性和安全性成为各类教育应用&#xff08;如远程课堂、智慧校园和教育云平台&#xff09;的核心支撑。然而&#xff0c;传统的 MPLS 专线方案成本高、扩展性差&#xff0c;而纯 SD-WAN 的方案在极高可靠性要求的…

[黑马头条]-文章列表加载

目录 1.1)需求分析 1.2)表结构分析 ap_article 文章基本信息表 ap_article_config 文章配置表 ap_article_content 文章内容表 导入文章数据库 实现思路 接口定义 功能实现 定义接口 编写mapper文件 编写业务层代码 实现类&#xff1a; 定义常量类 编写控制器代码 …

使用TIANAI-CAPTCHA进行行为验证码的生成和缓存的二次校验

1.导入依赖&#xff1a;<dependency><groupId>cloud.tianai.captcha</groupId><artifactId>tianai-captcha-springboot-starter</artifactId><version>1.5.2</version> </dependency>2.在application.yml中配置验证码相关配置…

db.refresh()的重复使用和db.rollback()

db.refresh()在 SQLAlchemy 中&#xff0c;db.refresh() 用于从数据库中重新加载对象的状态&#xff0c;确保对象属性与数据库中的实际数据保持一致。下面详细介绍其使用场景和作用&#xff1a;1.获取数据库生成的值当数据库自动生成字段&#xff08;如自增 ID、默认值、触发器…

《Web安全之机器学习入门》读书笔记总结

目录 一、案例总结 1、基础知识 &#xff08;1&#xff09;第1章 通向智能安全的旅程 &#xff08;2&#xff09;第2章 打造机器学习工具箱 &#xff08;3&#xff09;第3章 机器学习概述 &#xff08;4&#xff09;第4章 Web安全基础 2、安全案例 &#xff08;1&#…