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

This function calculates the saddlepoint cumulative distribution for the difference in means, \bar{x} - \bar{y} of gamma random variables under the null of equal distributions. Only valid for \bar{x} - \bar{y} \ne 0.

1	pgammaDifSaddle(q, nx, ny, alpha, lambda, lower.tail = TRUE)

`q`	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
q <- seq(-0.2, 0.2, 0.001)
plot(q, pgammaDifSaddle(q = q, nx, ny, alpha, lambda),
     type = "l", ylab = "F(q)")