similar to: Rprof

I'm having difficulty with Rprof. I have a documentation example test.qss that runs fine without profiling, but under Rprof, > Rprof() > source("test.qss") Error in standardGeneric("model.matrix") : UseMethod used in an inappropriate fashion Luke wrote about a similar circumstance last summer: # From: Luke Tierney (luke@stat.umn.edu) # Date: Fri Jul

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?

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 discovered recently that the phrase "Oja median" produces no hits in Jonathan Baron's very valuable R search engine. I found this surprising since I've long regarded this idea as one of the more interesting notions in the multivariate robustness literature. To begin to remedy this oversight I wrote a bivariate version and then decided that writing a general p-variate version

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

I have a vector of logLik values that I'd like to return and it would be nice if the the print method didn't run them together. Could I make a plea for using sep = " ", rather than sep = "" in print.logLik? url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu Department of Economics vox: 217-333-4558

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

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

> 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

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

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)

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

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

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

We have a solaris/sparc machine that has been running an old version of R-devel: Version 2.2.0 Under development (unstable) (2005-06-04 r34577) which was built as m64 from sources. Attempting to upgrade to 2.4.0 the configure step goes ok, but I'm getting early on from make: > gcc -m64 -L/opt/sfw/lib/sparcv9 -L/usr/lib/sparcv9 > -L/usr/openwin/lib/sparcv9 -L/usr/local/lib -o

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:

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

Thierry, I have this: if (require(MatrixModels) && require(Matrix)) { X <- model.Matrix(Terms, m, contrasts, sparse = TRUE) in my function rqss() I've tried variants of requireNamespace too without success. If I understand properly model.Matrix is from MatrixModels but it calls sparse.model.matrix which is part of Matrix, and it is the latter function that I'm not