Я работаю над вопросом в R для курса статистики и наткнулся на этот вопрос:
Используя функцию
power.anova.test в пакете basi c R stats, обнаружите значительную разницу между средствами обработки на основе случайных выборок.
сбалансированный план с 20 наблюдениями в каждой группе: общее среднее значение (8 + 8,5 + 9) / 3 = 8,5, между группами дисперсия = 5, внутригрупповая дисперсия - 7,29.
Определите размер образца, который потребуется для достижения степени 0,80.
Повторите шаг 9 для ряда лечебных средств. Сохраняйте среднее значение контрольной группы таким же (8), но используйте ряд возможных значений для средних двух групп лечения, скажем, между 8,1 и 9,5 с шагом 0,1. Каждый раз рассчитывайте размер выборки, необходимый в каждой группе, чтобы достичь мощности 80%. Изучите, как изменяется размер выборки по мере того, как лечение (t1 и t2) означает изменение. Подсказка: рассмотрите возможность реализации «двойного цикла» для выполнения вычислений для каждой комбинации «Treat1» и «Treat2» в диапазоне от 8,1 до 9,5 с шагом 0,1. Например:
n.min <- matrix(NA, 15, 15) # store minimum sample sizes for
different
# treatment mean combinations
mean.c <- 8
mean.t1 <- seq(8.1, 9.5, 0.1)
mean.t2 <- seq(8.1, 9.5, 0.1)
for (i in 1:15) {
for (j in 1:15){
mean.g <- ...
var.b <- ...
n.min[i,j] <- power.anova.test(groups=3,
between.var=var.b,
within.var=var.w, sig.level=0.05, power=0.80)$n
}
}
Это последний фрагмент кода, который я не совсем понимаю. Я знаю последовательную функцию и использовал ее в прошлом, но не понимаю следующего:
for (i in 1:15){
for (j in 1:15){
mean.g<-...
var.b<-...
n.min[i.j] <-power.anova.test(groups=3,
between.var=var.b,
within.var=var.w, sig.level=0.05, power=0.80)$n
Если кто-нибудь может мне помочь, я был бы очень признателен! Даже если это другой способ зацикливания power.anova.test, чтобы ответить на поставленный выше вопрос. Спасибо!