getL1Strata <-
function(treatment, dataset, match.on, breaks){
# Remove dropped covs
dataset <- dataset[,match.on]
if(length(match.on) == 1){
dataset <- as.matrix(dataset)
colnames(dataset) <- match.on
}
## stuff borrowed from cem.main to add user defined breaks
vnames <- colnames(dataset)
nv <- dim(dataset)[2]
mycut <- vector(nv, mode="list")
names(mycut) <- vnames
for (i in 1:nv) {
tmp <- reduceVar(dataset[[i]], breaks[[vnames[i]]])
dataset[[i]] <- tmp$x
mycut[[vnames[i]]] <- tmp$breaks
}
# Calculate strata
strata <- stratify(dataset)
return(list(strata=strata, mycut=mycut))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.