summary
method for class "ru"
.
print
method for an object object
of class
"summary.ru"
.
For summary.lm
: a list of the following components from
object
:
information about the ratio-of-uniforms bounding box, i.e.,
object$box
an estimate of the probability of acceptance, i.e.,
object$pa
a summary of the simulated values, via
summary(object$sim_vals)
For print.summary.ru
: the argument x
, invisibly.
# one-dimensional standard normal ----------------
x <- ru(logf = function(x) -x ^ 2 / 2, d = 1, n = 1000, init = 0)
summary(x)
#> ru bounding box:
#> box vals1 conv
#> a 1.000000 0.00000 0
#> b1minus -1.050542 -1.73205 0
#> b1plus 1.050542 1.73205 0
#>
#> estimated probability of acceptance:
#> [1] 0.8090615
#>
#> sample summary
#> V1
#> Min. :-3.40328
#> 1st Qu.:-0.66900
#> Median : 0.01794
#> Mean : 0.01114
#> 3rd Qu.: 0.67393
#> Max. : 3.10759
# two-dimensional normal with positive association ----------------
rho <- 0.9
covmat <- matrix(c(1, rho, rho, 1), 2, 2)
log_dmvnorm <- function(x, mean = rep(0, d), sigma = diag(d)) {
x <- matrix(x, ncol = length(x))
d <- ncol(x)
- 0.5 * (x - mean) %*% solve(sigma) %*% t(x - mean)
}
x <- ru(logf = log_dmvnorm, sigma = covmat, d = 2, n = 1000, init = c(0, 0))
summary(x)
#> ru bounding box:
#> box vals1 vals2 conv
#> a 1.0000000 0.000000e+00 0.000000e+00 0
#> b1minus -0.8008868 -1.320440e+00 -7.414894e-09 0
#> b2minus -0.8008868 -9.922367e-10 -1.320440e+00 0
#> b1plus 0.8008868 1.320440e+00 -2.548085e-09 0
#> b2plus 0.8008868 -9.851831e-10 1.320440e+00 0
#>
#> estimated probability of acceptance:
#> [1] 0.5307856
#>
#> sample summary
#> V1 V2
#> Min. :-3.72460 Min. :-3.77076
#> 1st Qu.:-0.68457 1st Qu.:-0.64997
#> Median :-0.05360 Median :-0.06968
#> Mean :-0.02358 Mean :-0.01355
#> 3rd Qu.: 0.64936 3rd Qu.: 0.64148
#> Max. : 3.67461 Max. : 4.24154