Это будет табулировать этот вектор с использованием установленных вами границ:
table(cut(pop, c(-Inf,337,343,Inf) ))
(-Inf,337] (337,343] (343, Inf]
87 645 268
Таким образом, доля значений (что также является вероятностью) составляет:
table(cut(pop, c(-Inf,337,343,Inf) ))[2]/length(pop)
(337,343]
0.645
Чтобы сделать это воспроизводимымвы бы использовали set.seed()
.
. Чтобы уточнить оценку, если ее задают как теоретический вопрос, вы можете либо смоделировать ее с помощью replicate
, либо использовать:
pnorm(343, 341.08, 3.07)-pnorm(337, 341.08, 3.07)
[1] 0.6422225
Первый метод требует только данных.Два других метода потребовали бы , зная , что данные получены из нормального распределения с указанными параметрами.