similar to: persp3d and rgl.viewpoint for rotating 3D plots

Dear R-users, I would like to place two 3D plots side-by-side in a rgl-setting. It would nice to have something like "par(mfrow=c(1,2))" for basic plots, or an array framework for wireframe(lattice) (see example below). I only managed to overlap two persp3d plots. My final idea would be to animate both surfaces using play3d(rgl). Thanks in advance for any help. Best, Carlo Giovanni

Dear R-users, I am going to try to be as clearer as possible, showing also an example. 1) I have a function (in my real case it's much more complex) 2) I use "optim" to minimize 3) I want to use as method L-BFGS-B for several reasons 4) I know I could use "trace=6" (in "control") in order to see "live" the procedure 5) I would like to see separately the

Dear R-users, I am working with some grid-search optimization for 13 values of an object function. At glance one may compute the object function for each possible combination of set of parameters, but in my case would not be feasible: taking for example 13 possible values for each parameters like, in logarithm scale, 10^seq(-3,3,.5) will lead to 13^13 combinations of results. As a second trial

Dear R users, I have been trying to draw the following 3d graphs: The solid region bounded above by the paraboloid z = 9 - x2 - y2 and below by the unit circle in the xy-plane. I wanted to visualize the solid region bounded by those two graphs. I could draw those two 3d graphs separately, but I could not get those two to work in one single plot. Here's my code for the paraboloid: x <-

Dear all, I have tried to plot a triangular matrix with the function persp3d(rgl). for example z=rbind(c(1,NA,NA,NA),c(5,3,NA,NA),c(4,2,9,NA),c(8,6,5,11)) x=1:4 y=1:4 persp3d(x,y,z,color="gray") The two extreme points are not plotted (value=1 and value=10). It seems because the half of the matrix have 'NA' and perp3d need planar quadrilateral face. So I decided to use the

Hi, I want to format my axis in my persp3d plot. With my data, which I attached I created a persp3d plot with the following code, which I summarized from different code snippets I found: library(rugarch)library(rgl)library(fGarch)fd <-as.data.frame(modelfit,which ='density')color <-rgb(85,141,85,maxColorValue=255)x <-seq(-0.2,0.2,length=100)y <-c(1:2318)f

G'day everyone, I've been trying to get an interactive OpenGL plot to work from within Python 2.6 using Rpy2 and the persp3d function. The problem is that the plot seems to freeze upon activation. All interactivity is lost even though it works fine when run from within the R console (version 2.11.1). Is this a known issue? Any ideas as to what might be causing this, or possibly how to fix

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 guys, I hit on a problem when I use rgl. Could you try to run the code here in this link and see why the first persp3d gives a red bounding box and the second shows black? http://rafb.net/p/g1i7ur33.html (sorry for not pasting the code here directly but my previous email got filtered by this list so I suspect my code looks like spam to the spam filter) I'm expecting black color to be

Dear R-users, I'm struggling in R in order to "squeeze" a matrix without using a for-loop. Although my case is a bit more complex, the following example should help you to understand what I would like to do, but without the slow for-loop. Thanks in advance, Carlo Giovanni Camarda A <- matrix(1:54, ncol=6) # my original matrix A.new <- matrix(nrow=3, ncol=6) # a new

Dear R-users, I would like to know whether there is a way to extract a diagonal of an inverse matrix without computing the inverse of the matrix itself. The size of my matrices are really huge and, also using sparse matrix, computing the inverse leads to storage problems and low speed. In other words, given a square matrix A, I aim to know diag(B), where B=solve(A), without computing solve(A).

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060728/957eb23c/attachment.pl

Dear R-friends, I use optim(par=c(mystartingpoints), fn=myloglikelihoodfunction, gr=NULL, method=c("L-BFGS-B"), ## I would like to do not use any bounds control=list(trace=6, ## just to see what it's going on maxit=c(20000)), ## to be sure the it doesn't stop reaching the max iterations

Hola a todos: Estoy intentando visualizar unos resultados de un programa en R de forma dinámica. Dichos resultados se pueden representar, básicamente, con un plano que gira sobre uno de sus ejes: para dar sensación de movimiento he utilizado la función rgl.clear(), ya que no encuentro los equivalentes en rgl de los parámetros add o new. El problema es que al limpiar la pantalla con rgl.clear la

Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function "par" and then struggling around with "axes" and labels etc. Luckily, I would always have "rows + columns = even number" and the same "ylim". The next one could be a

I don't understand what is happening! I have a (large) object sim1, an matrix list with dim c(101,101) where each element is an 3*3 matrix. I am subsetting that with a matrix coo, of dim c(100,2), of unique indices, but the resulting object has length 99, not 100 as expected. Code reproducing the problem follows: library(RandomFields) set.seed(123) sim0 <- GaussRF(x=seq(0, 100, by=1),

Hi list, I'm trying to limit the z axis plotted using persp3d {rgl} to a given range. I'm trying the following statement persp3d(g0,pa,D, xlim=range(g0),ylim=range(pa), zlim=range (zbounds),col="lightblue") where g0 and pa are both a range of 20 numbers from 0.05 to 1, and D is a 20 x 20 matrix with values ranging from 0 to 30000. I'm trying to set a cap on the Z

I'm using the command below on an open3d() object to create a shaded cube. Changes to myScalingFactor do not effect changes in the size of the cube. What is the correct approach? Mark shade3d(translate3d(cube3d(identityMatrix() * myScalingFactor),-6,1,-1),col="green", alpha = 0.2) -- Mark W. Kimpel MD ** Neuroinformatics ** Dept. of Psychiatry Indiana University School of Medicine

Hello, I'd like to draw a graphic using different color line segments to represent microarray data. It simply likes this: Data=cbind(c(1, 2, -1, 3, -2.2), c(1, 3, -1, -1.8, 4)) The first col No represents y'axis position (negative down and positive up). The second col No represents a length of segment (negative segment should be on left of y, and positive on right). Also the

Hello all, I use persp3d from the rgl-package to plot a sruface. The typical call is persp3d(x, y, z) With cooridinate-vectros x, y and a function-values matrix z. Now I have different z's, say z_1,...,z_n Question: Is it possible to generate an animation from a sequence of such calls, for different z's? I would like to see how the surface is beeing changed in the time. Thank you