Длительность т-дис

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:

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)
}

Вы снова делитесь на сигмы