data(meaudret)
wit1 <- withinpca(meaudret$env, meaudret$design$season, scan = FALSE,
scal = "partial")
kta1 <- ktab.within(wit1, colnames = rep(c("S1","S2","S3","S4","S5"), 4))
kta2 <- t(kta1)
pta1 <- pta(kta2, scann = FALSE)
pta1
plot(pta1)
Run the code above in your browser using DataLab