Histogramme pour le CAI :
par(xaxs = "i", yaxs = "i") hist(cai, xlim = c(0, 1), xlab="CAI")
Graphe des valeurs centrées et cumulées de G+C3 et CAI :
gc3c <- gc3 - mean(gc3) caic <- cai - mean(cai) plot(cumsum(gc3c), type = "p", ann = FALSE, col = "red") par(new = TRUE) plot(cumsum(caic), type = "p", ann = FALSE, col = "green")
Modification de EXP$KD pour le code génétique de M. genitalium :
KD2 <- EXP$KD KD2[57] <- KD2[59]