similar to: p-values from a resampled reference distribution

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

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

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

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

For a paper on "Reproducible Econometric Research," Achim Zeileis and I would be interested in hearing from anyone with (positive) experience and/or suggestions about GUI-based data analysis from a reproducibility perspective. We omitted coverage of this topic in an initial draft on the well-known principle: if you can't say something nice, say nothing at all,"

I've installed the rgl package on a Suse x86-64 machine (further details below) and it produces nice screen images. Unfortunately, rgl.snapshot attempts to make png files produces only the response "failed". For other graphics png() works fine, and capabilities indicates that it is there. If anyone has a suggestion of what might be explored at this point I'd be very

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

I was experimenting yesterday with a binomial make.link option for estimating student t binary response models, tentatively called gossit, and I noticed eventually that the R qt function doesn't like df < 1. Vaguely recalling that Splus didn't seem to mind such weirdness, I checked on our soon to be defunct Splus6.2 and sure enough, it produced plausible answers instead of R's