用于逆高斯分布的R语言程序

Question

逆高斯分布的累积密度函数可写为： F（t;μ，λ）=Φ[√（（λ）/ t）（t/μ-1）] +〖exp〗^（2λ/μ）Φ[ - √（λ/ t）（t / μ+ 1）]，其中Φ表示标准正态分布的累积分布函数。请告诉我如何在R

中编写此函数

Answer 1

正态分布包含在R中。有关如何计算正态分布的分位数的详细信息，请参阅?qnorm。它还将向您展示其他有用的功能。

Answer 2

可从CRAN获得的statmod包包括用于逆高斯分布的累积密度函数的函数pinvgauss()，以及该分布的其他基本概率函数。随附文章：

Giner，G和Smyth，GK（2016）。 statmod: probability calculations for the inverse Gaussian distribution。 R Journal 8（1），339-351。