Maximiza uma funcao de log-verossimilhanca via stats::optim() e estima
erros-padrao pela inversa da informacao de Fisher observada (Hessiana).
Arguments
- log_veross
Funcao
function(theta) ...que devolve a log-verossimilhanca (escalar) para o vetor de parametrostheta.- inicio
Vetor de valores iniciais para
theta.- nomes
Vetor opcional de nomes dos parametros.
- metodo
Metodo de
stats::optim().- digits
Inteiro.
Value
Uma lista com estimativas (tibble: parametro, estimativa,
erro_padrao, z, ic_inf, ic_sup) e ajuste (tibble:
log_veross, aic, bic... convergiu).
See also
Other inferencia:
rnp_bayes_conjugada(),
rnp_bootstrap(),
rnp_bootstrap_parametrico(),
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
set.seed(1); x <- rnorm(200, 5, 2)
ll <- function(th) sum(dnorm(x, th[1], th[2], log = TRUE))
rnp_emv(ll, inicio = c(0, 1), nomes = c("media", "dp"))$estimativas
#> # A tibble: 2 × 6
#> parametro estimativa erro_padrao z ic_inf ic_sup
#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 media 5.07 0.131 38.7 4.81 5.33
#> 2 dp 1.85 0.0927 20.0 1.67 2.04