Я использовал stan и pymc3, но некоторые функции сложно реализовать эффективно, используя встроенные функции stan или theano. Поэтому мне бы очень хотелось низкоуровневую обертку, где я могу указать пользовательскую вероятность и градиент как функцию R.