varcov1 <-
function (eitbg, center, scores, v12inv, Dmat, wts)
{
numc <- max(center)
n <- length(eitbg)
u <- (1:n)/(n + 1)
sc <- getScores(scores, u)
rs <- sc[rank(eitbg, ties.method = "first")]
brrp <- matrix(rep(0, n * n), ncol = n)
veeinv <- v12inv %*% diag(wts) %*% v12inv
part1 <- t(Dmat) %*% veeinv %*% Dmat
i1 <- 1
i2 <- 1
for (i in 1:numc) {
rsc <- rs[center == i]
ni <- length(rsc)
rscbar <- mean(rsc)
rsc <- rsc - rscbar
rrp <- rsc %*% t(rsc)
i2 <- i1 + ni - 1
brrp[i1:i2, i1:i2] <- rrp
i1 <- i1 + ni
}
mid <- t(Dmat) %*% v12inv %*% brrp %*% v12inv %*% Dmat
varcov1 <- solve(part1) %*% mid %*% solve(part1)
return(varcov1)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.