FAQ/unpinR - CBU statistics Wiki

Please enter your password of your account at the remote wiki below.
/!\ You should trust both wikis because the password could be read by the particular administrators.

Clear message
location: FAQ / unpinR

Equal samples t test for a given power in R

For a two-tailed test with type I error, alp, with inputted mean and average group sd will have power equal to pow.

[TYPE INTO R THE DESIRED INPUTS ALPHA, MEAN, SD AND N USING VALUES IN FORM BELOW].

The R syntax below

pow <- 0.8
beta <- 0.05
mdiff <- 2
sd <- 4

the the below will give the equal samples in nout.

fn <- function(n) {
(pow - (1 - pt(abs(qt(beta/2,(2*n)-2)),(2*n)-2,mdiff/sqrt(2*sd*sd/n))) )^2
 }

nout<- nlm(fn,2)
nout<- trunc(nout$estimate+1)
print(nout)

This can be run in SPSS 16 using R integrator plug-in (see here).