dhga (version 0.1)

HubStatus: Hub status of the genes in a gene co-expression network

Description

The function returns hub status of each gene in a gene co-expression network (GCN)

Usage

HubStatus(pvalue.stress, pvalue.control, alpha)

Arguments

pvalue.stress

pvalue.stress is a vector of p-values for genes in the GCN under stress condition. pvalue.stress can be obtained from the output of hub.stress function.

pvalue.control

pvalue.control is a vector of p-values for genes in the GCN under normal or control condition. pvalue.control can be obtained from the output of hub.control function.

alpha

alpha is a scalar representing statistical level of significance. Default is alpha=0.0001

Value

Examples

data(rice_salt)
pval.stres <- pvalue.hub(rice_salt, beta=6, m=18, s=80, plot=FALSE)
p1 <- pval.stres[, 2]

data(rice_normal)
pval.control <- pvalue.hub(rice_normal, beta=6, m=18, s=80, plot=FALSE)
p2 <- pval.control[, 2]
HubStatus(p1,p2,alpha=0.0001)