cond_power <- function(t1,n_1,n_2,n_2neu,alpha_loc){ w_1 <- sqrt(n_1) w_2 <- sqrt(n_2) cp <- 1-pnorm(qnorm(1-alpha_loc)*sqrt(w_1^2+w_2^2)/w_2-t1*(w_1/w_2+sqrt(n_2neu)/sqrt(n_1))) }