Skip to content

change susie_suff_stat to susierss #210

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
beckse wants to merge 1 commit into
base: master
Choose a base branch
from
Open

change susie_suff_stat to susierss #210

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 4 additions & 4 deletions finemapper.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -832,9 +832,9 @@ def finemap(self, locus_start, locus_end, num_causal_snps, use_prior_causal_prob
assert self.df_ld.notnull().all().all()
if residual_var is not None: residual_var_init = residual_var

if hasattr(self.susieR, 'susie_suff_stat'):
logging.info('Using susieR::susie_suff_stat()')
susie_obj = self.susieR.susie_suff_stat(
if hasattr(self.susieR, 'susierss'):
logging.info('Using susieR::susierss()')
susie_obj = self.susieR.susierss(
bhat=bhat.reshape((m,1)),
shat=np.ones((m,1)),
R=self.df_ld.values,
Expand Down Expand Up @@ -865,7 +865,7 @@ def finemap(self, locus_start, locus_end, num_causal_snps, use_prior_causal_prob
prior_weights=(prior_weights.reshape((m,1)) if use_prior_causal_prob else self.R_null)
)
else:
raise NotImplementedError('Only susie_suff_stat() and susie_bhat() are supported. Check your version of susieR')
raise NotImplementedError('Only susierss() and susie_bhat() are supported. Check your version of susieR')
susie_time = time.time()-t0
logging.info('Done in %0.2f seconds'%(susie_time))

Expand Down