Necesito obtener los percentiles 5 % y 95% en la salida de la instrucción bugs en R2openBUGS, para poder hacer un contraste de hipótesis unilateral, con un intervalo de credibilidad bayesiano del 95%. Les adjunto el código en r y openBUGS. El problema es que no encuentro la instrucción dentro de Bugs en R2openBUGS para poder indicar o setear más percentiles o cuantiles. ¿Pueden ustedes por favor ayudarme en esta tarea? Muchas gracias por la ayuda.
library("R2OpenBUGS")
library("coda")
library("lattice")
n4 = 10000
n5 = 1000
n6 = 1
n7 = 1
N= 10
r = 8
a = 1
b = 1
mod<-"mbp2022.txt"
data <- list( N = N,
r = r,
a = a,
b = b
)
sink(mod)
cat("
model{
pi1 ~ dbeta(a,b)
r ~ dbin(pi1, N)
}
",fill=TRUE)
sink()
modelfile <- mod
vars <-list("pi1")
inits <-function(){list(pi1=0.5)}
out <-bugs(
model.file=file.path(modelfile),
data=data,
inits=inits,
parameters.to.save=vars,
n.iter=n4,
n.burnin=n5,
n.chains=n6,
n.thin=n7,
summary.only = F,
debug=F,
DIC=F,
working.directory = getwd(),
)
cat("\n")
print (out, digits = 03)