African Estuary Fishes
This data set is a list containing relations between sites and fish species linked to dates.
data(baran95)
This list contains the following objects:
is a data frame 95 seinings and 33 fish species.
is a data frame 2 factors : date and site. The date has 6 levels (april 1993, june 1993,
august 1993, october 1993, december 1993 and february 1994) and the sites are defined by 4 distances to the
Atlantic Ocean (km03, km17, km33 and km46).
is a vector of species latin names.
Baran, E. (1995) Dynamique spatio-temporelle des peuplements de Poissons estuariens en Guinée (Afrique de l'Ouest). Thèse de Doctorat, Université de Bretagne Occidentale. Data collected by net fishing sampling in the Fatala river estuary.
See a data description at http://pbil.univ-lyon1.fr/R/pdf/pps027.pdf (in French).
data(baran95)
w <- dudi.pca(log(baran95$fau + 1), scal = FALSE, scann = FALSE,
nf = 3)
w1 <- wca(w, baran95$plan$date, scann = FALSE)
fatala <- ktab.within(w1)
stat1 <- statis(fatala, scan = FALSE, nf = 3)
mfa1 <- mfa(fatala, scan = FALSE, nf = 3)
if(adegraphicsLoaded()) {
g1 <- s.class(stat1$C.Co, baran95$plan$site, facets = baran95$plan$date,
pellipses.axes.draw = FALSE, ppoints.cex = 0.5, plot = FALSE)
n1 <- length(g1@ADEglist)
g2 <- ADEgS(lapply(1:n1, function(i) s.label(stat1$C.Co, plabels.cex = 0,
ppoints.cex = 0.5, plot = FALSE)), positions = g1@positions, plot = FALSE)
G1 <- superpose(g2, g1, plot = TRUE)
G2 <- kplot(stat1, arrow = FALSE, traject = FALSE, class = baran95$plan$site,
col.plabels.cex = 0, ppoints.cex = 0.5)
g3 <- s.class(mfa1$co, baran95$plan$site, facets = baran95$plan$date,
pellipses.axes.draw = FALSE, ppoints.cex = 0.5, plot = FALSE)
n2 <- length(g3@ADEglist)
g4 <- ADEgS(lapply(1:n2, function(i) s.label(mfa1$co, plabels.cex = 0,
ppoints.cex = 0.5, plot = FALSE)), positions = g3@positions, plot = FALSE)
G3 <- superpose(g4, g3, plot = TRUE)
} else {
par(mfrow = c(3, 2))
w2 <- split(stat1$C.Co, baran95$plan$date)
w3 <- split(baran95$plan$site, baran95$plan$date)
for (j in 1:6) {
s.label(stat1$C.Co[,1:2], clab = 0, sub = tab.names(fatala)[j], csub = 3)
s.class(w2[[j]][, 1:2], w3[[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
par(mfrow = c(3, 2))
w4 <- split(mfa1$co, baran95$plan$date)
for (j in 1:6) {
s.label(mfa1$co[, 1:2], clab = 0, sub = tab.names(fatala)[j], csub = 3)
s.class(w4[[j]][, 1:2], w3[[j]], clab = 2, axese = FALSE, add.plot = TRUE)
}
par(mfrow = c(1, 1))
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.