R rig 生成逆高斯随机偏差

说明

生成逆高斯随机偏差。

用法

rig(n,mean,scale)

参数

n	所需的偏差数量。如果长度 > 1，则长度被视为所需的偏差数。
mean	平均值向量。
scale	尺度参数值向量(lambda，见下文)

细节

如果 x 是返回的向量，则 E(x) = mean 而 var(x) = scale*mean^3 。有关密度和分布函数，请参阅statmod 包。使用的算法是 Gentle (2003) 的算法 5.7，基于 Michael 等人。 (1976)。注意，这里的scale是GLM意义上的尺度参数，它是通常的‘lambda’参数的倒数。

值

逆高斯随机偏差的向量。

例子

require(mgcv)
set.seed(7)
## An inverse.gaussian GAM example, by modify `gamSim' output... 
dat <- gamSim(1,n=400,dist="normal",scale=1)
dat$f <- dat$f/4 ## true linear predictor 
Ey <- exp(dat$f);scale <- .5 ## mean and GLM scale parameter
## simulate inverse Gaussian response...
dat$y <- rig(Ey,mean=Ey,scale=.2)
big <- gam(y~ s(x0)+ s(x1)+s(x2)+s(x3),family=inverse.gaussian(link=log),
          data=dat,method="REML")
plot(big,pages=1)
gam.check(big)
summary(big)

作者

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

参考

Gentle, J.E. (2003) Random Number Generation and Monte Carlo Methods (2nd ed.) Springer.

Michael, J.R., W.R. Schucany & R.W. Hass (1976) Generating random variates using transformations with multiple roots. The American Statistician 30, 88-90.

https://www.maths.ed.ac.uk/~swood34/