Extract the shape parameters (low, up, k, and phi)

Usage

get_shape(
  fit,
  conf.int = TRUE,
  conf.level = 0.95,
  method = c("wald", "profile")
)

Arguments

fit: A profile_tls fit from fit_tls().
conf.int: Logical; include confidence-interval columns (default TRUE).
conf.level: Confidence level for the interval (default 0.95).
method: Either "wald" (default) or "profile".

Value

A tibble of the shape rows (low, up, k, and phi for the beta-binomial family) from tidy_parameters().

Details

With method = "profile", a shape parameter receives a profile interval only where tidy_parameters() supports that coordinate; up and unsupported shape coordinates retain their documented Wald route.

Examples

d <- simulate_tls(family = "binomial", CTmax = 36, z = 4, seed = 1)
fit <- fit_tls(d, y = survived, n = total, time = duration, temp = temp,
               family = "binomial", tref = 1)
get_shape(fit)
#> # A tibble: 3 × 8
#>   parameter group estimate std.error conf.low conf.high interval_type scale   
#>   <chr>     <chr>    <dbl>     <dbl>    <dbl>     <dbl> <chr>         <chr>   
#> 1 low       NA      0.0199   0.00552   0.0115    0.0345 wald          logit   
#> 2 up        NA      0.977    0.00797   0.962     0.993  wald          identity
#> 3 k         NA      4.89     0.413     4.14      5.78   wald          log