similar to: Question about function DrawDensity3D {VecStatGraphs3D}

Hi, When plotting the density estimation with DrawDensity3D-function (in package VecStatGraphs3D) it often happens that the end product ie. layers are not intact , see figure in the link: https://www.dropbox.com/s/vzyaiu0vso8hjw2/20623868_10213333234554853_1371245353_n.png?dl=0 Can we somehow effect on that an get intact layers? -- Atte Tenkanen

Dear all, The documentation on the function open3d in the rgl package, reads: open3d opens a new rgl device, and sets the parameters as requested. I want the new rgl device (window) to be sized and positioned in a specific place on my screen. So, I try to set the "windowRect" parameter as follows: > library("rgl") > open3d(windowRect=c(100,100,500,500))

Hi, I am applying the following code to map pca loadings onto a 3d grid, my problem is this - the output only plots the spheres in the requested color (in this case "red") for the first argument. The sphere from the second argument appear as flat dark circles. Also the text3d() command only seems to work for a couple of the positions, with no text added in most cases. Could anyone offer

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

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

Hello, Using the rgl package, I can set the device window to any dimension (that I have tested): par3d(windowRect=c(1,1,700,700)) With windows I can't get the window to span from the top to the bottom of the monitor. In the following, no matter how large the ypinch value gets it stops, leaving about 2 inches of space at the bottom of my screen: windows(record=TRUE, ypinch=1100, xpinch=10,

Dear expeRts, I am a new user of rgl, below is my first trial to plot a simple function in 3d. I managed to put the axes in the right locations, but: (1) The xlab, ylab, and zlab arguments are ignored; how can I put in axes labels? (2) Since I removed the axes in persp3d() the viewport is too small; is it possible to keep the size of the viewport? (3) The box is not correctly drawn, there

Hi, I try to execute the seven lines of code below to plot a graph. But I am failing as the messages below show. Where am I going wrong? install.packages("rgl") library(rgl) y_hat = X%*%B_hat open3d(windowRect = c(100,100,900,900),family = "serif") color = rainbow(length(y_hat))[rank(y_hat)] plot3d(educ,exper,wage,col = color,type = "s",size = 0.5,xlim =

Hello all This is my first posting for some years. I am back using R again and must say I do like the language (regarding scripting, I also use matlab, perl, and bash). My question involves plotting a Pareto frontier in three dimensions. This is strictly a exercise in visualization, I make no attempt to extract the Pareto set (aka dominating subset) first. EXAMPLE PLOTS For some example

Hello everybody, I would like to solve the problem described below in R (if possible). I am very new to R and have therefore no idea what to even look for in the docu. Therefore, if someone could tell me if this is possible in R and where I might find the respective docu then I would already be very thankful. I have a large set (a few million) of points in 3D. I would like to create a

Dear all, I've been learning biplot (Gabriel, 1971) and some days ago I sent for this list a procedural function with invitation for a collaborative package. Jari Oksanen made some suggestions and I agree with all. So, I reworked the function under the object-oriented programming (OOP/S3). I think it is now a good frame for more resources. Below it is the function and a small script to

Hi, I'm trying to make sense of the image3d plots in misc3d. Unfortunately I can't find any authors on the help pages of the function to e-mail directly. In the following toy example, I'm puzzled as to why (1) the positions of the image plot on the axes include negative values when I think I've specified all to be 1:3 and (2) the logic that needs to be adhered to in order to

Hi everyone, I want to plot a 3D interpolation of the concentration of aquatic organisms. My goal would be to have the result represented as clouds with a density proportional to the abundance of organisms, so that I could fly (well, swim actually ;) ) through the scene and see the patches here and there. Basically, I want to do something like this: http://www.youtube.com/watch?v=27mo_Y-aU-c but

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

Hi, I guess I have a naive question. I use kde2d function in a standard way to estimate kernel densities of x and y (x and y are vectors) and plot them using image(). f1=kde2d(x,y) image(f1) But what if I want to see kernel estimates of three variables, x, y and z (a vector) plotted together ? Something in which x<->y is plotted and colored according to the corresponding value of z ?

Dears, I've been learning biplot (Gabriel, 1971) and I found the function 'biplot', inside of the package 'stats', useful but, a bit limited. So, I'm thinking to start a colaborative package to enhance this methods to other multivariate methods. In this way, I would like to start it, making public a new function (biplot.pca, still in development, but running) that make

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

Hello, I'm writing a pass using CodeExtractor to pull out a outmost loop as standalone function. But the 'Module Verifier' issues warning "function-local metadata used in wrong function". After looking through the CodeExtractor code, I found nothing about debugger intrinsics. Is this expected behavior? Thanks, yuanfang

Hi, I have two different imaging modalities (for the identification of areas of infarcted myocardium) that I need to compare regarding agreement and consistency. However, I don't think that methods like Cohen's Kappa, PCC, Bland-Altmann or ICC are sufficient here as there is not just a pairwise but also a spatial relationship between measured data points. For example if the results of the

I've written a very simple bit of code to plot a trajectory using rgl: x <- (c(0,-5.947,-11.496,-16.665,-21.474,-25.947,-30.116,-34.017,-37.684,-41.148,-44.435,-47.568,-50.567,-53.448,-56.223,-58.906)) y <- (c(0,33.729,65.198,94.514,121.784,147.151,170.797,192.920,213.717,233.360,252.003,269.774,286.781,303.117,318.858,334.071)) z <-