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


R varConstPower 常數加冪方差函數


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

說明

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

用法

varConstPower(const, power, form, fixed)

參數

const , power

可選的數值向量或數值列表,分別包含常數項和冪項的係數。除非在 form 中指定了分組因子,否則兩個參數的長度都必須為 1。如果任一參數的長度大於 1,則它必須具有將其元素標識為 form 中定義的分組因子級別的名稱。如果 form 中存在分組因子並且參數長度為 1,則其值將分配給所有分組級別。 const 隻允許使用正值。默認值為 numeric(0) ,這會導致在初始化 object 時將適當長度的零向量分配給係數(對應於常數方差等於 1)。

form

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

fixed

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

表示常數加冪方差函數結構的 varConstPower 對象,也繼承自類 varFunc

例子

vf1 <- varConstPower(1.2, 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.varConstPower

相關用法


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