-
Notifications
You must be signed in to change notification settings - Fork 5
/
README.Rmd
75 lines (57 loc) · 3.16 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
---
output: github_document
---
<!-- README.md is generated from README.Rmd. Please edit that file -->
```{r, include = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
message = FALSE,
fig.path = "man/figures/README-",
out.width = "100%"
)
```
```{r library, echo=FALSE}
library(badger)
```
# compareRhythms <img src='man/figures/logo.png' align="right" height="138.5" />
<!-- badges: start -->
```{r, echo = FALSE, results='asis'}
cat(
badge_github_version("compareRhythms", "green"),
badge_license(),
badge_codecov(),
badge_lifecycle("stable"),
badge_doi("10.1111/febs.16095", "yellow")
)
```
[](https://zenodo.org/badge/latestdoi/314899899)
<!-- badges: end -->
The goal of *compareRhythms* is to find features with altered circadian rhythm parameters (*amplitude* and *phase*) between the control and experimental groups. A fit of the function $A\cos(2\pi t/24 - \phi)$ to the data defines the rhythm parameters amplitude ($A$) and phase ($\phi$).
## Installation
You can install the current version of *compareRhythms* from GitHub with:
``` r
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install(version = "3.19") # This installs bioconductor, but version 3.15 onwards are also ok
BiocManager::install(c("SummarizedExperiment", "DESeq2", "edgeR", "limma", "rain")) # Packages needed by compareRhythms
install.packages("devtools") # if it is not already installed
devtools::install_github("cran/npsm") # Package archived by CRAN
devtools::install_github("cran/DODR") # Package archived by CRAN
devtools::install_github("bharathananth/compareRhythms", build_vignettes = TRUE, dependencies = TRUE)
```
## Usage summary
The analysis is run using a single function `compareRhythms()`. To execute this function, the three necessary ingredients are the timeseries data, the experimental design and parameters to choose and tune the method. The output of the function is a *data.frame* with the IDs of the differentially rhythmic features, the category it is classified under and optionally the rhythm parameters of the feature in the two groups. The differential rhythmicity categories are **gain** of, **loss** of, **change** of, or **same** rhythms (with respect to the reference/control group).
For complete examples, please refer to the vignette titled **User Guide** by running
```{r userguide, eval=FALSE}
library(compareRhythms)
vignette("UserGuide", package="compareRhythms")
```
or view the [documentation](https://bharathananth.github.io/compareRhythms/articles/UserGuide.html).
## Citation
Please cite *compareRhythms* in publications as:
- Software v1.0.0 https://doi.org/10.5281/zenodo.7699722
- Venn diagram analysis overestimates the extent of circadian rhythm reprogramming. FEBS J, 289: 6605-6621. https://doi.org/10.1111/febs.16095
The R code to perform all the analyses in this publication (using *compareRhythms*) can be found in the repository https://github.com/bharathananth/FEBSJ-VDA-overestimates
## Getting help
If you encounter a bug, please file a minimal reproducible example on [github](https://github.com/bharathananth/compareRhythms/issues).