Exemples de code R pour les exercices

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]