bbbtrim <-
function(J,K,L,data,tr=.2,grp=c(1:p),alpha=.05,p=J*K*L,nboot=600,pr=FALSE){
#
# Perform three-way anova, independent groups, based on trimmed means
#
# That is, there are three factors with a total of JKL independent groups.
#
# A bootstrap-t method is used to perform multiple comparisons
# The variable data is assumed to contain the raw
# data stored in list mode. data[[1]] contains the data
# for the first level of all three factors: level 1,1,1.
# data[[2]] is assumed to contain the data for level 1 of the
# first two factors and level 2 of the third factor: level 1,1,2
# data[[L]] is the data for level 1,1,L
# data[[L+1]] is the data for level 1,2,1. data[[2L]] is level 1,2,L.
# data[[KL+1]] is level 2,1,1, etc.
#
# The default amount of trimming is tr=.2
#
# It is assumed that data has length JKL, the total number of
# groups being tested.
#
if(is.list(data))data=listm(elimna(matl(data)))
if(is.matrix(data))data=listm(elimna(data))
if(!is.list(data))stop("Data are not stored in list mode or a matrix")
if(p!=length(data)){
print("The total number of groups, based on the specified levels, is")
print(p)
print("The number of groups in data is")
print(length(data))
print("Warning: These two values are not equal")
}
x=data
temp=con3way(J,K,L)
conA<-temp$conA
conB<-temp$conB
conC<-temp$conC
conAB<-temp$conAB
conAC<-temp$conAC
conBC<-temp$conBC
conABC=temp$conABC
Factor.A<-linconb(x,con=conA,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
Factor.B<-linconb(x,con=conB,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
Factor.C<-linconb(x,con=conC,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
Factor.AB<-linconb(x,con=conAB,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
Factor.AC<-linconb(x,con=conAC,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
Factor.BC<-linconb(x,con=conBC,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
Factor.ABC<-linconb(x,con=conABC,tr=tr,alpha=alpha,nboot=nboot,pr=pr)
list(Factor.A=Factor.A,Factor.B=Factor.B,Factor.C=Factor.C,
Factor.AB=Factor.AB,Factor.AC=Factor.AC,Factor.BC=Factor.BC,
Factor.ABC=Factor.ABC,pr=pr)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.