bulik / ldsc

LD Score Regression (LDSC)
GNU General Public License v3.0
652 stars 345 forks source link

Error while partitioning heritability with the baseline model #77

Open hhautakangas opened 7 years ago

hhautakangas commented 7 years ago

Hi,

does anyone know, what could be wrong when I receive this error message:

Reading annot matrix from ----/Filtered_baseline.[1-22] ... Traceback (most recent call last): File ".../ldsc/ldsc.py", line 644, in sumstats.estimate_h2(args, log) File ".../ldsc/ldscore/sumstats.py", line 372, in estimate_h2 df_results = hsqhat._overlap_output(ref_ld_cnames, overlap_matrix, M_annot, M_tot, args.print_coefficients) File ".../ldsc/ldscore/regressions.py", line 398, in _overlap_output overlap_matrix_prop[i, :] = overlap_matrix[i, :] / M_annot ValueError: could not broadcast input array from shape (53) into shape (1)

Many thanks, Heidi

rkwalters commented 7 years ago

Hi Heidi, The error looks like some sort of mismatch in the expected number of annotations. Is it possible there’s something unusual about your input files for partitioning? (Empty annotations, extra/unused columns, mismatched M values, etc)

Would probably need to know more about your full command line call and the structure of “Filtered_baseline” to debug further.

Cheers, Raymond

On Mar 14, 2017, at 8:24 AM, hhautakangas notifications@github.com wrote:

Hi,

does anyone know, what could be wrong when I receive this error message:

Reading annot matrix from ----/Filtered_baseline.[1-22] ... Traceback (most recent call last): File ".../ldsc/ldsc.py", line 644, in sumstats.estimate_h2(args, log) File ".../ldsc/ldscore/sumstats.py", line 372, in estimate_h2 df_results = hsqhat._overlap_output(ref_ld_cnames, overlap_matrix, M_annot, M_tot, args.print_coefficients) File ".../ldsc/ldscore/regressions.py", line 398, in _overlap_output overlap_matrix_prop[i, :] = overlap_matrix[i, :] / M_annot ValueError: could not broadcast input array from shape (53) into shape (1)

Many thanks, Heidi

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/bulik/ldsc/issues/77, or mute the thread https://github.com/notifications/unsubscribe-auth/AILEvb2RqHAK_tCJeFXHA1jNVqnuOSUZks5rlocPgaJpZM4McgNC.