Nothing
VE.Hajek.Total.NHT <- function(VecY.s, VecPk.s)
{
if(! is.vector(VecY.s) ){stop("VecY.s must be a vector.") }
if(! is.vector(VecPk.s) ){stop("VecPk.s must be a vector.") }
n <- length(VecY.s)
if(n != length(VecPk.s) ){stop("The lengths of VecY.s and VecPk.s are different.") }
if(anyNA(VecPk.s) ){stop("There are missing values in VecPk.s.") }
if(min(VecPk.s)<=0|max(VecPk.s)>1 ){stop("There are invalid values in VecPk.s.") }
if(anyNA(VecY.s) ){stop("There are missing values in VecY.s.") }
Doublen <- as.double(n)
OUTPUT <- .C("VE_Hajek_Total_NHT",
as.double(VecY.s),
as.double(VecPk.s),
n,
Doublen,
VarEst = double(1),
PACKAGE = "samplingVarEst")$VarEst
if(OUTPUT<0 ){warning("The variance estimate contains negative values.")}
OUTPUT
}
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.