Farlie-Gumbel-Morgenstern's Bivariate Distribution

Density, distribution function, and random generation for the (one parameter) bivariate Farlie-Gumbel-Morgenstern's distribution.

Usage

dbifgmcop(x1, x2, apar, log = FALSE)
pbifgmcop(q1, q2, apar)
rbifgmcop(n, apar)

Arguments

x1, x2, q1, q2

vector of quantiles.

n

number of observations. Same as in runif.

apar

the association parameter.

log

Logical. If TRUE then the logarithm is returned.

Value

dbifgmcop gives the density, pbifgmcop gives the distribution function, and rbifgmcop generates random deviates (a two-column matrix).

Author

T. W. Yee

Details

See bifgmcop, the VGAM family functions for estimating the parameter by maximum likelihood estimation, for the formula of the cumulative distribution function and other details.

See also

bifgmcop.

Examples

if (FALSE)  N <- 101; x <- seq(0.0, 1.0, len = N); apar <- 0.7
#> Error: object 'N' not found
ox <- expand.grid(x, x)
#> Error: object 'x' not found
zedd <- dbifgmcop(ox[, 1], ox[, 2], apar = apar)
#> Error: object 'apar' not found
contour(x, x, matrix(zedd, N, N), col = "blue")
#> Error: object 'x' not found
zedd <- pbifgmcop(ox[, 1], ox[, 2], apar = apar)
#> Error: object 'ox' not found
contour(x, x, matrix(zedd, N, N), col = "blue")
#> Error: object 'x' not found

plot(r <- rbifgmcop(n = 3000, apar = apar), col = "blue")
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'plot': object 'apar' not found
par(mfrow = c(1, 2))
hist(r[, 1])  # Should be uniform
#> Error: object 'r' not found
hist(r[, 2])  # Should be uniform
#> Error: object 'r' not found
 # \dontrun{}