mlhighHet: mlhighHet Performs heterogeneity analysis in gene expression

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/mlhighHet.R

Description

This function extracts features based on ML method, finds optimal cut-off values of features using sequencial Cox PH model and obtain the most consistent level according to the cut-offs.

Usage

1
mlhighHet(cols, idSurv, idEvent, idFrail, num, fold = 3, data)

Arguments

cols

A numeric vector of column numbers indicating the features for which the log Loss functions are to be computed

idSurv

The name of the survival time variable

idEvent

The name of the survival event variable

idFrail

The name of the frailty variable

num

Number of features to be selected

fold

An integer denoting number of folds in cross validation, default value 3

data

A data frame that contains the survival and covariate information for the subjects

Details

This function extracts features based on minimum log-Loss function using Cox proportional hazard model as learner method on a high dimensional survival data. For those selected genes, we obtain optimal cutoff values using minimum p-value in a Cox PH model. The Cox PH model is used sequencially for each combination of genes and all possible gene combinations are tested to obtain best possible combination with minimum BIC value. The subjects are classified according to different levels of those genes. Using a Cox PH frailty model, we obtain the most consistent level for which the frailty variance is minimum. The data is splited using cross validation technique. The performance measure is considered as logarithmic loss function. It is defined as,

L(f,t)=-log(f(t))

The CoxPH frailty model is defined as,

λ(t)=λ 0(t)ν exp{X'β}

where ν is called the frailty. The variance of the frailty term is considered as the heterogeneity among the subjects or patients. Gaussian distribution with mean 0 is considered for the distribution of frailty component.

Value

dataframes containing optimal gene cutoff values and most consistent level according to those cut-offs with frailty variance.

Author(s)

Atanu Bhattacharjee, Gajendra K. Vishwakarma & Souvik Banerjee

See Also

mlhighCox, mlhighFrail

Examples

1
2
data(hnscc)
mlhighHet(cols=c(27:32), idSurv="OS", idEvent="Death", idFrail="ID", num=2, fold = 3, data=hnscc)

highMLR documentation built on May 11, 2021, 5:08 p.m.

Related to mlhighHet in highMLR...