Блог

надеюсь, у вас все хорошо. Итак, меня попросили разработать и провести имитационное исследование, чтобы оценить вероятность того, что наблюдаемая максимальная длина пробега (пробег — это последовательность последовательных орлов или решек) для эксперимента по подбрасыванию честной монеты находится в интервале [9, 11] при размере выборки n = 1000, это мойпопытка пока

 result <-replicate(10000,{  #replicate 10000 times
experiment <- sample(c("T","H"),size=1000,replace=TRUE) #1000 flips
expe_run <- rle(experiment)     #find the run
expe_val <- expe_run$values      #values of run
expe_length <- expe_run$lengths       #length of run
as<-list(expe_length,expe_val)     #make a list for sapply function
max_run <-sapply(as, FUN=max)       #apply max function through out for both
head_run <-expe_length[which(expe_val=='H')]     # show the head run
tail_run <-expe_length[which(expe_val=='T')]      #show the tail run
max_run
})

probability <-table(result)/10000  #probability for run
probability
  

Проблема в том, что я не знаю, как закончить вопрос, который оценивает вероятность того, что наблюдаемая максимальная длина пробега для эксперимента с подбрасыванием честной монеты находится в интервале [9, 11], хотя я получил таблицу для каждой возможной вероятности. Не могли бы вы мне помочь? Спасибо

Вы ищете вероятность того, что максимальное значение in в течение 10 000 испытаний expe_length содержится в интервале [9,11] ?

Если да, то что-то вроде:

 result <-replicate(10000, { 
  ...

  expe_length <- expe_run$lengths #length of run

  max(expe_length) %in% 9:11
})
  

должно предоставить вам result в виде вектора TRUE / FALSE значений, указывающих, имела ли эта пробная версия максимальную длину прогона в этом интервале.

После этого,

 sum(result) / length(result)
  

даст вам пропорцию, которую вы ищете.