# NOT RUN {
## Phoneme data
library(fdcov)
library(fds)
# Create data set
data(aa); data(ao); data(dcl);data(iy);data(sh)
dat=cbind(aa$y[,1:20],ao$y[,1:20],dcl$y[,1:20],iy$y[,1:20],sh$y[,1:20])
dat=t(dat)
grp=c(rep(1,20),rep(2,20),rep(3,20),rep(4,20),rep(5,20))
# Test the equality of the covariance operators
p=ksample.perm(dat,grp,iter=100,only.glob=FALSE)
# Plot partial p-values
perm.plot(p,5, lab=c('aa','ao','dcl','iy','sh'))
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab