Я новичок в r и использую пакет elrm для точной логистической регрессии (ELR). Я хотел бы, чтобы можно было воспроизвести мои результаты, поэтому, судя по моим исследованиям, мне нужно заложить основу для симуляции.
Мой набор данных связан с повторными преступлениями, связанными с предыдущими характеристиками, уровнем риска по стандартной оценке и тем, как наша программа влияет на уровень повторных преступлений. Итак, ЛЕЧЕНИЕ — это моя основная независимая интересующая переменная, а уровень рецидивизма (RECID_sum/case_sum) — моя зависимая переменная.
Я пробовал это несколькими разными способами.
Сначала я побежал
set.seed(66672)
После этого я запускаю приведенный ниже код для ELR:
m.treatment <- elrm(formula = RECID_sum / case_sum ~ TREATMENT + FDRUGVSP +
FPROPVSP + FELEVENTVSP + RISK + street_time + TREATMENT, interest =~ TREATMENT,
iter=100000, dataset=final, burnIn=100)
И каждый раз я получал разные результаты, даже после того, как обязательно снова установил начальное значение перед запуском кода во второй раз.
Я также пробовал следующее:
m.treatment <- elrm(formula = RECID_sum / case_sum ~ TREATMENT + FDRUGVSP +
FPROPVSP + FELEVENTVSP + RISK + street_time + TREATMENT, interest =~ TREATMENT,
iter = 100000, dataset = final, burnIn = 100); set.seed(66672)
но каждый раз получал разные результаты. Я не получил никаких ошибок при запуске этого кода.
Я провел некоторое исследование по установке начальных значений для симуляций Монте-Карло в целом, но этот подход не кажется применимым здесь. Кажется, нет информации об установке начального числа для ELR в пакете elrm. Я могу пропустить это из-за моей относительной неопытности с r, и я могу не знать точно, что я ищу. Обычно я использую SPSS, который настолько неуклюж и, к сожалению, не может сделать ELR разумным образом.
Любая помощь, которую вы можете предложить, очень ценится, и извините меня за новизну, если это простой ответ.