createWeights: Internal weights function

In jon-mellon/weightBES: Raking weights for BES Internet panel

Usage

createWeights(current.targets, bes)

Arguments

current.targets
bes

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (current.targets, bes) 
{
    not.there <- names(current.targets)[!names(current.targets) %in% 
        names(bes)]
    if (length(not.there) > 0) {
        stop(paste("Var not present in BES", not.there, collapse = ""))
    }
    vars <- bes[complete.cases(bes[, names(current.targets)]), 
        ]
    if (any(!sapply(vars[, names(current.targets)], is.factor))) {
        stop("Following vars are not factor:", names(which(!sapply(vars[, 
            names(current.targets)], is.factor))))
    }
    testVar <- function(name) {
        if (all(sort(levels(vars[[name]])) == sort(names(current.targets[[name]])))) {
            return(NULL)
        }
        else {
            stop(paste0("Failed on ", name))
        }
    }
    sapply(names(current.targets), testVar)
    library(anesrake)
    output <- anesrake(inputter = current.targets, dataframe = vars, 
        caseid = vars$id, verbose = TRUE, cap = 15, choosemethod = "total", 
        type = "nolim", pctlim = 0.005, nlim = 5, iterate = TRUE, 
        force1 = TRUE)
    vars$test_weight <- output$weightvec
    vars <- vars[, c("id", "test_weight", "region.weighted")]
    return(vars)
  }

jon-mellon/weightBES documentation built on May 19, 2019, 7:26 p.m.