Comparison of the compound GB2 and kernel densities by group

Share:

Description

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

Usage

1
2
3
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

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.