similar to: R-help Digest, Vol 4, Issue 27 ( -Reply)

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, I am drawing four contour plots on a 2x2 layout. I need to downsize the contour line labels. cex and labcex do not seem to work. Any suggestions? Thanks, John. -- -------------------------------------------------------------------------- Dr. John Janmaat Department of Economics, Acadia University, Wolfville, NS, B4P 2R6 E-mail: jjanmaat at acadiau.ca Web:

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

Bonjour, J'ai t?l?charg? SciViews Insider que je trouve tr?s convivial. Par contre, je n'arrive pas ? comprendre comment enregistrer un script R en type de fichier R justement. Mes programmes fonctionnent tr?s bien, mais SciViews me propose uniquement de les enregistrer au format txt sous un type de fichier "bloc notes". Comment les enregistrer avec l'extension .R comme le

This only happens if 'file' is a text connection, and is the expected behaviour in that case: you cannot capture an incomplete line to a text connection. There seems no reason to break the documented behaviour in other cases to change something that you consider to a bug when file=NULL and the user does not produce complete output. It would be possible to make use of isIncomplete()

I used the packaged "MinGW-2.0.0-3.exe" exactly as specified on http://www.stats.ox.ac.uk/pub/Rtools/ - in fact I used these recommendations throughout. According to the release notes MinGW version 2.0.0 contains the following list of packages: GCC-3.2-core-20020817-1 binutils-2.13-20020903-1 mingw-runtime-2.2 w32api-2.0 gdb-5.1.1-1 make-3.79.1-20010722 (binary renamed as mingw32-make)

Hello All, Is there an analogue to \underbar or the AMS math \underline in graphical math expressions? Thanks, John. -- ===================================================================================== Dr. John Janmaat Department of Economics Acadia University Wolfville, Nova Scotia, Canada B4P 2R6 TEL: 902-585-1461 WWW: http://ace.acadiau.ca/~jjanmaat/ EMAIL: jjanmaat at acadiau.ca June

Hello, I'm playing around with cluster analysis, and am looking for methods to select the number of clusters. I am aware of methods based on a 'pseudo F' or a 'pseudo T^2'. Are there packages in R that will generate these statistics, and/or other statistics to aid in cluster number selection? Thanks, John. --

Hello All, Is there a fast way to find the index(row and column) of a point in a matrix? Thanks, John. -- -------------------------------------------------------------------------- Dr. John Janmaat Department of Economics, Acadia University, Wolfville, NS, B4P 2R6 E-mail: jjanmaat at acadiau.ca Web: http://ace.acadiau.ca/~jjanmaat Tel: 902-585-1461 Fax: 902-585-1070

Hello, Thanks to Douglas Bates, there is now a new benchmark suite (version 2.3) which is compatible with R 1.9.0 and the recent Matrix library (0.8-1 or above). You find it at http://www.sciviews.org/other/benchmark.htm. It compares R 1.9.0 under Windows with: S-PLUS 6.5, Matlab 6.0, O-Matrix 5.6, Octave 2.1.42, Scilab 2.7 and Ox 3.30. In short, R in its version 1.9.0 and with the new Matrix

Hello All, I'm looking to test a variable in a logit model (glm(..., binomial(link="logit"))) for exogeneity (endogeneity). At this point I am planning to try implementing Jeffery Grogger's "A Simple Test for Exogeneity in Probit, Logit, and Poisson Regression Models", Economic Letters, 1990. To do this, I need to be able to do an instrumental variables NLS

I have both the CRAN binary and my own compiled version of R v1.7.1 on the same machine (Dell Pentium III 800 MHz running NT 4) Using the example provided earlier: > phi <- 1.6180339887498949 > a <- floor(runif(750000)*1000) > system.time(b <- (phi^a - (-phi)^(-a))/sqrt(5))[3] I get 10.99 secs on the CRAN binary and 2.09 secs on my own compiled version. I hope this helps

Dear all, I am trying to get an estimate of the intercept for a linear model. In this case, I know the slope of the model, can anyone tell me how to constrain the formula in lm() so that it only estimates the intercept not the slope? Many thanks in advance, Sincerely, Liqing Zhang Dept. of Eco. Evol. Biol. Univ. of CA, Irvine email: lzhang at uci.edu >From VM Mon Apr 30 08:18:45 2001

Well.. plans are there from a long time to rewrite SciViews completely and make it platform independent (to work on Linux/Unix and MacOS X, as well as Windows). I have done some work in this direction when time permitted, but I am pretty busy with other work. During the holidays, I will continue to work in this direction. I will try to package a first running version of SciViews compatible with

Hi all, After being a long time off the R and R-SIG-GUI lists, because I was working hard on SciViews, I am happy to deliver the first public version of this program. It is still v. 0.1, that is, an incomplete alpha version. It currently installs only on Windows 2000 or XP (NOT on Windows 9X/ME/Millenium, nor any other platform!). You can download it at http://www.sciviews.org. Please, note that

Sorry for entering late in the thread, but I was absent last week. I am developping an alternate GUI for R (and also Splus, Octave, Scilab, Matlab, Ox & Mathematica). The project started 3 years ago, but it is progressing slowly, since I am currently the only developer on the project and had to do it during my extra time. Now, it is going a little faster because I got some time and money

This reminds me of a question I've never asked yet: Is there a trick to make this work when a username authentication is required for using the proxy? I need a user name and a password to connect to the Internet, and so I get: > version _ platform i386-pc-mingw32 arch x86 os Win32 system x86, Win32 status major

Hi, For those who are interested, I have update my R benchmark (to version 1.5.0) and also to Splus 6. They are available at: http://www.sciviews.org/other/benchmark.htm. A comparison is made between Matlab (5.3 & 6), R 1.5.0, Splus 6 rel 2, O-Matrix 5.1, Octave 2.1.31, Scilab 2.6, Rlab 2.1 and OX 3.00 under Windows 2000 pro. Overall, R is not the fastest package, but it is one of the

Hi Jason, I suppose you installed the Matrix library, and it is working on your computer? If yes, may be det.Matrix() was removed, or renamed in the Matrix library you have (I cannot check this for the latest version, because I am away of the office until August 1st), but I will do that next week. In the meantime, you can replace 'det.Matrix' by 'det.default', and it should run.

> From: Philippe Grosjean > > Duncan, > > I agree totally with you on all points, now that we clarified our > respective ideas. I am afraid I probably agree also with your last > point, from a theoretical point-of-view ("I still think we need more > glue and am working on that while we continue to experiment with the > design of GUIs for the next 5 years