get_robust_ucos: Recalibrate and summarize robust uncolocalized events.

In colocboost: Multi-Context Colocalization Analysis for QTL and GWAS Studies

Recalibrate and summarize robust uncolocalized events.

Description

get_robust_ucos get the uncolocalized events by discarding the weaker signals if "ucos_details" with "output_level = 2" exist.

Usage

get_robust_ucos(cb_output, npc_outcome_cutoff = 0.2, pvalue_cutoff = NULL)

Arguments

cb_output	Output object from colocboost analysis
npc_outcome_cutoff	Minimum threshold of normalized probability of colocalized traits in each CoS.
pvalue_cutoff	Maximum threshold of marginal p-values of colocalized variants on colocalized traits in each CoS.

Value

A "colocboost" object with some or all of the following elements:

cos_summary	A summary table for colocalization events.
vcp	The variable colocalized probability for each variable.
cos_details	A object with all information for colocalization results.
data_info	A object with detailed information from input data
model_info	A object with detailed information for colocboost model
ucos_details	A object with information for trait-specific effects if exists after removing weaker signals.

Source

See detailed instructions in our tutorial portal: https://statfungen.github.io/colocboost/articles/Interpret_ColocBoost_Output.html

See Also

Other colocboost_inference: get_ambiguous_colocalization(), get_colocboost_summary(), get_robust_colocalization()

Examples

# colocboost example
set.seed(1)
N <- 1000
P <- 100
# Generate X with LD structure
sigma <- 0.9^abs(outer(1:P, 1:P, "-"))
X <- MASS::mvrnorm(N, rep(0, P), sigma)
colnames(X) <- paste0("SNP", 1:P)
L <- 3
true_beta <- matrix(0, P, L)
true_beta[10, 1] <- 0.5 # SNP10 affects trait 1
true_beta[10, 2] <- 0.4 # SNP10 also affects trait 2 (colocalized)
true_beta[50, 2] <- 0.3 # SNP50 only affects trait 2
true_beta[80, 3] <- 0.6 # SNP80 only affects trait 3
Y <- matrix(0, N, L)
for (l in 1:L) {
  Y[, l] <- X %*% true_beta[, l] + rnorm(N, 0, 1)
}
res <- colocboost(X = X, Y = Y, output_level = 2)
res$ucos_details$ucos$ucos_index
filter_res <- get_robust_ucos(res, npc_outcome_cutoff = 0.2, pvalue_cutoff = 1e-6)
filter_res$ucos_details$ucos$ucos_index

colocboost documentation built on June 7, 2026, 9:08 a.m.