Доклиника, распределение не нормально. Выборка?, _помогите посчитать выборку_ Опции

[himik]

День добрый.

Планируем делать доклинику. Пилот не делали, но из литературы знаем, что распределение ненормально. Подскажите, пожалуйста, как рассчитать выборку для сравнения препаратов. Нулевых гипотез 2 вида: а) препарат А лучше препарата Б, б) препарат А не хуже препарата В.

Буду благодарен за помощь.

[himik]

А почему вы выбрали именно это распределение?

[p2004r]

А почему вы выбрали именно это распределение?

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

[himik]

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

Насколько я понимаю, при тех параметрах, при которых вы писали код на стр 1 обсуждения, Вейбулл вырождается в экспоненциальное, не?

[p2004r]

Насколько я понимаю, при тех параметрах, при которых вы писали код на стр 1 обсуждения, Вейбулл вырождается в экспоненциальное, не?

там написан код который фитит распределение в условия

Код

my.w.par <- function() {
    optim(c(14,87),
          function(x) {
              x1<-x[1];
              x2<-x[2];
              d <- rweibull(60000, shape=x1, scale = x2)
              (86.4 - mean(d))^2 +
                  15*(7.25 - sd(d))^2
          },  
          method = "BFGS")$par
}

plot(t(replicate(30, my.w.par())))

ну а рисунок как то так

Код

res <- sapply(c(5,8,10,20,30,60,90), function(i) quantile(replicate(100000, mean(rweibull(i, shape=14, scale = 88))), probs = c(0.25, 0.975)) )

plot(c(5, 8, 10,20,30,60,90), res[1,], ylim = c(70,100), type = "b")
lines(c(5, 8, 10,20,30,60,90), res[2,], type = "b")