R rTweedie 生成 Tweedie 随机偏差

说明

为 1 到 2 之间的幂生成 Tweedie 随机偏差。

用法

rTweedie(mu,p=1.5,phi=1)

参数

mu	要生成的偏差的期望值向量。为 mu 的每个元素生成一个偏差。
p	偏差的方差与其平均值成正比，即 mu 的幂 p 。 p 必须介于 1 和 2 之间。1 类似于泊松(如果 phi=1 则完全是泊松)，2 是伽玛。
phi	尺度参数。偏差的方差由 phi*mu^p 给出。

细节

1<p<2 的 Tweedie 随机变量是 N gamma 随机变量的总和，其中 N 具有泊松分布，平均值为 mu^(2-p)/((2-p)*phi) 。求和的 Gamma 随机变量具有形状参数 (2-p)/(p-1) 和尺度参数 phi*(p-1)*mu^(p-1)(请注意，此尺度参数与具有 Gamma 误差的 GLM 的尺度参数不同)。

这是 tweedie 包中的 rtweedie 的受限制但速度更快的版本。

值

随机向量偏离 Tweedie 分布，期望值向量 mu ，方差向量 phi*mu^p 。

例子

 library(mgcv)
 f2 <- function(x) 0.2 * x^11 * (10 * (1 - x))^6 + 10 *
            (10 * x)^3 * (1 - x)^10
 n <- 300
 x <- runif(n)
 mu <- exp(f2(x)/3+.1);x <- x*10 - 4
 y <- rTweedie(mu,p=1.5,phi=1.3)
 b <- gam(y~s(x,k=20),family=Tweedie(p=1.5))
 b
 plot(b) 

作者

Simon N. Wood simon.wood@r-project.org

参考

Peter K Dunn (2009). tweedie: Tweedie exponential family models. R package version 2.0.2. https://cran.r-project.org/package=tweedie

也可以看看

ldTweedie , Tweedie