Estimar valores-p a partir de simulações

Um funcionário em um aquário monitora quanto suas lontras do mar comem. A quantidade de comida que uma determinada lontra come é aproximadamente distribuída de forma normal com uma média de $17$17 quilos e um desvio-padrão de $1$1 quilo. Suspeitando que essa lontra não estava comendo o suficiente, ele tomou uma amostra aleatória de $n=10$n, equals, 10 dias e observou uma média amostral de $\bar x=16{,}5$x, with, \bar, on top, equals, 16, comma, 5 quilos de comida por dia.

Para verificar qual era a probabilidade de que uma amostra como essa acontecesse por acaso, o funcionário realizou uma simulação. Ele simulou $40$40 amostras de $n=10$n, equals, 10 valores de uma população normal com uma média de $17$17 quilos e um desvio-padrão de $1$1 quilo. Ele registrou a média dos valores em cada amostra. Veja as médias amostrais das $40$40 amostras:

O funcionário quer testar $H_0: \mu=17 \text{ kg}$H, start subscript, 0, end subscript, colon, mu, equals, 17, start text, space, k, g, end text versus $H_\text{a}: \mu<17 \text{ kg}$H, start subscript, start text, a, end text, end subscript, colon, mu, is less than, 17, start text, space, k, g, end text, sendo $\mu$mu a média real de comida por dia.