R predict.lme 來自 lme 對象的預測

說明

i 級別的預測是通過將總體預測(僅基於固定效應估計)和隨機效應對小於或等於 i 的分組級別的預測的估計貢獻相加而獲得的。結果值估計 i 級別的最佳線性無偏預測 (BLUP)。如果 newdata 中存在未包含在原始分組因子中的組值，則對於大於或等於未知組出現的級別的級別，相應的預測將設置為 NA。

用法

## S3 method for class 'lme'
predict(object, newdata, level = Q, asList = FALSE,
        na.action = na.fail, ...)

參數

object	繼承自類 "lme" 的對象，表示擬合的線性混合效果模型。
newdata	用於獲得預測的可選數據幀。固定效應模型和隨機效應模型中使用的所有變量以及分組因子都必須出現在 DataFrame 中。如果缺失，則返回擬合值。
level	一個可選的整數向量，給出用於獲得預測的分組級別。級別值從最外層到最內層分組遞增，級別零對應於總體預測。默認為最高或最內層的分組。
asList	可選的邏輯值。如果 TRUE 並在 level 中給出單個值，則返回的對象是一個列表，其中預測按組劃分；否則，根據 level 的長度，返回值是向量或數據幀。
na.action	一個函數，指示當 newdata 包含 NA 時應該發生什麽。默認操作 ( na.fail ) 會導致函數打印一條錯誤消息，並在存在任何不完整的觀察結果時終止。
...	該泛型的某些方法需要額外的參數。此方法中沒有使用任何內容。

值

如果在 level 中指定了單個分組級別，則返回值是包含按組分割的預測的列表 ( asList = TRUE ) 或包含預測的向量 ( asList = FALSE )；否則，當在 level 中指定多個分組級別時，返回的對象是一個數據幀，其中的列由不同級別的預測和分組因子給出。

例子

fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject)
newOrth <- data.frame(Sex = c("Male","Male","Female","Female","Male","Male"),
                      age = c(15, 20, 10, 12, 2, 4),
                      Subject = c("M01","M01","F30","F30","M04","M04"))
## The 'Orthodont' data has *no* 'F30', so predict  NA  at level 1 :
predict(fm1, newOrth, level = 0:1)

作者

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

也可以看看

lme , fitted.lme