[R] sample size 2-sample proportion test
Berta
ibanez at bioef.org
Tue Apr 17 14:10:47 CEST 2007
Hi R-users,
I want to calculate the sample size needed to carry out a 2-sample
proprotion test (with alfa=0.05, beta=0.8)
1.- probability of success in subpopulation A: 0.8
2.- probability of success in subpopulation B: 0.05
3.- percentage of population in subpopulation A = 5%,
4.- percentage of population in subpopulation B= 95% (the rest, A and B
complementary)
(Note that it is dificult to find people from subpopulation A)
Using library(Hmisc) (or others if exist?), I would like to know which one
of these three options (if any) should I use:
# FIRST:
bsamsize(p1=0.8, p2=0.05, fraction=0.05, alpha=.05, power=.8)
# n1 n2
# 1.46 27.82
# SECOND:
samplesize.bin(alpha=0.05, beta=0.8, pit=0.8, pic=0.05, rho=0.05)
# n= 42
# THIRD: To calculate sample size assuming fraction=0.5 and then multiplying
by 0.05:
bsamsize(p1=0.8, p2=0.05, fraction=0.5, alpha=.05, power=.8)
# n1 n2
# 5.47 5.47
# so that n=6*100/5=120
Taking into account the difficulty of obtaining the samples, which method
would be recomended? Or in other words, how many people should I find,
around 30-40 or 120?
Thanks a lot in advance,
Berta
More information about the R-help
mailing list