Analyse de la sensibilité de la sortie d'un automate cellulaire
Analyse de la sensibilité d'un automate cellulaire # Fonction qui calcule le nombre de cases occupées (=1) en fonction d'une variable X # qui est un data.frame qui contient les probabilités de mourir à chaque pas de temps (p) et # le nombre de cases (/9) pour obtenir un décès d'étouffement. # La variable t donne les temps pour lesquels on retourne la taille de la population # Si t est une valeur, on retourne juste la valeur pour ce temps, si t est un vecteur, # on retourne toutes les valeurs pour les temps t ac <- function(X, t=t) { Nfinal <- NULL for (grandtotal in 1:nrow(X)) { penc <- X[grandtotal, 1] v <- X[grandtotal, 2]*9 m <- matrix(data = 0, ncol=12, nrow=12) for (i in 1:10) { col <- sample(1:10, 1)+1 row <- sample(1:10, 1)+1 m[row, col] <- 1 } N <- NULL for (tps in 1:max(t)) { ...