Plot Saturation Curve for NBMM

Saturation curve for the effective number of independent measurements per donor for a negative binomial mixed model

Usage

plotSaturation(fit, method = c("counts", "measurements"))

Arguments

fit

model fit from fastlmm() or fastglmm()

method

use asymptotic values for "counts" or "measurements"

Details

See powerNBMM()

See also

meff() powerNBMM()

Examples

data(PsychAD)

# regression formula
form <- PTPRG ~ (1|SubID) + offset(log(libSize))

# fit NBMM in on PTPRG expression
fit <- fastglmm.nb(form, PsychAD)

# Saturation curve as mu increases
plotSaturation(fit, method="counts")


# Saturation curve as m increases
plotSaturation(fit, method="measurements")