Nothing
dvolkov <- function(x, theta, m, J, order=96, log=FALSE){
if(length(theta) > 1 | length(m) > 1 | length(J) > 1) stop ("Vectorization of parameters not implemented")
if(!is.finite(theta) | theta <= 0 | !is.finite(J) | J <= 0 | m <= 0 | m >= 1) return (rep(NaN, length(x)))
Cvolkov <- .C(C_volkov, res=as.double(rep(0,J)), theta0=as.double(theta), m0=as.double(m), J0=as.integer(J), N0=as.integer(order), total=as.double(0))
y <- c(0, Cvolkov$res) # index 0 should receive 0
if (any(!is.wholenumber(x))) warning("non integer values in x")
x[ !is.wholenumber(x) | x < 0 | x > J ] <- 0
y <- y[x+1]
if (any(is.nan(y))) warning ("NaNs produced")
if(log) return(log(y))
else return(y)
}
######### INTERNAL USE ONLY! #########
# Calculates the expected S from a Volkov distribution (used in rsad)
Svolkov <- function(theta, m, J, order=96) {
warning("The use of this function outside sads core code is deprecated and will be removed in a future version!")
if(!is.finite(theta) | theta <= 0 | !is.finite(J) | J <= 0 | m <= 0 | m >= 1) return (NaN)
Cvolkov <- .C(C_volkov, res=as.double(rep(0,J)), theta0=as.double(theta), m0=as.double(m), J0=as.integer(J), N0=as.integer(order), total=as.double(0))
Cvolkov$total
}
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.