benchmark/01a_run_kronos.R
In thomazbastiaanssen/kronos: Microbiome Oriented Circadian Rhythm Analysis Toolkit
library(tidyverse)
library(kronos)
data = species
# row.names(data) = data[, 1]
# data = data[,-1]
data = data[,metadata$woltka_ID]
out_list = fw_kronos(data, formula = ~ time(Timepoint), metadata = metadata, pairwise = F)
res_kronos = do.call(rbind, lapply(FUN = getKronos_groupwise, X = out_list))[,-1]
rm(list=setdiff(ls(), c("species", "metadata","res_jtk", "res_kronos", "res_limo", "res_cosinor2")))
thomazbastiaanssen/kronos documentation built on May 1, 2023, 6:45 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
library(tidyverse)
library(kronos)
data = species
# row.names(data) = data[, 1]
# data = data[,-1]
data = data[,metadata$woltka_ID]
out_list = fw_kronos(data, formula = ~ time(Timepoint), metadata = metadata, pairwise = F)
res_kronos = do.call(rbind, lapply(FUN = getKronos_groupwise, X = out_list))[,-1]
rm(list=setdiff(ls(), c("species", "metadata","res_jtk", "res_kronos", "res_limo", "res_cosinor2")))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.