April 1, 2012
By

Hadley Wickham (@hadleywickham) this week mentioned on Twitter his preference for saveRDS() over the more familiar save(). Being a new function to me, I thought I’d take a look… save() and load() will be familiar to many R users. They … Continue reading →

April 1, 2012
By

Hadley Wickham (@hadleywickham) this week mentioned on Twitter his preference for saveRDS() over the more familiar save(). Being a new function to me, I thought I’d take a look…

## Julia, I Love You

March 31, 2012
By

Julia is a new language for scientific computing that is winning praise from a slew of very smart people, including Harlan Harris, Chris Fonnesbeck, Douglas Bates, Vince Buffalo and Shane Conway. As a language, it has lofty design goals, which, if attained, will make it noticeably superior to Matlab, R and Python for scientific programming.

## Back to Blogging

March 31, 2012
By

If you’re subscribed to this blog, you’ve surely noticed the very long hiatus I’ve taken from writing over the last six months. I wish I’d kept up with blogging more faithfully this year, but, in my defense, I’ve been busy doing a few big things: I wrote a book with Drew Conway called Machine Learning

March 31, 2012
By

I've been doing more analysis of the Philadelphia Homicide data that the Philadelphia Inquirer has published, and presented some of it at the Philadelphia UseR group yesterday. My slides and source are on github.I should be clear tha...

## Draw Your Breast with CloudStat – A R Apps (for fun)

March 31, 2012
By

This is a simple apps, called “Draw Your Breast with R“  created with R to generate Breast alike graphics. With this Draw Your Breast with R apps, you can change 4 parameters which are Theta, Phi, Expand and Color to generate graphics like...

## Ggplot2, PubMed citation frequency and DSM-IV Axis I disorders by year

March 31, 2012
By

I searched PubMed for several major DSM-IV disorders and downloaded the hits. Using ggplot2 I plotted the number of publications each year for each disorder.

## Playing with XML-Package: Get No. of Google Search Hits with R

March 30, 2012
By

GoogleHits <- function(input) { require(XML) require(stringr) require(RCurl) url

## GBLUP example in R

March 30, 2012
By

Shirin Amiri was asking about GBLUP (genomic BLUP) and based on her example I set up the following R script to show how GBLUP works. Note that this is the so called marker model, where we estimate allele substitution effects of the markers and not individual based model, where genomic breeding values are inferred directly. The code:library(package="MatrixModels") dat <- data.frame(...

March 30, 2012
By

## See the wind

March 30, 2012
By

The image below isn't a bearskin rug in the shape of the USA. It's fact, it's a visualization of the wind flowing over the United States, as of 4PM EDT today, March 30. You can click through to see the current wind conditions, based on latest data from the National Digital Forecast Database. But more importantly, as long as...

## By: Super Nerdy Cool » Build 32 bit R on 64 bit Ubuntu by utilizing chroot

route, it will suffice for now. With the chroot method, first compile the 64 bit version of R the usual way. For the 32 bit version of R,

## R 2.15.0 "Easter Beagle" is released

March 30, 2012
By

The R core group today announced the availability of R 2.15.0, codenamed "Easter Beagle". If you build R yourself, the new source distribution (including packages for Debian Linux) is available now, and binaries for Windows, MacOS and Linux will be available from your local CRAN mirror over the next couple of days. (As of this writing, Windows binaries are...

## Build 32 bit R on 64 bit Ubuntu by utilizing chroot

March 30, 2012
By

In the past, I’ve described how one could build multiarch (64 bit and 32 bit) versions of R on a 64 bit Ubuntu machine. The method based on this thread no longer works as of R 2.13 or 2.14 I believe. I received advice from someone on #R over on freenode (forgot who) a few... Read more »

## Comparing Spectra with different math treatments

March 30, 2012
By

It is amazing the quantity of graphics you can develop with R, and how you can show and manage  these graphics. Here in the same plot we compare the raw demo spectra, treated with MSC, SNV and with the first derivative (differences between consecu...

## Slow R graphics via SSH

March 30, 2012
By

I do most of my R work on a server via SSH.  Viewing graphics (i.e. plot()) can be excruciatingly slow through the SSH pipe, but it's easy to make it fast.  Simply run:        X11.options(type="Xlib")once in the R session.  Tags: linuxgraphicsR

## Slow R graphics via SSH

March 30, 2012
By

I do most of my R work on a server via SSH.  Viewing graphics (i.e. plot()) can be excruciatingly slow through the SSH pipe, but it's easy to make it fast.  Simply run:        X11.options(type="Xlib")once in the R session....

## Maths can be cool (to impress your kids)

March 30, 2012
By

Just imagine that your kids need some help, to prepare fishes for April 1st, like Her: "please, Daddy, help us to draw some fishes" Me: "Sure, Daddy is a champion, actually, I do that everyday at work: drawing fishes - and more generally nice stuf...

## Mapped: British, Spanish and Dutch Shipping 1750-1800

March 30, 2012
By

I recently stumbled upon a fascinating dataset which co

## R 2.15.0 is released

March 30, 2012
By

Bellow is the announcement made by Peter Dalgaard: The build system rolled up R-2.15.0.tar.gz (codename “Easter Beagle”) at 9:00 this morning. This is the first release of the 2.15 series and contains several new features and changes; see the list below for details. You can get the source code from http://cran.r-project.org/src/base/R-2/R-2.15.0.tar.gz or wait for it to be mirrored at...

## R 2.15.0 is released!

March 30, 2012
By

The new R 2.15.0 is out! You know the drill! Get the source code from here.

## BoRdeaux (France): the first French R useR conference

March 30, 2012
By

The first French R useR conference will be held in Bordeaux on July 2-3. The Rencontres R are conceived as a place to present and share ideas on using the R statistical software. This meeting is designed to be a nationwide event where various topics belong, such as graphical tools, applied statistics, biostatistics, bayesian statistics,

## Bootstrap example

March 30, 2012
By

Bootstrap your way into robust inference. Wow, that was fun to write.. Introduction Say you made a simple regression, now you have your . You wish to know if it is significantly different from (say) zero. In general, people look … Continue reading →

## Slow R graphics via SSH

March 30, 2012
By

I do most of my R work on a server via SSH.  Viewing graphics (i.e. plot()) can be excruciatingly slow through the SSH pipe, but it's easy to make it fast.  Simply run:        X11.options(type="Xlib") once in the R session.  Tags:  R

## UPDATE: Install R 2.14.2 or R 2.15 in Debian Squeeze

March 29, 2012
By

Since not so long, I upgraded  the version of R that is included in the Debian repositories (2.11.1 in Squeeze) and I was very happy with it since 2010. But there was some packages that I could not install and then, I decided to upgrade R despite ...

## R-bloggers

March 29, 2012
By

OK, like the nightclub Kenan and Kel founded in Rigby’s shop, we all know this blog is ‘duh bomb’ (or will be as I add more than my solitary post). However, there are other good R blogs. I thoroughly recommend signing up to the R-b...

## See how Deloitte uses R for actuarial analysis

March 29, 2012
By

Many thanks to Jim Guszcza (Predictive Analytics lead at Deloitte Consulting and Assistant Professor at UW-Madison) who gave a great webinar presentation yesterday on actuarial analysis with R. Jim's demo (starting at the 20 minute mark in the video replay below) is a great way to get a sense of how R is used for exploratory data analysis and...

## Dividing the Sample Set in two (Validation & Training)

March 29, 2012
By

We have in the Demo sample set “66” samples.  In this post we´ll see one way to divide the set in two parts: one for “Validation” and another for Training or Calibration.The selection will be random. And we are going to use the command: