Nothing
Teststat1 <-
function(trt, k=2, est, cens)
{
alpha = t(est$alpha)
m = dim(alpha)[2]
rho = 1-alpha
rho[cens!=3,] = 0
index = est$ppairs
f = est$pf
denom = alpha %*% f
dj = 1:m
for(j in 1:m)
dj[j] = sum(alpha[which(cens!=3),j]*f[j]/denom[which(cens!=3)])
djl = matrix(0, nrow=k, ncol=m)
for(i in 1:k)
for(j in 1:m)
djl[i,j] = sum(alpha[which(cens!=3 & trt==i-1),j]*f[j]/denom[which(cens!=3 & trt==i-1)])
nj = c(sum(dj), sum(dj)-cumsum(dj)[-m]) + apply(rho, 2, sum)
njl = matrix(0, nrow=k, ncol=m)
for(i in 1:k)
njl[i,] = c(sum(djl[i,]), sum(djl[i,])-cumsum(djl[i,])[-m]) + apply(rho[which(trt==i-1),], 2, sum)
u = 1:k
for(i in 1:k)
u[i] = sum(ifelse(nj > 0, djl[i,] - njl[i,] * dj / nj, 0))
u
}
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.