Skip to main content
100%

Gist b989c26735fde8214c19

✓ Published0🌍 Public
Aakeshavan
Last edited Sep 4, 2014
Created on Sep 4, 2014

This example demonstrates statistical power calculations for several genetic association models (DRAB1_sGMV, HLAB_ucca, HLAB_sGM, NFkb_WM) using R’s noncentral F distribution via the rpy2 interface, with power derived from effect sizes and sample sizes listed in a pandas DataFrame. It imports the `stats` package and uses `statsmodels.api` alongside scipy.stats, defining a custom `pf` function that calls R’s `qf` and `pf` functions to compute critical values and noncentral power. The notebook also loads precomputed intraclass correlation coefficients from a pickled file and displays the power table as an HTML-rendered pandas DataFrame.

AI-generated description

Similar vizzes