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$‍ quilos e um desvio-padrão de $1$‍ quilo. Suspeitando que essa lontra não estava comendo o suficiente, ele tomou uma amostra aleatória de $n=10$‍ dias e observou uma média amostral de $\overline{x}=16,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$‍ amostras de $n=10$‍ valores de uma população normal com uma média de $17$‍ quilos e um desvio-padrão de $1$‍ quilo. Ele registrou a média dos valores em cada amostra. Veja as médias amostrais das $40$‍ amostras:

O funcionário quer testar $H_0:\mu =17\text{ kg}$‍ versus $H_{\text{a}}:\mu <17\text{ kg}$‍, sendo $\mu$‍ a média real de comida por dia.