Python h2_obs_to_liabの例

コード例 #1

ファイル: sumstats.py プロジェクト: prabindm/ldsc

def _get_rg_table(rg_paths, RG, args):
    '''Print a table of genetic correlations.'''
    t = lambda attr: lambda obj: getattr(obj, attr, 'NA')
    x = pd.DataFrame()
    x['p1'] = [rg_paths[0] for i in xrange(1, len(rg_paths))]
    x['p2'] = rg_paths[1:len(rg_paths)]
    x['rg'] = map(t('rg_ratio'), RG)
    x['se'] = map(t('rg_se'), RG)
    x['z'] = map(t('z'), RG)
    x['p'] = map(t('p'), RG)
    if args.samp_prev is not None and args.pop_prev is not None and\
            all((i is not None for i in args.samp_prev)) and all((i is not None for it in args.pop_prev)):
        c = reg.h2_obs_to_liab(1, args.samp_prev[1], args.pop_prev[1])
        x['h2_liab'] = map(lambda x: c * x, map(t('tot'), map(t('hsq2'), RG)))
        x['h2_liab_se'] = map(
            lambda x: c * x, map(t('tot_se'), map(t('hsq2'), RG)))
    else:
        x['h2_obs'] = map(t('tot'), map(t('hsq2'), RG))
        x['h2_obs_se'] = map(t('tot_se'), map(t('hsq2'), RG))

    x['h2_int'] = map(t('intercept'), map(t('hsq2'), RG))
    x['h2_int_se'] = map(t('intercept_se'), map(t('hsq2'), RG))
    x['gcov_int'] = map(t('intercept'), map(t('gencov'), RG))
    x['gcov_int_se'] = map(t('intercept_se'), map(t('gencov'), RG))
    return x.to_string(header=True, index=False) + '\n'

コード例 #2

ファイル: sumstats.py プロジェクト: EugeniaRadulescu/ldsc

def _get_rg_table(rg_paths, RG, args):
    '''Print a table of genetic correlations.'''
    t = lambda attr: lambda obj: getattr(obj, attr, 'NA')
    x = pd.DataFrame()
    x['p1'] = [rg_paths[0] for i in xrange(1, len(rg_paths))]
    x['p2'] = rg_paths[1:len(rg_paths)]
    x['rg'] = map(t('rg_ratio'), RG)
    x['se'] = map(t('rg_se'), RG)
    x['z'] = map(t('z'), RG)
    x['p'] = map(t('p'), RG)
    if args.samp_prev is not None and args.pop_prev is not None and\
            all((i is not None for i in args.samp_prev)) and all((i is not None for it in args.pop_prev)):
        c = reg.h2_obs_to_liab(1, args.samp_prev[1], args.pop_prev[1])
        x['h2_liab'] = map(lambda x: c * x, map(t('tot'), map(t('hsq2'), RG)))
        x['h2_liab_se'] = map(lambda x: c * x,
                              map(t('tot_se'), map(t('hsq2'), RG)))
    else:
        x['h2_obs'] = map(t('tot'), map(t('hsq2'), RG))
        x['h2_obs_se'] = map(t('tot_se'), map(t('hsq2'), RG))

    x['h2_int'] = map(t('intercept'), map(t('hsq2'), RG))
    x['h2_int_se'] = map(t('intercept_se'), map(t('hsq2'), RG))
    x['gcov_int'] = map(t('intercept'), map(t('gencov'), RG))
    x['gcov_int_se'] = map(t('intercept_se'), map(t('gencov'), RG))
    return x.to_string(header=True, index=False) + '\n'