bd1way1 <-
function(isub,xcen,est,misran,...){
#
# 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<-vector("numeric")
for (j in 1:ncol(xcen))val[j]<-est(xcen[isub,j],na.rm=misran,...)
bd1way1<-(length(val)-1)*var(val)
bd1way1
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.