VGAM (version 1.0-2)

biplackettcop: Plackett's Bivariate Copula Family Function

Description

Estimate the association parameter of Plackett's bivariate distribution (copula) by maximum likelihood estimation.

Usage

biplackettcop(link = "loge", ioratio = NULL, imethod = 1, nsimEIM = 200)

Arguments

link

Link function applied to the (positive) odds ratio $psi$. See Links for more choices and information.

ioratio

Numeric. Optional initial value for $psi$. If a convergence failure occurs try assigning a value or a different value.

imethod, nsimEIM

See CommonVGAMffArguments.

Value

An object of class "vglmff" (see vglmff-class). The object is used by modelling functions such as vglm and vgam.

Details

The defining equation is $$\psi = H \times (1-y_1-y_2+H) / ((y_1-H) \times (y_2-H))$$ where $P(Y1 <= 2="" 4="" y1,="" y2="" <="y2)=" h(y1,y2)$="" is="" the="" cumulative="" distribution="" function.="" density="" function="" $h(y1,y2)="$" $$\psi="" [1="" +="" (\psi-1)(y_1="" y_2="" -="" y_1="" y_2)="" ]="" \left(="" ]^2="" \psi="" (\psi-1)="" \right)^{3="" 2}$$="" for="" $psi=""> 0$. Some writers call $psi$ the cross product ratio but it is called the odds ratio here. The support of the function is the unit square. The marginal distributions here are the standard uniform although it is commonly generalized to other distributions.

If $psi=1$ then $h(y1,y2) = y1*y2$, i.e., independence. As the odds ratio tends to infinity one has $y1=y2$. As the odds ratio tends to 0 one has $y2=1-y1$.

Fisher scoring is implemented using rbiplackcop. Convergence is often quite slow.

References

Plackett, R. L. (1965) A class of bivariate distributions. Journal of the American Statistical Association, 60, 516--522.

See Also

rbiplackcop, bifrankcop.

Examples

## Not run: 
# ymat <- rbiplackcop(n = 2000, oratio = exp(2))
# plot(ymat, col = "blue")
# fit <- vglm(ymat ~ 1, fam = biplackettcop, trace = TRUE)
# coef(fit, matrix = TRUE)
# Coef(fit)
# vcov(fit)
# head(fitted(fit))
# summary(fit)
# ## End(Not run)