diff --git a/.github/workflows/Documenter.yml b/.github/workflows/Documenter.yml
index edf6c42..a234e87 100644
--- a/.github/workflows/Documenter.yml
+++ b/.github/workflows/Documenter.yml
@@ -24,7 +24,7 @@ jobs:
       - uses: actions/checkout@v2
       - uses: julia-actions/setup-julia@latest
         with:
-          version: '1.4'
+          version: '1.6'
       - name: Install dependencies
         run: julia --project=docs/ -e 'using Pkg; Pkg.develop(PackageSpec(path=pwd())); Pkg.instantiate()'
       - name: Build and deploy
@@ -32,4 +32,3 @@ jobs:
           GITHUB_TOKEN: ${{ secrets.GITHUB_TOKEN }} # For authentication with GitHub Actions token
           DOCUMENTER_KEY: ${{ secrets.DOCUMENTER_KEY }} # For authentication with SSH deploy key
         run: julia --project=docs/ docs/make.jl
-
diff --git a/.github/workflows/Tier1.yml b/.github/workflows/Tier1.yml
index c1e3a2a..690387c 100644
--- a/.github/workflows/Tier1.yml
+++ b/.github/workflows/Tier1.yml
@@ -17,7 +17,7 @@ jobs:
     runs-on: ${{ matrix.os }}
     strategy:
       matrix:
-        julia-version: [1.4,1.5]
+        julia-version: [1.4, 1.5, 1.6]
         julia-arch: [x64]
         os: [ubuntu-18.04, macos-10.15, windows-2019]
     steps:
diff --git a/Project.toml b/Project.toml
index 172c033..dbf15ec 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "ClinicalTrialUtilities"
 uuid = "535c2557-d7d0-564d-8ff9-4ae146c18cfe"
 authors = ["Vladimir Arnautov (mail@pharmcat.net)"]
-version = "0.5.0"
+version = "0.5.1"
 
 [deps]
 DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
@@ -9,21 +9,19 @@ Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 QuadGK = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 RecipesBase = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
-Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 Roots = "f2b01f46-fcfa-551c-844a-d8ac1e96c665"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 
 [compat]
-DataFrames = "0.19, 0.20, 0.21, 0.22, 1.0"
-Distributions = "0.16, 0.17, 0.18, 0.19, 0.20, 0.21, 0.22, 0.23, 0.24"
-QuadGK = "2.0, 2.1, 2.2, 2.3, 2.4"
-RecipesBase = "0.7, 0.8, 1"
-Reexport = "0.1, 0.2, 1.0"
-Roots = "0.7, 0.8, 1"
-SpecialFunctions = "0.8, 0.9, 0.10, 1"
-StatsBase = "0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.30, 0.31, 0.32, 0.33"
-julia = "1.0, 1.1, 1.2, 1.3, 1.4, 1.5"
+DataFrames = "0.22, 1"
+Distributions = "0.20, 0.21, 0.22, 0.23, 0.24"
+QuadGK = "2"
+RecipesBase = "1"
+Roots = "1"
+SpecialFunctions = "1"
+StatsBase = "0.30, 0.31, 0.32, 0.33"
+julia = "1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6"
 
 [extras]
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
diff --git a/cange.log b/cange.log
index f448e2b..e61ba2e 100644
--- a/cange.log
+++ b/cange.log
@@ -1,3 +1,6 @@
+v0.5.1
+
+
 v0.5.0
   - StableRNGs test
   - rng keyword for randomtable, randomseq
diff --git a/docs/Project.toml b/docs/Project.toml
index dca4fa8..6ab3c40 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -6,7 +6,7 @@ Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 
 [compat]
 
-Documenter = "0.23, 0.24"
-DataFrames = "0.19, 0.20"
-CSV = "0.5, 0.6"
-Plots = "1.2, 1.3"
+Documenter = "0.23, 0.24, 0.25, 0.26"
+DataFrames = "0.22, 1"
+CSV = "0.6, 0.7, 0.8"
+Plots = "1"
diff --git a/src/ClinicalTrialUtilities.jl b/src/ClinicalTrialUtilities.jl
index 1cd7231..7dbcbd9 100644
--- a/src/ClinicalTrialUtilities.jl
+++ b/src/ClinicalTrialUtilities.jl
@@ -14,21 +14,17 @@
 __precompile__(true)
 module ClinicalTrialUtilities
 
-using Distributions, Random, Roots, QuadGK, RecipesBase, Reexport
+using Distributions, Random, Roots, QuadGK, RecipesBase
 
-@reexport using StatsBase
+using StatsBase
 
 import SpecialFunctions
 import Base: show, findfirst, getproperty, showerror, getindex, length, in, iterate, eltype, deleteat!, findall
 import StatsBase.confint
 import DataFrames: DataFrame, DataFrames, unstack
 
-try
-    methods(SpecialFunctions.logabsgamma)
-    global lgamma(x) = SpecialFunctions.logabsgamma(x)[1]
-catch
-    global lgamma(x) = SpecialFunctions.lgamma(x)
-end
+
+lgamma(x) = SpecialFunctions.logabsgamma(x)[1]
 
 const ZDIST  = Normal()
 const LOG2   = log(2)