Среднее число выборок для плана двойной выборки с сокращением

Question

Я пытаюсь вычислить среднее число выборок для плана двойной выборки с сокращением. Я очень плохо знаком с циклом внутри R, поэтому я ищу некоторую помощь. Я пришел из инженерного образования без особого обучения программированию, поэтому мне это сложно. Сюжет, который я пытаюсь создать, сравнивает ASN с сокращением и без него.

Я пытался поместить цикл for в разные места без особой удачи.

asnWOcurt <- n4_1 + n4_2*(pbinom(c4_2,n4_1,p_vec) - 
pbinom(c4_1,n4_1,p_vec))
summ_m <- 0
summ_i <- 0

for (i in (c4_1+1):c4_2) {
  for (m in (c4_2-i+1):n4_2) {
    summ_m <- summ_m + m*pnbinom(c4_2-i+1,m,p_vec)
  }
  summ_i <- summ_i + (n4_2*dbinom(c4_2-i, n4_2, p_vec) + 
summ_m)*pbinom(i,n4_1,p_vec)
}
asnWcurt <- n4_1 + summ_i 

 plot(p_vec, asnWOcurt, type="o", col="blue", pch="o", lty=1, 
xlim=c(0,0.2), ylim=c(80,200) )
  points(p_vec, asnWcurt, col="red", pch="*")
  lines(p_vec, asnWcurt, col="red",lty=2)