Title: | Biomarker Optimal Segmentation System |
---|---|
Description: | The Biomarker Optimal Segmentation System R package, 'bossR', is designed for precision medicine, helping to identify individual traits using biomarkers. It focuses on determining the most effective cutoff value for a continuous biomarker, which is crucial for categorizing patients into two groups with distinctly different clinical outcomes. The package simultaneously finds the optimal cutoff from given candidate values and tests its significance. Simulation studies demonstrate that 'bossR' offers statistical power and false positive control non-inferior to the permutation approach (considered the gold standard in this field), while being hundreds of times faster. |
Authors: | Liuyi Lan [aut], Xing Li [aut]
|
Maintainer: | Xuekui Zhang <[email protected]> |
License: | GPL-3 |
Version: | 1.0.4 |
Built: | 2025-02-09 02:43:18 UTC |
Source: | https://github.com/cran/bossR |
Computes the regression coefficient of certain regression based on certain cutoff.
getbeta(data, point, type = 2)
getbeta(data, point, type = 2)
data |
A data frame which contains 3 columns for cox regression : survival time, censor status, Biomarker values. 2 columns for linear regression : Y, X. |
point |
A number to cut biomarker or X. |
type |
A number; if =1, will perform linear regression; if =2(default) will perform cox regression. |
An object with 3 class: Coefficient beta, number of samples of which the biomarker is greater than the point, standard error of coefficient estimation.
Given a set of data including survival time ,censor status and Biomarker values, return the Optimal cutoff for the biomarker.
getOC(data, cutoff, type = 2)
getOC(data, cutoff, type = 2)
data |
A data frame which contains 3 columns for cox regression : survival time, censor status, Biomarker values. 2 columns for linear regression : Y, X. |
cutoff |
Numeric vector of candidate cutoffs. |
type |
A number; if =1, will perform linear regression;if =2(default) will perform cox regerssion. |
Optimal cutoff for the biomarker, the FWER of the model
BOSS - Biomarker Optimal Segmentation System
cutoff=c(56,112,167,223,278,334,389,445) data(myGene) getOC(data=myGene,cutoff)
cutoff=c(56,112,167,223,278,334,389,445) data(myGene) getOC(data=myGene,cutoff)
Computes the distribution function of the multivariate normal distribution.
getpvalue(threshold, mu, n, Sigma)
getpvalue(threshold, mu, n, Sigma)
threshold |
A number. |
mu |
The mean vector of length n. |
n |
A number indicates dimension. |
Sigma |
The correlation matrix of dimension n. |
The evaluated distribution function
This data set gives the expression levels of gene data, overall survival time, and censoring status from 500 lung adenocarcinoma cases.
myGene
myGene
A dataframe containing 500 observations of 3 variables.
raw survival data come from https://tau.cmmt.ubc.ca/cSurvival/project_data/TCGA-LUAD/df_survival_o.csv and raw gene expression data come from https://tau.cmmt.ubc.ca/cSurvival/project_data/TCGA-LUAD/df_gene.csv
Xuanjin Cheng, Yongxing Liu, Jiahe Wang, Yujie Chen, Andrew Gordon Robertson, Xuekui Zhang, Steven J M Jones, and Stefan Taubert. (2022) cSurvival: a web resource for biomarker interactions in cancer outcomes and in cell lines. Briefings in Bioinformatics.