f2: Numeric or complex vector with splitplot levels
f3: Numeric or complex vector with splitsplitplot levels
block: Numeric or complex vector with blocks
resp: Numeric vector with responses
alpha.f: Level of significance of the F test (default is 0.05)
alpha.t: Significance level of the multiple comparison test (default is 0.05)
norm: Error normality test (default is Shapiro-Wilk)
homog: Homogeneity test of variances (default is Bartlett)
mcomp: Multiple comparison test (Tukey (default), LSD and Duncan)
Returns
Analysis of variance of fixed effects and multiple comparison test of Tukey, Scott-Knott, LSD or Duncan.
Examples
f1=rep(c("PD","PDE","C"), e =40);f1=factor(f1,unique(f1))f2=rep(c(300,400), e =20,3);f2=factor(f2,unique(f2))f3=rep(c("c1","c2","c3","c4"), e =5,6);f3=factor(f3,unique(f3))bloco=rep(paste("B",1:5),24); bloco=factor(bloco,unique(bloco))set.seed(10)resp=rnorm(120,50,5)PSUBFAT2DBC(f1,f2,f3,bloco,resp,alpha.f =0.5)# force triple interactionPSUBFAT2DBC(f1,f2,f3,bloco,resp,alpha.f =0.4)# force double interaction