當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文


R lme.lmList LME 適合 lmList 對象


R語言 lme.lmList 位於 nlme 包(package)。

說明

如果 random 中定義的隨機效應名稱是 lmList 對象係數名稱的子集,則獲得隨機效應協方差矩陣的初始估計(覆蓋 random 中給出的任何值)。用於獲取 fixed 的調用序列中的 formula(fixed)data 參數作為 fixeddata 參數以及函數調用中的任何其他附加參數傳遞給 lme.formula 。有關該函數的說明,請參閱 lme.formula 的文檔。

用法

## S3 method for class 'lmList'
lme(fixed, data, random, correlation, weights, subset, method,
    na.action, control, contrasts, keep.data)

參數

fixed

繼承自類 "lmList." 的對象,表示符合通用模型的 lm 列表。

data

包含此參數是為了與通用函數保持一致。在此方法函數中它被忽略。

random

不帶條件表達式的可選單側線性公式,或具有 formula 屬性的 pdMat 對象。此方法函數不允許進行多級分組。默認為由 formula(fixed) 右側組成的公式。

correlation

說明組內相關結構的可選 corStruct 對象。有關可用 corStruct 類的說明,請參閱 corClasses 的文檔。默認為 NULL ,對應於沒有組內相關性。

weights

可選的 varFunc 對象或說明組內異方差結構的單邊公式。如果作為公式給出,則它將用作 varFixed 的參數,對應於固定方差權重。有關可用 varFunc 類的說明,請參閱 varClasses 的文檔。默認為 NULL ,對應於同方差組內誤差。

subset

一個可選表達式,指示應在擬合中使用的 data 行的子集。這可以是邏輯向量,或者指示要包括哪些觀察編號的數值向量,或者要包括的行名稱的字符向量。默認情況下包括所有觀察結果。

method

一個字符串。如果"REML" 則通過最大化受限對數似然來擬合模型。如果"ML" 則對數似然最大化。默認為 "REML"

na.action

一個函數,指示當數據包含 NA 時應該發生什麽。默認操作 ( na.fail ) 會導致 lme 打印錯誤消息並在存在任何不完整的觀察結果時終止。

control

估計算法的控製值列表,用於替換函數 lmeControl 返回的默認值。默認為空列表。

contrasts

可選列表。請參閱 model.matrix.defaultcontrasts.arg

keep.data

邏輯:data 參數(如果提供和 DataFrame )是否應該保存為模型對象的一部分?

代表線性混合效應模型擬合的 lme 類的對象。 printplotsummary 等通用函數具有顯示擬合結果的方法。有關配合的組成部分,請參閱lmeObject。函數 residcoeffittedfixed.effectsrandom.effects 可用於提取其某些組件。

例子

fm1 <- lmList(Orthodont)
fm2 <- lme(fm1)
summary(fm1)
summary(fm2)

作者

José Pinheiro and Douglas Bates bates@stat.wisc.edu

參考

The computational methods follow the general framework of Lindstrom and Bates (1988). The model formulation is described in Laird and Ware (1982). The variance-covariance parametrizations are described in Pinheiro and Bates (1996). The different correlation structures available for the correlation argument are described in Box, Jenkins and Reinse (1994), Littel et al (1996), and Venables and Ripley, (2002). The use of variance functions for linear and nonlinear mixed effects models is presented in detail in Davidian and Giltinan (1995).

Box, G.E.P., Jenkins, G.M., and Reinsel G.C. (1994) "Time Series Analysis: Forecasting and Control", 3rd Edition, Holden-Day.

Davidian, M. and Giltinan, D.M. (1995) "Nonlinear Mixed Effects Models for Repeated Measurement Data", Chapman and Hall.

Laird, N.M. and Ware, J.H. (1982) "Random-Effects Models for Longitudinal Data", Biometrics, 38, 963-974.

Lindstrom, M.J. and Bates, D.M. (1988) "Newton-Raphson and EM Algorithms for Linear Mixed-Effects Models for Repeated-Measures Data", Journal of the American Statistical Association, 83, 1014-1022.

Littel, R.C., Milliken, G.A., Stroup, W.W., and Wolfinger, R.D. (1996) "SAS Systems for Mixed Models", SAS Institute.

Pinheiro, J.C. and Bates., D.M. (1996) "Unconstrained Parametrizations for Variance-Covariance Matrices", Statistics and Computing, 6, 289-296.

Venables, W.N. and Ripley, B.D. (2002) "Modern Applied Statistics with S", 4th Edition, Springer-Verlag.

也可以看看

lme , lmList , lmeObject

相關用法


注:本文由純淨天空篩選整理自R-devel大神的英文原創作品 LME fit from lmList Object。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。