RI нужна помощь в понимании кода R моего профессора, чтобы вычислить точный доверительный интервал

Question

Мне трудно понять конкретные строки в примере моего профессора для вычисления точного доверительного интервала (примечание: это не использует допущение нормальности)

# R code for Exact CI

n<- 19; k<- 3
alp<- 0.05; mar<- 1.0e-4

B <- 10^5; x <- matrix(0, B, 1)

# Calculate the left endpoint of the exact CI
for (i in 0:B) {
    p <- i/B
    x[i] <- pbinom((k-1), n, p)
}
x[(abs(x - (1-alp/2)) <= mar)] 
pi_L <- mean(which(abs(x - (1-alp/2)) <= mar)/B)

# Calculate the right endpoint of the exact CI
for (i in 0:B) {
    p <- i/B
    x[i] <- pbinom(k, n, p)
}
x[(abs(x - alp/2) <= mar)] 
pi_U <- mean(which(abs(x - alp/2) <= mar)/B)

# The obtained exact CI
pi_L; pi_U

Части, которые я не понимаю:

x[(abs(x - (1-alp/2)) <= mar)] 
pi_L <- mean(which(abs(x - (1-alp/2)) <= mar)/B)
#####
x[(abs(x - alp/2) <= mar)] 
pi_U <- mean(which(abs(x - alp/2) <= mar)/B)

Я не знаю, какую роль в этом играет «<= mar» или «which (abs (x - alp / 2)»). Любая помощь будет признательна, спасибо. </p>