dgammaDifSaddle: Saddlepoint density for difference in means, \bar{x} -...
In bdsegal/gammaDist: Distribution functions and MLE for gamma random variables

This function calculates the saddlepoint density for difference in means, \bar{x} - \bar{y} of gamma random variables under the null of equal distributions.

1	dgammaDifSaddle(z, nx, ny, alpha, lambda)

`z`	scalar or vector of quantiles
`nx`	sample size of first group
`ny`	sample size of second group
`alpha`	common parameter value
`lambda`	common parameter value

nx <- 60
ny <- 60
alpha <- 1
lambda <- 4
z <- seq(-0.2, 0.2, 0.001)
plot(z, dgammaDifSaddle(z = z, nx, ny, alpha, lambda),
     type = "l", ylab = "f(q)")