diff --git a/README.md b/README.md
index d0e86a3..d97b9c9 100644
--- a/README.md
+++ b/README.md
@@ -6,6 +6,7 @@
 [![Codecov](https://codecov.io/gh/JuliaGaussianProcesses/GPLikelihoods.jl/branch/master/graph/badge.svg)](https://codecov.io/gh/JuliaGaussianProcesses/GPLikelihoods.jl)
 [![Code Style: Blue](https://img.shields.io/badge/code%20style-blue-4495d1.svg)](https://github.com/invenia/BlueStyle)
 [![ColPrac: Contributor's Guide on Collaborative Practices for Community Packages](https://img.shields.io/badge/ColPrac-Contributor's%20Guide-blueviolet)](https://github.com/SciML/ColPrac)
+[![Aqua QA](https://raw.githubusercontent.com/JuliaTesting/Aqua.jl/master/badge.svg)](https://github.com/JuliaTesting/Aqua.jl)
 
 GPLikelihoods.jl provides a collection of likelihoods to be used as building
 blocks for defining non-Gaussian problems.  It is intended to be mainly
diff --git a/test/Project.toml b/test/Project.toml
index 6e2c775..feeb76e 100644
--- a/test/Project.toml
+++ b/test/Project.toml
@@ -1,4 +1,5 @@
 [deps]
+Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 Functors = "d9f16b24-f501-4c13-a1f2-28368ffc5196"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
@@ -7,6 +8,7 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [compat]
+Aqua = "0.8"
 Distributions = "0.19, 0.20, 0.21, 0.22, 0.23, 0.24, 0.25"
 Functors = "0.1, 0.2, 0.3, 0.4"
 StatsFuns = "0.9, 1"
diff --git a/test/runtests.jl b/test/runtests.jl
index 1552a99..5b25cba 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -1,11 +1,13 @@
 using GPLikelihoods
 using GPLikelihoods: GaussHermiteExpectation, MonteCarloExpectation
 using GPLikelihoods.TestInterface: test_interface
-using Test
-using Random
-using Functors
+
+using Aqua
 using Distributions
+using Functors
+using Random
 using StatsFuns
+using Test
 using Zygote
 
 @testset "GPLikelihoods.jl" begin
@@ -20,4 +22,8 @@ using Zygote
         include("likelihoods/negativebinomial.jl")
     end
     include("expectations.jl")
+    @testset "Code quality (Aqua.jl)" begin
+        Aqua.test_all(GPLikelihoods; ambiguities = false)
+        Aqua.test_ambiguities([GPLikelihoods, Base, Core]; recursive=false)
+    end
 end