[R] Problem with mclust surfacePlot function

About this list Date view Thread view Subject view Author view Attachment view

From: Mike White (mikewhite.diu@tiscali.co.uk)
Date: Fri 07 May 2004 - 00:03:25 EST


Message-id: <004601c43372$e66d2f20$835de150@FSSFQCV7BGDVED>

I am trying to follow the mclust examples in "MCLUST: Software for Model
Based Clustering, Density Estimation and Disriminant Analysis" by Chris
Fraley and Adrian Raftery, but I cannot reproduce the density and
uncertainty surfaces for the Lansing Woods maples. I am using R 1.8.1 with
the code below. The same code works fine in S-Plus 6.2

Am I missing something or is this a bug?

Thanks
Mike White

library(mclust)
data(lansing) # R only
maples<-lansing[as.character(lansing[,"species"]) == "maple", -3]
maplesBIC <- EMclust(maples)
maplesModel <- summary(maplesBIC, maples)
plotMaples2 <- function(type, what, transformation)
{
 out <- do.call("surfacePlot", c(maplesModel, list(data=maples, type=type,
what=what,transformation=transformation)))
 invisible()
}
par(pty="s", mfrow=c(1,2))
plotMaples2(type="contour", what="density", transformation="log")
par(pty="s")
plotMaples2(type="contour", what="uncertainty", transformation = "log")

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


About this list Date view Thread view Subject view Author view Attachment view

This archive was generated by hypermail 2.1.3 : Mon 31 May 2004 - 23:05:08 EST