Bernoulli Family for mcp
Examples
# Fit a binary-response model with a probit link
data = data.frame(time = 1:6, y = c(0, 0, 0, 1, 1, 1))
fit = mcp(list(y ~ 1), data, family = bernoulli(link = "probit"), par_x = "time", sample = FALSE)
mcp_pars(fit) # Show the parameters of the fitted Bernoulli model
#> # A tibble: 1 × 9
#> name part scope role segment dpar order group_col population_name
#> <chr> <chr> <chr> <chr> <int> <chr> <int> <chr> <chr>
#> 1 Intercept_1 predict… popu… fixe… 1 mu NA NA NA
