964 search results for "Latex"

Combining R and LaTeX with Sweave

January 6, 2009
By

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

Read more »

Combining R and LaTeX with Sweave

January 6, 2009
By

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

Read more »

Side-by-Side Box Plots with Patterns From Data Sets Stacked by reshape2 and melt() in R

Side-by-Side Box Plots with Patterns From Data Sets Stacked by reshape2 and melt() in R

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

Read more »

R 3.1.0 is released!

April 10, 2014
By
R 3.1.0 is released!

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

Read more »

categoryCompare Paper Finally Out!

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

Read more »

Stationarity of ARCH processes

April 6, 2014
By
Stationarity of ARCH processes

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

Read more »

Choose Your Own Data Adventure

April 5, 2014
By
Choose Your Own Data Adventure

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

Read more »

Regressions with Multiple Fixed Effects – Comparing Stata and R

April 5, 2014
By
Regressions with Multiple Fixed Effects – Comparing Stata and R

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

Read more »

Two R tutorials for beginners

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

Read more »

Inference for ARCH processes

April 2, 2014
By
Inference for ARCH processes

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

Read more »