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。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。