Saturday, September 7, 2013

Trying out Sublime Text 3

I've been a huge fan of TextWrangler for years. I use it for all of my coding (including with R), for taking notes during meetings, for taking notes on articles, and more. It's my most frequently used application. But, I'm ready for a change and am going to give Sublime Text 3 a try. It seems very powerful, relatively light, and incredibly extensible. 

I recently installed Sublime Text 3, added Package Control, and installed some R packages to try out. I'll give it a try for the next three weeks and see what I think. It might be time to say farewell to TextWrangler.

Monday, May 20, 2013

Recurrence Quantification Analysis

This page has some wonderful resources for recurrence analysis. One particularly useful resource on this site is the listing of software options for conducting recurrence analysis. After a fair amount of searching, I couldn't find an R package that computed the metrics from a recurrence quantification analysis. The tseriesChaos package provides a function for producing recurrence plots; but, I didn't see anything for quantifying these plots.

After digging through the different software options listed on this site, I tried out and really like the Commandline Recurrence Plots script offered by Norbert Marwan himself.

The script was very easy to setup on my Mac and, by using Rscript it was easy to combine with R code to (a) draw specific chunks of data for different individuals in my dataset; (b) compute and output the recurrence quantification metrics; (c) output the recurrence plot dataset for creating the actual plot; and, (d) producing the plot and creating a dataset of metrics.

I'll clean up, comment, and post the code that I used as soon as I can come up for air.

Monday, April 15, 2013

R Graphics Parameters -- Rows and Columns

For some reason I always forget the code for setting R's graphics parameters. And, I always need this same line. So, now I shan't forget it.


quartz(type="pdf",file="figure_NUM.pdf")
par(mfrow=c(3,2), cex=1, mar=c(2,2,2,2))
dev.off()