similar to: multiple persp plots

Any ideas on adding text to persp plots? The perspp function from S-plus is not implemented. I'm using R1.3.0 on Win2000. Thanks, Sundar -- Sundar Dorai-Raj, Ph.D. Assistant Professor Discrete & Statistical Sciences Office: (334) 844-3647 231 Allison Lab Fax: (334) 844-3611 Auburn University Email: dorairaj at dms.auburn.edu Auburn AL 36849-5307

Hi, I have submitted a package to CRAN in ftp://ftp.tuwien.ac.at/incoming/pastecs_1.0.tar.gz. The R extension manual advises to send a mail to cran at ci.tuwien.ac.at about it, but it fails. I have also tried WWWadmin at ci.tuwien.ac.at. Are these email adresses still valid? Best, Philippe Grosjean ...........]<(({?<...............<?}))><............................... ) ) ) )

Hi, I have 16 tables of 39 variables with around 1500 to 3500 cases in each table. These 16 tables are in 16 different data frames. I want to merge them in a single large data frame, so: All <- rbind(T1, T2, ...) The problem: it is VERY slow, and takes a lot of memory (> 500 Mb) on my Athlon XP1800+, 1.75 Go memory, Win 2000 Pro SP1, R 1.4.1. Any suggestion to improve speed, and/or memory

>>>>> "PG" == Philippe Grosjean <phgrosje at ulb.ac.be> writes: PG> I am interested by your list of students wishes, since I am PG> developing a GPL user interface for math/stat calculation PG> engines (among others: R and Splus, but also Octave, Scilab PG> and Mathematica,...). It will run under Windows only (for the PG> moment),

It's two days I try to answer to a mail of Achim Zeileis, the maintainer of strucchange package, using the mail address: zeileis at ci.tuwien.ac.at. I got the following error: ----- The following addresses had permanent fatal errors ----- <zeileis at ci.tuwien.ac.at> (reason: 550 5.7.1 <zeileis at ci.tuwien.ac.at>... Access denied) Can anyone tell me how to reach him?

using: identify(...., plot=FALSE) allows identifying points in a graph without printing the corresponding label in the graph. This works fine in R for Windows. However, when I copy the graph as metafile into Word, the corresponding label appears in black font in the copied picture. Obviously, identify() prints the label anyway, in "invisible" font or so, in the graph, and this

Hello, I am doing several self-starting growth functions for using with nls(). When I list the self-starting functions included in nls library, for instance, SSlogis, there is: > SSlogis function (input, Asym, xmid, scal) ... <environment: 03476D20> attr(,"class") [1] "selfStart" What is this <environment: 03476D20> instruction? By using deriv() and then

Generalised Inverse: The Moore-Penrose Generalisied Inverse is probably better defined as a pseudo-Inverse that arises in solving least squares problems. Another well known pseudo-Inverse is the so-called Drazin pseudo-Inverse. If memory serves (and it's been 10-12 years!) it can be obtained via a diagonalisation. Anyway, I dare say Prof. Ripley (among others) probably has "all the

Hello all, Since I got no answer to my first mail, I suspect I did not formulate my question adequately, so I try again. In R, one can debug or trace a function 'foo' with debug(foo) or trace(foo), respectively. This leads to a special treatment of the function 'foo' until one enter undebug(foo) or untrace(foo). I would like to know if there is a convenient way to know at any

Hi, Under Windows, I am exporting graphs from R (in metafile -emf- format), and want to combine them using a graphic package. It happens that the background of the graph is represented by a filled rectangle (with the color given by bg=...). This is to be expected. However, I find no options to disable this and produce a graph with transparent background. Does such an option exist? If not, then

Hello everybody, I am making a modified version of "Stephan Steinhaus' benchmark test for number crunching, v. 2, (see http://www.scinetificweb.com/ncrunch/ncrunch.pdf for the original version), comparing several functions of some math/stat software. R is not performing bad at all... except for the sorting of a 1,100,000 random vector (test #3) which is the worst of all (see cell F3 in

Hello everybody, I am making a modified version of "Stephan Steinhaus' benchmark test for number crunching, v. 2, (see http://www.scinetificweb.com/ncrunch/ncrunch.pdf for the original version), comparing several functions of some math/stat software. R is not performing bad at all... except for the sorting of a 1,100,000 random vector (test #3) which is the worst of all (see cell F3 in

Hi, I'm new to R (and math ;) Would somebody please be so kind as to direct me in plotting a 3D sphere? I tried something in the lines of: #### y <- x <- seq(-pi, pi, length=pi*10) f <- function(x,y) { z <- sqrt(pi - x^2 - y^2) #z[is.na(z)] <- 0 z } z <- outer(x, y, f) persp(x, y, z, theta = 120, phi = 30) #### I've also tried: .... make.surface.grid(...) ..

Hello, nls library provides 6 self-starting models, among them: SSfp, a four parameters logistic function. Its self-starting procedure involves several steps. One of these steps is: pars <- as.vector(coef(nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xydata, start = list(lscal = 0), algorithm = "plinear"))) which assumes an initial value of lscal equal to 0. If lscal

Thorsten is right. There is a direct formula for computing the Moore-Penrose inverse using the singular value composition of a matrix. This is incorporated in the following: mpinv <- function(A, eps = 1e-13) { s <- svd(A) e <- s$d e[e > eps] <- 1/e[e > eps] return(s$v %*% diag(e) %*% t(s$u)) } Hope it helps. Dietrich

I am currently working on surface plots (with R 2.0.1) and I want to improve the aspect of them. Having a typical 'cow-boy hat' surface [source: Spector (1994), p. 206]: pts <- seq(from = -8, to = 8, length = 50) cow <- function(x, y) { z <- sqrt(x^2+y^2) sin(z)/z } out <- outer(pts, pts, cow) sur <- persp(x=pts, y=pts, z=out, theta=45, phi=30) What I exactly wanted is:

Hello list. I very often need 3d scatterplots, and use scatterplot3D quite a lot. I am trying to modify persp() to plot scatterplots, and make use of the theta and phi arguments that persp() offers. I am having some difficulty passing the correct arguments to persp(). Here is my function so far. Much of it is copied from the persp() manpage. points3d <- function(x,y,z,

Hello list. I very often need 3d scatterplots, and use scatterplot3D quite a lot. I am trying to modify persp() to plot scatterplots, and make use of the theta and phi arguments that persp() offers. I am having some difficulty passing the correct arguments to persp(). Here is my function so far. Much of it is copied from the persp() manpage. points3d <- function(x,y,z,

I have this data set that both x & y are ordered vectors of length 600 & 700 respectively; z is a 600 by 700 matrix whose entry z[i,j] is either a missing value (indicated by 'NaN') or a real number between 0 and 1. The contour function contour(x,y,z) gives me a blank picture. I guess the reason is that most of z-entries are missing, only less than 1% are non missing.

Hi folks, I'm experiencing some hassle to add some text on a persp plot... Here is the code to generate the persp plot: x <- seq(-1.95, 1.95, length = 30) y <- seq(-1.95, 1.95, length = 30) z <- outer(x, y, function(a,b) a*b^2) persp(x, y, z,xlim=c(-2,2),ylim=c(-2,2),zlim=c(-8,8), phi=30, theta=-30,nticks=8,ticktype="detailed") I've tried the text() function :