R/t1waybtsqrk.R

In musto101/wilcox_R: What the package does (short line)

t1waybtsqrk <-
function(x,alpha=.05,nboot=599,B=nboot,tr=0.2){
#
#  One-way ANOVA for trimmed means, independent groups.
#  Uses a method studied by Ozdemir et al. 
#  
if(is.matrix(x))x=listm(x)
x=lapply(x,elimna)
T<-btsqrk(x,alpha=alpha,tr=tr)$teststat
means<-c()
ylist<-list(0)
TT<-c()
b<-floor((1-alpha)*B)
means<-sapply(x,mean,tr)
k<-length(x)
for (i in 1:B)
{
	for (j in 1:k)
	{ylist[[j]]<-(sample(x[[j]],length(x[[j]]),replace=T)-means[j])}
	TT<-c(TT,btsqrk(ylist,alpha,tr)$teststat)
}
TT=sort(TT)
pval<-mean(T<=TT,na.rm=TRUE)
list(test.stat=T,crit.value=TT[b],p.value=pval)
}