Today I did some experiments to learn how to combine R and LaTeX to create reproducible research reports. Here are my first results:First Demo ...

R news and tutorials contributed by (552) R bloggers

Today I did some experiments to learn how to combine R and LaTeX to create reproducible research reports. Here are my first results:First Demo ...

Today I did some experiments to learn how to combine R and LaTeX to create reproducible research reports. Here are my first results:First Demo ...

Introduction A while ago, one of my co-workers asked me to group box plots by plotting them side-by-side within each group, and he wanted to use patterns rather than colours to distinguish between the box plots within a group; the publication that will display his plots prints in black-and-white only. I gladly investigated how to

R 3.1.0 (codename “Spring Dance“) was released today! You can get the source code from http://cran.r-project.org/src/base/R-3/R-3.1.0.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. The full list of new features and bug fixes is provided below. Upgrading to R 3.1.0 You can

categoryCompare Paper Finally Out!I can finally say that the publication on my Bioconductor package categoryCompare is finally published in the Bioinformatics and Computational Biology section of Frontiers in Genetics. This has been a long time coming, and I wanted to give some background on the inspiration and development of the method and software.TL;DRThe software package has been...

In the context of AR(1) processes, we spent some time to explain what happens when is close to 1. if the process is stationary, if the process is a random walk if the process will explode Again, random walks are extremely interesting processes, with puzzling properties. For instance, as , and the process will cross the x-axis an infinite number...

The question is: can we automate scientific discovery, and what might an interface to such a tool look like. I’ve been experimenting with automating simple and complex data analysis and report generation tasks for biological data and mostly using R and LATEX. You can see some of my progress and challenges encountered in the presentation

In my paper on the impact of the recent fracking boom on local economic outcomes, I am estimating models with multiple fixed effects. These fixed effects are useful, because they take out, e.g. industry specific heterogeneity at the county level - or state specific time shocks. The models can take the form: where is

Two R tutorials for beginnersI am currently in the process of rescuing some of the pages from my now defunct datajujitsu.co.uk blogger blog and moving to this Github/Clojure/Bootstrap version. I also today gave a tutorial to the University of Manche...

Consider some ARCH() process, say ARCH(), where with a Gaussian (strong) white noise . > n=500 > a1=0.8 > a2=0.0 > w= 0.2 > set.seed(1) > eta=rnorm(n) > epsilon=rnorm(n) > sigma2=rep(w,n) > for(t in 3:n){ + sigma2=w+a1*epsilon^2+a2*epsilon^2 + epsilon=eta*sqrt(sigma2) + } > par(mfrow=c(1,1)) > plot(epsilon,type="l",ylim=c(min(epsilon)-.5,max(epsilon))) > lines(min(epsilon)-1+sqrt(sigma2),col="red") (the red line is the conditional variance process). > par(mfrow=c(1,2)) > acf(epsilon,lag=50,lwd=2)...