add compat and fix deprecations (#29)

* add compat and fix deprecations

* change to logbeta and loggamma

Author: raphaelsaavedra

Project.toml

@@ -16,3 +16,9 @@ HypothesisTests = "09f84164-cd44-5f33-b23f-e6b0d136a0d5"
 
 [targets]
 test = ["Test", "Random", "HypothesisTests"]
+
+[compat]
+julia = "~1"
+Distributions = "~0.21"
+Optim = "~0.19"
+SpecialFunctions = "~0.8"

src/distributions/beta.jl

@@ -28,7 +28,7 @@ Proof somewhere
 function log_likelihood(::Type{Beta}, y::Vector{T}, param::Vector{Vector{T}}, n::Int) where T
     loglik = 0.0
     for i in 1:n
-        loglik += (param[i][1] - 1)*log(y[i]) + (param[i][2] - 1)*log(1 - y[i]) - lbeta(param[i][1], param[i][2])
+        loglik += (param[i][1] - 1)*log(y[i]) + (param[i][2] - 1)*log(1 - y[i]) - logbeta(param[i][1], param[i][2])
     end
     return -loglik
 end

src/distributions/gamma.jl

@@ -25,7 +25,7 @@ Proof somewhere
 function log_likelihood(::Type{Gamma}, y::Vector{T}, param::Vector{Vector{T}}, n::Int) where T
     loglik = 0.0
     for i in 1:n
-        loglik += (param[i][1] - 1)*log(y[i]) - y[i]/param[i][2] - lgamma(param[i][1]) - param[i][1]*log(param[i][2])
+        loglik += (param[i][1] - 1)*log(y[i]) - y[i]/param[i][2] - loggamma(param[i][1]) - param[i][1]*log(param[i][2])
     end
     return -loglik
 end

src/distributions/poisson.jl

@@ -18,7 +18,7 @@ Proof somewhere
 function log_likelihood(::Type{Poisson}, y::Vector{Int}, param::Vector{Vector{T}}, n::Int) where T
     loglik = zero(T)
     for i in 1:n
-        loglik += y[i]*log(param[i][1]) - param[i][1] - lfactorial(y[i])
+        loglik += y[i]*log(param[i][1]) - param[i][1] - logfactorial(y[i])
     end
     return -loglik
 end