Publications by Douglas Bates

The name

11.03.2012

When I bought the first computer for use in our Statistics Department – a Vax 11/750 that cost about a quarter of a million dollars in 1983 – I was considered extravagant because I purchased and installed a second megabyte of memory for the machine. Related To leave a comment for the author, please follow the link and comment o...

673 sym

An RcppEigen example

11.03.2012

R is an Open Source project providing an interactive language and environment for statistical computing.  It has become the lingua franca for research in statistical methods.  Because R is an interpreted language it is comparatively easy to develop applications (called packages) for it.  However, the resulting code can sometimes be slow, wh...

5528 sym

A Julia version of the multinomial sampler

12.03.2012

In the previous post on RcppEigen I described an example of sampling from collection of multinomial distributions represented by a matrix of probabilities.  In the timing example the matrix was 100000 by 5 with each of the 100000 rows summing to 1.  The objective is to create a vector of 100000 1-based indices representing a sample...

2169 sym R (586 sym/2 pcs)

Julia functions for the Rmath library

26.03.2012

1.1 Signatures of the d-p-q-r functions in libRmath Users of R are familiar with the functions for evaluating properties of and for sampling from probability distributions, the so-called d-p-q-r functions. The designation d-p-q-r reflects the fact that, for each distribution, there are up to 4 such functions, each beginning with one of these l...

10242 sym R (6653 sym/14 pcs)

PK models in R and in Julia

29.03.2012

1.1 Pharmacokinetic models with an analytic solution Pharmacokinetics is the study of the absorption and elimination of drugs and their metabolites in the body. As described in the wikipedia article there are several parameters, such rate constants, volumes of distribution and clearances, that we wish to estimate from pharmacokinetic data, whic...

5323 sym R (3525 sym/8 pcs)

An R programmer looks at Julia

07.04.2012

In January of this year I first saw mention of the Julia language in the release notes for LLVM. I mentioned this to Dirk Eddelbuettel and later we got in contact with Viral Shahregarding a Debian package for Julia. There are many aspects of Julia that are quite intriguing to an R programmer. I am interested in programming languages for “Comp...

11082 sym R (6766 sym/18 pcs)

The Simple Gibbs example in Julia

21.05.2012

The Gibbs sampler discussed on Darren Wilkinson’s blog and also on Dirk Eddelbuettel’s blog has been implemented in several languages, the first of which was R. In preparation for a session at useR!2012 on “What other languages should R users know about?”, Dirk, Chris Fonnesbeck and I have considered implementations of thi...

4780 sym R (2847 sym/14 pcs)