dMX: Evaluate the d (pdf) and p (cdf) functions from GAMLSS...

Description Usage Arguments Value Author(s) Examples

View source: R/gamlssMX-10-11-07.R

Description

The functions dMX() and pMX() can be used to evaluated the pdf (p function) and the cdf (p function) receptively from a gamlss.family mixture. The functions getpdfMX() and getpdfNP() can be used to evaluate the fitted d function at a specified observation and therefore for plotting the fitted distribution of a fitted model at this observation.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
dMX(y, mu = list(mu1 = 1, mu2 = 5), sigma = list(sigma1 = 1, sigma2 = 1), 
       nu = list(nu1 = 1, nu2 = 1), tau = list(tau1 = 1, tau2 = 1), 
       pi = list(pi1 = 0.2, pi2 = 0.8), family = list(fam1 = "NO", fam2 = "NO"), 
       log = FALSE, ...)
pMX(q, mu = list(mu1 = 1, mu2 = 5), sigma = list(sigma1 = 1, sigma2 = 1), 
       nu = list(nu1 = 1, nu2 = 1), tau = list(tau1 = 1, tau2 = 1), 
       pi = list(pi1 = 0.2, pi2 = 0.8), family = list(fam1 = "NO", fam2 = "NO"), 
       log = FALSE, ...)
getpdfMX(object = NULL, observation = 1)
getpdfNP(object = NULL, observation = 1)

Arguments

y,q

vector of quantiles

mu

a vector of mu's

sigma

a vector of sigma's

nu

a vector of nu's

tau

a vector of tau's

pi

a vector of pi's

family

a vector of GAMLSS family's

log

whether the log of the function or not

object

a fitted gamlssMX object

observation

the observation number in which we want to plot the fitted mixture

...

for extra arguments

Value

Returns values or pdf or cdf.

Author(s)

Mikis Stasinopoulos

Examples

1
2
3
4
5
6
fyNO<-dMX(y=seq(0,3,.01), mu=list(1.253, 0.1876), sigma=list(exp(-0.6665 ), exp(-2.573 )),
                  pi=list(0.4079609, 0.5920391 ), family=list("NO","NO") )
plot(fyNO~seq(0,3,.01), type="l")                  
FyNO<-pMX(q=seq(0,3,.01), mu=list(1.253, 0.1876), sigma=list(exp(-0.6665 ), exp(-2.573 )),
                  pi=list(0.4079609, 0.5920391 ), family=list("NO","NO") )
plot(FyNO~seq(0,3,.01), type="l")  

gamlss.mx documentation built on May 29, 2017, 12:21 p.m.

Search within the gamlss.mx package
Search all R packages, documentation and source code