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


R varExp 指數方差函數


R語言 varExp 位於 nlme 包(package)。

說明

該函數是varExp 類的構造函數,表示指數方差函數結構。令 表示方差協變量, 表示在 處計算的方差函數,指數方差函數定義為 ,其中 是方差函數係數。當存在分組因子時,每個因子級別使用不同的

用法

varExp(value, form, fixed)

參數

value

帶有方差函數係數的可選數值向量或數值列表。 Value 的長度必須為 1,除非在 form 中指定了分組因子。如果 value 的長度大於 1,則它必須具有將其元素標識為 form 中定義的分組因子級別的名稱。如果 form 中存在分組因子並且 value 的長度為 1,則其值將分配給所有分組級別。默認值為 numeric(0) ,這會導致在初始化 object 時將適當長度的零向量分配給係數(對應於常數方差等於 1)。

form

~ v~ v | g 形式的可選單邊公式,指定方差協變量 v 以及可選的係數分組因子 g。方差協變量必須計算為數值向量,並且可能涉及使用 "." 的表達式,表示擬合模型對象,可以從中提取擬合值 ( fitted(.) ) 和殘差 ( resid(.) )(這允許方差協變量為在目標函數優化期間更新)。當 form 中存在分組因子時,其每個級別使用不同的係數值。可以同時指定多個分組變量,並用 * 運算符分隔,如 ~ v | g1 * g2 * g3 中所示。在這種情況下,每個分組變量的水平將粘貼在一起,並使用生成的因子對觀測值進行分組。默認為~ fitted(.),表示由擬合模型對象的擬合值給出的方差協變量,無分組因子。

fixed

可選的數值向量或數值列表,指定方差函數中的部分或全部係數應固定的值。如果在 form 中指定了分組因子,則 fixed 必須具有標識要固定哪些係數的名稱。 fixed 中包含的係數在目標函數優化期間不允許發生變化。默認為 NULL ,對應於沒有固定係數。

表示指數方差函數結構的 varExp 對象,也繼承自類 varFunc

例子

vf1 <- varExp(0.2, form = ~age|Sex)

作者

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.

也可以看看

varClasses , varWeights.varFunc , coef.varExp

相關用法


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