https://github.com/cran/aster
Tip revision: aa47935123bfca8a22cbc8345d658d0c1713a289 authored by Charles J. Geyer on 14 December 2023, 15:20:02 UTC
version 1.1-3
version 1.1-3
Tip revision: aa47935
rktp.R
rktp <- function(n, k, mu, xpred = 1) {
if (length(n) > 1)
n <- length(n)
stopifnot(n == as.integer(n))
stopifnot(n > 0)
stopifnot(is.numeric(xpred))
stopifnot(length(xpred) > 0)
stopifnot(all(xpred == as.integer(xpred)))
stopifnot(is.numeric(mu))
stopifnot(length(mu) > 0)
stopifnot(all(mu > 0))
stopifnot(is.numeric(k))
stopifnot(length(k) >= 0)
stopifnot(all(k == as.integer(k)))
.C(C_aster_rktp,
n = as.integer(n),
lenxp = length(xpred),
lenmu = length(mu),
lenk = length(k),
xpred = as.double(xpred),
mu = as.double(mu),
k = as.integer(k),
result = double(n))$result
}
rnzp <- function(n, mu, xpred = 1) rktp(n, 0, mu, xpred)