augPred
位於 nlme
包(package)。 說明
預測值是在 primary
的指定值處獲得的。如果 object
具有分組結構(即 getGroups(object)
不是 NULL
),則為每個組獲取預測值。如果level
具有多個元素,則針對max(level)
分組因子的每一級別獲得預測。如果在預測模型中使用除 primary
之外的其他協變量,則使用它們的平均值(數值協變量)或最頻繁值(分類協變量)來獲取預測值。原始觀察結果也包含在返回的對象中。
用法
augPred(object, primary, minimum, maximum, length.out, ...)
## S3 method for class 'lme'
augPred(object, primary = NULL,
minimum = min(primary), maximum = max(primary),
length.out = 51, level = Q, ...)
參數
object |
使用 |
primary |
一個可選的單邊公式,指定用於生成增強預測的主要協變量。默認情況下,如果可以從用於生成 |
minimum |
主要協變量的可選下限。默認為 |
maximum |
主要協變量的可選上限。默認為 |
length.out |
一個可選整數,其中包含用於評估預測的主要協變量值的數量。默認為 51。 |
level |
指定所需預測級別的可選整數向量。級別從最外層到最內層分組遞增,級別 0 代表總體(固定效應)預測。默認為最內層。 |
... |
一些泛型方法可能需要額外的參數。 |
值
具有四列的 DataFrame ,分別表示主協變量的值、組(如果 object
沒有分組結構,則所有元素將為 1
)、預測值或觀察值以及類型第三列中的值:original
表示觀察值,predicted
(單個或無分組因子)或 predict.groupVar
(多級分組),其中 groupVar
替換為實際分組變量名稱 (fixed
用於人口預測)。返回的對象繼承自類 "augPred"
。
注意
這個函數是通用的;可以編寫方法函數來處理特定類的對象。已經具有此函數方法的類包括: gls
、 lme
和 lmList
。
例子
fm1 <- lme(Orthodont, random = ~1)
augPred(fm1, length.out = 2, level = c(0,1))
作者
José Pinheiro and Douglas Bates bates@stat.wisc.edu
參考
Pinheiro, J. C. and Bates, D. M. (2000), Mixed-Effects Models in S and S-PLUS, Springer, New York.
也可以看看
相關用法
- R as.matrix.corStruct corStruct 對象的矩陣
- R anova.lme 比較擬合對象的可能性
- R asTable 將 groupedData 轉換為矩陣
- R as.matrix.reStruct reStruct 對象的矩陣
- R anova.gls 比較擬合對象的可能性
- R asOneFormula 組合一組對象的公式
- R as.matrix.pdMat pdMat 對象的矩陣
- R allCoef 從一組對象中提取係數
- 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 nlmeStruct 非線性混合效應結構
- R predict.nlme 來自 nlme 對象的預測
- R corSymm 一般相關結構
- R qqnorm.gls gls 對象殘差的正態圖
- R pdCompSymm 具有複合對稱結構的正定矩陣
注:本文由純淨天空篩選整理自R-devel大神的英文原創作品 Augmented Predictions。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。