Skip to contents

Reads inst/extdata/reference.csv: one row per published number, with its source, the R call that recomputes it with this package (call), the expression that extracts the number from the result (extract, in terms of x), the printed value and the digits to compare. The same file drives the test of published values and the vignette that reproduces the published examples. Each distinct call is evaluated once.

Usage

stab_reference(source = NULL, evaluate = FALSE, level = 1)

Arguments

source

optional regular expression; keep rows whose source matches.

evaluate

logical; evaluate the rows and add columns computed and agrees (default FALSE).

level

largest level evaluated (1 = fast rows, 2 = slow rows; default 1); rows above it get NA.

Value

A data.table with columns id, source, quantity, call, extract, value, digits, level, note (and computed, agrees).

See also

Other stability utilities: stab_notation(), stab_save(), stab_theme()

Author

Muhammad Yaseen (myaseen208@gmail.com)

Examples

ref <-
  stab_reference(
      source   = NULL
    , evaluate = FALSE
    , level    = 1
    )
ref
#>                     id                                    source
#>                 <char>                                    <char>
#>  1:  eskridge1990_EV_1 Eskridge (1990), Crop Science 30, Table 4
#>  2:  eskridge1990_EV_2 Eskridge (1990), Crop Science 30, Table 4
#>  3:  eskridge1990_EV_3 Eskridge (1990), Crop Science 30, Table 4
#>  4:  eskridge1990_EV_4 Eskridge (1990), Crop Science 30, Table 4
#>  5:  eskridge1990_EV_5 Eskridge (1990), Crop Science 30, Table 4
#>  6:  eskridge1990_EV_6 Eskridge (1990), Crop Science 30, Table 4
#>  7:  eskridge1990_EV_7 Eskridge (1990), Crop Science 30, Table 4
#>  8:  eskridge1990_EV_8 Eskridge (1990), Crop Science 30, Table 4
#>  9:  eskridge1990_EV_9 Eskridge (1990), Crop Science 30, Table 4
#> 10: eskridge1990_EV_10 Eskridge (1990), Crop Science 30, Table 4
#> 11: eskridge1990_EV_11 Eskridge (1990), Crop Science 30, Table 4
#> 12: eskridge1990_EV_12 Eskridge (1990), Crop Science 30, Table 4
#> 13:  eskridge1990_FW_1 Eskridge (1990), Crop Science 30, Table 4
#> 14:  eskridge1990_FW_2 Eskridge (1990), Crop Science 30, Table 4
#> 15:  eskridge1990_FW_3 Eskridge (1990), Crop Science 30, Table 4
#> 16:  eskridge1990_FW_4 Eskridge (1990), Crop Science 30, Table 4
#> 17:  eskridge1990_FW_5 Eskridge (1990), Crop Science 30, Table 4
#> 18:  eskridge1990_FW_6 Eskridge (1990), Crop Science 30, Table 4
#> 19:  eskridge1990_FW_7 Eskridge (1990), Crop Science 30, Table 4
#> 20:  eskridge1990_FW_8 Eskridge (1990), Crop Science 30, Table 4
#> 21:  eskridge1990_FW_9 Eskridge (1990), Crop Science 30, Table 4
#> 22: eskridge1990_FW_10 Eskridge (1990), Crop Science 30, Table 4
#> 23: eskridge1990_FW_11 Eskridge (1990), Crop Science 30, Table 4
#> 24: eskridge1990_FW_12 Eskridge (1990), Crop Science 30, Table 4
#> 25:  eskridge1990_SH_1 Eskridge (1990), Crop Science 30, Table 4
#> 26:  eskridge1990_SH_2 Eskridge (1990), Crop Science 30, Table 4
#> 27:  eskridge1990_SH_3 Eskridge (1990), Crop Science 30, Table 4
#> 28:  eskridge1990_SH_4 Eskridge (1990), Crop Science 30, Table 4
#> 29:  eskridge1990_SH_5 Eskridge (1990), Crop Science 30, Table 4
#> 30:  eskridge1990_SH_6 Eskridge (1990), Crop Science 30, Table 4
#> 31:  eskridge1990_SH_7 Eskridge (1990), Crop Science 30, Table 4
#> 32:  eskridge1990_SH_8 Eskridge (1990), Crop Science 30, Table 4
#> 33:  eskridge1990_SH_9 Eskridge (1990), Crop Science 30, Table 4
#> 34: eskridge1990_SH_10 Eskridge (1990), Crop Science 30, Table 4
#> 35: eskridge1990_SH_11 Eskridge (1990), Crop Science 30, Table 4
#> 36: eskridge1990_SH_12 Eskridge (1990), Crop Science 30, Table 4
#> 37:  eskridge1990_ER_1 Eskridge (1990), Crop Science 30, Table 4
#> 38:  eskridge1990_ER_2 Eskridge (1990), Crop Science 30, Table 4
#> 39:  eskridge1990_ER_3 Eskridge (1990), Crop Science 30, Table 4
#> 40:  eskridge1990_ER_4 Eskridge (1990), Crop Science 30, Table 4
#> 41:  eskridge1990_ER_5 Eskridge (1990), Crop Science 30, Table 4
#> 42:  eskridge1990_ER_6 Eskridge (1990), Crop Science 30, Table 4
#> 43:  eskridge1990_ER_7 Eskridge (1990), Crop Science 30, Table 4
#> 44:  eskridge1990_ER_8 Eskridge (1990), Crop Science 30, Table 4
#> 45:  eskridge1990_ER_9 Eskridge (1990), Crop Science 30, Table 4
#> 46: eskridge1990_ER_10 Eskridge (1990), Crop Science 30, Table 4
#> 47: eskridge1990_ER_11 Eskridge (1990), Crop Science 30, Table 4
#> 48: eskridge1990_ER_12 Eskridge (1990), Crop Science 30, Table 4
#>                     id                                    source
#>                 <char>                                    <char>
#>                                                        quantity
#>                                                          <char>
#>  1:  EV of entry 1 (alpha = 0.05) from the summaries of Table 3
#>  2:  EV of entry 2 (alpha = 0.05) from the summaries of Table 3
#>  3:  EV of entry 3 (alpha = 0.05) from the summaries of Table 3
#>  4:  EV of entry 4 (alpha = 0.05) from the summaries of Table 3
#>  5:  EV of entry 5 (alpha = 0.05) from the summaries of Table 3
#>  6:  EV of entry 6 (alpha = 0.05) from the summaries of Table 3
#>  7:  EV of entry 7 (alpha = 0.05) from the summaries of Table 3
#>  8:  EV of entry 8 (alpha = 0.05) from the summaries of Table 3
#>  9:  EV of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 10: EV of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 11: EV of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 12: EV of entry 12 (alpha = 0.05) from the summaries of Table 3
#> 13:  FW of entry 1 (alpha = 0.05) from the summaries of Table 3
#> 14:  FW of entry 2 (alpha = 0.05) from the summaries of Table 3
#> 15:  FW of entry 3 (alpha = 0.05) from the summaries of Table 3
#> 16:  FW of entry 4 (alpha = 0.05) from the summaries of Table 3
#> 17:  FW of entry 5 (alpha = 0.05) from the summaries of Table 3
#> 18:  FW of entry 6 (alpha = 0.05) from the summaries of Table 3
#> 19:  FW of entry 7 (alpha = 0.05) from the summaries of Table 3
#> 20:  FW of entry 8 (alpha = 0.05) from the summaries of Table 3
#> 21:  FW of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 22: FW of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 23: FW of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 24: FW of entry 12 (alpha = 0.05) from the summaries of Table 3
#> 25:  SH of entry 1 (alpha = 0.05) from the summaries of Table 3
#> 26:  SH of entry 2 (alpha = 0.05) from the summaries of Table 3
#> 27:  SH of entry 3 (alpha = 0.05) from the summaries of Table 3
#> 28:  SH of entry 4 (alpha = 0.05) from the summaries of Table 3
#> 29:  SH of entry 5 (alpha = 0.05) from the summaries of Table 3
#> 30:  SH of entry 6 (alpha = 0.05) from the summaries of Table 3
#> 31:  SH of entry 7 (alpha = 0.05) from the summaries of Table 3
#> 32:  SH of entry 8 (alpha = 0.05) from the summaries of Table 3
#> 33:  SH of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 34: SH of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 35: SH of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 36: SH of entry 12 (alpha = 0.05) from the summaries of Table 3
#> 37:  ER of entry 1 (alpha = 0.05) from the summaries of Table 3
#> 38:  ER of entry 2 (alpha = 0.05) from the summaries of Table 3
#> 39:  ER of entry 3 (alpha = 0.05) from the summaries of Table 3
#> 40:  ER of entry 4 (alpha = 0.05) from the summaries of Table 3
#> 41:  ER of entry 5 (alpha = 0.05) from the summaries of Table 3
#> 42:  ER of entry 6 (alpha = 0.05) from the summaries of Table 3
#> 43:  ER of entry 7 (alpha = 0.05) from the summaries of Table 3
#> 44:  ER of entry 8 (alpha = 0.05) from the summaries of Table 3
#> 45:  ER of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 46: ER of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 47: ER of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 48: ER of entry 12 (alpha = 0.05) from the summaries of Table 3
#>                                                        quantity
#>                                                          <char>
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                               call
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                             <char>
#>  1: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  2: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  3: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  4: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  5: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  6: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  7: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  8: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  9: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 10: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 11: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 12: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 13: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 14: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 15: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 16: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 17: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 18: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 19: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 20: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 21: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 22: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 23: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 24: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 25: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 26: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 27: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 28: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 29: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 30: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 31: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 32: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 33: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 34: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 35: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 36: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 37: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 38: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 39: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 40: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 41: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 42: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 43: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 44: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 45: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 46: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 47: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 48: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                               call
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                             <char>
#>      extract value digits level
#>       <char> <num>  <int> <int>
#>  1:  x$EV[1]  2.24      2     1
#>  2:  x$EV[2]  2.04      2     1
#>  3:  x$EV[3]  1.72      2     1
#>  4:  x$EV[4]  2.14      2     1
#>  5:  x$EV[5]  1.46      2     1
#>  6:  x$EV[6]  1.90      2     1
#>  7:  x$EV[7]  1.67      2     1
#>  8:  x$EV[8]  1.74      2     1
#>  9:  x$EV[9]  2.83      2     1
#> 10: x$EV[10]  2.36      2     1
#> 11: x$EV[11]  2.26      2     1
#> 12: x$EV[12]  2.49      2     1
#> 13:  x$FW[1]  4.52      2     1
#> 14:  x$FW[2]  4.38      2     1
#> 15:  x$FW[3]  3.76      2     1
#> 16:  x$FW[4]  3.81      2     1
#> 17:  x$FW[5]  3.87      2     1
#> 18:  x$FW[6]  4.28      2     1
#> 19:  x$FW[7]  4.03      2     1
#> 20:  x$FW[8]  4.17      2     1
#> 21:  x$FW[9]  4.16      2     1
#> 22: x$FW[10]  4.33      2     1
#> 23: x$FW[11]  3.10      2     1
#> 24: x$FW[12]  4.60      2     1
#> 25:  x$SH[1]  2.47      2     1
#> 26:  x$SH[2]  2.20      2     1
#> 27:  x$SH[3]  1.58      2     1
#> 28:  x$SH[4]  1.81      2     1
#> 29:  x$SH[5]  1.75      2     1
#> 30:  x$SH[6]  2.18      2     1
#> 31:  x$SH[7]  2.20      2     1
#> 32:  x$SH[8]  2.17      2     1
#> 33:  x$SH[9]  2.28      2     1
#> 34: x$SH[10]  2.19      2     1
#> 35: x$SH[11]  1.41      2     1
#> 36: x$SH[12]  2.24      2     1
#> 37:  x$ER[1]  4.34      2     1
#> 38:  x$ER[2]  3.85      2     1
#> 39:  x$ER[3]  3.12      2     1
#> 40:  x$ER[4]  3.35      2     1
#> 41:  x$ER[5]  3.18      2     1
#> 42:  x$ER[6]  3.70      2     1
#> 43:  x$ER[7]  3.67      2     1
#> 44:  x$ER[8]  3.55      2     1
#> 45:  x$ER[9]  3.78      2     1
#> 46: x$ER[10]  3.76      2     1
#> 47: x$ER[11]  2.86      2     1
#> 48: x$ER[12]  3.73      2     1
#>      extract value digits level
#>       <char> <num>  <int> <int>
#>                                                                                                                                        note
#>                                                                                                                                      <char>
#>  1:                                                                                                                                        
#>  2:                                                                                                                                        
#>  3:                                                                                                                                        
#>  4:                                                                                                                                        
#>  5:                                                                                                                                        
#>  6:                                                                                                                                        
#>  7:                                                                                                                                        
#>  8:                                                                                                                                        
#>  9:                                                                                                                                        
#> 10:                                                                                                                                        
#> 11:                                                                           inputs of Table 3 are rounded to two decimals; computed 2.265
#> 12:                                                                                                                                        
#> 13:                                                                                                                                        
#> 14:                                                                                                                                        
#> 15:                                                                                                                                        
#> 16:                                                                                                                                        
#> 17:                                                                                                                                        
#> 18:                                                                                                                                        
#> 19:                                                                                                                                        
#> 20: Table 3 prints b = 1.27 for entry 8 (as for entry 7); the published FW and ER fit b = 1.22, so the printed slope is taken as a misprint
#> 21:                                                                                                                                        
#> 22:                                                                                                                                        
#> 23:                                                                                                                                        
#> 24:                                                                                                                                        
#> 25:                                                                                                                                        
#> 26:                                                                                                                                        
#> 27:                                                                                                                                        
#> 28:                                                                                                                                        
#> 29:                                                                                                                                        
#> 30:                                                                                                                                        
#> 31:                                                                                                                                        
#> 32:                                                                                                                                        
#> 33:                                                                                                                                        
#> 34:                                                                                                                                        
#> 35:                                                                                                                                        
#> 36:                                                                                                                                        
#> 37:                                                                           inputs of Table 3 are rounded to two decimals; computed 4.356
#> 38:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.844
#> 39:                                                                                                                                        
#> 40:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.363
#> 41:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.174
#> 42:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.709
#> 43:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.679
#> 44:                                                                        see the note on FW of entry 8: the published value fits b = 1.22
#> 45:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.763
#> 46:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.751
#> 47:                                                                           inputs of Table 3 are rounded to two decimals; computed 2.854
#> 48:                                                                                                                                        
#>                                                                                                                                        note
#>                                                                                                                                      <char>
# \donttest{
stab_reference(
    source   = "."
  , evaluate = TRUE
  , level    = 1
  )
#>                     id                                    source
#>                 <char>                                    <char>
#>  1:  eskridge1990_EV_1 Eskridge (1990), Crop Science 30, Table 4
#>  2:  eskridge1990_EV_2 Eskridge (1990), Crop Science 30, Table 4
#>  3:  eskridge1990_EV_3 Eskridge (1990), Crop Science 30, Table 4
#>  4:  eskridge1990_EV_4 Eskridge (1990), Crop Science 30, Table 4
#>  5:  eskridge1990_EV_5 Eskridge (1990), Crop Science 30, Table 4
#>  6:  eskridge1990_EV_6 Eskridge (1990), Crop Science 30, Table 4
#>  7:  eskridge1990_EV_7 Eskridge (1990), Crop Science 30, Table 4
#>  8:  eskridge1990_EV_8 Eskridge (1990), Crop Science 30, Table 4
#>  9:  eskridge1990_EV_9 Eskridge (1990), Crop Science 30, Table 4
#> 10: eskridge1990_EV_10 Eskridge (1990), Crop Science 30, Table 4
#> 11: eskridge1990_EV_11 Eskridge (1990), Crop Science 30, Table 4
#> 12: eskridge1990_EV_12 Eskridge (1990), Crop Science 30, Table 4
#> 13:  eskridge1990_FW_1 Eskridge (1990), Crop Science 30, Table 4
#> 14:  eskridge1990_FW_2 Eskridge (1990), Crop Science 30, Table 4
#> 15:  eskridge1990_FW_3 Eskridge (1990), Crop Science 30, Table 4
#> 16:  eskridge1990_FW_4 Eskridge (1990), Crop Science 30, Table 4
#> 17:  eskridge1990_FW_5 Eskridge (1990), Crop Science 30, Table 4
#> 18:  eskridge1990_FW_6 Eskridge (1990), Crop Science 30, Table 4
#> 19:  eskridge1990_FW_7 Eskridge (1990), Crop Science 30, Table 4
#> 20:  eskridge1990_FW_8 Eskridge (1990), Crop Science 30, Table 4
#> 21:  eskridge1990_FW_9 Eskridge (1990), Crop Science 30, Table 4
#> 22: eskridge1990_FW_10 Eskridge (1990), Crop Science 30, Table 4
#> 23: eskridge1990_FW_11 Eskridge (1990), Crop Science 30, Table 4
#> 24: eskridge1990_FW_12 Eskridge (1990), Crop Science 30, Table 4
#> 25:  eskridge1990_SH_1 Eskridge (1990), Crop Science 30, Table 4
#> 26:  eskridge1990_SH_2 Eskridge (1990), Crop Science 30, Table 4
#> 27:  eskridge1990_SH_3 Eskridge (1990), Crop Science 30, Table 4
#> 28:  eskridge1990_SH_4 Eskridge (1990), Crop Science 30, Table 4
#> 29:  eskridge1990_SH_5 Eskridge (1990), Crop Science 30, Table 4
#> 30:  eskridge1990_SH_6 Eskridge (1990), Crop Science 30, Table 4
#> 31:  eskridge1990_SH_7 Eskridge (1990), Crop Science 30, Table 4
#> 32:  eskridge1990_SH_8 Eskridge (1990), Crop Science 30, Table 4
#> 33:  eskridge1990_SH_9 Eskridge (1990), Crop Science 30, Table 4
#> 34: eskridge1990_SH_10 Eskridge (1990), Crop Science 30, Table 4
#> 35: eskridge1990_SH_11 Eskridge (1990), Crop Science 30, Table 4
#> 36: eskridge1990_SH_12 Eskridge (1990), Crop Science 30, Table 4
#> 37:  eskridge1990_ER_1 Eskridge (1990), Crop Science 30, Table 4
#> 38:  eskridge1990_ER_2 Eskridge (1990), Crop Science 30, Table 4
#> 39:  eskridge1990_ER_3 Eskridge (1990), Crop Science 30, Table 4
#> 40:  eskridge1990_ER_4 Eskridge (1990), Crop Science 30, Table 4
#> 41:  eskridge1990_ER_5 Eskridge (1990), Crop Science 30, Table 4
#> 42:  eskridge1990_ER_6 Eskridge (1990), Crop Science 30, Table 4
#> 43:  eskridge1990_ER_7 Eskridge (1990), Crop Science 30, Table 4
#> 44:  eskridge1990_ER_8 Eskridge (1990), Crop Science 30, Table 4
#> 45:  eskridge1990_ER_9 Eskridge (1990), Crop Science 30, Table 4
#> 46: eskridge1990_ER_10 Eskridge (1990), Crop Science 30, Table 4
#> 47: eskridge1990_ER_11 Eskridge (1990), Crop Science 30, Table 4
#> 48: eskridge1990_ER_12 Eskridge (1990), Crop Science 30, Table 4
#>                     id                                    source
#>                 <char>                                    <char>
#>                                                        quantity
#>                                                          <char>
#>  1:  EV of entry 1 (alpha = 0.05) from the summaries of Table 3
#>  2:  EV of entry 2 (alpha = 0.05) from the summaries of Table 3
#>  3:  EV of entry 3 (alpha = 0.05) from the summaries of Table 3
#>  4:  EV of entry 4 (alpha = 0.05) from the summaries of Table 3
#>  5:  EV of entry 5 (alpha = 0.05) from the summaries of Table 3
#>  6:  EV of entry 6 (alpha = 0.05) from the summaries of Table 3
#>  7:  EV of entry 7 (alpha = 0.05) from the summaries of Table 3
#>  8:  EV of entry 8 (alpha = 0.05) from the summaries of Table 3
#>  9:  EV of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 10: EV of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 11: EV of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 12: EV of entry 12 (alpha = 0.05) from the summaries of Table 3
#> 13:  FW of entry 1 (alpha = 0.05) from the summaries of Table 3
#> 14:  FW of entry 2 (alpha = 0.05) from the summaries of Table 3
#> 15:  FW of entry 3 (alpha = 0.05) from the summaries of Table 3
#> 16:  FW of entry 4 (alpha = 0.05) from the summaries of Table 3
#> 17:  FW of entry 5 (alpha = 0.05) from the summaries of Table 3
#> 18:  FW of entry 6 (alpha = 0.05) from the summaries of Table 3
#> 19:  FW of entry 7 (alpha = 0.05) from the summaries of Table 3
#> 20:  FW of entry 8 (alpha = 0.05) from the summaries of Table 3
#> 21:  FW of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 22: FW of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 23: FW of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 24: FW of entry 12 (alpha = 0.05) from the summaries of Table 3
#> 25:  SH of entry 1 (alpha = 0.05) from the summaries of Table 3
#> 26:  SH of entry 2 (alpha = 0.05) from the summaries of Table 3
#> 27:  SH of entry 3 (alpha = 0.05) from the summaries of Table 3
#> 28:  SH of entry 4 (alpha = 0.05) from the summaries of Table 3
#> 29:  SH of entry 5 (alpha = 0.05) from the summaries of Table 3
#> 30:  SH of entry 6 (alpha = 0.05) from the summaries of Table 3
#> 31:  SH of entry 7 (alpha = 0.05) from the summaries of Table 3
#> 32:  SH of entry 8 (alpha = 0.05) from the summaries of Table 3
#> 33:  SH of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 34: SH of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 35: SH of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 36: SH of entry 12 (alpha = 0.05) from the summaries of Table 3
#> 37:  ER of entry 1 (alpha = 0.05) from the summaries of Table 3
#> 38:  ER of entry 2 (alpha = 0.05) from the summaries of Table 3
#> 39:  ER of entry 3 (alpha = 0.05) from the summaries of Table 3
#> 40:  ER of entry 4 (alpha = 0.05) from the summaries of Table 3
#> 41:  ER of entry 5 (alpha = 0.05) from the summaries of Table 3
#> 42:  ER of entry 6 (alpha = 0.05) from the summaries of Table 3
#> 43:  ER of entry 7 (alpha = 0.05) from the summaries of Table 3
#> 44:  ER of entry 8 (alpha = 0.05) from the summaries of Table 3
#> 45:  ER of entry 9 (alpha = 0.05) from the summaries of Table 3
#> 46: ER of entry 10 (alpha = 0.05) from the summaries of Table 3
#> 47: ER of entry 11 (alpha = 0.05) from the summaries of Table 3
#> 48: ER of entry 12 (alpha = 0.05) from the summaries of Table 3
#>                                                        quantity
#>                                                          <char>
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                               call
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                             <char>
#>  1: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  2: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  3: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  4: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  5: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  6: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  7: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  8: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>  9: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 10: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 11: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 12: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 13: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 14: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 15: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 16: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 17: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 18: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 19: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 20: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 21: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 22: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 23: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 24: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 25: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 26: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 27: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 28: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 29: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 30: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 31: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 32: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 33: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 34: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 35: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 36: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 37: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 38: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 39: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 40: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 41: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 42: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 43: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 44: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 45: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 46: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 47: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#> 48: .stab_safety_core(m = c(4.69, 4.49, 3.93, 4.16, 4.20, 4.54, 4.62, 4.65, 4.68, 4.53, 3.86, 4.60), s2 = c(2.22, 2.21, 1.81, 1.51, 2.77, 2.57, 3.22, 3.13, 1.27, 1.74, 0.94, 1.65), b = c(1.08, 1.05, 0.92, 0.84, 1.15, 1.12, 1.27, 1.27, 0.76, 0.91, 0.65, 1.00), sh = c(0.02, 0.14, 0.24, 0.24, 0.41, 0.26, 0.37, 0.48, 0.33, 0.22, 0.42, 0.26), s2d = c(0.03, 0.15, 0.23, 0.19, 0.35, 0.23, 0.20, 0.36, 0.21, 0.21, 0.16, 0.28), sy2 = 1.87, se2 = 1.80, E = 15, alpha = 0.05)
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                               call
#>                                                                                                                                                                                                                                                                                                                                                                                                                                                                             <char>
#>      extract value digits level
#>       <char> <num>  <int> <int>
#>  1:  x$EV[1]  2.24      2     1
#>  2:  x$EV[2]  2.04      2     1
#>  3:  x$EV[3]  1.72      2     1
#>  4:  x$EV[4]  2.14      2     1
#>  5:  x$EV[5]  1.46      2     1
#>  6:  x$EV[6]  1.90      2     1
#>  7:  x$EV[7]  1.67      2     1
#>  8:  x$EV[8]  1.74      2     1
#>  9:  x$EV[9]  2.83      2     1
#> 10: x$EV[10]  2.36      2     1
#> 11: x$EV[11]  2.26      2     1
#> 12: x$EV[12]  2.49      2     1
#> 13:  x$FW[1]  4.52      2     1
#> 14:  x$FW[2]  4.38      2     1
#> 15:  x$FW[3]  3.76      2     1
#> 16:  x$FW[4]  3.81      2     1
#> 17:  x$FW[5]  3.87      2     1
#> 18:  x$FW[6]  4.28      2     1
#> 19:  x$FW[7]  4.03      2     1
#> 20:  x$FW[8]  4.17      2     1
#> 21:  x$FW[9]  4.16      2     1
#> 22: x$FW[10]  4.33      2     1
#> 23: x$FW[11]  3.10      2     1
#> 24: x$FW[12]  4.60      2     1
#> 25:  x$SH[1]  2.47      2     1
#> 26:  x$SH[2]  2.20      2     1
#> 27:  x$SH[3]  1.58      2     1
#> 28:  x$SH[4]  1.81      2     1
#> 29:  x$SH[5]  1.75      2     1
#> 30:  x$SH[6]  2.18      2     1
#> 31:  x$SH[7]  2.20      2     1
#> 32:  x$SH[8]  2.17      2     1
#> 33:  x$SH[9]  2.28      2     1
#> 34: x$SH[10]  2.19      2     1
#> 35: x$SH[11]  1.41      2     1
#> 36: x$SH[12]  2.24      2     1
#> 37:  x$ER[1]  4.34      2     1
#> 38:  x$ER[2]  3.85      2     1
#> 39:  x$ER[3]  3.12      2     1
#> 40:  x$ER[4]  3.35      2     1
#> 41:  x$ER[5]  3.18      2     1
#> 42:  x$ER[6]  3.70      2     1
#> 43:  x$ER[7]  3.67      2     1
#> 44:  x$ER[8]  3.55      2     1
#> 45:  x$ER[9]  3.78      2     1
#> 46: x$ER[10]  3.76      2     1
#> 47: x$ER[11]  2.86      2     1
#> 48: x$ER[12]  3.73      2     1
#>      extract value digits level
#>       <char> <num>  <int> <int>
#>                                                                                                                                        note
#>                                                                                                                                      <char>
#>  1:                                                                                                                                        
#>  2:                                                                                                                                        
#>  3:                                                                                                                                        
#>  4:                                                                                                                                        
#>  5:                                                                                                                                        
#>  6:                                                                                                                                        
#>  7:                                                                                                                                        
#>  8:                                                                                                                                        
#>  9:                                                                                                                                        
#> 10:                                                                                                                                        
#> 11:                                                                           inputs of Table 3 are rounded to two decimals; computed 2.265
#> 12:                                                                                                                                        
#> 13:                                                                                                                                        
#> 14:                                                                                                                                        
#> 15:                                                                                                                                        
#> 16:                                                                                                                                        
#> 17:                                                                                                                                        
#> 18:                                                                                                                                        
#> 19:                                                                                                                                        
#> 20: Table 3 prints b = 1.27 for entry 8 (as for entry 7); the published FW and ER fit b = 1.22, so the printed slope is taken as a misprint
#> 21:                                                                                                                                        
#> 22:                                                                                                                                        
#> 23:                                                                                                                                        
#> 24:                                                                                                                                        
#> 25:                                                                                                                                        
#> 26:                                                                                                                                        
#> 27:                                                                                                                                        
#> 28:                                                                                                                                        
#> 29:                                                                                                                                        
#> 30:                                                                                                                                        
#> 31:                                                                                                                                        
#> 32:                                                                                                                                        
#> 33:                                                                                                                                        
#> 34:                                                                                                                                        
#> 35:                                                                                                                                        
#> 36:                                                                                                                                        
#> 37:                                                                           inputs of Table 3 are rounded to two decimals; computed 4.356
#> 38:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.844
#> 39:                                                                                                                                        
#> 40:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.363
#> 41:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.174
#> 42:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.709
#> 43:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.679
#> 44:                                                                        see the note on FW of entry 8: the published value fits b = 1.22
#> 45:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.763
#> 46:                                                                           inputs of Table 3 are rounded to two decimals; computed 3.751
#> 47:                                                                           inputs of Table 3 are rounded to two decimals; computed 2.854
#> 48:                                                                                                                                        
#>                                                                                                                                        note
#>                                                                                                                                      <char>
#>     computed agrees
#>        <num> <lgcl>
#>  1: 2.239223   TRUE
#>  2: 2.044749   TRUE
#>  3: 1.717076   TRUE
#>  4: 2.138770   TRUE
#>  5: 1.462418   TRUE
#>  6: 1.903099   TRUE
#>  7: 1.668416   TRUE
#>  8: 1.739957   TRUE
#>  9: 2.826344   TRUE
#> 10: 2.360289   TRUE
#> 11: 2.265255  FALSE
#> 12: 2.487147   TRUE
#> 13: 4.516157   TRUE
#> 14: 4.381348   TRUE
#> 15: 3.756157   TRUE
#> 16: 3.812315   TRUE
#> 17: 3.874045   TRUE
#> 18: 4.279236   TRUE
#> 19: 4.033281   TRUE
#> 20: 4.063281  FALSE
#> 21: 4.158472   TRUE
#> 22: 4.334427   TRUE
#> 23: 3.099438   TRUE
#> 24: 4.600000   TRUE
#> 25: 2.470971   TRUE
#> 26: 2.198984   TRUE
#> 27: 1.580679   TRUE
#> 28: 1.810679   TRUE
#> 29: 1.754749   TRUE
#> 30: 2.179191   TRUE
#> 31: 2.196979   TRUE
#> 32: 2.166325   TRUE
#> 33: 2.279415   TRUE
#> 34: 2.192224   TRUE
#> 35: 1.409223   TRUE
#> 36: 2.239191   TRUE
#> 37: 4.356252  FALSE
#> 38: 3.843752  FALSE
#> 39: 3.122228   TRUE
#> 40: 3.363170  FALSE
#> 41: 3.173751  FALSE
#> 42: 3.709173  FALSE
#> 43: 3.679071  FALSE
#> 44: 3.501856  FALSE
#> 45: 3.763400  FALSE
#> 46: 3.751275  FALSE
#> 47: 2.854345  FALSE
#> 48: 3.729625   TRUE
#>     computed agrees
#>        <num> <lgcl>
# }