Nothing
test_that("metrics calculations work in 1 patch/strata", {
f <- 0.3
q <- 0.9
g <- 1/10
tau <- 11
p <- 1
Omega <- diag(g, p)
M <- 500
W <- 800
vc <- (f*q/g) * exp(-g*tau) * (f*q*M)/W
calV <- metric_calV(f = f, q = q, Omega = Omega, tau = tau, M = M, W = W)
expect_true(approx_equal(as.vector(calV), vc))
Psi <- matrix(1,1,1)
beta <- t(Psi) %*% diag_inverse(W)
b <- 0.55
c <- 0.15
r <- 1/200
D <- c/r
calD <- metric_calD(W = W, beta = beta, b = b, D = D, H = W)
expect_true(approx_equal(b*c/r, as.vector(calD)))
calR <- metric_calR(b = b, beta = beta, calV = calV, W = W, D = D, H = W)
r0 <- (f*q/g) * exp(-g*tau) * (f*q*M)/W * (b*c/r)
expect_true(approx_equal(r0, as.vector(calR)))
calZ <- metric_calZ(Omega = Omega, tau = tau, f = f, q = q, M = M, W = W, calD = calD)
expect_true(approx_equal(r0, as.vector(calZ)))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.