Как найти логарифмическую вероятность для обобщенного лямбда-распределения?

Question

У меня следующий r-код для оценки максимального правдоподобия обобщенного лямбда-распределения. У меня вопрос как вычислить

library(gld)
library(GLDEX)
example.data <- rgl(200,c(3,1,.4,-0.1),param="fkml")
example.fit <- fit.fkml(example.data,"ML",return.data=TRUE)
print(example.fit)
summary(example.fit)
plot(example.fit,one.page=FALSE)

Спасибо большое.