Fits smoothing spline regression models using scalable algorithms designed for large samples. Seven marginal spline types are supported: linear, cubic, different cubic, cubic periodic, cubic thin-plate, ordinal, and nominal. Random effects and parametric effects are also supported. Response can be Gaussian or non-Gaussian: Binomial, Poisson, Gamma, Inverse Gaussian, or Negative Binomial.

Easy-to-use, efficient, flexible and scalable statistical tools.
  Package bigstatsr provides and uses Filebacked Big Matrices via memory-mapping.
  It provides for instance matrix operations, Principal Component Analysis,
  sparse linear supervised models, utility functions and more
  <doi:10.1093/bioinformatics/bty185>.

Selecting linear and generalized linear models for large data sets
    using modified stepwise procedure and modern selection criteria (like
    modifications of Bayesian Information Criterion). Selection can be 
    performed on data which exceed RAM capacity. Bogdan et al., (2004)
    <doi:10.1534/genetics.103.021683>.

Fits Cox model via stochastic gradient descent. This implementation avoids computational instability of the standard Cox Model when dealing large datasets. Furthermore, it scales up with large datasets that do not fit the memory. It also handles large sparse datasets using proximal stochastic gradient descent algorithm. For more details about the method, please see Aliasghar Tarkhan and Noah Simon (2020) <arXiv:2003.00116v2>.

Extend the bigmemory package with 'table', 'tapply', and 'split'
    support for 'big.matrix' objects. The functions may also be used with native R
    matrices for improving speed and memory-efficiency.

For studying recurrent disease and death with competing
    risks, comparisons based on the well-known cumulative incidence function
    can be confounded by different prevalence rates of the competing events.
    Alternatively, comparisons of the conditional distribution of the survival
    time given the failure event type are more relevant for investigating the
    prognosis of different patterns of recurrence disease. This package implements
    a nonparametric estimator for the conditional cumulative incidence function
    and a nonparametric conditional bivariate cumulative incidence function for the
    bivariate gap times proposed in Huang et al. (2016) <doi:10.1111/biom.12494>.