Log-verossimilhança do modelo beta ordinal com dispersão fixa — betaregscale_log

Esta função calcula a log-verossimilhança de um conjunto de dados para uma variável beta ordinal, aplicando diferentes funções de ligação.

Usage

betaregscale_log_vero(
  param,
  formula,
  dados,
  link = "logit",
  link_phi = "logit",
  acumulada = TRUE,
  ncuts = 100,
  type = "m",
  lim = 0.5,
  repar = "2"
)

Arguments

param: Vetor numérico de parâmetros, incluindo coeficientes de regressão e o parâmetro phi.
formula: Fórmula que indica a relação entre a variável dependente e as variáveis independentes.
dados: Um conjunto de dados que contém a variável dependente e as variáveis independentes especificadas na fórmula.
link: Nome da função de ligação a ser usada. Pode ser uma das seguintes: "logit", "probit", "cauchit", "cloglog" ou "identity". O padrão é "logit".
link_phi: Função de ligação para phi. Uma dentre "log","sqrt","identity"
acumulada: Um valor lógico indicando se a log-verossimilhança acumulada deve ser calculada. O padrão é TRUE.
ncuts: Número de cortes para a variável ordinal. O padrão é 100.
type: Tipo de intervalo. "m" = meio; "l" = esquerda e "r" = direita.
lim: Região de incerteza da medida. Padrão 0.5.
repar: Tipo de reparametrização sendo 0, 1 ou 2. Padrão 2.

Value

Retorna a soma da log-verossimilhança dos dados.

Examples

# Criar um conjunto de dados de exemplo
set.seed(421)
dados <- data.frame(x1 = rnorm(100), x2 = rnorm(100))
param <- c(0, 0.5,-0.2, 1 / 5)
phi <- 1 / 5
formula <- y ~ x1 + x2
dados_simulados <- betaregscale_simula_dados(
 formula = ~ x1 + x2,
 dados = dados,
 betas = c(0, 0.5,-0.2),
 phi = phi,
 link = "logit",
 link_phi = "logit",
 ncuts = 100,
 type = "m"
 )
log_verossimilhanca <- betaregscale_log_vero(param, formula, dados_simulados)
print(log_verossimilhanca)
#> [1] -390.8556