# NOT RUN {
x <- model.matrix(~ factor(wool) + factor(tension), warpbreaks)
y <- warpbreaks$breaks
dprob <- make.dpois(x = x, y = y)
params <- list(rep(1, 4))
update_probabilities(dprob, params)
# }
Run the code above in your browser using DataLab