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