时间戳

首页 > 解决方案 > 用 R 最大化对数似然函数

r - 用 R 最大化对数似然函数

问题描述

所以我想最大化似然函数，以便估计四个变量（a，b，c，sigma^2），R但它不起作用我不知道为什么。在这里，我估计年龄 30,60 和 90。对于每个年龄，我将其截断为 30，因此您可以看到 30 岁时 x 变为 0，60 岁时变为 30，依此类推。但是优化不起作用我不知道为什么；有人可以帮助我吗？

d30 <- 2975.1
d60 <- 11456.38
d90 <- 2977.08
r30 <- 1531956.05
r60 <- 650404.58
r90 <- 9728.47
Logg <- function(a,b,c,s)
{
  (d30*log((((a*exp(b*0))/(1+((s*a/b)*(exp(b*0)-1))))+c)*r30)-log(factorial(d30))-((((a*exp(b*0))/(1+((s*a/b)*(exp(b*0)-1))))+c)*r30))+
    (d60*log((((a*exp(b*30))/(1+((s*a/b)*(exp(b*30)-1))))+c)*r60)-log(factorial(d60))-((((a*exp(b*30))/(1+((s*a/b)*(exp(b*30)-1))))+c)*r60))+
    (d90*log((((a*exp(b*60))/(1+((s*a/b)*(exp(b*60)-1))))+c)*r90)-log(factorial(d90))-((((a*exp(b*60))/(1+((s*a/b)*(exp(b*60)-1))))+c)*r90))
}
Logg(1,1,1,1)
optim(c(0,0,0,0),Logg,control = list(fnscale=-1))

对数似然函数

对数似然函数

拉姆达定义

拉姆达定义

亩定义

亩定义

标签： roptimizationestimationmle

解决方案

也许这是因为你试图用数值计算非常高的阶乘？对于如此高的数字，R 只会返回无穷大：

> d30<-2975.1
> factorial(d30)
[1] Inf

推荐阅读