使用R语言做极大似然估计实例

简介: 使用R语言做极大似然估计实例

在普遍的理解中,最大似然估计是使用已知的样本结果信息来反向推断最有可能导致这些样本结果的模型参数值!

换句话说,最大似然估计提供了一种在给定观测数据的情况下评估模型参数的方法,即“模型已确定且参数未知”。

在所有双射函数的意义上,极大似然估计是不变的   ,如果   的极大似然估计     。

让   ,   等于   中的似然函数。由于   是的最大似然估计  

因此,   的最大似然估计  。

例如,伯努利分布为   ,  

给定样本   ,概率是

则对数似然

与ICI

因此,一阶条件

何时满足   。为了说明,考虑以下数据

> X

\[1\] 0 0 1 1 0 1 1 1 1 0 0 0 1 0 1

(负)对数似然

> loglik=function(p){

+ -sum(log(dbinom(X,size=1,prob=p)))

+ }

我们可以在下面看到

> plot(u,v,type="l",xlab="",ylab="")

根据以上计算,我们知道的极大似然估计  是

> mean(X)

\[1\] 0.53

数值为

$par

\[1\] 0.53


$value

\[1\] 10.36


$counts

function gradient

20 NA


$convergence

\[1\] 0


$message

NULL

我们没有说优化是在区间内 。但是,我们的概率估计值属于 。为了确保最优值在 ,我们可以考虑一些约束优化程序

ui=matrix(c(1,-1),2,1), ci=c(0,-1)

$par

\[1\] 0.53


$value

\[1\] 10.36


$counts

function gradient

20 NA


$convergence

\[1\] 0


$message

NULL


$outer.iterations

\[1\] 2


$barrier.value

\[1\] 6.91e-05

在上一张图中,我们达到了对数似然的最大值

> abline(v=opt$par,col="red")

另一种方法是考虑   (如指数分布)。则对数似然

这里

因此,一阶条件

满足

从数值角度来看,我们有相同的最优值

(opt=optim(0,loglik))

$par

\[1\] 0.13


$value

\[1\] 10.36


$counts

function gradient

20 NA


$convergence

\[1\] 0


$message

NULL


> exp(opt$par)/(1+exp(opt$par))

\[1\] 0.53


相关文章
|
3天前
|
数据挖掘
R语言中的动态线性模型
【4月更文挑战第27天】本文探讨了R语言中动态线性模型(DLMS)在处理自相关时间序列数据的应用。DLMs基于状态空间模型,包含观测和状态方程,能适应新信息并进行预测。使用`dlm`包可构建和估计模型,通过实例展示了如何预测股票价格。模型解释与验证涉及拟合优度、预测准确性和模型诊断。R还支持多变量、非线性及贝叶斯DLMs等高级主题,扩展了时间序列分析的能力。`dlm`包与其他工具一起,使R成为动态线性模型分析的强大平台。
|
8天前
【视频】线性回归中的贝叶斯推断与R语言预测工人工资数据|数据分享-3
【视频】线性回归中的贝叶斯推断与R语言预测工人工资数据|数据分享
17 0
【视频】线性回归中的贝叶斯推断与R语言预测工人工资数据|数据分享-3
|
4天前
|
算法 vr&ar Python
【视频】什么是梯度下降?用线性回归解释和R语言估计GARCH实例
【视频】什么是梯度下降?用线性回归解释和R语言估计GARCH实例
|
4天前
|
算法 数据挖掘
R语言中的贝叶斯统计方法
【4月更文挑战第26天】R语言在贝叶斯统计中发挥着重要作用,提供如"BUGS"、"Stan"、"JAGS"等包来处理复杂模型和数值计算。贝叶斯方法基于概率论,涉及先验分布、似然函数、后验分布和MCMC模拟。"BUGS"适用于复杂层次模型,"Stan"则在大规模数据和复杂模型上有优势。
|
5天前
|
机器学习/深度学习 前端开发 自动驾驶
【视频】什么是Bootstrap自抽样及应用R语言线性回归预测置信区间实例|数据分享
【视频】什么是Bootstrap自抽样及应用R语言线性回归预测置信区间实例|数据分享
11 2
|
8天前
|
数据挖掘
【视频】线性回归中的贝叶斯推断与R语言预测工人工资数据|数据分享-2
【视频】线性回归中的贝叶斯推断与R语言预测工人工资数据|数据分享
16 0
【视频】线性回归中的贝叶斯推断与R语言预测工人工资数据|数据分享-2
|
14天前
在R语言中使用概率分布:dnorm,pnorm,qnorm和rnorm
在R语言中使用概率分布:dnorm,pnorm,qnorm和rnorm
16 0
|
14天前
|
算法 Windows
R语言通过WinBUGS对MGARCH和MSV模型进行贝叶斯估计和比较
R语言通过WinBUGS对MGARCH和MSV模型进行贝叶斯估计和比较
16 0
|
14天前
R语言对回归模型进行协方差分析
R语言对回归模型进行协方差分析
23 0
|
13天前
|
数据可视化
R语言最优化问题中的共轭函数
R语言最优化问题中的共轭函数

热门文章

最新文章