similar to: export/import matrix

Hi Folks, I'm trying to get rgl.Sweave to produce plots with transparency. However, it just seems to produce opaque plots when pdf is the output type. Perhaps this is a known issue? I'll just use .png in the meantime, but wanted to see about this, as I didn't see it in the documentation (though it's possible I missed it). Thanks, Allie \documentclass{article} \title {rgl

Hi R users and experts, I am interested in learning more about the use of 3D plots. Specifically, I want to add points and lines to a surface plot. And get the axes and labels plotted also. Here is what I have tried with an example data set : library(rgl) vol2 <- 2*volcano # Exaggerate the relief library(reshape) mvol2 <- melt(vol2) str(mvol2) # First, persp and persp3d plots do not succeed

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 list I'm having trouble with inverting color gradation. As seen in the "rgl" example, library(rgl) example(rgl.surface) I understand that I can assign colors to heights for each point. col <- colorlut[ y-ylim[1]+1 ] # assign colors to heights for each point However, I am now trying to assign inverted colors. (top colored green, and bottom colored gray). Is there any

Running: R Version 2.3.1 (2006-06-01), debian stable I'm testing an application that uses rgl (installed with apt-get install r-cran-rgl) but fails, at least partly because that application uses a more current version. There is also a mysterious error msg from Xlib: > library(car) > library(rgl) > iris.mod <- lm(cbind(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) ~

Dear List, Some independent variable were missing in calculation using lm and glm (gaussian). (X= Y1+Y2+…..+Y16, Independent number: 16 variable) However, those variables did work well in cor(X, Y) respectively. str(dataframe) was also run to ensure that the variables were all numbers. Moreover, the missing variables were different in lm and glm. In lm, 3 factors were not taken into

Dear List, I tried to do cca based on species data and environmental variables (formula instead of community data). However, there was an error saying row sums must be >0. I searched the previous related messages but found few solutions. Please kindly help and thank you in advance. code This is vegan 1.17-3 Warning message: package 'vegan' was built under R version 2.10.1

Hello, I am using package MuMIn to calculate AIC for a full model with 10 explanatory variables. Thanks in advance in sharing your experience. Q1 In the AIC list of all models, each model is differentiated by model number. Please kindly advise if it is possible to find the corresponding explanatory variable(s) for the model number. Q2 error message I tried to display sub-model with only

hi, all I am looking at R package RGL to draw a colored mesh/surface plot like this one (from matlab). http://www.mathworks.com/access/helpdesk/help/techdoc/visualize/cbar.gif The key features I am looking for is surfaced with grid and color, but not the terrain-like gradient. but I didn't come even close to it after browsing through rgl help file. have anyone drawn something like this

In R, I'm imported a data frame of 2,321,123 by 4 called "dataF". I converted the data frame "dataF" to a matrix dataM <- as.matrix(dataF) Does R have an efficient routine to treat the special elements that contain "inf" in them. For example, can you separate the rows that have "inf" elements from the matrix into a separate matrix without

Dear all, I would like to put two perspective plots into one plot. The help page for ?persp shows how one can add points and lines but not another perspective plot. data(volcano) z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) ## Don't draw the grid lines : border = NA persp(x,

?rev On Jun 16, 2009, at 5:13 AM, Naoki Irie wrote: > Dear list > > > I'm having trouble with inverting color gradation. > > As seen in the "rgl" example, > > library(rgl) > example(rgl.surface) > > I understand that I can assign colors to heights for each point. > > col <- colorlut[ y-ylim[1]+1 ] # assign colors to heights for each >

Hi all, I need a little help for construct an state's map on R. The first problem is to get the data. I have a datafile of longitude and latitude in the follow format: trajectory latitude longtude T -22.045618 -51.287056 T -22.067078 -51.265888 T -22.067039 -51.207249 T -22.059690 -48.089695 T -22.075529 -48.074608 T -22.072460 -48.044472 T -22.062767 -48.298473 T -22.077349

Hello I have a table describing butterfly range traits. It is composed of three columns as below Species name range size (X) latitude of range midpoint (Y) There are 11 kinds of butterflies. Each has its range size, and the latitude of each range midpoint ranges from 9 to 19. I would like to have the average range size of every degree of latitude. For example, the average range

Hi, I'm trying to save heatmaps as PDFs. However, the PDF version of the heatmaps (Heatmap_CAFvsTNF_run2.pdf) is blurred when compared to its counterpart, which was saved manually by using the software "Grab" (Heatmap_CAFvsTNF_run2.tiff). -----R code-------- sample_output <- "stroma_run2" filename <-

Hi this is more of an OpenLDAP issue, but since it is mostly used with Samba, some experience would be helpful. We have a SLES9 PDC, running samba-3.0.14a, openldap 2.2.6 Sometimes, the OpenLDAP gets corrupted, no ideas why. It's a bad thing. We have 100+ computers, 300+ users, which shall get 350 computers and 800 users. What would you suggest in order to be reliable ? Reliability is

Among others, datam contains the columns: logconc, tm, dose, subj, bilirubin. None of these are factor variables. The following compartment models work (the first still has not converged after 100 interations): res1 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))* (exp(-exp(p2)*tm)-exp(-exp(p1)*tm))),start=list(fixed=c(5,-2,-0.1)), fixed=list(p1+p2+p3~1),control=list(maxIter=100),

On Thu, Dec 14, 2017 at 1:58 PM, Berend Hasselman <bhh at xs4all.nl> wrote: > >> On 14 Dec 2017, at 19:36, lily li <chocold12 at gmail.com> wrote: >> >> Hi R users, >> >> I have a question about reading from text files. The file has the structure >> below: >> >> Time Column1 Column2 >>

Hi - How can I add different notes in different panels? My data looks like ID Dose Visit Value 1 1 0 0.5 -6.5802e-02 2 1 0 1.0 2.4085e-01 3 1 0 1.5 -2.2907e-01 4 1 0 2.0 2.4074e-01 ... ... ... ... ... 270 45 30 3.0 -8.1316e-01 271 46 30 0.5 -2.0786e-01 272 46 30 1.0 -2.9336e-01 273 46 30 1.5 -5.5657e-01 274 46 30 2.0 -2.1659e-01

Hi the list, I use rgl to produce a 3D graph. I would like to "show" this graph to some collaborator. Is there a way to save it and send it to someone else? Christophe Genolini