outlier <- rbinom(150000,1,0.1)

outlier.mat <- matrix(outlier,ncol=30)

data <- rnorm(150000)

data.mat <- matrix(data,ncol=30)

data.mat <- data.mat*(1+4*outlier.mat)

avgs <- apply(data.mat,1,mean)

sds <- apply(data.mat,1,sd)

tstats <- sqrt(30)*avgs/sds

hist(tstats,main="t Statistics")

sum(abs(tstats)>qt(0.975,29))/5000