Skip to contents

Cubic Ecological Data

bf88.Rd

bf88 is a list of 6 data frames corresponding to 6 stages of vegetation.
Each data frame gives some bird species informations for 4 counties.

Usage

data(bf88)

Format

A list of six data frames with 79 rows (bird species) and 4 columns (counties).
The 6 arrays (S1 to S6) are the 6 stages of vegetation.
The attribut 'nomesp' of this list is a vector of species French names.

Source

Blondel, J. and Farre, H. (1988) The convergent trajectories of bird communities along ecological successions in european forests. Oecologia (Berlin), 75, 83–93.

Examples

data(bf88)
fou1 <- foucart(bf88, scann = FALSE, nf = 3)
fou1
#> Foucart's  COA
#> class: foucart coa dudi
#> $call: foucart(X = bf88, scannf = FALSE, nf = 3)
#> table  number: 6 
#> 
#> $nf: 3 axis-components saved
#> $rank: 3
#> eigen values: 0.5278 0.3591 0.3235
#> blo    vector       6      blocks
#>  vector length mode    content       
#>  $cw    4      numeric column weights
#>  $lw    79     numeric row weights   
#>  $eig   3      numeric eigen values  
#> 
#>  data.frame nrow ncol content             
#>  $tab       79   4    modified array      
#>  $li        79   3    row coordinates     
#>  $l1        79   3    row normed scores   
#>  $co        4    3    column coordinates  
#>  $c1        4    3    column normed scores
#> 
#>      **** Intrastructure ****
#> 
#>  data.frame nrow ncol content                     
#>  $Tli       474  3    row coordinates (each table)
#>  $Tco       24   3    col coordinates (each table)
#>  $TL        474  2    factors for Tli             
#>  $TC        24   2    factors for Tco             
#> 

if(adegraphicsLoaded()) {
  g1 <- scatter(fou1, plot = FALSE)
  g2 <- s.traject(fou1$Tco, fou1$TC[, 1], plines.lty = 1:length(levels(fou1$TC[, 1])), plot = FALSE)
  g3 <- s.traject(fou1$Tco, fou1$TC[, 2], plines.lty = 1:length(levels(fou1$TC[, 2])), plot = FALSE)
  g41 <- s.label(fou1$Tco, plot = FALSE)
  g42 <- s.label(fou1$co, plab.cex = 2, plot = FALSE)
  g4 <- superpose(g41, g42, plot = FALSE)
  G1 <- ADEgS(list(g1, g2, g3, g4), layout = c(2, 2))
  
  G2 <- kplot(fou1, row.plab.cex = 0, psub.cex = 2)

} else {
  par(mfrow = c(2,2))
  scatter(fou1)
  s.traject(fou1$Tco, fou1$TC[, 1])
  s.traject(fou1$Tco, fou1$TC[, 2])
  s.label(fou1$Tco)
  s.label(fou1$co, add.p = TRUE, clab = 2)
  par(mfrow = c(1, 1))
  kplot(fou1, clab.c = 2, clab.r = 0, csub = 3)
}
#> Error in s.label(dfxy = fou1$li, xax = 1, yax = 2, plot = FALSE, storeData = TRUE,     pos = -3, plabels = list(cex = 0.75), xlim = c(-1.58905122803571,     1.42885036980866), ylim = c(-1.41491882964118, 1.60298276820319    )): non convenient selection for dfxy (can not be converted to dataframe)