Convirtamolo en función
genera_datos <- function(sitios=15, visitas=4, probabilidad=0.5) {
datos <- matrix(NA, sitios, visitas) # matriz vacia donde vamos a poner los datos
for (i in 1:sitios){ # variable i
y <- rbinom(visitas, 1, probabilidad) # dados guardados en y
datos [i,] <- y # y pasa a la fila i de la tabla de datos
}
return(datos) # muestra los datos al finalizar el bucle
}
# llamemos la funcion
genera_datos() #con los valores por defecto
## [,1] [,2] [,3] [,4]
## [1,] 1 1 0 0
## [2,] 1 0 1 1
## [3,] 1 0 0 0
## [4,] 0 0 1 0
## [5,] 0 0 0 1
## [6,] 0 0 0 1
## [7,] 1 1 1 0
## [8,] 1 1 1 1
## [9,] 1 1 1 0
## [10,] 0 1 0 0
## [11,] 1 0 0 1
## [12,] 1 0 1 1
## [13,] 1 0 1 0
## [14,] 1 0 0 1
## [15,] 1 1 1 1
genera_datos(30, 6, 0.5) # con 30 sitios, 6 visitas repetidas y probabilidad 0.5
## [,1] [,2] [,3] [,4] [,5] [,6]
## [1,] 1 1 0 1 0 1
## [2,] 1 1 1 1 0 0
## [3,] 0 0 1 1 1 0
## [4,] 0 0 1 0 0 1
## [5,] 1 1 1 1 0 1
## [6,] 1 0 1 1 0 0
## [7,] 0 1 0 1 0 1
## [8,] 0 0 1 0 0 1
## [9,] 1 0 0 1 1 1
## [10,] 1 1 0 1 1 0
## [11,] 1 0 1 1 1 0
## [12,] 1 1 1 0 1 0
## [13,] 0 0 1 0 1 0
## [14,] 1 1 0 0 0 0
## [15,] 0 0 1 1 1 0
## [16,] 0 1 0 0 0 0
## [17,] 1 1 0 1 1 1
## [18,] 1 0 1 0 0 0
## [19,] 1 1 0 0 0 1
## [20,] 0 1 1 1 1 0
## [21,] 1 0 1 1 1 1
## [22,] 1 0 0 1 1 0
## [23,] 1 1 0 0 0 0
## [24,] 0 1 0 0 1 0
## [25,] 1 1 0 0 1 1
## [26,] 1 0 1 1 0 0
## [27,] 1 0 1 0 0 1
## [28,] 0 0 0 0 1 0
## [29,] 1 1 1 0 0 0
## [30,] 1 1 0 0 1 1