A Base R, simple implementation of the No-Underrun Sampler. This package aims to mostly directly implement the algorithm as described by th
Changed to vapply
+1
-1
+1
-1
R/NURS.R
+1
-1
R/NURS.R
···
77
77
seq_len(n_ext),
78
78
\(i) (if (B_k) right_point else left_point) + i * (2 * B_k - 1) * h * rho
79
79
)
80
-
log_ext <- sapply(orbit_ext, logpdf)
80
+
log_ext <- vapply(orbit_ext, logpdf, numeric(1))
81
81
82
82
# Recursive sub-stopping criterion
83
83
if (NURS_sub_stop(log_ext, log_eps_h)) break