Delta is the ratio of group 2 to group 1, alpha is (two-tailed) type I error, group means, mean1 and mean2, with standard deviations, sd1 and sd2, and n1 is the number in group 1. Power is outputted.
[COPY AND PASTE INTO R]
fn <- function (alpha,delta,mean1,mean2,sd1,sd2,n1) { mdiff <- mean1 - mean2 sd <- sqrt((n1*sd1*sd1+n1*delta*sd2*sd2)/(n1+delta*n1-2)) pow <- 1 - pt(abs(qt(beta/2, n1+(delta*n1)-2)), n1+(delta*n1)-2, mdiff/sqrt((sd*sd)/n1+(sd*sd)/(n1*delta))) cat("Power for (two-tailed) unequal sized groups unpaired t-test =") print(pow) }
[THEN ENTER INPUTS IN FORM BELOW]
delta <- 2 alpha <- 0.05 mean1 <- 2 mean2 <- 0 sd1 <- 4 sd2 <- 5 n1 <- 45
AND RUN
fn(alpha,delta,mean1,mean2,sd1,sd2,n1)