R语言-点估计与极大似然估计

这里写自定义目录标题


提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档


前言

提示:这里可以添加本文要记录的大概内容:
例如:随着人工智能的不断发展,机器学习这门技术也越来越重要,很多人都开启了学习机器学习,本文就介绍了机器学习的基础内容。


提示:以下是本篇文章正文内容,下面案例可供参考

一、点估计

对于随机变量来说,矩是其最广泛,最常用的数字特征,母体的各阶矩一般与的分布中所含的未知参数有关,有的甚至就等于未知参数。由辛钦大数定律知,简单随机子样的子样原点矩依概率收敛到相应的母体原点矩。这就启发我们想到用子样矩替换母体矩,进而找出未知参数的估计,基于这种思想求估计量的方法称为矩法。用矩法求得的估计称为矩法估计,简称矩估计。它是由英国统计学家皮尔逊Pearson于1894年提出的。
例如我们来看看正态分布N(0,1)的矩估计效果。

x<-rnorm(100) #产生N(0,1)的100个随机数

mu<-mean(x) #对N(mu,sigma)中的mu做矩估计

sigma<-var(x) #这里的var并不是样本方差的计算函数,而是修正的样本方差,其实也就是x的总体方差

mu

[1] -0.1595923

sigma

[1] 1.092255

二、极大似然估计

极大似然估计的基本思想是:基于样本的信息参数的合理估计量是产生获得样本的最大概率的参数值。
设X1,X2,…,X100为来自P(lambda)的独立同分布样本,那么似然函数为:

L(lambda,x)=lambda^(x1+x2+…+x100)exp(10lambda)/(gamma(x1+1)…gamma(x100+1))

x<-rpois(100,2)

sum(x)

[1] 215

ga(x)#这是一个求解gamma(x1+1)…gamma(x100+1)的函数,用gamma函数求阶乘是为了提高计算效率

[1] 1.580298e+51


版权声明:本文为weixin_44986915原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明。