Structural multivariate-univariate linear mixed model solver for estimation of multiple random effects with unknown variance-covariance structures (e.g., heterogeneous and unstructured) and known covariance among levels of random effects (e.g., pedigree and genomic relationship matrices) (Covarrubias-Pazaran, 2016 <doi:10.1371/journal.pone.0156744>; Maier et al., 2015 <doi:10.1016/j.ajhg.2014.12.006>; Jensen et al., 1997). REML estimates can be obtained using the Direct-Inversion Newton-Raphson and Direct-Inversion Average Information algorithms for the problems r x r (r being the number of records) or using the Henderson-based average information algorithm for the problem c x c (c being the number of coefficients to estimate). Spatial models can also be fitted using the two-dimensional spline functionality available.

Integrates the 13C nuclear magnetic resonance spectra using different integration ranges. Output depends on the method chosen. For the Molecular Mixing Model, a measurement of the fitting quality is given by its R-factor. For more details see: <doi:10.5281/zenodo.10137768>.

Formulas for calculating sound velocity, water pressure, depth, 
    density, absorption and sonar equations.

Simulates the cultural evolution of quantitative traits of bird song. 'SongEvo' is an individual- (agent-) based model. 'SongEvo' is spatially-explicit and can be parameterized with, and tested against, measured song data. Functions are available for model implementation, sensitivity analyses, parameter optimization, model validation, and hypothesis testing.

Estimate the abundance of cell clones from the
	     distribution of lengths of DNA fragments (as created by
	     sonication, whence `sonicLength').  The algorithm in
	     "Estimating abundances of retroviral insertion sites from
	     DNA fragment length data" by Berry CC, Gillet NA, Melamed
	     A, Gormley N, Bangham CR, Bushman FD. Bioinformatics;
	     2012 Mar 15;28(6):755-62 is implemented.  The
	     experimental setting and estimation details are described
	     in detail there. Briefly, integration of new DNA in a
	     host genome (due to retroviral infection or gene therapy)
	     can be tracked using DNA sequencing, potentially allowing
	     characterization of the abundance of individual cell
	     clones bearing distinct integration sites. The locations
	     of integration sites can be determined by fragmenting the
	     host DNA (via sonication or fragmentase), breaking the
	     newly integrated DNA at a known sequence, amplifying the
	     fragments containing both host and integrated DNA,
	     sequencing those amplicons, then mapping the host
	     sequences to positions on the reference genome. The
	     relative number of fragments containing a given position
	     in the host genome estimates the relative abundance of
	     cells hosting the corresponding integration site, but
	     that number is not available and the count of amplicons
	     per fragment varies widely.  However, the expected number
	     of distinct fragment lengths is a function of the
	     abundance of cells hosting an integration site at a given
	     position and a certain nuisance parameter. The algorithm
	     implicitly estimates that function to estimate the
	     relative abundance.