I haven't looked at the source so I don't know exactly what is going on, but I
think I have a work around. While running your example I noticed that ellipse()
does not just add the ellipse to the plot produced by plot(), it replots the
figure. However, just running ellipse() without plot() generates an error
"Error in if (coorgeo == "acomp") { : argument is of length zero" so ellipse
needs the plot environment produced by plot(). Moving the pdf() file works on
my Windows machine:
> plot(winters.acomp, main="Winters Creek", cex=0.5)
> pdf("winters-pdf.pdf")
> ellipses(mean=mn, var=vr, r=r, steps=72, thinRatio=NULL, aspanel=FALSE,
+ col='red', lwd=2)
> dev.off()
-------------------------------------
David L Carlson
Department of Anthropology
Texas A&M University
College Station, TX 77840-4352
-----Original Message-----
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Rich Shepard
Sent: Tuesday, October 14, 2014 4:21 PM
To: r-help@r-project.org
Subject: [R] Ternary Plots Do Not Display Ellipses in PDF
A rather strange situation here and I've not found the source of the
problem.
The point is to print a ternary plot matrix of compositional data with
ellipses enclosing 95% of the variance in each plot. The ellipses display on
the monitor, dev = x11cairo (see attached winters-x11cairo.pdf), but not when
sent directly to a file, dev = pdf (see attached winters-pdf.pdf).
Here's winters.acomp:
structure(c(0.0666666666666667, 0.0612244897959184, 0.0434782608695652,
0.043956043956044, 0.05, 0.0161290322580645, 0.6, 0.571428571428571,
0.623188405797101, 0.593406593406593, 0.433333333333333, 0.629032258064516,
0.0666666666666667, 0.0612244897959184, 0.101449275362319, 0.0659340659340659,
0.0666666666666667, 0.032258064516129, 0.244444444444444, 0.26530612244898,
0.217391304347826, 0.263736263736264, 0.366666666666667, 0.290322580645161,
0.0222222222222222, 0.0408163265306122, 0.0144927536231884, 0.032967032967033,
0.0833333333333333, 0.032258064516129), .Dim = c(6L, 5L), .Dimnames = list(
NULL, c("filter", "gather", "graze", "predate", "shred")), class = "acomp")
And this is the command sequence:
> library(compositions)
> plot(winters.acomp, main="Winters Creek", cex=0.5)
> r <- sqrt(qchisq(p=0.95, df=4))
> mn <- mean(winters.acomp)
> vr <- var(winters.acomp)
> plot(winters.acomp, main="Winters Creek", cex=0.5)
> ellipses(mean=mn, var=vr, r=r, steps=72, thinRatio=NULL, aspanel=FALSE,
col='red', lwd=2)
# monitor plot window is manually closed.
> pdf("winters-pdf.pdf")
> plot(winters.acomp, main="Winters Creek", cex=0.5)
> ellipses(mean=mn, var=vr, r=r, steps=72, thinRatio=NULL, aspanel=FALSE,
col='red', lwd=2)
> dev.off()
What am I not seeing here that causes the different outputs?
Rich
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
