View source: R/grin.assoc.lsn.outcome.R
grin.assoc.lsn.outcome | R Documentation |
Performs statistical association analysis between binary gene-lesion events and clinical outcomes of interest, including binary outcomes (e.g., Minimal Residual Disease (MRD)) and time-to-event outcomes (e.g., Event-Free Survival (EFS), and Overall Survival (OS)). The function supports both univariate and covariate-adjusted logistic regression and Cox proportional hazards models.
grin.assoc.lsn.outcome(
lsn.mtx,
clin.data,
annotation.data,
clinvars,
covariate = NULL
)
lsn.mtx |
A binary lesion matrix where each row represents a unique gene-lesion pair (e.g., |
clin.data |
A clinical data |
annotation.data |
A gene annotation |
clinvars |
A character vector of clinical outcome variables to analyze. Binary variables (e.g., MRD) should be coded as |
covariate |
Optional. A character vector specifying covariates to include as model adjustments (e.g., risk group, age, gender, etc...). |
For each gene-lesion pair in the binary lesion matrix, the function can performs:
Logistic regression for binary outcomes (e.g., MRD), producing odds ratios (OR), 95_confidence intervals (CI), p-values, and FDR-adjusted q-values.
Cox proportional hazards models for survival outcomes (e.g., EFS, OS), producing hazard ratios (HR), 95\
Models can optionally be adjusted for covariates such as clinical or demographic factors. Summary counts of patients with and without lesions, stratified by outcome status, are also included in the output.
A results data.frame
containing gene annotation and association statistics for each gene-lesion pair across the specified clinical outcomes. The output includes:
Odds ratio (OR), lower and upper 95CI, p-value, and q-value (FDR) for logistic regression models.
Hazard ratio (HR), lower and upper 95CI, p-value, and q-value for Cox proportional hazards models.
Patient counts for those with and without lesions, and corresponding outcome event statuses.
Abdelrahman Elsayed abdelrahman.elsayed@stjude.org and Stanley Pounds stanley.pounds@stjude.org
Andersen, P. K., & Gill, R. D. (1982). Cox's regression model for counting processes: A large sample study.
Therneau, T. M., & Grambsch, P. M. (2000). Modeling Survival Data: Extending the Cox Model.
Dobson, A. J. (1990). An Introduction to Generalized Linear Models.
prep.binary.lsn.mtx
, coxph
, glm
data(lesion_data)
data(hg38_gene_annotation)
data(clin_data)
# Step 1: Prepare gene-lesion overlap
gene.lsn <- prep.gene.lsn.data(lesion_data, hg38_gene_annotation)
gene.lsn.overlap <- find.gene.lsn.overlaps(gene.lsn)
# Step 2: Create a binary lesion matrix (minimum 5 patients per lesion)
lsn.binary.mtx <- prep.binary.lsn.mtx(gene.lsn.overlap, min.ngrp = 5)
# Step 3: Create survival objects and add to clinical data
library(survival)
clin_data$EFS <- Surv(clin_data$efs.time, clin_data$efs.censor)
clin_data$OS <- Surv(clin_data$os.time, clin_data$os.censor)
# Step 4: Specify outcomes of interest
clinvars <- c("MRD.binary", "EFS", "OS")
# Step 5: Run association analysis
assc.outcomes <- grin.assoc.lsn.outcome(lsn.binary.mtx,
clin_data,
hg38_gene_annotation,
clinvars)
# Optional: Adjust for covariates using the 'covariate' argument
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.