All notable changes to this project will be documented in this file.
The format is based on Keep a Changelog and this project adheres to Semantic Versioning.
-
New functions:
doseRiderLMM
anddoseRiderGAMM
.- Description: Introducing support for Linear Mixed Models (LMM) and Generalized Additive Mixed Models (GAMM) through
doseRiderLMM
anddoseRiderGAMM
functions respectively.
- Description: Introducing support for Linear Mixed Models (LMM) and Generalized Additive Mixed Models (GAMM) through
-
New feature in
lmm.R
:- Functions:
fit_lmm
andcreate_lmm_formula
.- Description: The
fit_lmm
function allows fitting of linear mixed models whilecreate_lmm_formula
assists in formula creation for the LMM.
- Description: The
- New Parameter:
omic
.- Description: The parameter lets the user choose between a Gaussian distribution and a negative binomial model, enhancing the flexibility of the model.
- Functions:
-
New function:
smooth_pathway_trend
.- Description: Utilizes the
predict()
function withre.form = NA
to eliminate the fixed effect, providing a smoothed trend for pathway analysis.
- Description: Utilizes the
-
Initial creation of
bmd.R
:- Purpose: Development of functionalities to calculate the Benchmark Dose (BMD). Note: Still in development phase.
-
New option:
center_values
inplot_smooth
.- Description: This option allows for centering and scaling the expression values for each gene.
- Usage: When
center_values
is enabledT
, the expression values of each gene are centered around their mean and scaled by dividing by their standard deviation.
-
New function:
DoseRiderParallel
.- Description: This function allows parallel processing of gene sets using multiple cores. It takes the same arguments as the original DoseRider function but runs the gene set analysis in parallel to improve computation speed.
- Usage:
DoseRiderParallel(se, gmt, dose_col = "dose", sample_col = "sample", covariate = "", omic = "rnaseq", minGSsize = 5, maxGSsize = 300, method = "fdr", num_cores = 5)
- Calculation of
p_value_cubic
inDoseRider
for the comparison between linear and cubic splines.- Description: Previously, only the linear base vs cubic
p.value
was calculated during the comparison. Now, we also calculate thep.value
for the cubic spline. - Reason: To obtain the most non-linear significant gene sets.
- Description: Previously, only the linear base vs cubic