t-dis的LLhood

Question

以下代码中t分布的对数似然函数有什么错误：

Rdpack::viewRd()

Answer 1

由于cdf，您只是错过了除以sigma的功能：

loglike <- function(p,x){
  mu <- p[1]
  sgm <- p[2]
  df <- p[3}
 ll <- dt((x-mu)/sgm, df=d)\sgm
 sum(log(ll)
}

与要计算正态分布的LL相同：

loglike <- function(p,x){
  mu <- p[1]
  sgm <- p[2]
 ll <- dnorm((x-mu)/sgm)/sgm
sum(log(ll)
}

您再次除以sigma