diff --git a/fine-mapping/finemapping_results_wrangler.ipynb b/fine-mapping/finemapping_results_wrangler.ipynb index 8bfa97f..51a05c3 100644 --- a/fine-mapping/finemapping_results_wrangler.ipynb +++ b/fine-mapping/finemapping_results_wrangler.ipynb @@ -325,8 +325,8 @@ " idx = 0\n", " for (f in files) {\n", " idx = idx + 1\n", - " rds_file = dscrutils::read_dsc(f)\n", " unit = dirname(f)\n", + " rds_file = dscrutils::dscread(unit, tools::file_path_sans_ext(basename(f)))\n", " processed_dat[[idx]] = get_${source}_output(unit, rds_file)\n", " cat(\"We are at unit #\", idx, \":\", unit, \"\\n\")\n", " }\n", diff --git a/fine-mapping/summary_statistics_finemapping.ipynb b/fine-mapping/summary_statistics_finemapping.ipynb index e77c526..8a507ff 100644 --- a/fine-mapping/summary_statistics_finemapping.ipynb +++ b/fine-mapping/summary_statistics_finemapping.ipynb @@ -376,7 +376,7 @@ " }\n", " library(susieR)\n", " L = min(${L}, length(z_scores))\n", - " s = susie_z(z_scores, r, r_tol = 5e-3, L = L, min_abs_corr = 0.1,\n", + " s = susie_rss(z_scores, r, r_tol = 5e-3, L = L, min_abs_corr = 0.1,\n", " prior_weights=prior, track_fit=${\"T\" if track_fit else \"F\"})\n", " s$prior = prior\n", " s$z = z_scores\n",