similar to: list of functions flagged with debug() or trace()

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

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

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

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, I would like to plot two superposed surfaces (original data, and smoothed models). Basically, I would like to create a usual persp plot with box, and then superpose a second surface (without box and axes) on the top plane of the box defined by the first graph. I don't know if it is clear :-( Is it possible to draw such a second surface in an existing persp plot in R? Another (related)

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 all, I have a bit of a newbie question here that I hope y'all can help with. I've run a cluster analysis using hclust on about 500 objects (using R1.3 under Win 2000). The problem is that the tips of the dendrogram are so close together on the plot that the labels overlap and are unreadable. I've used "cex" to reduce the label sizes but this isn't sufficient with so

I stumbled across this link and was wondering if this has any relation to the R project. http://www.fas.umontreal.ca/BIOL/Casgrain/en/labo/R/index.html -jennifer ----------------------------------- J. Steinbachs, Ph.D. Computational Biologist http://compbiology.org ----------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

Hi, I make a glm analysis with binomial errors and its OK. Now I need plot the ajusted logist line in my points. How do I make it? | * * * * -> I need plot an ajusted | * * line in this points. | * * | * | * * | * * |* * * * |_________________________ Exist an abline like function for ajust others formulas? Exist in R

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Is it possible to create a time series from existing data frame? One column of that frame contains dates (like "01/01/01") and the second one the values. The period is non-regular (it is not strictly daily nor working-days-daily). Thank you. lukas - -- Lukas Kubin lukas.kubin at permonik.com phone: 00420603836180 -----BEGIN PGP

Dear R and S users, I've been working on fitting finite mixture of negative exponential distributions using maximum likelihood based on the example given in MASS. So far I had much success in fitting two components. The problem started when I tried to extend the procedure to fit three components. More specifically, likelihood = sum( ln(c1*exp(-x/lambda1)/lambda1 + c2*exp(-x/lambda2)/lambda2

Hi, I need help regarding estimating a linear model where restrictions are imposed on the coefficients. An example is as follows: Y_{t+2}=a1Y_{t+1} + a2 Y_t + b x_t + e_t restriction a1+ a2 =1 Is there a function or a package that can estimate the coefficient of a model like this? I want to estimate the coefficients rather than test them. Thank you for your help Ahmad Abu Hammour --------------

Colleagues ---------------------------------- System info: Version 1.2.3 (2001-04-26) on NT ESS v. 5.1.18 using emacs ver. 20.4 ---------------------------------- I am wondering if there is a problem with the function rep(). Both the commands using rep() below were expected to produce 8 twos, but only the second did. x <- rep(2,40*(1-0.8)) length(x) y <- rep(2,40*0.2) length(y)

Hi, Is there a simple command to pause a script until the user hits return? Alternatively, is there a command which prompts the user for some input and assigns the input into a variable? dave +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ |David Richmond It works on a | + Dept. of Sociology complex scientific + |Saint Mary's College