snipar.correlate module

snipar.correlate.compute_corr(z1, z2, v1, v2, r, l)[source]
snipar.correlate.estimate_cov(z1, z2, v1, v2, r, w)[source]
snipar.correlate.estimate_cov_weights(v1, v2, r, var_1, var_2, c, l)[source]
snipar.correlate.estimate_var(z, v, w)[source]
snipar.correlate.estimate_weights(v_est, v, l)[source]
snipar.correlate.jacknife(z1, z2, v1, v2, r, l, n_blocks, block_size)[source]
snipar.correlate.jacknife_est(z1, z2, v1, v2, r, l, n_blocks)[source]
snipar.correlate.read_sumstats_file(sumstats_file, chrom)[source]
snipar.correlate.read_sumstats_files(sumstats_files, chroms)[source]
snipar.correlate.reweighted_estimate_cov(z1, z2, v1, v2, r, l, var_1, var_2, tol=1e-08, max_iter=1000)[source]
snipar.correlate.reweighted_estimate_var(z, v, l, tol=1e-08, max_iter=1000)[source]
class snipar.correlate.sumstats(chrom, sid, pos, A1, A2, freqs, direct, direct_SE, avg_NTC, avg_NTC_SE, population, population_SE, r_direct_avg_NTC, r_direct_pop, ldscores=None, map=None)[source]

Bases: object

compute_ld_scores(bedfiles, chroms, ld_wind, ld_out=None)[source]
concatenate(s2)[source]
cor_direct_avg_NTC(n_blocks)[source]
cor_direct_pop(n_blocks)[source]
filter(filter_pass)[source]
filter_NAs()[source]
filter_corrs(max_Z)[source]
filter_maf(min_maf)[source]
scores_from_ldsc(ld_files)[source]