One of key questions in experimental design is “How many measurements should be included in the sample?”
In order to answer this, 4 parameters are required:
- type I error, value selected for α (a), the risk of rejecting a true hypothesis (risk of false positive), normally leading to the level of confidential interval, eg we normally confidential with 95% of confidential interval. Under a normal distribution, the upper critical value for this 1.96.
- type II error, value of β (b), the risk of accepting a false null hypothesis when a particular value of the alternative hypothesis is true (risk of false negative). In practice, we accept 10% of such risk. This translated as the upper critical value is 1.282.
- value of the population standard deviation (δ, sd). This value can be obtained from a pilot study or literatures.
- margin of error, m for example, at 5% (standard value of 0.05). this can be explained as, the true mean is μ, we expect to estimate the mean (ý) with 5% error, the deviation from the true or margin m is 0.05u. The expression is
ý – m < µ < ý + m
The standard deviation or variance can be calculated from a pilot study. For a binary trait, the variance is obtained as p(1-p), p is the probability of one of the event, for example, disease. For a continuous traits, such body weight or height, this is initial assessment. It should be standardised for simplicity. The simple formula for this two-side test is
N = [(zα/2 + zβ )/(m/sd)]^2
here is a simple R function
pow<-function (a,b,sd,m){
za<-qnorm(1-a/2.0)
zb<-qnorm(1-b)
pow<-((za+zb)/(m/sd))^2
return(pow)}
pow(0.05, 0.10, 10,3)
> pow(0.05, 0.1, 10,3)
[1] 116.7491
to be continued
This is a interesting and useful topic. to be followed
http://www.statsoft.com/textbook/power-analysis/#power_doe4
http://en.wikipedia.org/wiki/Sample_size
Test proportion
http://www.itl.nist.gov/div898/handbook/prc/section2/prc242.htm
test mean
http://www.itl.nist.gov/div898/handbook/prc/section2/prc222.htm
Wonderful blog! I truly love how it? s easy on my eyes as well as the data are well written. I am wondering how I can be notified whenever a new post has been made. I have subscribed to your rss feed which need to do the trick! Have a nice day!