metricMCB.cv: Calculation of model AUC for Methylation Correlation Blocks...

Description Usage Arguments Value Author(s) References Examples

View source: R/metricMCB_cv.R

Description

To enable quantitative analysis of the methylation patterns within individual Methylation Correlation Blocks across many samples, a single metric to define the methylated pattern of multiple CpG sites within each block. Compound scores which calculated all CpGs within individual Methylation Correlation Blocks by SVM model were used as the compound methylation values of Methylation Correlation Blocks.

Usage

1
metricMCB.cv(MCBset,data_set,Surv,nfold,Method,seed,silent)

Arguments

MCBset

Methylation Correlation Block information returned by the IndentifyMCB function.

data_set

methylation matrix used for training the model in the analysis.

Surv

Survival function contain the survival information for training.

nfold

fold used in the cross validation precedure.

Method

model used to calculate the compound values for multiple Methylation correlation blocks. Options include "svm" "cox" and "lasso". The default option is SVM method.

seed

seed int for cross validation sampling.

silent

Ture indicates that processing information and progress bar will be shown.

Value

Object of class list with elements (XXX will be replaced with the model name you choose):

MCB_matrix Prediction results of model.
auc_results AUC results for each model.

Author(s)

Xin Yu

References

Xin Yu et al. 2019 Predicting disease progression in lung adenocarcinoma patients based on methylation correlated blocks using ensemble machine learning classifiers (under review)

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
#import datasets
data(demo_survival_data)
datamatrix<-create_demo()
data(demo_MCBinformation)
#select MCB with at least 3 CpGs.
demo_MCBinformation<-demo_MCBinformation[demo_MCBinformation[,"CpGs_num"]>2,]

trainingset<-colnames(datamatrix) %in% sample(colnames(datamatrix),0.6*length(colnames(datamatrix)))
testingset<-!trainingset
#create the results using Cox regression. 
mcb_cox_res<-metricMCB.cv(MCBset = demo_MCBinformation,
               data_set = datamatrix,
               Surv = demo_survival_data,
               Method = "cox")

EnMCB documentation built on Dec. 22, 2020, 2 a.m.