Summarizing EM mixture model fits
summary.mixEM.Rdsummary method for class mixEM.
Usage
# S3 method for class 'mixEM'
summary(object, digits=6, ...)Arguments
- object
an object of class
mixEMsuch as a result of a call tonormalmixEM- digits
Significant digits for printing values
- ...
further arguments passed to
printmethod.
Details
summary.mixEM prints parameter estimates for
each component of a fitted mixture model.
The estimates printed vary with the type of model.
Value
The function summary.mixEM prints the final loglikelihood
value at the solution as well as a matrix of values for each component
that could include:
- lambda
The estimated mixing weights
- mu
The estimated mean parameters
- sigma
The estimated standard deviations
- theta
The estimated multinomial parameters
- beta
The estimated regression parameters
Examples
data(faithful)
attach(faithful)
#> The following objects are masked from faithful (pos = 10):
#>
#> eruptions, waiting
#> The following objects are masked from faithful (pos = 11):
#>
#> eruptions, waiting
#> The following objects are masked from faithful (pos = 12):
#>
#> eruptions, waiting
set.seed(100)
out <- normalmixEM(waiting, mu=c(50,80), sigma=c(5,5), lambda=c(.5,.5))
#> number of iterations= 15
summary(out)
#> summary of normalmixEM object:
#> comp 1 comp 2
#> lambda 0.360885 0.639115
#> mu 54.614805 80.091037
#> sigma 5.871184 5.867761
#> loglik at estimate: -1034.002