CompoundAuxDensPlot {GB2}

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