Learn R Programming

⚠️There's a newer version (3.3) of this package.Take me there.

saemix (version 3.2)

Stochastic Approximation Expectation Maximization (SAEM) Algorithm

Description

The 'saemix' package implements the Stochastic Approximation EM algorithm for parameter estimation in (non)linear mixed effects models. The SAEM algorithm (i) computes the maximum likelihood estimator of the population parameters, without any approximation of the model (linearisation, quadrature approximation,...), using the Stochastic Approximation Expectation Maximization (SAEM) algorithm, (ii) provides standard errors for the maximum likelihood estimator (iii) estimates the conditional modes, the conditional means and the conditional standard deviations of the individual parameters, using the Hastings-Metropolis algorithm (see Comets et al. (2017) ). Many applications of SAEM in agronomy, animal breeding and PKPD analysis have been published by members of the Monolix group. The full PDF documentation for the package including references about the algorithm and examples can be downloaded on the github of the IAME research institute for 'saemix': .

Copy Link

Version

Install

install.packages('saemix')

Monthly Downloads

556

Version

3.2

License

GPL (>= 2)

Maintainer

Emmanuelle Comets

Last Published

June 27th, 2023

Functions in saemix (3.2)

conddist.saemix

Estimate conditional mean and variance of individual parameters using the MCMC algorithm
coef.saemix

Extract coefficients from a saemix fit
dataGen.case

Bootstrap datasets
SaemixObject-class

Class "SaemixObject"
PD1.saemix

Data simulated according to an Emax response model, in SAEM format
SaemixRes-class

Class "SaemixRes"
SaemixModel-class

Class "SaemixModel"
backward.procedure

Backward procedure for joint selection of covariates and random effects
SaemixData-class

Class "SaemixData"
compare.saemix

Model comparison with information criteria (AIC, BIC).
fim.saemix

Computes the Fisher Information Matrix by linearisation
createSaemixObject

Create saemix objects with only data filled in
cow.saemix

Evolution of the weight of 560 cows, in SAEM format
fitted.saemix

Extract Model Predictions
discreteVPCTTE

VPC for time-to-event models
epilepsy.saemix

Epilepsy count data
default.saemix.plots

Wrapper functions to produce certain sets of default plots
[

Get/set methods for SaemixData object
discreteVPC

VPC for non Gaussian data models
forward.procedure

Backward procedure for joint selection of covariates and random effects
logLik

Extract likelihood from a saemixObject resulting from a call to saemix
llgq.saemix

Log-likelihood using Gaussian Quadrature
initialize-methods

Methods for Function initialize
llis.saemix

Log-likelihood using Importance Sampling
mydiag

Matrix diagonal
map.saemix

Estimates of the individual parameters (conditional mode)
knee.saemix

Knee pain data
plot.SaemixData

Plot of longitudinal data
plot,SaemixModel,SaemixData-method

Plot model predictions for a new dataset. If the dataset is large, only the first 20 subjects (id's) will be shown.
predict.SaemixModel

Predictions for a new dataset
print-methods

Methods for Function print
lung.saemix

NCCTG Lung Cancer Data, in SAEM format
plot-methods

Methods for Function plot
predict-methods

Methods for Function predict
npdeSaemix

Create an npdeObject from an saemixObject
plotDiscreteData

Plot non Gaussian data
oxboys.saemix

Heights of Boys in Oxford
plot,SaemixModel,ANY-method

Plot model predictions using an SaemixModel object
psi-methods

Functions to extract the individual estimates of the parameters and random effects
plot,SaemixObject,ANY-method

General plot function from SAEM
readSaemix-methods

Methods for Function read
replaceData

Replace the data element in a SaemixObject object
rapi.saemix

Rutgers Alcohol Problem Index
saemix.plot.data

Functions implementing each type of plot in SAEM
saemix.plot.select

Plots of the results obtained by SAEM
readSaemix,SaemixData-method

Create a longitudinal data structure from a file or a dataframe Helper function not intended to be called by the user
saemix.bootstrap

Bootstrap for saemix fits
saemix.internal

Internal saemix objects
saemix

Stochastic Approximation Expectation Maximization (SAEM) algorithm
resid.saemix

Extract Model Residuals
saemixPredictNewdata

Predictions for a new dataset
show-methods

Methods for Function show
saemix.plot.setoptions

Function setting the default options for the plots in SAEM
showall-methods

Methods for Function showall
saemixModel

Function to create a SaemixModel object
saemix.predict

Compute model predictions after an saemix fit
saemixControl

List of options for running the algorithm SAEM
saemixData

Function to create a SaemixData object
summary-methods

Methods for Function summary
simulateDiscreteSaemix

Perform simulations under the model for an saemixObject object defined by its log-likelihood
subset

Data subsetting
[,SaemixRes-method

Get/set methods for SaemixRes object
simulate.SaemixObject

Perform simulations under the model for an saemixObject object
[,SaemixModel-method

Get/set methods for SaemixModel object
[,SaemixObject-method

Get/set methods for SaemixObject object
step.saemix

Stepwise procedure for joint selection of covariates and random effects
stepwise.procedure

Stepwise procedure for joint selection of covariates and random effects
theo.saemix

Pharmacokinetics of theophylline
xbinning

Internal functions used to produce prediction intervals (from the npde package)
yield.saemix

Wheat yield in crops treated with fertiliser, in SAEM format
toenail.saemix

Toenail data
vcov

Extracts the Variance-Covariance Matrix for a Fitted Model Object
transform

Transform covariates
transformCatCov

Transform covariates
validate.names

Name validation (## )Helper function not intended to be called by the user)
validate.covariance.model

Validate the structure of the covariance model
transformContCov

Transform covariates
testnpde

Tests for normalised prediction distribution errors