Probabilidad empírica en R con x1+x2>2*x3
Tengo una pregunta: ¿cómo generar 100 veces un vector de 3 números aleatorios (x1, x2, x3) del 1 al 10 para poder encontrar la probabilidad empírica de que x1+x2>2*x3?
Sé cómo generar 3 números aleatorios del 1 al 10 con sample(1:10, 3), ¿pero cómo hacer esto n veces?
Gracias de antemano.
davy.ai
Para generar 100 vectores de 3 números aleatorios del 1 al 10, puedes usar el siguiente código:
“`{r}
set.seed(123) # para reproducibilidad
n <- 100
vec <- matrix(sample(1:10, n*3, replace = TRUE), ncol = 3)
La función
apply
aplica la funciónsum(x[1:2]) > 2*x[3]
a cada fila devec
, lo cual devuelve un valor lógico que indica six1 + x2 > 2*x3
. La funciónmean
luego calcula la media de estos valores lógicos, lo cual representa la probabilidad empírica dex1 + x2 > 2*x3
sobre los 100 vectores.