similar to: rgl: plot3d and ellipse3d

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

Dear all, I write to the community to know if other people expierenced the same problem. I want to save a rgl scene where I plotted a 3d cloud of spheres whose radius is set proportional to a given variable and coloured by another group variable and with 95% ellipsoids for distribution of any group. I had no problem at all for coding but I had two problem when saving the scene (in any format)

Dear list, I'm trying to visualise some ellipsoidal shapes in 3D. Their position, axes, and angular orientation can be arbitrary. I saw an ellipse3d function in rgl; however it is heavily oriented towards the statistical concept of ellipse of confidence, whilst I am just concerned with the geometrical object. Below is my current implementation. It is quite slow when creating many shapes, so

I have a 3xN matrix of parameters obtained from N regressions where the 3 parameters are jointly statistically significant. I would like to reproduce a 3D confidence ellipsoid projecting 2D ellipses onto the XY plane as in Figure 5.2 in this

Hello, I have created a 3d scatter plot using plot3d. Is it possible to export the resulting graph so that it can be placed in a pdf, powerpoint, or something of the sort and still be able to rotate it? Thanks [[alternative HTML version deleted]]

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

For a paper dealing with generalized ellipsoids, I want to illustrate in 3D an ellipsoid that is unbounded in one dimension, having the shape of an infinite cylinder along, say, z, but whose cross-section in (x,y) is an ellipse, say, given by the 2x2 matrix cov(x,y). I've looked at rgl:::cylinder3d, but don't see any way to make it accomplish this. Does anyone have any ideas? thx,

Dear RGL experts, I haven't been able to add greek letters to my rgl plot3d. I have tried "expression" with no success. Here is the interested bit: > library(rgl) > cb <- cube3d() > plot3d(cb,xlab=expression(alpha),ylab="",zlab="",box=FALSE,alpha=0.5) The expression(alpha) appears as "alpha", rather than as a greek symbol. I suspect greek

Dear Users: I wish to save 3d scatter plots I have generated using the plot3d command in the rgl package. Any image format such as .tiff would work. Can anyone tell me how to do this? Thank you [[alternative HTML version deleted]]

Buenas tardes. ¿Sabe alguno de vosotros cómo indicar la variable con la que identificar los puntos con plot3d? library(rgl) plot3d(Data$RTLML,Data$JD,Data$SB) Las muestras pertenecen a una de dos categorías, según la variable "Family", pero no sé cómo hacer que me las represente de diferente color. Gracias, Manuel . -- Dr Manuel

Dear R Users, i' wondering, if there is any possibility to render interactive stereoscopic 3d plots in R and use the functionality of those new 3d displays (e.g. with poarisation technology) like this one: > http://shop.fujitsu.com/de/displays/display-p23t-6p-fpr-3d Does anyone know if there are developments toward this feature? Or perhaps someone know how to implement such a

Dear List, I don't get scatter3d to color the sheres according to the '|' argument. library(car) scatter3d(prestige ~ income + education|type, data=Prestige) The spheres on my screen are all colored the same and they are not conditional on Prestige$type. On the other hand: Fit3d and Ellipse3d are colored according to the group argument. rgl_0.92.879 car_2.0-12 R version 2.15.0

*Summary*: The latest Windows binary version of rgl_081.708 from R-Forge has some problem that causes rgl.snapshot() to fail, at least on my system. Thereafter, *all* rgl 3D graphics are rendered without any text labels. The last problem remains even after (a) removing rgl and re-installing from CRAN, (b) re-starting with a fresh installation of all of R 2.7.2. *System*: Win XP Pro A vexing

I have a general style question about R coding. Suppose I'm writing a function (foo1) that calls other functions (foo2, foo3, ...) which have complicated argument lists (e.g. optim(), plot()), _and_ I may be calling several different functions in the body of foo1. Since foo2 and foo3 have different sets of arguments, I can't just use "..." ; I did write some code a while ago

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,

[Using R 2.7.2 on Windows XP] After re-building our heplots package, I've begun to get the following error from sessionInfo(), even though it passes R CMD check and builds without errors: > sessionInfo() Error in x$Priority : $ operator is invalid for atomic vectors In addition: Warning message: In FUN(c("MASS", "heplots", "car", "rgl",

Dear all, I have a question about the 'sdev' value returned by the princomp function (which does principal components analysis). On the help page for princomp it says 'sdev' is 'the standard deviations of the principal components'. However, when I calculate the principal components for the USArrests data set, I don't find this to be the case: Here is how I

Brian, I just tried your version rgl_0.70-2 with R-2.5.0 on FreeBSD 7.0-CURRENT (i386) and it works! Thank you very much for this competently and very fast help. Next I will try if it compiles on amd64, too ... Rainer Prof Brian Ripley schrieb: > This looks to me like a problem in your OpenGL, I am afraid. > > I've made available a (completely unofficial) revised tarball at

On p. 48 of "Statistics Complements" to the 3rd MASS edition, http://www.stats.ox.ac.uk/pub/MASS3/VR3stat.pdf I read that the orthogonal rotations of Z Lambda^-1 remain uncorrelated, where Z is the PC and Lambda is the diag matrix of singular values. However, the example below that text is > A <- loadings(ir.pca) %*% diag(ir.pca$sdev) If ir.pca$sdev are the singular values,