Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- tstat <- matrix(0,1000,11)
- for (j in 1:nrow(tstat)) {
- x <- rpearsonIV(66, m=8.06, nu=8.7, location=6.0087, scale=9.752)
- for (i in 1:ncol(tstat)) {
- delta <- (i-1)/10
- tstat[j, i] <- t.test(x+delta, alternative="greater")$p.value
- }
- }
- tpower <- colMeans(tstat < 0.05)
Add Comment
Please, Sign In to add comment