我试图熟悉R中的mixtools包,但我遇到了一个奇怪的问题。该包应该能够“排队”绘图,以便绘图函数首先绘制一个绘图,然后绘制另一个绘图等,直到达到绘图输入中的绘图数量。
然而,这在我的R版本(Mac上的R 3.1.2 GUI 1.65 Mavericks版本(6833))中不起作用。
运行此代码(来自
http://www.inside-r.org/packages/cran/mixtools/docs/plot.mixEM
)给我一个情节。
data(faithful)
attach(faithful)
set.seed(100)
out <- normalmixEM(waiting, arbvar = FALSE, verb = TRUE,
epsilon = 1e-04)
plot(out, density = TRUE, w = 1.1)
但下面这张图,应该是两幅图,只给了我两个空白的石英窗:
##Fitting randomly generated data with a 2-component location mixture of bivariate normals.
x.1 <- rmvnorm(40, c(0, 0))
x.2 <- rmvnorm(60, c(3, 4))
X.1 <- rbind(x.1, x.2)
out.1 <- mvnormalmixEM(X.1, arbvar = FALSE, verb = TRUE,
epsilon = 1e-03)
plot(out.1, density = TRUE, alpha = c(0.01, 0.05, 0.10),
marginal = TRUE)
有人知道该怎么办吗?
谢谢