Analyze RNA-seq data with variance component score test accounting
for data heteroscedasticity through precision weights. Perform both gene-wise and gene set
analyses, and can deal with longitudinal data. Method is detailed in: Agniel & Hejblum
(2017)
tcgsaseq is a package for analyzing RNA-seq data. The 2 main functions of the package are varseq and tcgsa_seq:
varseq.tcgsa_seq.The method implemented in this package is detailed in the following article:
The easiest way to get tcgsaseq is to install it from CRAN:
install.packages("tcgsaseq")
Or to get the development version from GitHub:
devtools::install_github("denisagniel/tcgsaseq")
-- Denis Agniel and Boris Hejblum
tcgsaseq R package$pvals output are back to being the same regardless of the test
chosen being asympotic or permutation$pvals output when which_test = permutationspreprocessed = FALSE) in sp_weightssp_weights when facing NA or NaNvaluesWARNING instead of ERROR when a gene set with no measured genes is testedNA support through na.rm_... logical argumentspreprocessed = FALSE)varseq and tcgsa_seq wrapper functions, including
an option to compute heteroscedasticity weights without conditioning on the variable(s)
to be testedvarseq for gene-wise testing