est_q: Encontra o quantil da distribuição Qui-Quadrado inf.

In prdm0/LikRatioTest: Teste da Razao de Verossimilhanca Generalizada

Description

Encontra o quantil da distribuição Qui-Quadrado inf.

Usage

est_q(fn, alpha = 0.05, bilateral = FALSE, step = 1e-04, c, k)

Arguments

fn	Recebe a distribuição Qui-Quadrado inf.
alpha	Nível de significância adotado.
bilateral	Se TRUE, retorna os quantis para um teste bilateral. O padrão considera bilateral = FALSE.
step	Tamanho do passo da integral numeérica responsável pela obtenção dos quantis da Qui-Quadrado inf. O padrão considera step = 1e-4.
c	Parâmetro da distribuição Qui-Quadrado inf.
k	Parâmetro da distribuição Qui-Quadrado inf.

Details

A função recebe como argumento uma função densidade de probabilidade ou uma função de probabilidade que é passada como argumento à f. O objeto passado como argumento de f deverá ser implementado conforme os exemplos. Note que se a função for implementada segundo os exemplos, não haverá a necessidade de implementar a função sob a hipótese nula, visto que são as mesmas funções. Para especificar os parâmetros que serão fixados, i.e, para especificar a distribuição sob a hipótese nula, utiliza-se o argumento par0 que receberá uma lista formada por dois vetores. O primeiro vetor da lista de verá ser um vetor de strings com os nomes das variáveis que deseja-se fixar e o segundo vetor deverá conter os valores que serão atribuidos à cada uma das variáveis passada ao primeiro vetor.

Author(s)

Pedro Rafael D. Marinho

Examples

fdp_chisq_inf <- function(par, x) {
  k <- par[1]
  c <- par[2]
  dchisq(x = x, df = k) * (1 - (1 - pchisq(q = x, df = k)) ^ c + c * pchisq(q = x, df = k) *
                        (1 - pchisq(q = x, df = k)) ^ (c - 1))
}
est_q(fn = fdp_chisq_inf, alpha = 0.05, bilateral = FALSE, c = 1, k = 1)

prdm0/LikRatioTest documentation built on Dec. 18, 2019, 6:33 a.m.