nlme
位于 nlme
包(package)。 说明
该通用函数符合 Lindstrom 和 Bates (1990) 中说明的公式中的非线性混合效应模型,但允许嵌套随机效应。组内误差允许相关和/或具有不等方差。
用法
nlme(model, data, fixed, random, groups, start, correlation, weights,
subset, method, na.action, naPattern, control, verbose)
## S3 method for class 'formula'
nlme(model, data, fixed, random, groups, start, correlation, weights,
subset, method, na.action, naPattern, control, verbose)
参数
model |
非线性模型公式,响应位于 |
data |
包含 |
fixed |
|
random |
可选地,以下任何一个:(i)形式为 |
groups |
|
start |
可选的数值向量,或固定效应和随机效应的初始估计列表。如果声明为数值向量,它将在内部转换为具有单个组件 |
correlation |
说明组内相关结构的可选 |
weights |
可选的 |
subset |
一个可选表达式,指示应在拟合中使用的 |
method |
一个字符串。如果 |
na.action |
一个函数,指示当数据包含 |
naPattern |
表达式或公式对象,指定哪些返回值被视为缺失。 |
control |
估计算法的控制值列表,用于替换函数 |
verbose |
可选的逻辑值。如果打印 |
值
代表非线性混合效应模型拟合的 nlme
类对象。 print
、 plot
和 summary
等通用函数具有显示拟合结果的方法。有关配合的组成部分,请参阅nlmeObject
。函数 resid
、 coef
、 fitted
、 fixed.effects
和 random.effects
可用于提取其某些组件。
注意
该函数不会在内部进行任何缩放:当响应缩放时优化效果最好,因此其方差约为 1。
例子
fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc),
data = Loblolly,
fixed = Asym + R0 + lrc ~ 1,
random = Asym ~ 1,
start = c(Asym = 103, R0 = -8.5, lrc = -3.3))
summary(fm1)
fm2 <- update(fm1, random = pdDiag(Asym + lrc ~ 1))
summary(fm2)
作者
José Pinheiro and Douglas Bates bates@stat.wisc.edu
参考
The model formulation and computational methods are described in Lindstrom, M.J. and Bates, D.M. (1990). The variance-covariance parametrizations are described in Pinheiro and Bates (1996).
Lindstrom, M.J. and Bates, D.M. (1990) "Nonlinear Mixed Effects Models for Repeated Measures Data", Biometrics, 46, 673-687.
Pinheiro, J.C. and Bates., D.M. (1996) "Unconstrained Parametrizations for Variance-Covariance Matrices", Statistics and Computing, 6, 289-296.
For the different correlation structures, variance functions and links,
see ‘References’ in lme
.
也可以看看
nlmeControl
, nlme.nlsList
, nlmeObject
, nlsList
, nlmeStruct
, pdClasses
, reStruct
, varFunc
, corClasses
, varClasses
相关用法
- R nlmeStruct 非线性混合效应结构
- R nlme.nlsList 从 nlsList 对象进行 NLME 拟合
- R nlmeControl nlme Fit 的控制值
- R nlsList.selfStart nlsList 通过 selfStart 函数拟合
- R nlsList 具有通用模型的 nls 对象列表
- R needUpdate.modelStruct 检查 modelStruct 对象是否需要更新
- R needUpdate 检查是否需要更新
- R Pixel X 射线像素强度随时间的变化
- R corARMA ARMA(p,q) 相关结构
- R getGroupsFormula 提取分组公式
- R corRatio 有理二次相关结构
- R logLik.glsStruct glsStruct 对象的对数似然
- R intervals.lmList lmList 系数的置信区间
- R corLin 线性相关结构
- R plot.augPred 绘制 augPred 对象
- R print.varFunc 打印 varFunc 对象
- R recalc 重新计算压缩线性模型对象
- R Variogram.corSpher 计算 corSpher 对象的半变异函数
- R getGroups.lme 提取 lme 对象组
- R predict.nlme 来自 nlme 对象的预测
- R corSymm 一般相关结构
- R qqnorm.gls gls 对象残差的正态图
- R pdCompSymm 具有复合对称结构的正定矩阵
- R [.pdMat 下标 pdMat 对象
- R pdConstruct.pdBlocked 构造 pdBlocked 对象
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Nonlinear Mixed-Effects Models。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。