On 28/08/2015 4:00 PM, Atte Tenkanen wrote:
> Hi,
>
> DrawDensity3D-function in package VecStatGraphs3D utilizes
> rgl.points-function {rgl}:
You shouldn't be using rgl.points. Use points3d.
The pointsize is a material property (see ?material3d), and isn't a
vector, so you only get one size per call to points3d.
Duncan Murdoch
>
> function (vectors, Div = 40, Layers = 3, DrawAxes = FALSE)
> {
> open3d(windowRect = c(100, 100, 800, 800))
> bg3d("white")
> Cx = vectors[, 1]
> Cy = vectors[, 2]
> Cz = vectors[, 3]
> Cr <- kde3d(x = Cx, y = Cy, z = Cz, n = Div)
> th <- seq(min(Cr$d), max(Cr$d), len = Layers + 2)
> ramp <- colorRamp(c("white", "yellow", "red"))
> colo <- rgb(ramp(seq(0, 1, length = Layers)), maxColorValue = 255)
> al <- seq(0.1, 0.6, len = Layers)
> module = sqrt(Cx * Cx + Cy * Cy + Cz * Cz)
> spheres3d(0, 0, 0, radius = max(module), color = "black",
> front = "line", back = "line", lwd = 1, smooth = TRUE,
> lit = TRUE, line_antialias = FALSE, alpha = 0.2)
> x <- c(0, max(module), 0, 0)
> y <- c(0, 0, max(module), 0)
> z <- c(0, 0, 0, max(module))
> labels <- c("", "X", "Y", "Z")
> i <- c(1, 2, 1, 3, 1, 4)
> text3d(x, y, z, labels, adj = 0.8, cex = 1.5, font = 2, color =
> "black")
> segments3d(x[i], y[i], z[i], lwd = 3)
> rgl.points(x = Cx, y = Cy, z = Cz, size = 3, color = "black")
> contour3d(Cr$d, level = th[c(-1, -(Layers + 2))], x = Cr$x,
> y = Cr$y, z = Cr$z, alpha = al, color = colo, add = TRUE,
> engine = "rgl", fill = TRUE, smooth = 2, material = "shiny")
> if (DrawAxes == TRUE) {
> axes3d()
> }
> }
>
> Is it somehow possible to define the sizes of the points all separately?
>
> I tried by adding ”Psize” to function arguments and changing
>
> rgl.points(x = Cx, y = Cy, z = Cz, size = Psize, color = "black”),
>
> then giving individual point size to each point but this does not work.
>
> This does’t work either:
>
> for(i in 1:length(Cx))
> {
> rgl.points(x=Cx[i], y=Cz[i], z=Cz[i], size=PSize[i], col= Colors[i])
> }
>
> Atte Tenkanen
>
> ______________________________________________
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.
>
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
