Error using Turing + ReverseDiff with PyJulia

I am getting an error when I use Turing and Reverse Diff with Pyjulia, shown below. Does anyone have any experience getting this to work? Also created a MWE here:

Error


Stacktrace:

 [1] setadbackend(::Symbol) at /Users/ada/.julia/packages/Turing/fbY6B/src/core/ad.jl:7

 [2] top-level scope at /Users/ada/vyasa/Code/tapestry/dd915c6f27e8b6b4b2d1549d0a2a8371/500d611a28e02c18a31923e427138e5b/rdiff_test.jl:4

 [3] include(::String) at ./client.jl:439

 [4] #invokelatest#1 at ./essentials.jl:712 [inlined]

 [5] invokelatest(::Any, ::Any) at ./essentials.jl:711

 [6] _pyjlwrap_call(::Function, ::Ptr{PyCall.PyObject_struct}, ::Ptr{PyCall.PyObject_struct}) at /Users/ada/.julia/packages/PyCall/zqDXB/src/callback.jl:28

 [7] pyjlwrap_call(::Ptr{PyCall.PyObject_struct}, ::Ptr{PyCall.PyObject_struct}, ::Ptr{PyCall.PyObject_struct}) at /Users/ada/.julia/packages/PyCall/zqDXB/src/callback.jl:49

in expression starting at /Users/ada/vyasa/Code/tapestry/dd915c6f27e8b6b4b2d1549d0a2a8371/500d611a28e02c18a31923e427138e5b/rdiff_test.jl:4?