| HubStatus {dhga} | R Documentation |
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
The function returns a list with two components. First component returns the list of genes along with their hub status. Second component gives a table containing number of genes under different categories of hubs viz. housekeeping, unique to stress and unique to normal
Author(s)
Samarendra Das and Baidya Nath Mandal
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)