Nothing
VE.EB.HT.Ratio <- function(VecY.s, VecX.s, VecPk.s, MatPkl.s, VecAlpha.s = rep.int(1, length(VecPk.s)))
{
if(! is.vector(VecY.s) ){stop("VecY.s must be a vector.") }
if(! is.vector(VecX.s) ){stop("VecX.s must be a vector.") }
if(! is.vector(VecPk.s) ){stop("VecPk.s must be a vector.") }
if(! is.vector(VecAlpha.s) ){stop("VecAlpha.s must be a vector.") }
if(! is.matrix(MatPkl.s) ){stop("MatPkl.s must be a matrix.") }
DimMat <- dim(MatPkl.s)
DimMatR <- DimMat[1]
DimMatC <- DimMat[2]
if(DimMatR != DimMatC ){stop("MatPkl.s must be a square matrix. Number of rows and columns has to be equal.") }
n <- length(VecY.s)
if(n != length(VecPk.s) ){stop("The lengths of VecY.s and VecPk.s are different.") }
if(n != length(VecX.s) ){stop("The lengths of VecY.s and VecX.s are different.") }
if(n != length(VecAlpha.s) ){stop("The lengths of VecY.s and VecAlpha.s are different.") }
if(n != DimMatR ){stop("The lengths of VecY.s, VecPk.s and dimensions of MatPkl.s are different.") }
if(anyNA(VecPk.s) ){stop("There are missing values in VecPk.s.") }
if(anyNA(VecAlpha.s) ){stop("There are missing values in VecAlpha.s.") }
if(min(VecPk.s)<=0|max(VecPk.s)>1 ){stop("There are invalid values in VecPk.s.") }
if(min(VecAlpha.s)<0 ){stop("There are invalid values in VecAlpha.s.") }
if(any(VecAlpha.s==0) ){warning("There are zero values in VecAlpha.s. It is valid, but estimates get unstable.")}
if(anyNA(MatPkl.s) ){stop("There are missing values in MatPkl.s.") }
if(min(MatPkl.s)<=0|max(MatPkl.s)>1 ){stop("There are invalid values in MatPkl.s.") }
if(anyNA(VecY.s) ){stop("There are missing values in VecY.s.") }
if(anyNA(VecX.s) ){stop("There are missing values in VecX.s.") }
if(any(VecX.s==0) ){warning("There are zero values in the denominator variable VecX.s.") }
VecRhok.s <- (1/VecPk.s)^(1-VecAlpha.s)
VecEstTheta_k <- .C("Est_Ratio_Excluding_All_ElementsByRho",
as.double(VecY.s),
as.double(VecX.s),
as.double(VecPk.s),
VecRhok.s,
n,
VectVarEst = double(n),
PACKAGE = "samplingVarEst")$VectVarEst
EstTheta <- Est.Ratio(VecY.s, VecX.s, VecPk.s)
VecPseudo.s <- (EstTheta - VecEstTheta_k)/(VecPk.s*VecRhok.s)
OUTPUT <- .C("VE_HT_form",
VecPseudo.s,
as.double(VecPk.s),
as.double(c(MatPkl.s)),
n,
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.