From cfdc4f3a618e460b4a8c9333a9724e824f18bee4 Mon Sep 17 00:00:00 2001 From: Gao Wang Date: Thu, 26 Oct 2023 14:01:22 -0400 Subject: [PATCH] Update documentation --- _sources/code/fine_mapping/SuSiE/SuSiE.ipynb | 9 ++++----- code/fine_mapping/SuSiE/SuSiE.html | 9 ++++----- searchindex.js | 2 +- 3 files changed, 9 insertions(+), 11 deletions(-) diff --git a/_sources/code/fine_mapping/SuSiE/SuSiE.ipynb b/_sources/code/fine_mapping/SuSiE/SuSiE.ipynb index b892aab57..fe8a8ef97 100644 --- a/_sources/code/fine_mapping/SuSiE/SuSiE.ipynb +++ b/_sources/code/fine_mapping/SuSiE/SuSiE.ipynb @@ -385,7 +385,6 @@ " refine=TRUE,\n", " compute_univariate_zscore=FALSE,\n", " coverage=${coverage})\n", - " fitted[[r]]$sets_secondary = susie_get_cs(fitted[[r]], fdat$residual_X_scaled[[r]], coverage=${secondary_coverage})\n", " fitted[[r]]$analysis_time <- proc.time() - st\n", " fitted[[r]] <- post_process_susie(fitted[[r]], fdat, r, signal_cutoff = ${pip_cutoff})\n", " if (!is.null(fitted[[r]]$sets$cs)) {\n", @@ -400,13 +399,13 @@ " }\n", " # generate weights for TWAS using some alternative approaches --- this is not exactly fine-mapping but makes sense to do it here for data production\n", " fitted[[r]]$susie_weights = susie_weights(fdat$fitted[[r]])\n", - " fitted[[r]]$susie_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$susie_weights, fdat$residual_Y_scaled[[r]])\n", + " fitted[[r]]$susie_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$susie_weights, fdat$residual_Y_scaled[[r]])^2\n", " fitted[[r]]$enet_weights = glmnet_weights(fdat$residual_X_scaled[[r]], fdat$residual_Y_scaled[[r]])\n", - " fitted[[r]]$enet_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$enet_weights, fdat$residual_Y_scaled[[r]])\n", + " fitted[[r]]$enet_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$enet_weights, fdat$residual_Y_scaled[[r]])^2\n", " fitted[[r]]$lasso_weights = glmnet_weights(fdat$residual_X_scaled[[r]], fdat$residual_Y_scaled[[r]], alpha = 1)\n", - " fitted[[r]]$lasso_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$lasso_weights, fdat$residual_Y_scaled[[r]])\n", + " fitted[[r]]$lasso_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$lasso_weights, fdat$residual_Y_scaled[[r]])^2\n", " fitted[[r]]$mr_ash_weights = mr_ash_weights(fdat$residual_X_scaled[[r]], fdat$residual_Y_scaled[[r]], beta.init=fitted[[r]]$lasso_weights)\n", - " fitted[[r]]$mr_ash_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$mr_ash_weights, fdat$residual_Y_scaled[[r]])\n", + " fitted[[r]]$mr_ash_r2 = cor(fdat$residual_X_scaled[[r]] %*% fitted[[r]]$mr_ash_weights, fdat$residual_Y_scaled[[r]])^2\n", " }\n", " names(fitted) <- names(fdat$residual_Y_scaled)\n", " saveRDS(fitted, ${_output:ar}, compress='xz')" diff --git a/code/fine_mapping/SuSiE/SuSiE.html b/code/fine_mapping/SuSiE/SuSiE.html index 9bdb16afa..fd1a169c0 100644 --- a/code/fine_mapping/SuSiE/SuSiE.html +++ b/code/fine_mapping/SuSiE/SuSiE.html @@ -782,7 +782,6 @@

Univariate SuSiE