swh:1:snp:ffdd0a7d2c8ea15ad41d45b3b178f668bd942287
Tip revision: 713027f11850cd2b542fcd4f21527c74ee4042d8 authored by Derek Young on 30 January 2013, 00:00:00 UTC
version 0.4.6
version 0.4.6
Tip revision: 713027f
summary.mixEM.R
summary.mixEM <- function(object, digits=6, ...){
x <- object
o <- switch(x$ft,
"multmixEM" = rbind(x$lambda, t(x$theta)),
"normalmixEM" = rbind(x$lambda, x$mu, x$sigma),
"repnormmixEM" = rbind(x$lambda, x$mu, x$sigma),
"regmixEM" = rbind(x$lambda, x$sigma, x$beta),
"regmixEM.lambda" = rbind(x$lambda, x$sigma, x$beta),
"regmixEM.mixed" = rbind(x$lambda, x$sigma, x$beta),
"regmixEM.loc" = rbind(x$sigma, x$beta),
"regmixEM.chgpt" = rbind(x$lambda, x$sigma),
"logisregmixEM" = rbind(x$lambda, x$beta),
"poisregmixEM" = rbind(x$lambda, x$beta),
"mvnormalmixEM" = rbind(x$lambda, matrix(unlist(x$mu), byrow=TRUE, nrow=length(x$lambda))),
stop("Unknown mixEM object of type ", x$ft))
colnames(o) <- paste("comp",1:ncol(o))
rownames(o) <- switch(x$ft,
"multmixEM" = c("lambda", paste("theta", 1:ncol(x$theta), sep="")),
"normalmixEM" = c("lambda", "mu", "sigma"),
"repnormmixEM" = c("lambda", "mu", "sigma"),
"regmixEM" = c("lambda", "sigma", paste("beta", 1:nrow(x$beta), sep="")),
"regmixEM.lambda" = c("lambda", "sigma", paste("beta", 1:nrow(x$beta), sep="")),
"regmixEM.mixed" = c("lambda", "sigma", paste("beta", 1:nrow(x$beta), sep="")),
"regmixEM.loc" = c("sigma", paste("beta", 1:nrow(x$beta), sep="")),
"regmixEM.chgpt" = c("lambda", "sigma"),
"logisregmixEM" = c("lambda", paste("beta", 1:nrow(x$beta), sep="")),
"poisregmixEM" = c("lambda", paste("beta", 1:nrow(x$beta), sep="")),
"mvnormalmixEM" = c("lambda", paste("mu", 1:length(x$mu), sep="")))
cat("summary of", x$ft, "object:\n")
print(o, digits=digits)
cat("loglik at estimate: ", x$loglik, "\n")
}