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).

Cheers,
Roman



2010/1/18 Ben Bolker <bolker@ufl.edu>

>  I think it's easier than that ... This looks doable by setting
> par(lty=2) before you draw the plot, then setting lty.axis and lty.grid
> to override that setting ...
>
> >      data(dune.env)
> >      ord <- cca(dune ~ A1 + Moisture, dune.env)
> >      ordiplot3d(ord)
> > par(lty=2)
> >      ordiplot3d(ord,type="h")  ## oops, everything is dashed
> >      ordiplot3d(ord,type="h",lty.axis=1,lty.grid=1)
>
>  However, this also draws everything else --  the biplot arrows and
> coordinate axes that ordiplot3d() adds to the basic scatterplot3d plot
> -- with dashed lines as well, which you might not want.  To override
> this you would have to hack ordiplot3d in some way -- the easiest thing
> would be to add lty=1 to the pl$points3d() and arrows() calls, or to
> add par() calls before and after the scatterplot3d() call to set lty
> to 2 (or whatever) and then change it back to lty=1 again.
>
>  Ben Bolker
>
> Jari Oksanen wrote:
> > Roman,
> >
> >
> > On 18/01/2010 19:35, "Roman Luštrik" <roman.lustrik@gmail.com> wrote:
> >
> >> I'm trying to plot a CCA plot in 3D with ordiplot3d().
> >> Is there a way to plot vertical lines for sites (type="h") as dashed,
> >> instead of solid, which is the default?
> >>
> > I don't know a way to achieve this. The ordiplot3d() uses scatterplot3d()
> of
> > the scatterplot3d package, and passes all arguments to that underlying
> > function. Howevever, scatterplot3d() seems to have no option of setting
> > vertical segments in type = "h" and therefore ordiplot3d cannot set them.
> > The only way seems to be to edit the code of scatterplot3d and add an
> > argument for the purpose.
> >
> > Cheers, Jari Oksanen
> >
> >> Here is the code I use if of any help:
> >>
> >> cca.plot3d <- ordiplot3d(CCA, type="h", angle=40)
> >> text(cca.plot3d, "points", pos=c(4, 4, 2, 4, 3, 4, 2, 2, 4, 2))
> >> text(cca.plot3d, "arrows", cex=0.7)
> >> points(cca.plot3d, "points")
> >>
> >> Cheers,
> >> Roman
> >>
> >> [[alternative HTML version deleted]]
> >>
> >> _______________________________________________
> >> R-sig-ecology mailing list
> >> R-sig-ecology@r-project.org
> >> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
> >
> > _______________________________________________
> > R-sig-ecology mailing list
> > R-sig-ecology@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
>
>
> --
> Ben Bolker
> Associate professor, Biology Dep't, Univ. of Florida
> bolker@ufl.edu / people.biology.ufl.edu/bolker
> GPG key: people.biology.ufl.edu/bolker/benbolker-publickey.asc
>
>

	[[alternative HTML version deleted]]

