Proporcijos PI: 3 ar daugiau grupių

Daugiareikšmio (k ≥ 3) kintamojo reikšmių proporcijų vienu metu skaičiuojamų pasikliautinųjų intervalų (PI) skaičiavimo funkcija, kuri yra patobulintas DescTools::MultinomCI() variantas. Rezultatas – duomenų lentelė.

ci_multinom(
  x,
  method = "goodman",
  conf.level = 0.95,
  gr_colname = "group",
  ...
)

Arguments

x

Vektorius su grupių dydžiais. Geriausia, jei vektoriaus elementai turėtų prasmingus pavadinimus (žiūrėti pavyzdžius).

method

Skaičiavimo metodas ("goodman", "sisonglaz", "cplus1" ir kiti variantai, aprašyti DescTools::MultinomCI() dokumentacijoje).

conf.level

Pasikliovimo lygmuo. Numatytoji reikšmė – 0.95.

gr_colname

Stulpelio pavadinimas (kabutėse), kuriame bus parašyti grupių pavadinimai. Numatytoji reikšmė yra "group".

...

Kiti parametrai, kuriuos priima DescTools::MultinomCI(). Žiūrėti šios funkcijos dokumentaciją.

Value

Rezultatas – duomenų lentelė, kurios stulpeliai:

  • group arba kitas vartotojo pasirinktas pavadinimas stulpeliui su grupių pavadinimams, numatytoji reikšmė (<fct>).

  • est (<dbl>) – proporcijos įvertis.

  • lwr.ci, upr.ci (<dbl>) – (lower CI, upper CI) apatinė ir viršutinė proporcijos pasikliautinojo intervalo ribos.

  • x (<int>) – Grupės dydis.

  • n (<int>) – Įvykių skaičius iš viso. / Imties dydis.

Details

Ši funkcija naudojama taip pat, kaip DescTools::MultinomCI(), tik numatytasis metodas yra Goodman metodas, o rezultatas – duomenų lentelė, o ne vektorius. Todėl rezultatą galima patogiai braižyti naudojant ggplot2.

Examples

# Dažniai be pavadinimų
ci_multinom(c(20, 35, 54))
#> # A tibble: 3 × 6
#>   group   est lwr.ci upr.ci     x     n
#>   <fct> <dbl>  <dbl>  <dbl> <int> <int>
#> 1 1     0.183  0.111  0.287    20   109
#> 2 2     0.321  0.225  0.435    35   109
#> 3 3     0.495  0.384  0.607    54   109

# Nurodytas skaičiavimo metodas
ci_multinom(c(20, 35, 54), method = "goodman")
#> # A tibble: 3 × 6
#>   group   est lwr.ci upr.ci     x     n
#>   <fct> <dbl>  <dbl>  <dbl> <int> <int>
#> 1 1     0.183  0.111  0.287    20   109
#> 2 2     0.321  0.225  0.435    35   109
#> 3 3     0.495  0.384  0.607    54   109

# Dažniai su grupių pavadinimais
x <- c("dideli" = 20, "vidutiniai" = 35, "maži" = 54)
ci_multinom(x, method = "goodman")
#> # A tibble: 3 × 6
#>   group        est lwr.ci upr.ci     x     n
#>   <fct>      <dbl>  <dbl>  <dbl> <int> <int>
#> 1 dideli     0.183  0.111  0.287    20   109
#> 2 vidutiniai 0.321  0.225  0.435    35   109
#> 3 maži       0.495  0.384  0.607    54   109

# Dažniai su grupių pavadinimais ir jungimo operatorius
c("dideli" = 20, "vidutiniai" = 35, "maži" = 54) |>
  ci_multinom()
#> # A tibble: 3 × 6
#>   group        est lwr.ci upr.ci     x     n
#>   <fct>      <dbl>  <dbl>  <dbl> <int> <int>
#> 1 dideli     0.183  0.111  0.287    20   109
#> 2 vidutiniai 0.321  0.225  0.435    35   109
#> 3 maži       0.495  0.384  0.607    54   109

# Kitas metodas
c("dideli" = 33, "vidutiniai" = 35, "maži" = 30) |>
  ci_multinom(method = "sisonglaz")
#> # A tibble: 3 × 6
#>   group        est lwr.ci upr.ci     x     n
#>   <fct>      <dbl>  <dbl>  <dbl> <int> <int>
#> 1 dideli     0.337  0.235  0.447    33    98
#> 2 vidutiniai 0.357  0.255  0.468    35    98
#> 3 maži       0.306  0.204  0.417    30    98