similar to: CIS inquiries

I sent the message below to r-sig-mac yesterday, but having no reply I decided to explore a bit myself and found that editing: /Library/Frameworks/R.framework/Versions/2.3/Resources/share/make/ shlib.mk yzzy: diff shlib.mk shlib.mk~ 3c3 < include $(R_HOME)/etc/Makeconf --- > include $(R_HOME)/etc${R_ARCH}/Makeconf restored the functionality of R CMD INSTALL. Is this a known issue?

Motivated by Deepayan's recent inquiries about the efficiency of the R 'quantile' function: http://tolstoy.newcastle.edu.au/R/devel/05/11/3305.html http://tolstoy.newcastle.edu.au/R/devel/06/03/4358.html I decided to try to revive an old project to implement a version of the Floyd and Rivest (1975) algorithm for finding quantiles with O(n) comparisons. I used

----- Forwarded Message ---- From: Nnamdi <nnamdii at yahoo.com> To: roger koenker <rkoenker at uiuc.edu> Sent: Friday, April 18, 2008 1:08:58 PM Subject: Re: [R] efficiently replacing values in a matrix I tried the sparse matrix implementation, still there are issues: a <- matrix(nrow=10000,ncol=10000) > a.csr <- as.matrix.csr(a) Error in if (nnz == 0) { : missing value

Does anyone have (good) experience converting tables of tournament results into dendrogram-like graphics? Tables, for example, like this: read.table(url("http://www.econ.uiuc.edu/~roger/research/ncaa/NCAA.d")) Any pointers appreciated. RK url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu Department of Economics vox: 217-333-4558

> On Mar 26, 2015, at 8:50 AM, Martyn Plummer <plummerM at iarc.fr> wrote: > > On Wed, 2015-03-25 at 15:12 -0500, Roger Koenker wrote: >> Thierry, >> >> I have this: >> >> if (require(MatrixModels) && require(Matrix)) { >> X <- model.Matrix(Terms, m, contrasts, sparse = TRUE) > > You have this in the current release, which

I'm trying to install the current R-devel on a G4 powerbook. I'm using ./configure --enable-R-shlib --with-blas='-framework vecLib' --with-lapack with rudjer: gcc -v Reading specs from /usr/libexec/gcc/darwin/ppc/3.1/specs Thread model: posix Apple Computer, Inc. GCC version 1175, based on gcc version 3.1 20020420 (prerelease) rudjer: g77 -v Reading specs from

I'd like to mention that there is a new quantile regression package "nprq" on CRAN for additive nonparametric quantile regression estimation. Models are structured similarly to the gss package of Gu and the mgcv package of Wood. Formulae like y ~ qss(z1) + qss(z2) + X are interpreted as a partially linear model in the covariates of X, with nonparametric components defined as

I have a pdf scan of several pages of data from a quite famous old paper by C.S. Pierce (1873). I would like (what else?) to convert it into an R dataframe. Somewhat to my surprise the pdf seems to already be in a character recognized form, since I can search for numerical strings and they are nicely found. Of course, as is usual with such tables there are also headings and column

Negative indexing is often handy, but I'm in need of an appropriate idiom for handling cases in which the index set can be null: x <- rnorm(5) a <- 1:5 s <- rep(FALSE,5) y <- x[-a[s]] # I'd like y == x but instead one has x[-a[s]] == x[a[s]] == numeric(0), which is rather # unfortunate -- so far the best I have come up with is: as <- ifelse(length(a[s]),-a[s],TRUE)

My apologies, the data can now be found at: url <- "http://www.econ.uiuc.edu/~roger/research/ebayes/velo.d" x <- scan(url,skip = 1) If I could get each of the histograms to mimic what is produced by hist(x, 100, freq = FALSE) I?ve experimented with xlim, ylim, without success so far... url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu

Greetings, I'm organizing summer econometrics lunch meetings to discuss thesis work. The first meeting will be this friday July 18 12-1pm in the conference room on the third floor of Wohlers. The first talk will be by Lingjie Ma Control Variate Estimation of Structural Quantile Regression Models url: www.econ.uiuc.edu/~roger/my.html Roger Koenker email rkoenker@uiuc.edu Department of

An extremely modest proposal: It would be nice if packages could have a FAQ and if faq(package.name) would produce this faq. And if, by default faq() FAQ() would produce the admirable R faq... Apologies in advance if there is already a mechanism like this, but help.search() didn't reveal anything. url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at

I'm writing a new predict method and would like to be able to pass an argument called "se" via the "..." mechanism. However, predict has a "se.fit" argument that wants to interpret my specification of "se" as a partially matched version of se.fit. Surely there a standard treatment for this ailment, but I can't seem to find it. url:

What is the current state of the R-art for ordered probit models, and more esoterically is there any available R strategy for ordered cauchit models, i.e. ordered multinomial alternatives with a cauchy link function. MCMC is an option, obviously, but for a univariate latent variable model this seems to be overkill... standard mle methods should be preferable. (??) Googling reveals that spss

It being friday, I would like to call your attention to an innovative data analysis by Leslie Lamport available from: http://research.microsoft.com/users/lamport/pubs/hair.pdf I particularly liked the graphics. url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu Department of Economics vox: 217-333-4558 University of Illinois

In the last few days I've received couple of messages pointing out that our SparseM package fails to install on the patched version of 1.7.0. Laurent Gaultier kindly suggested that replacing: setGeneric("as.matrix.csr") by setGeneric("as.matrix.csr", function(x, nrow, ncol, eps) standardGeneric("as.matrix.csr")) was sufficient to fix the problem.

For the record: I had an inquiry about a problem with the function lm.fit.recursive in my quantreg package on a 64bit gentoo linux system. After some debugging the problem was traced to the fact that on this system g77 took the view that the statement: data one/1.0d0/ meant that the symbol one should take the value zero. Replacing such data statements with their parameter statement

I'm still having difficulties with methods...is the following behavior expected? (This is without loading any libraries) R : Copyright 2003, The R Development Core Team Version 1.7.0 Under development (unstable) (2003-03-31) > ls() character(0) > > diag(2) [,1] [,2] [1,] 1 0 [2,] 0 1 > setGeneric("diag") [1] "diag" > diag(2) Error in

In a private response to Tony Plate's suggestion to replace version() output with sessionInfo() in R-help requests, > roger koenker wrote: >> Thanks for this, it would seem useful to have version numbers for >> the packages too? and Tony replied, > > Sounds sensible to me! If I were you I'd send a message to R-devel > suggesting this. AFAIK, some changes to

Johannes, You can minimize an model expression by just putting the ~ on the left and everything else on the righthand side, but I don't think that this is really what you want. In the NLS expression this would ignore the jacobian of the transformation from errors to response, and in nlrq there is the same problem, however you can adjust for the jacobian by rescaling by the geometric mean of