Coleta estatística do ajuste

Coleta as estimativas e suas estatísticas para um objeto de ajuste do modelo beta ordinal com censura intervalar. Coleta também as estatísticas de bondade do ajuste como a log-verossimilhança, o AIC e o BIC do modelo

Usage

betaregscale_coef(fit, alpha = 0.05)

Arguments

fit

Objeto do sjuste retornado das funções betaregscale_fit e betaregscale_fit_z

alpha

Nível de significância do alpha para os intervalos de confiança. Padrão 0.05 bilateral.

Value

Lista contendo as estimativas coeficientes, suas estatísticas e a bondade do ajuste.

Examples

# Simulando dados
n <- 50
fx <- ~ x1 + x2
fz <- ~ z1
dados <- data.frame(
  x1 = rnorm(n),
  x2 = rnorm(n),
  z1 = runif(n),
  z2 = runif(n)
 )
 dados_simulados <- betaregscale_simula_dados_z(
   formula_x = fx,
   formula_z = fz,
   dados = dados,
   betas = c(0.2,-0.5, 0.3),
   zetas = c(1, 1.2),
   link = "logit",
   link_phi = "logit",
   ncuts = 100
   )

 fit_z <- betaregscale_fit_z(
  formula = y ~ x1 + x2,
  dados = dados_simulados,
  link = "logit",
  link_phi = "logit",
  num_hessiana = TRUE
  )

coe <- betaregscale_coef(fit_z)
coe$est
#>      variable   estimate    ci_lower  ci_upper        se    t_value
#> 1 (Intercept)  0.3392146 -0.04532723 0.7237565 0.1961984  1.7289365
#> 2          x1 -0.2523367 -0.76791961 0.2632462 0.2630573 -0.9592459
#> 3          x2  0.5046937  0.05678142 0.9526060 0.2285309  2.2084268
#> 4         phi  0.8237284  0.50488388 1.1425728 0.1626787  5.0635279
#>        p_value
#> 1 9.052750e-02
#> 2 3.424499e-01
#> 3 3.223702e-02
#> 4 7.100409e-06
coe$gof
#>      logLik      AIC      BIC
#> 1 -145.4782 300.9564 310.5165