cbmhd_subMC <-
function(data,cbmhd_subMC,x,y,q,q2,n1,n2){
np1=n1+1
nall=n1+n2
mb=outer(x[data[1:n1]],y[data[np1:nall]],"-")
est=hd(mb,q)+hd(mb,q2)
est
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.