This package provides a collection of utility functions for community ecology analyses, with emphasis on workflows using the phyloseq and vegan packages. Includes functions for normalizing OTU tables, computing alpha diversity via rarefaction (using a fast C++ implementation), differential abundance comparisons with compact letter displays, primer checking for amplicon sequencing, plotting QIIME 2/DADA2 generated transition stats and miscellaneous helpers for ordination plots and taxonomic name formatting.
Create a hexagon tile map display from spatial polygons. Each polygon is represented by a hexagon tile, placed as close to it's original centroid as possible, with a focus on maintaining spatial relationship to a focal point. Developed to aid visualisation and analysis of spatial distributions across Australia, which can be challenging due to the concentration of the population on the coast and wide open interior.
This package provides methods for metabolomics and spectral data analysis, including data import, preprocessing, visualization, univariate and multivariate analysis, machine learning, feature selection, and pathway analysis. The package supports analytical workflows for different data types used in metabolomics and spectroscopy. Some optional functionality uses the suggested packages cyjShiny and specmine.datasets'. The package specmine.datasets is maintained separately at <https://github.com/PedroFontao/specmine.datasets>.
This package provides predictive accuracy tools to evaluate time-to-event survival models. This includes calculating the concordance probability estimate that incorporates the follow-up time for a particular study developed by Devlin, Gonen, Heller (2020)<doi:10.1007/s10985-020-09503-3>. It also evaluates the concordance probability estimate for nested Cox proportional hazards models using a projection-based approach by Heller and Devlin (under review).
Process and summarize DAS data files. These files are typically, but do not have to be DAS <https://swfsc-publications.fisheries.noaa.gov/publications/TM/SWFSC/NOAA-TM-NMFS-SWFSC-305.PDF> data produced by the Southwest Fisheries Science Center (SWFSC) program WinCruz'. This package standardizes and streamlines basic DAS data processing, and includes a PDF with the DAS data format requirements expected by the package.
Tensor-train is a compact representation for higher-order tensors. Some algorithms for performing tensor-train decomposition are available such as TT-SVD, TT-WOPT, and TT-Cross. For the details of the algorithms, see I. V. Oseledets (2011) <doi:10.1137/090752286>, Yuan Longao, et al (2017) <doi:10.48550/arXiv.1709.02641>, I. V. Oseledets (2010) <doi:10.1016/j.laa.2009.07.024>.
This package provides tools for converting data from complex or irregular layouts to a columnar structure. For example, tables with multilevel column or row headers, or spreadsheets. Header and data cells are selected by their contents and position, as well as formatting and comments where available, and are associated with one other by their proximity in given directions. Functions for data frames and HTML tables are provided.
This package offers functions to process multiple ChIP-seq BAM files and detect allele-specific events. It computes allele counts at individual variants (SNPs/SNVs), implements extensive QC (quality control) steps to remove problematic variants, and utilizes a Bayesian framework to identify statistically significant allele-specific events. BaalChIP is able to account for copy number differences between the two alleles, a known phenotypical feature of cancer samples.
The first day of any MMWR week is Sunday. MMWR week numbering is sequential beginning with 1 and incrementing with each week to a maximum of 52 or 53. MMWR week #1 of an MMWR year is the first week of the year that has at least four days in the calendar year. This package provides functionality to convert dates to MMWR day, week, and year and the reverse.
Phylogenetic clustering (phyloclustering) is an evolutionary continuous time Markov Chain model-based approach to identify population structure from molecular data without assuming linkage equilibrium. The package phyclust provides a convenient implementation of phyloclustering for DNA and SNP data, capable of clustering individuals into subpopulations and identifying molecular sequences representative of those subpopulations. It is designed in C for performance and interfaced with R for visualization.
squallms is a Bioconductor R package that implements a "semi-labeled" approach to untargeted mass spectrometry data. It pulls in raw data from mass-spec files to calculate several metrics that are then used to label MS features in bulk as high or low quality. These metrics of peak quality are then passed to a simple logistic model that produces a fully-labeled dataset suitable for downstream analysis.
An R wrapper for agena.ai <https://www.agena.ai> which provides users capabilities to work with agena.ai using the R environment. Users can create Bayesian network models from scratch or import existing models in R and export to agena.ai cloud or local API for calculations. Note: running calculations requires a valid agena.ai API license (past the initial trial period of the local API).
It calculates the Air Pollution Tolerance Index (APTI) of plant species using biochemical parameters such as chlorophyll content, leaf extract pH, relative water content, and ascorbic acid content. It helps in identifying tolerant species for greenbelt development and pollution mitigation studies. It includes a shiny app for interactive APTI calculation and visualisation. For method details see, Sahu et al. (2020).<DOI:10.1007/s42452-020-3120-6>.
Fits joint Bayesian four-parameter logistic (4PL) models to thermal-tolerance proportion data, extracts the classical thermal load sensitivity quantities (z, CTmax at 1 hour, T_crit) with full posterior uncertainty, and predicts heat-injury accumulation and survival under fluctuating temperature regimes with optional Sharpe-Schoolfield repair. Models are fitted with Stan via the brms package. Implements the framework described in Noble, Arnold, Nakagawa and Pottier (in preparation).
The Bayesian Adjustment for Confounding (BAC) algorithm (Wang et al., 2012) can be used to estimate the causal effect of a continuous exposure on a continuous outcome. This package provides an approximate sensitivity analysis of BAC with regards to the hyperparameter omega. BACprior also provides functions to guide the user in their choice of an appropriate omega value. The method is based on Lefebvre, Atherton and Talbot (2014).
Selection of informative features like genes, transcripts, RNA seq, etc. using Bootstrap Maximum Relevance and Minimum Redundancy technique from a given high dimensional genomic dataset. Informative gene selection involves identification of relevant genes and removal of redundant genes as much as possible from a large gene space. Main applications in high-dimensional expression data analysis (e.g. microarray data, NGS expression data and other genomics and proteomics applications).
Conducts conditional random sampling on observed values in sparse matrices. Useful for training and test set splitting sparse matrices prior to model fitting in cross-validation procedures and estimating the predictive accuracy of data imputation methods, such as matrix factorization or singular value decomposition (SVD). Although designed for applications with sparse matrices, CRASSMAT can also be applied to complete matrices, as well as to those containing missing values.
One way to choose the number of principal components is via the reconstruction error. This package is designed mainly for this purpose. Graphical representation is also supported, plus some other principal component analysis related functions. References include: Jolliffe I.T. (2002). Principal Component Analysis. <doi:10.1007/b98835> and Mardia K.V., Kent J.T. and Bibby J.M. (1979). Multivariate Analysis. ISBN: 978-0124712522. London: Academic Press.
Calculates and visualises cumulative percent decay curves, which are typically calculated from metagenomic taxonomic profiles. These can be used to estimate the level of expected endogenous taxa at different abundance levels retrieved from metagenomic samples, when comparing to samples of known sampling site or source. Method described in Fellows Yates, J. A. et. al. (2021) Proceedings of the National Academy of Sciences USA <doi:10.1073/pnas.2021655118>.
This package provides S4 classes and methods for storing dense and sparse matrices in DuckDB databases. The package supports constructing database-backed matrices from base R and Matrix objects, extracting slices and summaries, performing arithmetic and selected linear algebra operations, and materializing results for larger-than-memory workflows. It integrates with dbProject to keep database paths, live connections, and lazy matrix tables synchronized across interactive analyses.
Automatically performs exploratory data analysis (EDA) for tabular datasets, including data summaries, missing value analysis, descriptive statistics, visualizations, correlation analysis, outlier detection, and automated report generation. The package provides a streamlined workflow for rapid data exploration and produces publication-ready tables and graphics. For methodological details see Tukey (1977, ISBN:9780201076165), Pearson (1895) <doi:10.1098/rspl.1895.0041>, and Wickham (2014) <doi:10.18637/jss.v059.i10>.
This package provides a syntax to create and combine Gaussian process kernels in greta'. You can then use these kernels to define either full rank or sparse Gaussian processes, and project a fitted process to new coordinates. This lets you include Gaussian process components in the statistical models you write with greta'. This is an extension to the greta software, Golding (2019) <doi:10.21105/joss.01601>.
This package implements hierarchical conformal prediction for clustered data with missing responses. The method uses repeated cluster-level splitting and within-cluster subsampling to accommodate dependence, and inverse-probability weighting to correct distribution shift induced by missingness. Conditional densities are estimated by inverting fitted conditional quantiles (linear quantile regression or quantile regression forests), and p-values are aggregated across resampling and splitting steps using the Cauchy combination test.
This package provides a multi-layered untargeted pipeline for high-throughput LC/HRMS data processing to extract signals of organic small molecules. The package performs ion pairing, peak detection, peak table alignment, retention time correction, aligned peak table gap filling, peak annotation and visualization of extracted ion chromatograms (EICs) and total ion chromatograms (TICs). The IDSL.IPA package was introduced in <doi:10.1021/acs.jproteome.2c00120> .