Skip to contents

Core functions

fastlmm()
Fast Linear Mixed Model with 1 Random Effect
fastglmm()
Fit Generalized Linear Mixed Model via PQL
fastglmm.nb()
Fit Negative Binomial Mixed Model via PQL

Example Data

PsychAD
Gene expession of PTPRG in microglia in 60k cells

Generics

anova(<fastlmm>)
ANOVA Tables
coef(<fastlmm>)
Extract Model Coefficients
cooks.distance(<fastlmm>)
Cook's Distance Metric
deviance(<fastlmm>)
Model Deviance
ddf()
Denominator Degrees of Freedom
df.residual(<fastlmm>)
Residual Degrees-of-Freedom
edf()
Effective Degrees-of-Freedom of Model Fit
extractAIC(<fastlmm>)
Extract AIC from a Fitted Model
family(<fastlmm>) family(<fastglmm>) family(<glmmPQL>)
Extract Model Family
fitted(<fastlmm>) fitted(<fastglmm>)
Extract Model Fitted Values
fixef(<fastlmm>)
Extract Fixed Effects
hatvalues(<fastlmm>)
Return Diagonals of Hat Matrix
linearHypothesis()
Test Linear Hypothesis
logLik(<fastlmm>)
Extract Log-Likelihood
model.frame(<fastlmm>)
Extracting the Model Frame from a Fit
nobs(<fastlmm>)
Extract the Number of Observations from a Fit
plot(<fastlmm>)
Diagnostic Plots for Model Fits
predict(<fastlmm>) predict(<fastglmm>)
Model Predictions
ranef(<fastlmm>)
Extract the modes of the random effect
residuals(<fastlmm>)
Extract Model Residuals
rstudent(<fastlmm>)
Extract Studentized Residuals
sigma(<fastlmm>)
Extract Residual Standard Deviation
simulate(<fastlmm>)
Simulate Responses
summary(<fastlmm>)
Object Summaries and Hypothesis Testing
terms(<fastlmm>)
Model Terms
varpart()
Variance Partitioning Analysis
vcov(<fastlmm>)
Calculate Variance-Covariance Matrix for a Fitted Model Object

Other

as.fastlmm()
Convert list to fastlmm class
getFamilyString()
Convert GLM family to string
indicator_decomp()
Spectral decomposition of factor indicator matrix
preprocess_indicator()
Create sparse indicator matrix
fastlmm.fit()
Fitter Function for Linear Mixed Model
nb_theta()
Estimate theta of the Negative Binomial
refitModel()
Refit fastlmm model with new delta
dispersion()
Overdispersion parameter phi for quasi-likelihood
varianceTerms()
Array of variance component estimates
negative.binomial()
Family function for Negative Binomial GLMs

Pure R Implementations (for intuition only)

fastlmm_R()
Fit linear mixed model using SVD of covariance
fastglmm_R()
Fit generalized linear mixed model via PQL
fastglmm_R.nb()
Fit negative binomial mixed model via PQL