Evaluanción de rate function (la derivada de la mean)

Usage

imusaokomoto(t, theta)

mmusaokomoto(t, theta)

log_likelihood_region_musaokomoto(t, tau_left, tau_right, theta)

log_prior_region_musaokomoto(theta, hyperparameters)

D_log_prior_region_musaokomoto(theta, hyperparameters)

D_log_likelihood_region_musaokomoto(
  vec_d_i,
  tau1,
  tau2,
  theta,
  nhpp_dist = "W"
)

Arguments

t: valor a evaluar (real)
theta: vector de parámetros de verosimilitud del NHPP
tau_left: valor del primer punto de cambio
tau_right: valor del segundo punto de cambio
hyperparameters: Hyperparameters
vec_d_i: vector de días en los que hubo revases entre el los puntos de cambio tau1 y tau2
tau1: valor del primer punto de cambio
tau2: valor del segundo punto de cambio
nhpp_dist: nombre de tasa de NHPP

Value

regresa la evaluación

Examples

t <- c(1.4, 2.8)
theta <- c(1.2, 2.1, 3.2)

imusaokomoto(t, theta)
#> [1] 0.8076923 0.5250000

t <- c(1.4, 2.8)
theta <- c(1.2, 2.1, 3.2)
mmusaokomoto(t, theta)
#> [1] 1.623699 2.528343

log_likelihood_region_musaokomoto(DataCPSim, 0, 575, theta = c(0.5, 0.5))
#> [1] -2576.488

hyper <- data.frame(
  hyperprior_shape = c(1, 3),
  hyperprior_scale = c(2, 1.2)
)
log_prior_region_musaokomoto(theta = c(0.5, 0.5), hyper)
#> [1] -2.831777