count_1_31: A set of functions to plot gamlss.family distributions

Description Usage Arguments Details Value Note Author(s) References See Also Examples

Description

Those functions are used in the distribution book of gamlss, see Rigby et. al 2019.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
binom_1_31(family = BI, mu = c(0.1, 0.5, 0.7), bd = NULL, miny = 0,
           maxy = 20, cex.y.axis = 1.2, cex.all = 1.5)
           
binom_2_33(family = BB, mu = c(0.1, 0.5, 0.8), sigma = c(0.5, 1, 2), 
           bd = NULL, miny = 0, maxy = 10, cex.y.axis = 1.5, 
           cex.all = 1.5)  
           
binom_3_33(family = ZIBB, mu = c(0.1, 0.5, 0.8), sigma = c(0.5, 1, 2),
           nu = c(0.01, 0.3), bd = NULL, miny = 0, maxy = 10, 
           cex.y.axis = 1.5, cex.all = 1.5, cols = c("darkgray", "black"), 
           spacing = 0.3)

contR_2_12(family = "NO", mu = c(0, -1, 1), sigma = c(1, 0.5, 2), 
          cols = c(1, 2, 3), ltype = c(1, 2, 3), maxy = 7, 
          no.points = 201, y.axis.lim = 1.1)
          
contR_3_11(family = "PE", mu = 0, sigma = 1, nu = c(1, 2, 3), 
          cols = c(1, 2, 3), maxy = 7, no.points = 201, 
          ltype = c(1, 2, 3), y.axis.lim = 1.1)          

contR_4_13(family = "SEP3", mu = 0, sigma = 1, nu = c(0.5, 1, 2), 
         tau = c(1, 2, 5), cols = c(1, 2, 3), maxy = 7, 
         no.points = 201, ltype = c(1, 2, 3), y.axis.lim = 1.1)
         
contRplus_2_11(family = GA, mu = 1, sigma = c(0.1, 0.6, 1), 
          cols = c(1, 2, 3), maxy = 4, no.points = 201, 
          y.axis.lim = 1.1, ltype = c(1, 2, 3))         

contRplus_3_13(family = "BCCG", mu = 1, sigma = c(0.15, 0.2, 0.5), 
          nu = c(-2, 0, 4), cols = c(1, 2, 3), maxy = 4, 
          ltype = c(1, 2, 3), no.points = 201, y.axis.lim = 1.1)
          
contRplus_4_33(family = BCT, mu = 1, sigma = c(0.15, 0.2, 0.5), 
          nu = c(-4, 0, 2), tau = c(100, 5, 1), cols = c(1, 2, 3), 
          maxy = 4, ltype = c(1, 2, 3), no.points = 201, 
          y.axis.lim = 1.1)

contR01_2_13(family = "BE", mu = c(0.2, 0.5, 0.8), 
            sigma = c(0.2, 0.5, 0.8), cols = c(1, 2, 3), 
            ltype = c(1, 2, 3), maxy = 7, no.points = 201, 
            y.axis.lim = 1.1, maxYlim = 10, cex.axis = 1.2)

contR01_4_33(family = GB1, mu = c(0.5), sigma = c(0.2, 0.5, 0.7),
            nu = c(1, 2, 5), tau = c(0.5, 1, 2), 
            cols = c(1, 2, 3, 4, 5, 6), maxy = 0.999, 
            ltype = c(1, 2, 3), no.points = 201, 
            y.axis.lim = 1.1, maxYlim = 10)

count_1_31(family = PO, mu = c(1, 2, 5), miny = 0, maxy = 10, 
           cex.y.axis = 1.2, cex.all = 1.5)
           
count_1_22(family = PO, mu = c(1, 2, 5, 10), miny = 0, 
           maxy = 20, cex.y.axis = 1.2, cex.all = 1.5) 
           
count_2_32(family = NBI, mu = c(0.5, 1, 5), sigma = c(0.1, 2), 
           miny = 0, maxy = 10, cex.y.axis = 1.5, cex.all = 1.5)

count_2_32R(family = NBI, mu = c(1, 2), sigma = c(0.1, 1, 2), 
           miny = 0, maxy = 10, cex.y.axis = 1.5, cex.all = 1.5) 
           
count_2_33(family = NBI, mu = c(0.1, 1, 2), sigma = c(0.5, 1, 2), 
          miny = 0, maxy = 10, cex.y.axis = 1.5, cex.all = 1.5) 
          
count_3_32(family = SICHEL, mu = c(1, 5, 10), sigma = c(0.5, 1), 
           nu = c(-0.5, 0.5), miny = 0, maxy = 10, cex.y.axis = 1.5, 
           cex.all = 1.5, cols = c("darkgray", "black"), spacing = 0.2)          
           
count_3_33(family = SICHEL, mu = c(1, 5, 10), sigma = c(0.5, 1, 2), 
           nu = c(-0.5, 0.5, 1), miny = 0, maxy = 10, cex.y.axis = 1.5,
           cex.all = 1.5, cols = c("darkgray", "black"), spacing = 0.3)

Arguments

family

a gamlss family distribution

mu

the mu parameter values

sigma

The sigma parameter values

nu

the nu parameter values

tau

the tau parameter values

bd

the binomial denominator

miny

minimal value for the y axis

maxy

maximal value for the y axis

cex.y.axis

the size of the letters in y axis

cex.all

the overall size

cols

colours

spacing

spacing between plots

ltype

The type of lines used

no.points

the number of points in the curve

y.axis.lim

the maximum value for the y axis

cex.axis

the size of the letters in the axis

maxYlim

the maximum permissible value for Y

Details

Th function plot different types of continuous and discrete distributions: i) contR: continuous distribution defined on minus infinity to plus infinity, ii) contRplus: continuous distribution defined from zero to plus infinity, iii) contR01: continuous distribution defined from zero to 1, iv) bimom binomial type discrete distributions, v) count count type discrete distributions.

The first number after the first underline in the name of the function indicates the number of parameters in the distribution. The two numbers after the second underline indicate how may rows and columns are in the plot.

Value

The result is a plot

Note

more notes

Author(s)

Mikis Stasinopoulos, Robert Rigby, Gillian Heller, Fernada De Bastiani

References

Rigby, R. A. and Stasinopoulos D. M. (2005). Generalized additive models for location, scale and shape,(with discussion), Appl. Statist., 54, part 3, pp 507-554.

Rigby R.A., Stasinopoulos D. M., Heller G. and De Bastiani F., (2019) Distributions for modelling location scale and shape: Using GAMLSS in R, Chapman and Hall/CRC (in press).

Stasinopoulos D. M. Rigby R.A. (2007) Generalized additive models for location scale and shape (GAMLSS) in R. Journal of Statistical Software, Vol. 23, Issue 7, Dec 2007, http://www.jstatsoft.org/v23/i07.

Stasinopoulos D. M., Rigby R.A., Heller G., Voudouris V., and De Bastiani F., (2017) Flexible Regression and Smoothing: Using GAMLSS in R, Chapman and Hall/CRC.

See Also

gamlss.family

Examples

1

Stan125/gamlss.dist documentation built on May 12, 2019, 7:38 a.m.