Prueba binomial en R

Question

Necesito hacer una prueba Bernoulli en R.

Digamos que tengo un conjunto de datos especificados para insertar (10 muestras aleatorias de 40). De esto hay cuatro opciones posibles, y deseo asignar la probabilidad de que se elija en base a sus medias y sd las muestras generadas aleatoriamente provienen de una distribución normal con media de 25 y sd de 4.

No estoy seguro de cómo codificar esto en R. Gracias por cualquier ayuda.

Answer 1

Si quiere comprobar que su muestra procede de una distribución Bernoulli con parámetro $p=0.3$ hay que contar el número de aciertos (los 1) y el tamaño de la muestra. He aquí un ejemplo

# Simulate a random sample of size 17 with p=.3
set.seed(123)
simulated_sample <- runif(17) < 0.3
x <- sum(simulated_sample == TRUE)    # (number of successes)
n <- length(simulated_sample) # (sample size)
# Do the test.
binom.test(x, n, p=0.3)
#   Exact binomial test
# 
# data:  x and n 
# number of successes = 4, number of trials = 17, p-value = 0.7919
# alternative hypothesis: true probability of success is not equal to 0.3 
# 95 percent confidence interval:
#  0.06810774 0.49899327 
# sample estimates:
# probability of success 
#              0.2352941