Nothing
Pkl.Hajek.U <- function(VecPk.U)
{
if(! is.vector(VecPk.U) ){stop("VecPk.U must be a vector.") }
if(anyNA(VecPk.U) ){stop("There are missing values in VecPk.U.") }
if(min(VecPk.U)<=0|max(VecPk.U)>1){stop("There are invalid values in VecPk.U.") }
N <- length(VecPk.U)
if(all(VecPk.U==1) ){OUTPUT <- 1 }
else
{
OUTPUT <- .C("Pkl_Hajek_U",
as.double(VecPk.U),
N,
VectMatProb = double(N*N),
PACKAGE = "samplingVarEst")$VectMatProb
}
if(min(OUTPUT) <= 0 ){warning("Some values of 2nd order (joint) inclusion probabilities are so tiny that they are assumed zero by R.")}
matrix(OUTPUT, nrow=N, ncol=N)
}
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.