Skip to contents

Estimacao bayesiana conjugada

Source: R/verossimilhanca.R
rnp_bayes_conjugada.Rd

Atualiza distribuicoes a priori conjugadas com os dados para tres modelos classicos: Beta-Binomial, Gama-Poisson e Normal-Normal (variancia conhecida).

Usage

rnp_bayes_conjugada(
  familia = c("beta_binomial", "gama_poisson", "normal_normal"),
  priori,
  dados,
  conf = 0.95,
  digits = 4L
)

Arguments

familia

String: "beta_binomial", "gama_poisson" ou "normal_normal".

priori

Vetor nomeado de hiperparametros a priori (beta_binomial: a, b; gama_poisson: a, b; normal_normal: mu0, sigma0).

dados

Lista com os dados (beta_binomial: sucessos, n; gama_poisson: soma, n; normal_normal: x e sigma conhecido).

conf

Nivel de credibilidade do intervalo.

digits

Inteiro.

Value

tibble com parametro, valor (hiperparametros a posteriori), alem de media_post, ic_inf, ic_sup.

See also

Other inferencia: rnp_bootstrap(), rnp_bootstrap_parametrico(), rnp_emv(), rnp_ic_bootstrap(), rnp_ic_verossimilhanca(), rnp_informacao_fisher(), rnp_jackknife(), rnp_log_verossimilhanca(), rnp_metodo_momentos(), rnp_poder_teste(), rnp_tamanho_amostra_teste(), rnp_teste_aderencia(), rnp_teste_binomial(), rnp_teste_grubbs(), rnp_teste_ks(), rnp_teste_normalidade(), rnp_teste_permutacao(), rnp_teste_proporcoes(), rnp_teste_qui_quadrado(), rnp_teste_razao_veross(), rnp_teste_runs(), rnp_teste_score(), rnp_teste_sinais(), rnp_teste_wald()

Examples

rnp_bayes_conjugada("beta_binomial", priori = c(a = 1, b = 1),
                    dados = list(sucessos = 8, n = 10))
#> # A tibble: 2 × 5
#>   parametro valor media_post ic_inf ic_sup
#>   <chr>     <dbl>      <dbl>  <dbl>  <dbl>
#> 1 a             9       0.75  0.482  0.940
#> 2 b             3       0.75  0.482  0.940