Как внедрить пакет адаптивной выборочной метрополии (ARMS) в R?
У меня есть модель dlm (динамическая линейная модель), которую я создал с использованием пакета dlm в R. Я хотел бы реализовать адаптивную выборку отклонения метрополии, чтобы ограничить модель эволюции состояния dlm.
Я прочитал виньетку dlm, в которой написано, что может быть способ сделать это с помощью метода выборочной адаптации в метрополии. В справочном пакете говорится, что для этой функции необходимо указать целевую плотность журнала. Какова плотность цели журнала? Определение целевой плотности журнала неясно в пакете справки. В пакете говорится, что это лог-плотность, однако есть несколько примеров, в которых он определяется иначе. Каково непротиворечивое определение логарифмической целевой плотности?
Ниже приведен пример из справочной функции пакета оружий, которая реализует функцию вооружений для экспоненциальной функции с параметром 1. Однако целевая плотность записи определяется как -x (второй параметр функции).
y <- arms(5, function(x) -x, function(x) (x>0)*(x<70), 5000)
summary(y); hist(y,prob=TRUE,main="Exponential(1)")
curve(exp(-x),0,8,add=TRUE)