similar to: name of a name

Dear list, I wish to visualise some 4D data as a kind of colour / translucent cloud in 3D. I haven't seen such plots in R (but perhaps I missed a feature of rgl). The easiest option I found would be to export the data in povray's df3 (density file) format and visualise it with povray. The format specification baffles me a little, http://www.povray.org/documentation/view/3.6.1/374/ ;

Dear list, I made the following example of a proto object that contains some data and a spline interpolation. I don't understand why test$predict() fails with this error message: Error: evaluation nested too deeply: infinite recursion / options(expressions=)? Best regards, baptiste test <- proto(source = data.frame(x=1:10, y=rnorm(10)), raw = function(.){

Hi list, This is a very basic question about lattice: I wish to add some vertical lines in each panel of a xyplot as demonstrated in this example: > library(lattice) > > xx <- seq(1, 10, length=100) > x <- rep(xx, 4) > y <- c(cos(xx), sin(xx), xx, xx^2/10) > fact <- factor(rep(c("cos", "sin", "id", "square"), each=100))

Dear list, I've been trying this for a few hours and I just don't understand how lattice works with groups and subscripts. Consider the following example, > xx <- seq(1, 10, length=100) > x <- rep(xx, 4) > y <- c(cos(xx), sin(xx), xx, xx^2/10) > fact <- factor(rep(c("cos", "sin", "id", "square"), each=100)) > fact2

Dear all, I am stumped at what should be a painfully easy task: predicting from an lm object. A toy example would be this: XX <- matrix(runif(8),ncol=2) yy <- runif(4) model <- lm(yy~XX) XX.pred <- data.frame(matrix(runif(6),ncol=2)) colnames(XX.pred) <- c("XX1","XX2") predict(model,newdata=XX.pred) I would have expected the last line to give me the

Hi, I wonder whether there is a way to generate a polygon (a triangle in my case) with color gradient using grid.polygon() in package grid? I tried something like library(grid) grid.polygon(x=c(0, 0.5, 1), y=c(0.5, 1, 0.5), gp=gpar(col=NA, fill=colorRampPalette(c("green", "lightgray"), space="Lab")(200))) But am only

On 09/24/2015 09:23 AM, Jean-Baptiste Vignaud wrote: > Hello; > > Some issue with virtual storage and fts-lucene: > > Sep 24 08:12:16 imap(jbv00 at x <jbv00 at d1.dc.xandmail.com>xx): Panic: file > virtual-storage.c: line 370 (virtual_backend_box_close): assertion failed: > (mbox->backends_open_count > 0) Fixed: http://hg.dovecot.org/dovecot-2.2/rev/4b5b8e6a4a32

Hello, I was asked to try the following code on R, gamma.mles function (xx,shape0,rate0) { n<- length(xx) xbar<- mean(xx) logxbar<- mean(log(xx)) theta<-c(shape0,rate0) repeat { theta0<- theta shape<- theta0[1] rate<- theta0[2] S<- n*matrix(c(log(rate)-digamma(shape)+logxbar,shape/rate-xbar),ncol=1) I<- n*matrix(c(trigamma(shape),-1/rate,-1/rate,shape/rate^2),ncol=2)

I have data in a format like this: name ssex sex view num rating rt ahl4 f m f 56 -108 2246 ahl4 f m f 74 85 1444 ahl4 f m f 52 151 1595 ahl4 f m f 85 1 1447 ahl4 f m f 53 46 1716 ahl4 f m f 37 145 1276 ahl4 f m f 50 98 1465 ahl4 f m f 51 -26 1322 ahl4 f m f 38 -97 1790 ahl4 f m f 14 -158 865 ... ahl4 f m p 43 -136 1669 ahl4 f m p 10 -59 808 ahl4 f m p 67 -111 1279 ahl4 f m p 85 -86 994 ahl4 f m

Hi, I used the write.table function to save data in txt file, and this is the output: V1 V2 V3 V4 1 YAL005C 21 14 11 2 YAL007C 2 1 4 3 YAL012W 8 16 3 4 YAL016W 24 23 23 5 YAL019W 3 3 2 6 YAL020C 2 4 2 7 YAL021C 7 5 5 8 YAL022C 3 1

The Subject line mostly says it. I'm designing it as a semester-long, 3 hours per week, course that takes in students who got the basics of R in stats classes, but don't really know how to program in it. Translation: if their own examples don't look enough like examples from previous work, they're stumped. Does anybody have a text for an intermediate R course (but not too

I have some Bernoulli data something like this: x<-sort(runif(100,1,20)) p<-pnorm(x,10,3) y<-as.numeric(runif(x)<p) plot(x,y) lines(x,p) This plot is not very satisfactory because the ogive does not visually fit the (0,1) points very well, and also because the points tend to fall on top of one another. The second problem can be eliminated by adding vertical jitter. However I was

Dear List, I have a 47 species * 83 samples matrix containing percentage abundance data. I have two cluster analyses one of the samples and one of the species, and have ordered the rows and columns of the species by samples matrix according to these two cluster analyses. So far so good! Now what I want to do is create a plot with the species dendrogram at the top of the plot, the samples

Dear R users, I have a data frame (test) including five columns of upper (numeric), lower (numeric), observed (numeric), best_sim (numeric) and stname (factor with 80 levels, each level with different length). Now I would like to write a short program to draw one graph as follow for each level of stname but I would like also to draw each time 12 graphs for the 12 levels of stname in the same

Hi All, I hope someone can help me. We use Rsync to transfer files to and from our Australian and US server and we are experiencing files being overwritten and a few files disappearing (their contents disappear). I have setup ntp on each machine pointing to the same ntp server and set the hwclock to utc. We have found that the files on the US server are x:xx:01 secs and the AU server is

Hello, I am new to R program, therefore, I am sorry if this is a really stupid question. I wrote a simple function and for some reason it doesn't work ReturnsGrid = function(x,y,m){ for (i in 1:m){ grid[i] <- x + (i-1)*(y-x)/m } grid } xx=ReturnsGrid(0,9,3) Thanks a lot! [[alternative HTML version deleted]]

Hi, I want to make the following array of numbers: -3 -3 -3 -3 -3 -3 -3 -2 -2 -2 -2 -2 -2 -2 ... 3 3 3 3 3 3 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 ... -3 -2 -1 0 1 2 3 (3 would be N, a painful example to type number). Here is my dirty attempt to do it, > N<-3 > > x<-c(-N:N) > > rj<-rbind(matrix(outer(matrix(1,1,2*N+1),x),nrow=1),rep(x,2*N+1)) > It

Hi list, i have a problem in one of my centos boxes. I installed CentOS 4.3. I addedd additionl interface card in this server and when i execute ifconfig command i see devNNN(where n is a random number). And when i restart the server the NNNN always change. dev9926 Link encap:Ethernet HWaddr xx:xx:xx:xx:xx:xx inet addr:xxx.xxx.xxx.xxx Bcast:xxx.xxx.xxx.xxx Mask: xxx.xxx.xxx.xxx

Dear all, Can anybody tell me how to export a 3d figure made with the plot3d function? I'm careless about whether it's still interactive or not in another format, as long I can get it out of R. Thanks! Alejandro Gonz?lez Departamento de Biodiversidad y Conservaci?n Real Jard?n Bot?nico Consejo Superior de Investigaciones Cient?ficas Claudio Moyano, 1 28014 Madrid, Spain Tel +0034

Why does R think these numbers ***are*** equal? In a somewhat bizarre set of circumstances I calculated x0 <- 0.03580067 x1 <- 0.03474075 y0 <- 0.4918823 y1 <- 0.4474461 dx <- x1 - x0 dy <- y1 - y0 xx <- (x0 + x1)/2 yy <- (y0 + y1)/2 chk <- yy*dx - xx*dy + x0*dy - y0*dx If you think about it ***very*** carefully ( :-) )