Как получить другой результат с симуляцией повторения? - PullRequest
Новая
       7

Как получить другой результат с симуляцией повторения?

0 голосов
/ 15 октября 2018

Я пытался проверить стандартную ошибку случайно выбранной выборки из 24 выбранных из 64 объектов, но когда я пытаюсь запустить ее 10 раз, она дает 10 одинаковых ответов.Как правильно проверить это с повторением?Это мой код и результаты. 

po=c(900,822,781,805,670,1238,573,634,578,487,442,451,459,464,400,366,364,317,328,302,288,291,253,291,308,272,284,255,270,214,195,260,209,183,163,253,232,260,201,147,292,164,143,169,139,170,150,143,113,115,123,154,140,119,130,127,100,107,114,111,163,116,122,134)

se_1 = sd(sample(po[1:64],24,replace = FALSE)) * sqrt(1/24)

rep(se_1, times=10)
#[1] 44.69907 44.69907 44.69907 44.69907 44.69907 44.69907 44.69907 44.69907
#[9] 44.69907 44.69907

1 Ответ

0 голосов
/ 15 октября 2018

Сделать se_1 функцией 

se_1 <- function() sd(sample(po[1:64],24,replace = FALSE)) * sqrt(1/24)

Тогда вы можете так 

replicate(se_1(), times=10)

, или вы можете просто сделать 

replicate(sd(sample(po[1:64],24,replace = FALSE)) * sqrt(1/24), times=10)

replicate().Выражение, которое вы хотите повторно вызывать, в то время как rep() просто повторяет существующее значение (или значения) определенное количество раз.Он не вычисляет «новые» значения.

...