power_test: Função para o cálculo do poder do teste.

In prdm0/LikRatioTest: Teste da Razao de Verossimilhanca Generalizada

Description

Função para o cálculo do poder do teste.

Usage

power_test(
  N = 1000L,
  B = 250L,
  n,
  f,
  sig,
  q,
  kicks,
  par0,
  ncores,
  p = 0.5,
  step = 0.001,
  ...
)

Arguments

N	Número de réplicas de Monte-Carlo.
B	Número de réplicas de bootstrap.
n	Tamanho da amostra.
sig	Nível de significância.
q	Função quantílica para a geração da amostra.
kicks	Chutes iniciais.
par0	Lista de parâmetros para serem fixados na hipótese nula.
ncores	Quantidade de cores utilizados na paralelização.
p	Valor do parâmetro p.
step	Passos para das integrais calculadas para obtenção do quantil da distribuiçao Chisq-Inf.
...	Parâmetros para geração da amostra verdadeira.

Examples

# pdf_ew <- function(par, x, var = NULL) {
#   alpha <- par[1]
#   sigma <- par[2]
#   theta <- par[3]
#
#   if (is.list(var))
#     eval(parse(text = paste(var[[1]], " <- ", unlist(var[[2]]), sep = "")))
#
#   alpha * theta / sigma * (1 - exp(-(x / sigma) ^ alpha)) ^ (theta - 1) *
#     exp(-(x / sigma) ^ alpha) * (x / sigma) ^ (alpha - 1)
# }
#
# rew <- function(n, alpha, sigma, theta){
#   u <- runif(n, 0, 1)
#   sigma * (-log(1 - u ^ (1 / theta))) ^ (1 / alpha)
# }
#
# set.seed(1L, kind = "L'Ecuyer-CMRG")
#
# tictoc::tic()
# power_test(N = 1e4L,
#            B = 250L,
#            n = 150L,
#            f = pdf_ew,
#            sig = 0.05,
#            q = rew,
#            kicks = c(1, 1, 1),
#            par0 = list(c("alpha", "theta"), c(1.5, 1.7)),
#            ncores = 8L,
#            p = 0.5,
#            step = 1e-3,
#            alpha = 1.5,
#            sigma = 1.5,
#            theta = 1.7
#
# )
# tictoc::toc()

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