使用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


相关文章
|
6月前
|
存储 算法 Windows
课程视频|R语言bnlearn包:贝叶斯网络的构造及参数学习的原理和实例(下)
课程视频|R语言bnlearn包:贝叶斯网络的构造及参数学习的原理和实例
|
6月前
|
算法 数据可视化 数据挖掘
课程视频|R语言bnlearn包:贝叶斯网络的构造及参数学习的原理和实例(上)
课程视频|R语言bnlearn包:贝叶斯网络的构造及参数学习的原理和实例
|
6月前
|
数据可视化
R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码2
R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码
|
6月前
|
数据可视化 数据挖掘
R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码1
R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码
|
6月前
|
机器学习/深度学习 监控 数据可视化
R语言SOM神经网络聚类、多层感知机MLP、PCA主成分分析可视化银行客户信用数据实例2
R语言SOM神经网络聚类、多层感知机MLP、PCA主成分分析可视化银行客户信用数据实例
|
6月前
|
数据可视化 数据挖掘 Python
【视频】因子分析简介及R语言应用实例:对地区经济研究分析重庆市经济指标
【视频】因子分析简介及R语言应用实例:对地区经济研究分析重庆市经济指标
|
6月前
|
数据可视化
R语言nlme、nlmer、lme4用(非)线性混合模型non-linear mixed model分析藻类数据实例(上)
R语言nlme、nlmer、lme4用(非)线性混合模型non-linear mixed model分析藻类数据实例
|
6月前
|
机器学习/深度学习 算法 搜索推荐
【视频】决策树模型原理和R语言预测心脏病实例|数据分享(上)
【视频】决策树模型原理和R语言预测心脏病实例|数据分享
|
6月前
|
资源调度 安全
R语言nlme、nlmer、lme4用(非)线性混合模型non-linear mixed model分析藻类数据实例(下)
R语言nlme、nlmer、lme4用(非)线性混合模型non-linear mixed model分析藻类数据实例
|
6月前
|
数据可视化
R语言GAMLSS模型对艾滋病病例、降雪量数据拟合、预测、置信区间实例可视化
R语言GAMLSS模型对艾滋病病例、降雪量数据拟合、预测、置信区间实例可视化
下一篇
无影云桌面