dGSEAwrite
is supposed to write out enrichment results.
dGSEAwrite(
eTerm,
which_content = c("gadjp", "adjp", "pvalue", "FWER", "FDR", "qvalue",
"nES", "ES"),
which_score = c("gadjp", "adjp", "FWER", "FDR", "qvalue", "nES"),
cutoff = 0.1,
filename = NULL,
keep.significance = T
)
an object of class "eTerm"
the content will be written out. It includes two categories: i) based on "adjp" for adjusted p value, "gadjp" for globally adjusted p value, "pvalue" for p value, "FWER" for family-wise error rate, "FDR" for false discovery rate, "qvalue" for q value; ii) based on "ES" for enrichment score, "nES" for normalised enrichment score. For the former, the content is : first -1*log10-transformed, and then multiplied by -1 if nES is negative.
which statistics/score will be used for declaring the significance. It can be "adjp" for adjusted p value, "gadjp" for globally adjusted p value, "FWER" for family-wise error rate, "FDR" for false discovery rate, "qvalue" for q value
a cutoff to declare the signficance. It should be used together with 'which_score'
a character string naming a filename
logical to indicate whether or not to mask those insignfiicant by NA. By default, it sets to true to mask those insignfiicant by NA
a data frame with following components:
setID
: term ID
setSize
: the number of genes in the set
name
: term name
namespace
: term namespace
distance
: term distance
sample names
: sample names in the next columns
# NOT RUN {
#output <- dGSEAwrite(eTerm, which_content="gadjp", which_score="gadjp", filename="eTerm.txt")
# }
Run the code above in your browser using DataLab