Skip to content

NameError: name 'mses' is not defined #12

@WeiCSong

Description

@WeiCSong

Hi,

Thanks for this great tool on pleiotropic selection. I ran J-PALM with
python ~/palm/palm.py --traitDir palmdir/ --metadata ~/ldsc/all_gwas/CDmeta --maxp 1 --B 1000 --traits ID602,ID21,ID58,ID285,JIA,ID630

and got the following error

Loading likelihoods...
/lustre/home/acct-bmelgn/bmelgn-3/palm/palm.py:104: RuntimeWarning: divide by zero encountered in true_divide ] 0%
if np.logical_not(np.any(stats.chi2.sf((np.array(loc_betas)/np.array(loc_ses))**2,df=1) < args.maxp)):
[=======================================================================================================================================================] 100%
Analyzing 324 loci...
Traceback (most recent call last):
File "/lustre/home/acct-bmelgn/bmelgn-3/palm/palm.py", line 279, in
_main(args)
File "/lustre/home/acct-bmelgn/bmelgn-3/palm/palm.py", line 256, in _main
omega,ses,margOmega,Mses,D,Dses = _T_inference(statistics,args)
File "/lustre/home/acct-bmelgn/bmelgn-3/palm/palm.py", line 239, in _T_inference
mstats = mcoeffs,mbetas,mses,mx0,mmults,mpvals,mtraitNames
NameError: name 'mses' is not defined

Could you have a look at this issue? Thanks for your help.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions