Diff for "FAQ/powprogs" - CBU statistics Wiki
location: Diff for "FAQ/powprogs"
Differences between revisions 1 and 25 (spanning 24 versions)
Revision 1 as of 2012-06-26 14:38:36
Size: 2909
Editor: PeterWatson
Comment:
Revision 25 as of 2017-06-02 12:14:00
Size: 4911
Editor: PeterWatson
Comment:
Deletions are marked like this. Additions are marked like this.
Line 5: Line 5:
For a theoretical background and details of specialist software have a look at graduate seminar on power at the [:StatsCourse2006:Graduate Statistics Programme October-December 2006]. There is also [http://homepages.gold.ac.uk/aphome/cc16work.doc a worked example] using "Method 2" on a t-test. For a theoretical background and details of specialist software have a look at graduate seminar on power at the [[StatsCourse2006|Graduate Statistics Programme October-December 2006]]. There is also [[http://homepages.gold.ac.uk/aphome/cc16work.doc|a worked example]] using "Method 2" on a t-test. [[attachment:vif.pdf | F. Y. Hsieh, Philip W. Lavori, Harvey J. Cohen and John R. Feussner (2003) An Overview of Variance Inflation Factors for Sample-Size Calculation ]] ''Eval Health Prof'' '''26''' 239-257 mentions various formulae for power calculations.

Note in some cases one needs to inflate the total sample size required if there is a natural clustering in the data such as patients being assessed by different exercise therapists. For example if there are b patients assessed by each exercise therapist with an intra-therapist correlation (ICC) then the design effect equals 1 + [(b-1)ICC]. The total sample needs to be multiplied by this design effect to give sample size adjusted for the clustering effect. A further adjustment for variations in cluster sizes can be made (measured by the coefficient of variation) can be incorporated into the formula giving DE = 1 + (b(1+cv^2)-1)ICC. This extra adjustment for differing cluster sizes is not needed for small cvs e.g. cv < 0.23 (see [[attachment:cvDE.pdf | page 26 of this presentation]]).
Line 9: Line 11:
 * [:FAQ/power/onesamp:One sample t-test]  * [[FAQ/power/onesamp|One sample t-test]]
Line 11: Line 13:
 * [:FAQ/power/unpaired:Unpaired t-tests (equal group sizes)]  * [[FAQ/power/unpaired|Unpaired t-tests (equal group sizes)]]
Line 13: Line 15:
 * [:FAQ/power/unpairedUneq:Unpaired t-tests (unequal group sizes)]  * [[FAQ/power/unpairedUneq|Unpaired t-tests (unequal group sizes)]]
Line 15: Line 17:
 * [:FAQ/power/pairt:Paired t-tests]  * [[FAQ/power/pairt|Paired t-tests]]
Line 17: Line 19:
 * [:FAQ/power/owAnovaN:Regression including One-Way ANOVA and ANCOVA]  * [[FAQ/power/owAnovaN|Regression including One-Way ANOVA and ANCOVA]]
Line 19: Line 21:
 * [:FAQ/power/pPowN: Comparing a single proportion with a constant (sign test)]  * [[FAQ/power/pPowN| Comparing a single proportion with a constant (sign test)]]
Line 21: Line 23:
 * [:FAQ/power/propsn: Comparing two independent proportions]  * [[FAQ/power/propsn| Comparing two independent proportions]]
Line 23: Line 25:
 * [:FAQ/power/prop1sn:Comparing three or more independent proportions]  * [[FAQ/power/prop1sn|Comparing three or more independent proportions]]
Line 25: Line 27:
 * [:FAQ/power/mcn: Comparing two related proportions (McNemar test)]  * [[FAQ/power/mcn| Comparing two related proportions (McNemar test)]]
Line 27: Line 29:
 * [:FAQ/power/FisherrN: Comparing two independent correlations from two different samples]  * [[FAQ/power/FisherrN| Comparing two independent correlations from two different samples]]
Line 29: Line 31:
 * [:FAQ/power/rmPowN: A term in any Anova (including repeated measures)]  * [[FAQ/power/rmPowN| A term in any Anova (including repeated measures)]]
Line 31: Line 33:
 * [:FAQ/power/llogN: A single predictor in a multiple binary logistic regression]  * [[FAQ/power/llogN| A single predictor in a multiple binary logistic regression]]
Line 33: Line 35:
 * [:FAQ/power/hazN: Survival analysis]  * [[FAQ/power/hazN| Survival analysis]]
Line 38: Line 40:
 * [:FAQ/power/onesampn:One sample t-test]  * [[FAQ/power/onesampn|One sample t-test]]
Line 40: Line 42:
 * [:FAQ/power/unpaireqn:Unpaired t-tests (equal group sizes)]  * [[FAQ/power/unpaireqn|Unpaired t-tests (equal group sizes)]]
Line 42: Line 44:
 * [:FAQ/power/unpairn:Unpaired t-tests (unequal group sizes)]  * [[FAQ/power/unpairn|Unpaired t-tests (unequal group sizes)]]
Line 44: Line 46:
 * [:FAQ/power/pairn:Paired t-tests]  * [[FAQ/power/pairn|Paired t-tests]]
Line 46: Line 48:
 * [:FAQ/power/owanova:Regression including One-Way ANOVA and ANCOVA]  * [[FAQ/power/owanova|Regression including One-Way ANOVA and ANCOVA]]
Line 48: Line 50:
 * [:FAQ/power/pPow: Comparing a single proportion with a constant (sign test)]  * [[FAQ/power/pPow| Comparing a single proportion with a constant (sign test)]]
Line 50: Line 52:
 * [:FAQ/power/props:Comparing two independent proportions]  * [[FAQ/power/props|Comparing two independent proportions]]
Line 52: Line 54:
 * [:FAQ/power/prop1s:Comparing three or more independent proportions]  * [[FAQ/power/prop1s|Comparing three or more independent proportions]]
Line 54: Line 56:
 * [:FAQ/power/mcnemarN: Comparing two related proportions (McNemar test)]  * [[FAQ/power/mcnemarN| Comparing two related proportions (McNemar test)]]
Line 56: Line 58:
 * [:FAQ/power/Fisherr: Comparing two independent correlations from two different samples]  * [[FAQ/power/Fisherr| Comparing two independent correlations from two different samples]]
Line 58: Line 60:
 * [:FAQ/power/rmPow:A term in any Anova (including repeated measures)]  * [[FAQ/power/rmPow|A term in any Anova (including repeated measures)]]
Line 60: Line 62:
 * [:FAQ/power/llogPow: A single predictor in a multiple binary logistic regression]  * [[FAQ/power/llogPow| A single predictor in a multiple binary logistic regression]]
Line 62: Line 64:
 * [:FAQ/power/haz: Survival analysis]  * [[FAQ/power/haz| Survival analysis]]

 * [[FAQ/power/roc| ROC Analysis]]
Line 65: Line 69:
Additional power freeware (including the popular G*POWER (currently version 3)) is available for download from [http://www.epibiostat.ucsf.edu/biostat/sampsize.html#PCSize here.] Some examples using G*POWER 3 are in Howell (2013). Additional power freeware (including the popular G*POWER (currently version 3)) is available for download from [[http://www.epibiostat.ucsf.edu/biostat/sampsize.html#PCSize|here.]] Some examples using G*POWER 3 are in Howell (2013). There are also some power calculators mentioned in the Power Grad talks and [[http://powerandsamplesize.com/Calculators | here including Survival Analysis power computations here]] and for Relative Risk [[https://www.stat.ubc.ca/~rollin/stats/ssize/caco.html | here]] where the calculations are the same as in __Comparing Proportions for Two Independent Samples__ setting p1=p0 (probability of adverse event in the control group) and p2= p0*RR/(1 + p0*(RR - 1)). See Schesselman, J. (1982), Case Control Studies, p. 145. Risk Ratios are ratios of group probabilities of a negative event where Odds Ratios are ratios of the group odds of a negative event [[http://www.theanalysisfactor.com/the-difference-between-relative-risk-and-odds-ratios/ | as described here.]]

Other power calculators [[http://www.ai-therapy.com/psychology-statistics/effect-size-calculator | here.]]

Power computations

Power computations can be performed in SPSS and R only using syntax. For SPSS users Chris Aberson has syntax for power calculations in his book. See reference below.

For a theoretical background and details of specialist software have a look at graduate seminar on power at the Graduate Statistics Programme October-December 2006. There is also a worked example using "Method 2" on a t-test. F. Y. Hsieh, Philip W. Lavori, Harvey J. Cohen and John R. Feussner (2003) An Overview of Variance Inflation Factors for Sample-Size Calculation Eval Health Prof 26 239-257 mentions various formulae for power calculations.

Note in some cases one needs to inflate the total sample size required if there is a natural clustering in the data such as patients being assessed by different exercise therapists. For example if there are b patients assessed by each exercise therapist with an intra-therapist correlation (ICC) then the design effect equals 1 + [(b-1)ICC]. The total sample needs to be multiplied by this design effect to give sample size adjusted for the clustering effect. A further adjustment for variations in cluster sizes can be made (measured by the coefficient of variation) can be incorporated into the formula giving DE = 1 + (b(1+cv^2)-1)ICC. This extra adjustment for differing cluster sizes is not needed for small cvs e.g. cv < 0.23 (see page 26 of this presentation).

Sample sizes required for a given power

Power required for given sample sizes

Additional power freeware (including the popular G*POWER (currently version 3)) is available for download from here. Some examples using G*POWER 3 are in Howell (2013). There are also some power calculators mentioned in the Power Grad talks and here including Survival Analysis power computations here and for Relative Risk here where the calculations are the same as in Comparing Proportions for Two Independent Samples setting p1=p0 (probability of adverse event in the control group) and p2= p0*RR/(1 + p0*(RR - 1)). See Schesselman, J. (1982), Case Control Studies, p. 145. Risk Ratios are ratios of group probabilities of a negative event where Odds Ratios are ratios of the group odds of a negative event as described here.

Other power calculators here.

References

Aberson CL (2010) Applied power analysis for the behavioral sciences. Routledge:London. This book contains examples of computing effect sizes and power using SPSS.

Howell DC (2013) Statistical methods for psychology. 8th Edition. International Edition. Wadsworth:Belmont,CA.

None: FAQ/powprogs (last edited 2017-08-25 12:19:05 by PeterWatson)