У меня проблема при выполнении кода:
library(openEBGM)
rm(list=ls())
#data <- caers
data <- caers[1:100,]
processed <- processRaw(data,zeroes = FALSE,stratify = FALSE)
theta_init <- data.frame(alpha1 = c(0.2, 0.1, 0.3, 0.5, 0.2),
beta1 = c(0.1, 0.1, 0.5, 0.3, 0.2),
alpha2 = c(2, 10, 6, 12, 5),
beta2 = c(4, 10, 6, 12, 5),
p = c(1/3, 0.2, 0.5, 0.8, 0.4))
system.time(
hyper_estimates_squashed <- autoHyper(data = processed, theta_init = theta_init, squashed = FALSE, zeroes = FALSE)
)
Функция autoHyper не сходится в подмножестве данных caers.В целом, все работает нормально.
Пожалуйста, помогите :) Спасибо, Питер