泊松分布的似然图和对数似然估计量

问题描述 投票:0回答:0

我正在尝试计算 R 中估计量向量的似然和对数似然。问题是,当我计算一个值的似然时,我一直得到 0。另一方面,对数似然很好。我猜这是因为阶乘在分母上太大了?

pois_lik <- function(x, lmin, lmax){
  numl <- 1000
  n <- length(x)
  lval <- numeric(numl)
  pl <- numeric(numl)
  lpl <- numeric(numl)
  lval <- seq(from=lmin, to=lmax, length.out=numl)
  for(k in 1:numll){
    pl[k] <- prod(dpois(x,lambda=lval[k]))
    lpl[k] <- sum(log(dpois(x,lambda=lval[k])))
  }
  pl_df <- data.frame(lval, pl, lpl)
  return(pl_df)
}


df <- pois_lik(sample, 1, 20)

这是我的代码。我的示例有 6115 个条目,范围为 {0,1,..,12}。

r function statistics factorial log-likelihood
© www.soinside.com 2019 - 2024. All rights reserved.