data(baran95)
w <- dudi.pca(log(baran95$fau+1), scal = FALSE, scann = FALSE,
nf = 3)
w1 <- within(w, baran95$plan$date, scann = FALSE)
fatala <- ktab.within(w1)
stat1 <- statis(fatala, scan = FALSE, nf = 3)
w1 <- split(stat1$C.Co, baran95$plan$date)
w2 <- split(baran95$plan$site, baran95$plan$date)
par(mfrow = c(3,2))
for (j in 1:6) {
s.label(stat1$C.Co[,1:2], clab = 0,
sub = tab.names(fatala)[j], csub = 3)
s.class(w1[[j]][,1:2], w2[[j]], clab = 2, axese = FALSE,
add.plot = TRUE)
}
par(mfrow = c(1,1))
kplot(stat1, arrow = FALSE, traj = FALSE, clab = 2, uni = TRUE,
class = baran95$plan$site) #simpler
mfa1 <- mfa(fatala, scan = FALSE, nf = 3)
w1 <- split(mfa1$co, baran95$plan$date)
w2 <- split(baran95$plan$site, baran95$plan$date)
par(mfrow = c(3,2))
for (j in 1:6) {
s.label(mfa1$co[,1:2], clab = 0,
sub = tab.names(fatala)[j], csub = 3)
s.class(w1[[j]][,1:2], w2[[j]], clab = 2, axese=FALSE,
add.plot = TRUE)
}
par(mfrow = c(1,1))
Run the code above in your browser using DataLab