summary.matchit.exact <- function(object, covariates = FALSE, ...) {
XX <- object$X
treat <- object$treat
qbins <- max(object$subclass,na.rm=TRUE)
if(!covariates){
q.table <- as.data.frame(matrix(0,qbins,3))
names(q.table) <- c("Treated","Control","Total")
for(i in 1:qbins){
qi <- object$subclass==i
q.table[i,] <- c(sum(treat[qi]==1, na.rm=T), sum(treat[qi]==0, na.rm=T),
length(treat[qi & !is.na(qi)]))
}
} else {
kk <- ncol(XX)
q.table <- as.data.frame(matrix(0,qbins,kk+3))
names(q.table) <- c("Treated","Control","Total",dimnames(XX)[[2]])
for(i in 1:qbins){
qi <- object$subclass==i & !is.na(object$subclass==i)
q.table[i,] <- c(sum(treat[qi]==1, na.rm=T), sum(treat[qi]==0, na.rm=T),
length(treat[qi & !is.na(qi)]),as.numeric(XX[qi,,drop=F][1,]))
}
}
ntab <- table(factor(!is.na(object$subclass),
levels=c("TRUE","FALSE")), treat)
nn <- rbind(table(treat),
ntab[c("TRUE","FALSE"),])
dimnames(nn) <- list(c("All","Matched","Discarded"),
c("Control","Treated"))
## output
res <- list(q.table = q.table, nn = nn, subclass = object$subclass,
treat = object$treat, call = object$call)
class(res) <- c("summary.matchit.exact", "summary.matchit")
return(res)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.