MixtureInf (version 1.1)

plotmix.norm: Histogram of the observations and the fitted density of a mixture of normals with unequal variance

Description

Plot the histogram of the observations and the fitted density of a mixture of normals with unequal variance.

Usage

plotmix.norm(x, theta, hist = 1, comp = TRUE, k = 20, h = 1, main = "", xlab = "Observations", ylab = "")

Arguments

x

data, can be either a vector or a matrix with the 1st column being the observed values and the 2nd column being the corresponding frequencies.

theta

parameter values, output of emtest.exp or pmle.exp, or a vector of parameter values, include mixing proportions, mixing means and mixing variances.

hist

style of histogram. hist = 0,1,2 are permitted. hist = 0 means no histogram, hist = 1,2 give two styles, default value: hist = 1.

comp

a parameter specifies the form of the component densities in the graph. comp = T switches on the component densities in the graph. If comp = F, no component densities will be included. default value: comp = T.

k

number of bins for the histogram, default value: k = 20.

h

a number between 0 and 1, specifying the height of the density function, default value: h = 1 (complete density function).

main

title of graph, default value: main = "".

xlab

label of x-axis, default value: xlab = "Observations".

ylab

label of y-axis, default value: ylab = "".

Value

Return the histogram of the observations and the plot of the fitted density.

See Also

emtest.norm, pmle.norm, rmix.norm

Examples

#generate a random sample from a 2 component normal mixture,
#plot the histogram of the observations and the density function with a set of 
#given parameters.
x <- rmix.norm(200,c(0.3,0.7),c(2,5),c(1,2))
theta <- c(0.25,0.75,1.8,5.2,1.1,2.1)
plotmix.norm(x,theta)