Использование диагностики pymc и последующего резюме на образцах не из pymc

Question

Я хотел бы использовать диагностические и итоговые функции pymc на образцах, которые я сгенерировал без использования pymc. В качестве одного примера я хотел бы использовать pcc's mc_error в моем собственном наборе образцов.

Некоторые из диагностических функций pymc могут принимать np.array образцов, но, если я правильно понимаю, другим, кажется, нужен объект трассировки. Как я могу преобразовать свой собственный набор образцов в объект трассировки, который может обрабатывать pymc? Я был бы признателен за любые указания.

Мои образцы в настоящее время хранятся в виде словаря, где каждый ключ является именем переменной I am sampling и указывает на np.array образцов.

Спасибо!

Answer 1

Вы говорите о функции batchsd? Это то, что используется для вычисления ошибки MC в PyMC 2. Его расположение в модуле pymc.database.base и может быть использовано на любом массиве.

Модуль pymc.diagnostics содержит все функции диагностики конвергенции и должен работать с массивами numpy.