bd1GLOB1 <-
function(isub,xcen,est,...){
#
# Compute test statistic for bd1way
#
# isub is a vector of length n,
# a bootstrap sample from the sequence of integers
# 1, 2, 3, ..., n
#
# xcen is an n by J matrix containing the input data
#
val<-NA
val=est(xcen[isub,],...)$center
val=as.vector(val)
test.stat<-(length(val)-1)*var(val)
test.stat
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.