bsqrmbt <-
function(x,y,alpha=.05,bend=1.28,nboot=599,SEED=TRUE){
#
# Goal: Test hypothesis that two independent groups have
# equal population M-measures of location.
# A bootstrap-t method is used.
# The method used was derived by F. Ozdemir
#
if(SEED)set.seed(2) # set seed of random number generator so that
# results can be duplicated.
T<-bsqrm(x,y,alpha,bend)$teststat
TT<-0
bsqrmbt<-numeric(2)
xone<-onestep(x,bend=bend)
yone<-onestep(y,bend=bend)
for(j in 1:nboot)
{
xx<-(sample(x,length(x),replace=TRUE)-xone)
yy<-(sample(y,length(y),replace=TRUE)-yone)
TT[j]<-bsqrm(xx,yy,alpha,bend)$teststat
}
TT<-sort(TT)
bott<-round(alpha*nboot)+1
bsqrmbt<-TT[nboot-bott]
pv=mean(T<=TT)
list(critval=bsqrmbt,teststat=TRUE,p.value=pv)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.