search for: segments3d

Last week I asked about data ellipses with rgl:::ellipse3d() with lines showing the principal axes. (The goal is a visual demonstration of PCA as a rotation of variable space to component space.) I was trying, unsuccessfully, to use princomp() to generate the PCA axes and plot them using segments3d: > > PC <- princomp(trees) > > sdev <- PC$sdev # component standard deviations > > sd <- sqrt(diag(cov)) # variable standard deviations > > > > # vectors in variable space of principal components > > vec <- matrix(mu,3,3, byrow=TRUE) + diag...

...eemed overly complex. The particular case I have now is a little simpler. foo2 (points3d) and foo3 (spheres3d) are both functions from the rgl package that pass arguments such as color, alpha, etc. along to an rgl.material() command with lots of possible arguments. (The function may also call segments3d or lines3d, but these all have the same arguments as points3d.) However, to change the size of points you use "size"; to change the size of spheres you use "radius". Do I (1) add "radius" to the argument list, cluttering up the argument list for one particular special...

...ncomp(trees) sdev <- PC$sdev # component standard deviations sd <- sqrt(diag(cov)) # variable standard deviations # vectors in variable space of principal components vec <- matrix(mu,3,3, byrow=TRUE) + diag(sd) %*% PC$loadings for (j in 1:3) { mat <- rbind(mu, vec[j,]) segments3d(mat, col="red") } Can someone help? thanks, -Michael -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Street http://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3...

Hi all, I am trying to do a 3d plot where the x,y,z axes intersects with the origin (0,0,0) using the plot3d() funtion in the rgl package without success. I looked back at the past archives on this subject and someone suggested using djmrgl package. I searched and found it, installed it but when I try to load it I get the error ... Error in inDL(x, as.logical(local), as.logical(now), ...) :

Dear folks, I have created a plot on RGL device :   x = 1:6 y = seq(-12, 5, by=1) z = matrix(0, length(y), length(x)) z[13,3] = 1; z[13,4] = 1.011765 surface3d(x, y, t(z), col=rainbow(1000)) grid3d(c("x-", "y-", "z"))   Now I want to draw 2 lines along x=3 & x=4, over the surface (with different colour). Could somebody help me how to draw that?   Thanks,

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On 10-11-08 08:05 AM, wolfgang.polasek at gmail.com wrote: > Thanks, what if a I want to add to the plot a second column of observations, like a type h line from the x,y plane? > Wolfgang see segments3d. rgl is somewhat more do-it-yourself/construct-your-own-geometrical objects than the various 2D graphic systems (a bit more than base graphics, certainly more than lattice or ggplot) [forwarding back to r-help] -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.10 (GNU/Linux) Comment: Using Gn...

Dear all, I have a question about the plotting of the dynamical orbit by using rgl. My code is as follows: open3d() par3d(cex=2) bg3d("white") for(i in 1:out.length){ ind<-which plot3d(x[1:i],y[1:i],out[,"z"][1:i],xlim=c(-10,10),ylim=c(-10,10),zlim=c(-10,10),col="red",cex=2) Sys.sleep(0.0001) } ,where x,y,out[,"z"] are the position of the object in

...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, mat...

...;rgl-Ex.R' failed. >> The error most likely occurred in: >> >>> ### * 3dobjects >>> >>> flush(stderr()); flush(stdout()) >>> >>> ### Name: points3d >>> ### Title: add primitive set shape >>> ### Aliases: points3d lines3d segments3d triangles3d quads3d >>> ### Keywords: dynamic >>> >>> ### ** Examples >>> >>> # Show 12 random vertices in various ways. >>> >>> M <- matrix(rnorm(36), 3, 12, dimnames=list(c('x','y','z'), >> +...

...gcv) size = max(g)/20 * sphere.factor if (clear3d) clear3d() if (spheres) spheres3d(g, col=col.obj, radius=size, alpha=.5) else text3d(g, texts=rownames(g), col=col.obj, alpha=.5) aspect3d(aspect3d) for(i in 1:nrow(hl)) { segments3d(rbind(matrix(0, nc=3), hl[i,]*var.red), col=col.var) } text3d(hl*var.red, texts=rownames(hl), col=col.var) if(simple.axes) { axes3d(c('x', 'y', 'z')) } e...

Hi all, I am considering writing a 3D visualisation package for R. Purposes: - Visualise 3 dimension data - Manipulate the data or series or any other object with the mouse as in any 3D software (3D Studio Max, Lightwave...) - Eventually display more than just 3D objects. Justification: - It's always interesting to plot data before performing any analysis - Manipulating (rotating,

...box) { require(rgl) size = max(g)/20 * sphere.factor if (clear3d) clear3d() if (spheres) spheres3d(g, col=col.obj, radius=size, alpha=.5) else text3d(g, texts=rownames(g), col=col.obj, alpha=.5) aspect3d(aspect3d) for(i in 1:nrow(hl)) { segments3d(rbind(matrix(0, nc=3), hl[i,]*var.factor), col=col.var) } text3d(hl*var.factor, texts=rownames(hl), col=col.var) if(simple.axes) { axes3d(c('x', 'y', 'z')) title3d(xlab=paste('PC', lambda.ini, sep...