CompoundAuxDensPlot: Comparison of the compound GB2 and kernel densities by group

In GB2: Generalized Beta Distribution of the Second Kind: Properties, Likelihood, Estimation

Comparison of the compound GB2 and kernel densities by group

Description

Function dplot.cavgb2 produces a plot in which the compound and kernel (Epanechnikov) densities are plotted by group.

Usage

dplot.cavgb2(group, x, shape1, scale, shape2, shape3, pl0, pl, w=rep(1,length(x)), 
xmax = max(x)*(2/3), ymax=2e-05, decomp="r", choicecol=1:length(levels(group)), 
xlab="") 

Arguments

group	numeric; a factor variable giving the group membership of each sampled unit.
x	numeric; can be a vector. The value(s) at which the density is calculated, used for the kernel estimate only. x is positive.
shape1, scale, shape2, shape3	numeric; positive parameters of the GB2 distribution. On the plot they are denotes as a, b, p, q and pl0 respectively.
pl0	numeric; a vector of initial proportions defining the number of components and the weight of each component density in the decomposition. Sums to one.
pl	numeric; a vector of fitted proportions (output of pkl.cavgb2). Sums to one. If pl is equal to pl0, we obtain the GB2 distribution.
w	numeric; weights.
xmax	numeric; scale on the horizontal axis. By default is equal to max(x)*(2/3).
ymax	numeric; scale on the vertical axis. By default is equal to 2e-05.
decomp	string; specifying if the decomposition of the GB2 is done with respect to the right tail ("r") or the left tail ("l") of the distribution. By default, decomp = "r" - right tail decomposition.
choicecol	numeric vector of length the number of groups; defines the color with which the density curves will be plotted for each group.
xlab	string; label for x. The default is " ".

Details

The legend is placed interactively.

Value

dplot.cavgb2 plots a graph with two curves - the GB2 density, the compound GB2 per group and the corresponding kernel estimate.

Author(s)

Monique Graf and Desislava Nedyalkova

GB2 documentation built on June 22, 2022, 9:07 a.m.