Skip to contents

Intervalo de confianca para um parametro escalar baseado na verossimilhanca perfilada: {theta : 2 (l(hat) - l(theta)) <= qchisq(conf, 1)}.

Usage

rnp_ic_verossimilhanca(log_veross, emv, intervalo, conf = 0.95, digits = 4L)

Arguments

log_veross

Funcao escalar function(theta) ....

emv

Estimativa de maxima verossimilhanca.

intervalo

Vetor c(min, max) de busca dos limites.

conf

Nivel de confianca.

digits

Inteiro.

Value

tibble com emv, limite_inferior, limite_superior, conf.

Examples

set.seed(1); x <- rpois(100, 3)
ll <- function(l) sum(dpois(x, l, log = TRUE))
rnp_ic_verossimilhanca(ll, emv = mean(x), intervalo = c(0.5, 8))
#> # A tibble: 1 × 4
#>     emv limite_inferior limite_superior  conf
#>   <dbl>           <dbl>           <dbl> <dbl>
#> 1  3.05            2.72            3.41  0.95