diff --git a/Project.toml b/Project.toml
index 2377a5d..f3a6620 100644
--- a/Project.toml
+++ b/Project.toml
@@ -7,6 +7,7 @@ version = "0.2.5"
 Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 DiffResults = "163ba53b-c6d8-5494-b064-1a9d43ac40c5"
 ElasticArrays = "fdbdab4c-e67f-52f5-8c3f-e7b388dad3d4"
+ElasticPDMats = "2904ab23-551e-5aed-883f-487f97af5226"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
 GaussianProcesses = "891a1506-143c-57d2-908e-e1f8e92e6de9"
 NLopt = "76087f3c-5699-56af-9a33-bf431cd00edd"
@@ -18,6 +19,7 @@ TimerOutputs = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
 [compat]
 DiffResults = "0.0.4, 1"
 ElasticArrays = "0.3, 0.4, 1"
+ElasticPDMats = "0.2.3"
 ForwardDiff = "0.9, 0.10"
 GaussianProcesses = "0.9, 0.10, 0.11, 0.12"
 NLopt = "0.4, 0.5, 0.6, 1"
diff --git a/src/models/gp.jl b/src/models/gp.jl
index 50f6ef4..b8d40e8 100644
--- a/src/models/gp.jl
+++ b/src/models/gp.jl
@@ -37,7 +37,7 @@ directly on the GP parameters, e.g.
 function MAPGPOptimizer(; every = 10, kwargs...)
     MAPGPOptimizer(0, every,
                    merge(defaultoptions(MAPGPOptimizer),
-                         kwargs.data))
+                         kwargs)) #Changed kwargs.data -> kwargs
 end
 function optimizemodel!(o::MAPGPOptimizer, model::GPBase)
     if o.i % o.every == 0