[R-sig-eco] ordiplot3d

Jari Oksanen jari.oksanen at oulu.fi
Tue Jan 19 09:47:32 CET 2010


On 19/01/10 10:17 AM, "Roman Luštrik" <roman.lustrik at gmail.com> wrote:

> Thank you Jari and Ben.
> I have tried something along the lines what Ben suggested but stopped and
> asked here if there's a better way, as I didn't want to invent the wheel
> (yet crude one).
> 
Roman,

I think the only sustainable alternative is to edit scatterplot3d(). This is
not very difficult. The following adds argument type.hplot for the line type
of vertical segments in type = "h" plots:

--- scatterplot3d.R.orig    2010-01-19 10:38:36.000000000 +0200
+++ scatterplot3d.R    2010-01-19 10:41:25.000000000 +0200
@@ -12,6 +12,7 @@
      font.axis = par("font.axis"), font.lab = par("font.lab"),
      lty.axis = par("lty"), lty.grid = par("lty"), lty.hide=NULL,
-     log = "", ...)
+     lty.hplot = par("lty"), log = "", ...)
      # log not yet implemented
+     ## lty.hplot sets the lty for line segments in type = "h"
 {
     ## Uwe Ligges <ligges at statistik.tu-dortmund.de>,
@@ -249,5 +250,5 @@
     if(type == "h") {
         z2 <- dat$y * yz.f + z.min
-        segments(x, z, x, z2, col = col, cex = cex.symbols, ...)
+        segments(x, z, x, z2, col = col, cex = cex.symbols, lty =
lty.hplot, ...)
         points(x, z, type = "p", col = col, pch = pch, cex = cex.symbols,
...)
     }

Cheers, Jari



More information about the R-sig-ecology mailing list