fit.norm: Internal function

In rauschenberger/semisup: Semi-Supervised Mixture Model

Description

This function fits the semi-supervised Gaussian mixture model. It is called by fit.wrap.

Usage

fit.norm(y, z, it.em, epsilon)

Arguments

y	observations: numeric vector of length n
z	class labels: integer vector of length n, with entries 0, 1 and NA
it.em	(maximum) number of iterations in the EM algorithm: positive integer (defaults to 100)
epsilon	convergence criterion for the EM algorithm: non-negative numeric (defaults to 1e-04)

Value

This function returns the parameter estimates, the posterior probabilities, and the likelihood.

See Also

This is an internal function. The user functions are mixtura and scrutor.

Examples

# data simulation
n <- 100
z <- rep(0:1,each=n/2)
y <- rnorm(n=n,mean=2*z,sd=1)
z[(n/4):n] <- NA

# model fitting
fit.norm(y,z,it.em=100,epsilon=1e-04)

rauschenberger/semisup documentation built on May 11, 2020, 4 p.m.