We introduce factor models designed to jointly analyze high-dimensional count data from multiple studies by extracting study-shared and specified factors. Our factor models account for heterogeneous noises and overdispersion among counts with augmented covariates. We propose an efficient and speedy variational estimation procedure for estimating model parameters, along with a novel criterion for selecting the optimal number of factors and the rank of regression coefficient matrix. More details can be referred to Liu et al. (2024) <doi:10.48550/arXiv.2402.15071>.

Evolutionary black box optimization algorithms building on the bbotk package. miesmuschel offers both ready-to-use optimization algorithms, as well as their fundamental building blocks that can be used to manually construct specialized optimization loops. The Mixed Integer Evolution Strategies as described by Li et al. (2013) <doi:10.1162/EVCO_a_00059> can be implemented, as well as the multi-objective optimization algorithms NSGA-II by Deb, Pratap, Agarwal, and Meyarivan (2002) <doi:10.1109/4235.996017>.

It provides functions to compute the values of different modifications of the Rand and Wallace indices. The indices are used to measure the stability or similarity of two partitions obtained on two different sets of units with a non-empty intercept. Splitting and merging of clusters can (depends on the selected index) have a different effect on the value of the indices. The indices are proposed in Cugmas and Ferligoj (2018) <http://ibmi.mf.uni-lj.si/mz/2018/no-1/Cugmas2018.pdf>.

This package provides a likelihood-based approach to modeling species distributions using presence-only data. In contrast to the popular software program MAXENT, this approach yields estimates of the probability of occurrence, which is a natural descriptor of a species distribution.

This package provides functions, which make matrix creation conciser (such as the core package's function m() for rowwise matrix definition or runifm() for random value matrices). Allows to set multiple matrix values at once, by using list of formulae. Provides additional matrix operators and dedicated plotting function.

Alternative implementation of the beautiful MissForest algorithm used to impute mixed-type data sets by chaining random forests, introduced by Stekhoven, D.J. and Buehlmann, P. (2012) <doi:10.1093/bioinformatics/btr597>. Under the hood, it uses the lightning fast random forest package ranger'. Between the iterative model fitting, we offer the option of using predictive mean matching. This firstly avoids imputation with values not already present in the original data (like a value 0.3334 in 0-1 coded variable). Secondly, predictive mean matching tries to raise the variance in the resulting conditional distributions to a realistic level. This would allow, e.g., to do multiple imputation when repeating the call to missRanger(). Out-of-sample application is supported as well.

This package offers three important components: (1) to construct a use-defined linear mixed model, (2) to employ one of linear mixed model approaches: minimum norm quadratic unbiased estimation (MINQUE) (Rao, 1971) for variance component estimation and random effect prediction; and (3) to employ a jackknife resampling technique to conduct various statistical tests. In addition, this package provides the function for model or data evaluations.This R package offers fast computations for large data sets analyses for various irregular data structures.

Perform sensitivity analysis on ordinary differential equation based models, including ad-hoc graphical analyses based on structured sequences of parameters as well as local sensitivity analysis. Functions are provided for creating inputs, simulating scenarios and plotting outputs.

Model based clustering using the multivariate multiple Scaled t (MST) and multivariate multiple scaled contaminated normal (MSCN) distributions. The MST is an extension of the multivariate Student-t distribution to include flexible tail behaviors, Forbes, F. & Wraith, D. (2014) <doi:10.1007/s11222-013-9414-4>. The MSCN represents a heavy-tailed generalization of the multivariate normal (MN) distribution to model elliptical contoured scatters in the presence of mild outliers (also referred to as "bad" points) and automatically detect bad points, Punzo, A. & Tortora, C. (2021) <doi:10.1177/1471082X19890935>.

Generates derived parameter(s) from Monte Carlo Markov Chain (MCMC) samples using R code. This allows Bayesian models to be fitted without the inclusion of derived parameters which add unnecessary clutter and slow model fitting. For more information on MCMC samples see Brooks et al. (2011) <isbn:978-1-4200-7941-8>.

This package provides tools for performing mathematical morphology operations, such as erosion and dilation, on data of arbitrary dimensionality. Can also be used for finding connected components, resampling, filtering, smoothing and other image processing-style operations.

Given a set of data points, a clustering is defined as a disjoint partition where each pair of sets in a partition has no overlapping elements. This package provides 25 methods that play a role somewhat similar to distance or metric that measures similarity of two clusterings - or partitions. For a more detailed description, see Meila, M. (2005) <doi:10.1145/1102351.1102424>.

Response Surface Designs (RSDs) involving factors not all at same levels are called Mixed Level RSDs (or Asymmetric RSDs). In many practical situations, RSDs with asymmetric levels will be more suitable as it explores more regions in the design space. (J.S. Mehta and M.N. Das (1968) <doi:10.2307/1267046>. "Asymmetric rotatable designs and orthogonal transformations").This package contains function named ATORDs_I() for generating asymmetric third order rotatable designs (ATORDs) based on third order designs given by Das and Narasimham (1962). Function ATORDs_II() generates asymmetric third order rotatable designs developed using t-design of unequal set sizes, which are smaller in size as compared to design generated by function ATORDs_I(). In general, third order rotatable designs can be classified into two classes viz., designs that are suitable for sequential experimentation and designs for non-sequential experimentation. The sequential experimentation approach involves conducting the trials step by step whereas, in the non-sequential experimentation approach, the entire runs are executed in one go (M. N. Das and V. Narasimham (1962) <doi:10.1214/AOMS/1177704374>. "Construction of Rotatable Designs through Balanced Incomplete Block Designs"). ATORDs_I() and ATORDs_II() functions generate non-sequential asymmetric third order designs. Function named SeqTORD() generates symmetric sequential third order design in blocks and also gives G-efficiency of the given design. Function named Asymseq() generates asymmetric sequential third order designs in blocks (M. Hemavathi, Eldho Varghese, Shashi Shekhar and Seema Jaggi (2020) <doi:10.1080/02664763.2020.1864817>. "Sequential asymmetric third order rotatable designs (SATORDs)"). In response surface design, situations may arise in which some of the factors are qualitative in nature (Jyoti Divecha and Bharat Tarapara (2017) <doi:10.1080/08982112.2016.1217338>. "Small, balanced, efficient, optimal, and near rotatable response surface designs for factorial experiments asymmetrical in some quantitative, qualitative factors"). The Function named QualRSD() generates second order design with qualitative factors along with their D-efficiency and G-efficiency. The function named RotatabilityQ() calculates a measure of rotatability (measure Q, 0 <= Q <= 1) given by Draper and Pukelshiem(1990) for given a design based on a second order model, (Norman R. Draper and Friedrich Pukelsheim(1990) <doi:10.1080/00401706.1990.10484635>. "Another look at rotatability").

This package provides the function mixor for fitting a mixed-effects ordinal and binary response models and associated methods for printing, summarizing, extracting estimated coefficients and variance-covariance matrix, and estimating contrasts for the fitted models.

This package provides a comprehensive set of tools for working with order statistics, including functions for simulating order statistics, censored samples (Type I and Type II), and record values from various continuous distributions. Additionally, it offers functions to compute moments (mean, variance, skewness, kurtosis) of order statistics for several continuous distributions. These tools assist researchers and statisticians in understanding and analyzing the properties of order statistics and related data. The methods and algorithms implemented in this package are based on several published works, including Ahsanullah et al (2013, ISBN:9789491216831), Arnold and Balakrishnan (2012, ISBN:1461236444), Harter and Balakrishnan (1996, ISBN:9780849394522), Balakrishnan and Sandhu (1995) <doi:10.1080/00031305.1995.10476150>, Genç (2012) <doi:10.1007/s00362-010-0320-y>, Makouei et al (2021) <doi:10.1016/j.cam.2021.113386> and Nagaraja (2013) <doi:10.1016/j.spl.2013.06.028>.

Models and predicts multiple output features in single random forest considering the linear relation among the output features, see details in Rahman et al (2017)<doi:10.1093/bioinformatics/btw765>.

Starting on the afternoon of July 17, 2026, Kevin Kruse fired off an astonishing array of BlueSky replies to an initial post of his featuring a certain government figure: <https://bsky.app/profile/did:plc:cnpe7qvcyjrhm6w7w7e4atur/post/3mqum4mxsuk2g>. This lasted a week and generated nearly seven hundred posts. A second wave started on August 12, 2026, with this post: <https://bsky.app/profile/kevinmkruse.bsky.social/post/3mstvbjpagc2a>. A third wave started on August 17, 2026, with <https://bsky.app/profile/kevinmkruse.bsky.social/post/3mtcpiw7gi22j>. A fourth wave started on August 27, 2026, with <https://bsky.app/profile/kevinmkruse.bsky.social/post/3mu3pugs2yk2f>. A fifth wave ran on August 30, 2026, beginning with <https://bsky.app/profile/kevinmkruse.bsky.social/post/3mudbzy5ksk25>. A sixth wave started September 4, 2026, with <https://bsky.app/profile/kevinmkruse.bsky.social/post/3muparqtdkk2w>. A seventh wave started September 12, 2026, with <https://bsky.app/profile/kevinmkruse.bsky.social/post/3mvdol6xu5k2s>. All of the over fourteen hundred posts from these series start with My man ... and make for excellent input to a fortunes'-like package. So this small package obliges and offers a random draw each time its myman() function is called. The overall package structure follows package fortunes', and atrrr was used to (bulk-)retrieve posts. Neither package is required to run this package to display random selections.

This package provides a variety of functions useful for data analysis, selection, manipulation, and graphics.

Facilitates creation and manipulation of metric graphs, such as street or river networks. Further facilitates operations and visualizations of data on metric graphs, and the creation of a large class of random fields and stochastic partial differential equations on such spaces. These random fields can be used for simulation, prediction and inference. In particular, linear mixed effects models including random field components can be fitted to data based on computationally efficient sparse matrix representations. Interfaces to the R packages INLA and inlabru are also provided, which facilitate working with Bayesian statistical models on metric graphs. The main references for the methods are Bolin, Simas and Wallin (2024) <doi:10.3150/23-BEJ1647>, Bolin, Kovacs, Kumar and Simas (2023) <doi:10.1090/mcom/3929> and Bolin, Simas and Wallin (2023) <doi:10.48550/arXiv.2304.03190> and <doi:10.48550/arXiv.2304.10372>.

Correlation coefficients for multivariate data, namely the squared correlation coefficient and the RV coefficient (multivariate generalization of the squared Pearson correlation coefficient). References include Mardia K.V., Kent J.T. and Bibby J.M. (1979). "Multivariate Analysis". ISBN: 978-0124712522. London: Academic Press.

The current version (0.3.0) streamlines the underlying code, adds a feasibility check to maicWt and maxessWt', extends exmWt.2ipd with options target.ipd and method for one-sided MAIC weighting between two IPDs, and adds wtTrtDiff for the weighted treatment-effect difference with a Wald confidence interval based on Section 5 of Glimm & Yau (2026). The second version (0.2.0) contains implementation for exact matching which is an alternative to propensity score matching (see Glimm & Yau (2026) <doi:10.1080/19466315.2025.2507378>). The initial version (0.1.2) contains a collection of easy-to-implement tools for checking whether a MAIC can be conducted, as well as an alternative way of calculating weights (see Glimm & Yau (2022) <doi:10.1002/pst.2210>.).

Estimation functions and diagnostic tools for mean length-based total mortality estimators based on Gedamke and Hoenig (2006) <doi:10.1577/T05-153.1>.

Conducts and visualizes propensity score analysis for multilevel, or clustered data. Bryer & Pruzek (2011) <doi:10.1080/00273171.2011.636693>.

Shiny web application to run meta-analyses. Essentially a graphical front-end to package meta for R. Can be useful as an educational tool, and for quickly analyzing and sharing meta-analyses. Provides output to quickly fill in GRADE (Grading of Recommendations, Assessment, Development and Evaluations) Summary-of-Findings tables. Importantly, it allows further processing of the results inside R, in case more specific analyses are needed.