MixtureDistribution: Mixture Distribution functions

Description Usage Arguments Value Author(s) See Also Examples

Description

pdf, cdf, inverse cdf, and random deviates of a user defined mixture distribution with equal weights.

Usage

1
2
3
4
dMixtureDistribution(x, functionList)
pMixtureDistribution(x, functionList)
qMixtureDistribution(p, functionList, xMin, xMax, nPoints=1000, logScale=FALSE)
rMixtureDistribution(n, functionList)

Arguments

x

The value at which the cdf or pdf are computed

functionList

A list of functions forming the mixture distribution. The functions must be a continuous CDF dependent on 'x'.

p

A vector of probabilities.

xMin

The minimum value for which probabilities and quantiles are calculated.

xMax

The maximum value for which probabilities and quantiles are calculated.

nPoints

The number of points at which the mixed CDF is calculated. Default=1000.

logScale

Are the CDF estimates to be performed on a logarithmic scale? Default=FALSE

n

The number of random deviates to compute

Value

A vector of densities, probabilities, quantiles, or random deviates from the mixture distribution.

Author(s)

Rob Carnell

See Also

pnorm

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
  qMixtureDistribution(c(0.5), list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"), 0, 4)
  pMixtureDistribution(2, list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  pMixtureDistribution(c(0.5,1,2,4,5,6), list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  x1 <- rMixtureDistribution(10000, list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  mean(x1)
  rMixtureDistribution(1, list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  rMixtureDistribution(3, list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  rMixtureDistribution(4, list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  integrate(dMixtureDistribution, 0, 4, functionList=list("pnorm(x, 1, 0.1)","pnorm(x, 2, 0.1)","pnorm(x, 3, 0.1)"))
  integrate(dMixtureDistribution, 0, Inf, functionList=list("pnorm(x, 1, 0.1)","plnorm(x, .2, 0.1)","pnorm(x, 3, 0.1)"))

MixtureDistribution documentation built on May 2, 2019, 5:25 p.m.