当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R varPower 功率方差函数


R语言 varPower 位于 nlme 包(package)。

说明

该函数是varPower类的构造函数,表示幂方差函数结构。令 表示方差协变量, 表示在 处计算的方差函数,幂方差函数定义为 ,其中 是方差函数系数。当存在分组因子时,每个因子级别使用不同的

用法

varPower(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 ,对应于没有固定系数。

表示幂方差函数结构的 varPower 对象,也继承自类 varFunc

例子

vf1 <- varPower(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.

也可以看看

varWeights.varFunc , coef.varPower

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Power Variance Function。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。