Расчет размера выборки с F-распределением в R-программе
Мне нужно рассчитать требуемый размер выборки населения в нецентральном F. Определил пилотные выборки и мощность. Я не уверен, что код правильный и хорошо представлен для этого.
- Дайте мне знать, если что-то не так или требуется какая-либо модификация.
- Мне нужно построить эффективную мощность для всех рассчитанных размеров выборки. Как я могу написать код, чтобы построить это?
Мне также нужно добавить абстракцию x=y к графику без определения x и y. Это возможно?
x <- c(1,3,6,7,15,10,12,8,40,5) y <- c(2,4,6,8,23,10,40,27,10,22) n0=10 alpha <- 0.01 beta <- 0.9 delta <- 1 mean(x) mean(y) sxx=sum((x-mean(x))^2) sxx sxy=sum((y-mean(y))*(x-mean(x))) sxy b1hat=sxy/sxx b0hat=((mean(y)-(b1hat*mean(x)))) b0hat yhat= b0hat + b1hat*x yhat sigma2=sum((y-yhat)^2)/(n0-2) sigma2 lambda1=n0*(delta^2)*((1+mean(x))^2)/sigma2 lambda1 betaeff1=(1-pf(qf((1-0.01),2,n0-2),2,n0-2,lambda1)) betaeff1 while(betaeff1 < beta){ n0=n0+1 lambda1=n0*(delta^2)*((1+mean(x))^2)/sigma2 betaeff1=(1-pf(qf((1-0.01),2,n0-2),2,n0-2,lambda1)) print(betaeff1) print(n0)}