R openEBGM с функцией автогипера - PullRequest
0 голосов
/ 19 февраля 2019

У меня проблема при выполнении кода:

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.В целом, все работает нормально.

Пожалуйста, помогите :) Спасибо, Питер

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...