t-dis的LLhood

Question

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

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

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