原文链接：http://tecdat.cn/?p=18970

原文出处：拓端数据部落公众号

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

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

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

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

因此，是的最大似然估计。

例如，伯努利分布为，

给定样本，概率是

则对数似然

与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="")

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

> 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

最受欢迎的见解

1.Matlab马尔可夫链蒙特卡罗法（MCMC）估计随机波动率（SV，Stochastic Volatility）模型

2.基于R语言的疾病制图中自适应核密度估计的阈值选择方法

3.WinBUGS对多元随机波动率模型：贝叶斯估计与模型比较

4.R语言回归中的hosmer-lemeshow拟合优度检验

5.matlab实现MCMC的马尔可夫切换ARMA – GARCH模型估计

6.R语言区间数据回归分析

7.R语言WALD检验 VS 似然比检验

8.python用线性回归预测股票价格

9.R语言如何在生存分析与Cox回归中计算IDI，NRI指标

拓端tecdat|使用R语言做极大似然估计实例相关推荐

使用R语言做极大似然估计实例
最近我们被客户要求撰写关于极大似然估计的研究报告,包括一些图形和统计输出. 在普遍的理解中,最大似然估计是使用已知的样本结果信息来反向推断最有可能导致这些样本结果的模型参数值! 换句话说,最大似然估计 ...
R语言自定义极大似然估计函数、假设检验函数
R语言自定义极大似然函数函数.假设检验函数目录一.自定义极大似然函数二.自定义假设检验函数目录一.自定义极大似然函数 1.求出似然函数,以正太分布.指数分布为例正太分布似然函数: 指数分布 ...
拓端tecdat荣获掘金社区入驻新人奖
2021年7月,由掘金发起了"入驻成长礼"颁奖活动.本次活动邀请到知名开发者.服务机构代表等业界人士. 据了解,掘金社区"新入驻创作者礼"主要对已经积累了一定历 ...
拓端tecdat荣获2022年度51CTO博主之星
相信技术,传递价值,这是51CTO每一个技术创作者的动力与信念,2022 年度,拓端tecdat 作为新锐的数据分析咨询公司,在51CTO平台上,不断的输出优质的技术文章,分享前沿创新技术,输出最佳生 ...
r语言数据变量分段_R数据分析：用R语言做meta分析
这里以我的一篇meta分析为例,详细描述meta分析的一般步骤,该例子实现的是效应量β的合并 R包:metafor或meta包,第一个例子以metafor包为例. 1.准备数据集 2.异质性检验 in ...
使用r语言做garch模型_使用GARCH估计货币波动率
使用r语言做garch模型 Asset prices have a high degree of stochastic trends inherent in the time series. In o ...
通过R语言做灰色预测
通过R语言做灰色预测 GM(1,1)模型的定义数据的检验与处理数据的生成级比检验 GM(1,1)建模生成累加数据和均值数据构造矩阵BBB及数据向量YYY,有计算: 建立模型,求解,并还原数 ...
R语言|做点儿好玩儿的东西
电影<少年班>中两次出现过这样一道题: 在电影的最后,出现了这道题的解答,是一条心形线. 我觉得还挺有意思的,就试着用R语言做了一下.代码如下. heart<-data.frame( ...
笔记—R语言做矩阵散点图
笔记-R语言做矩阵散点图 #install.packages("corrgram") #install.packages("GGally") install.p ...
R语言做主成分分析在我国城镇居民家庭平均每人全年消费性支出分析中的应用
统计数据建模课程论文题目: R语言做主成分分析在我国城镇居民家庭平均每人全年消费性支出分析中的应用姓名: 学号: 班级: 授课教师: 完成时间: 2020年07月07日摘要居民消 ...

拓端tecdat|使用R语言做极大似然估计实例

原文链接：http://tecdat.cn/?p=18970

原文出处：拓端数据部落公众号

拓端tecdat|使用R语言做极大似然估计实例相关推荐

最新文章

热门文章