Buen dia.
Quiero generar una función que permita modificar la estructura del vector "ADN" de tal forma que sustituya las bases con una probabilidad del 10%. Es decir que si en la primera entrada tenemos una "C" la función decida si es remplazada por alguna de las otras 3 bases incluida dentro del vector "bases"
bases <- c("A","G","C","T")
P <- runif(1,0,1) # Valor de probabilidad de que ocurra un base, único para
cada una.
ADN <- sample(bases,size = 21,replace = TRUE,prob = c(P,P,P,P))
ADN
[1] "C" "G" "T" "T" "G" "G" "G" "C" "C" "T" "C" "A" "A" "T" "G" "G" "A" "G" "T" "T" "T"
¿Pueden ayudarme con esta duda por favor?