plotmixt function

Plot for 1- to 3-dimensional normal and t-mixture density functions

Plot for 1- to 3-dimensional normal and t-mixture density functions.

plotmixt(mus, sigmas, Sigmas, props, dfs, dist="normal", draw=TRUE,
   deriv.order=0, which.deriv.ind=1, binned=TRUE, ...)

Arguments

• mus: (stacked) matrix of mean vectors
• sigmas: vector of standard deviations (1-d)
• Sigmas: (stacked) matrix of variance matrices (2-d, 3-d)
• props: vector of mixing proportions
• dfs: vector of degrees of freedom
• dist: "normal" - normal mixture, "t" - t-mixture
• draw: flag to draw plot. Default is TRUE.
• deriv.order: derivative order
• which.deriv.ind: index of which partial derivative to plot
• binned: flag for binned estimation of contour levels. Default is TRUE.
• ...: other graphics parameters, see plot.kde

Returns

If draw=TRUE, the 1-d, 2-d plot is sent to graphics window, 3-d plot to graphics/RGL window. If draw=FALSE, then a kdde-like object is returned.

Examples

## bivariate 
mus <- rbind(c(0,0), c(-1,1))
Sigma <- matrix(c(1, 0.7, 0.7, 1), nr=2, nc=2) 
Sigmas <- rbind(Sigma, Sigma)
props <- c(1/2, 1/2)
plotmixt(mus=mus, Sigmas=Sigmas, props=props, display="filled.contour", lwd=1)

## trivariate 
mus <- rbind(c(0,0,0), c(-1,0.5,1.5))
Sigma <- matrix(c(1, 0.7, 0.7, 0.7, 1, 0.7, 0.7, 0.7, 1), nr=3, nc=3) 
Sigmas <- rbind(Sigma, Sigma)
props <- c(1/2, 1/2)
plotmixt(mus=mus, Sigmas=Sigmas, props=props, dfs=c(11,8), dist="t")