我正在尝试计算 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}。