From 6af0a986c82b5a22e75001642121c8732f354988 Mon Sep 17 00:00:00 2001 From: xuewei cao <36172337+xueweic@users.noreply.github.com> Date: Tue, 15 Apr 2025 20:56:44 -0400 Subject: [PATCH] minor fix --- vignettes/ColocBoost_tutorial_advance.Rmd | 10 +++++----- vignettes/ColocBoost_tutorial_basic.Rmd | 17 +++++++++-------- ...locBoost_tutorial_strong_colocalization.Rmd} | 0 3 files changed, 14 insertions(+), 13 deletions(-) rename vignettes/{ColocBoost_tutorial_cbsummary.Rmd => ColocBoost_tutorial_strong_colocalization.Rmd} (100%) diff --git a/vignettes/ColocBoost_tutorial_advance.Rmd b/vignettes/ColocBoost_tutorial_advance.Rmd index 22c93cc..7ecfee9 100644 --- a/vignettes/ColocBoost_tutorial_advance.Rmd +++ b/vignettes/ColocBoost_tutorial_advance.Rmd @@ -34,7 +34,7 @@ data("Ind_5traits") X <- Ind_5traits$X[[1]] Y <- do.call(cbind, Ind_5traits$Y) res <- colocboost(X = X, Y = Y) -res$summary_table +res$cos_summary ``` **Example of using summary statistics** @@ -44,9 +44,9 @@ For summary statistics, similar optimization is applied using `Sumstat_5traits` ```{r oneLD} data("Sumstat_5traits") sumstat <- Sumstat_5traits$sumstat -LD <- Sumstat_5traits$LD +LD <- get_cormat(Ind_5traits$X[[1]]) res <- colocboost(sumstat = sumstat, LD = LD) -res$summary_table +res$cos_summary ``` **Example of combining individual-level data and summary statistics** @@ -60,9 +60,9 @@ X <- Ind_5traits$X[[1]] Y <- do.call(cbind, Ind_5traits$Y[1:3]) data("Sumstat_5traits") sumstat <- Sumstat_5traits$sumstat[4:5] -LD <- Sumstat_5traits$LD +LD <- get_cormat(Ind_5traits$X[[1]]) res <- colocboost(X = X, Y = Y, sumstat = sumstat, LD = LD) -res$summary_table +res$cos_summary ``` diff --git a/vignettes/ColocBoost_tutorial_basic.Rmd b/vignettes/ColocBoost_tutorial_basic.Rmd index 99c37fc..47a09be 100644 --- a/vignettes/ColocBoost_tutorial_basic.Rmd +++ b/vignettes/ColocBoost_tutorial_basic.Rmd @@ -39,9 +39,9 @@ This command initiates the colocalization analysis, applying the ColocBoost meth After running the analysis, you can explore the results to identify colocalized variants and review the summary statistics. This output will provide insights into which variants are colocalized across the different phenotypes and offer a comprehensive overview of the statistical results from the colocalization analysis. ```{r indResults} -res$coloc_results$coloc_csets$csets_index +res$cos_details$cos$cos_index -res$summary_table +res$cos_summary ``` @@ -54,7 +54,8 @@ This tutorial demonstrates how to analyze summary statistics data using `colocbo To get started, load the `Sumstat_5traits` dataset into your R session. Note: The `Sumstat_5traits` dataset includes only one LD matrix that applies to all traits. To demonstrate handling multiple traits, we replicate this single LD matrix for each trait as follows. To analyze the data using summary statistics, apply the colocboost function specifying the summary statistics and LD matrices. ```{r sumstat} data("Sumstat_5traits") -LD <- lapply(1:5, function(i) Sumstat_5traits$LD) +data("Ind_5traits") +LD <- get_cormat(Ind_5traits$X[[1]]) res <- colocboost(sumstat = Sumstat_5traits$sumstat, LD = LD) ``` *Note*: This step duplicates the single LD matrix into a list of five matrices, one for each trait. This is to mimic scenarios where different traits might have different LD structures. ColocBoost allows for the input of a single LD matrix if the LD across traits is consistent. For more advanced usage involving different LD matrices or more complex setups, please refer to the advanced tutorial (URL). @@ -64,9 +65,9 @@ res <- colocboost(sumstat = Sumstat_5traits$sumstat, LD = LD) After running the analysis, you can explore the results to identify colocalized variants and review the summary statistics. This output will provide insights into which variants are colocalized across the different phenotypes and offer a comprehensive overview of the statistical results from the colocalization analysis. ```{r sumstatResults} -res$coloc_results$coloc_csets$csets_index +res$cos_details$cos$cos_index -res$summary_table +res$cos_summary ``` This section of the analysis provides insights into which variants are colocalized across the different phenotypes and offers a comprehensive overview of the statistical results from the colocalization analysis. @@ -85,7 +86,7 @@ data("Sumstat_5traits") X <- Ind_5traits$X[1:3] Y <- Ind_5traits$Y[1:3] sumstat <- Sumstat_5traits$sumstat[4:5] -LD <- lapply(1:2, function(i) Sumstat_5traits$LD) +LD <- get_cormat(Ind_5traits$X[[1]]) ``` *Note*: This step duplicates the single LD matrix into a list of two matrices, one for each trait. This is to mimic scenarios where different traits might have different LD structures. ColocBoost allows for the input of a single LD matrix if the LD across traits is consistent. For more advanced usage involving different LD matrices or more complex setups, please refer to the advanced tutorial (URL). @@ -98,9 +99,9 @@ res <- colocboost(X = X, Y = Y, sumstat = sumstat, LD = LD) After running the analysis, you can explore the results to identify colocalized variants and review the summary statistics. This output will provide insights into which variants are colocalized across the different phenotypes and offer a comprehensive overview of the statistical results from the colocalization analysis. ```{r mixResults} -res$coloc_results$coloc_csets$csets_index +res$cos_details$cos$cos_index -res$summary_table +res$cos_summary ``` This section of the analysis provides insights into which variants are colocalized across the different phenotypes and offers a comprehensive overview of the statistical results from the colocalization analysis. diff --git a/vignettes/ColocBoost_tutorial_cbsummary.Rmd b/vignettes/ColocBoost_tutorial_strong_colocalization.Rmd similarity index 100% rename from vignettes/ColocBoost_tutorial_cbsummary.Rmd rename to vignettes/ColocBoost_tutorial_strong_colocalization.Rmd