## Ecological networks from abundance distributions

January 6, 2011
By

Another grad student and I tried recently to make a contribution to our understanding of the relationship between ecological network structure (e.g., nestedness) and community structure (e.g., evenness)......Alas, I had no luck making new insights. How...

## Graph gallery in R

January 6, 2011
By

R is sometime criticized for producing graphs not as elaborated as Matlab ones, or other softwares’. Here is a link to a graph gallery by Romain François to “enhance your data visualization with R”. The corresponding R code is given. Might be useful for ENSAE students for ‘statap’ projects. Below are four examples. The maps

## formatR update (0.1-6)

January 6, 2011
By

A new version of the formatR package is available on CRAN now (binary packages are still on the way). There are three major updates: the inline comments will also be preserved in most cases (in earlier versions, only single lines of comments are pres...

## Learning R — Documentation

January 6, 2011
By

We use R a lot.  R takes care of many our basic data management needs.  R is an awesome statistical analysis package.  R allows you to produce exceptional data graphics.  The only problem is … R has a wicked learning …   read ...

## Short review of the R book

January 5, 2011
By

David Scott wrote a review of Introducing Monte Carlo Methods with R in the International Statistical Review that is rather negative, since the main bulk reads as follows: I found some aspects of the book very disappointing. The first chapter (“Basic R Programming”) has some unfortunate mistakes and some statements, which are contentious at least

## My first R package: zipcode

January 5, 2011
By

My first package, zipcode, is now available on CRAN. It contains the CivicSpace database of 43,191 U.S. zip codes.

## sab-R-metrics: Introduction to R

January 5, 2011
By

In a recent post, I briefly mentioned that I may turn a majority of the focus of this blog to teaching R commands for use with sabermetric analysis. Only a few days later, Ricky Zanker began a new column at The Hardball Times doing just that. But that's okay. Hopefully both his and mine...

## sab-R-metrics: Introduction to R

January 5, 2011
By

In a recent post, I briefly mentioned that I may turn a majority of the focus of this blog to teaching R commands for use with sabermetric analysis. Only a few days later, Ricky Zanker began a new column at The Hardball Times doing just that. But that's okay. Hopefully both his and mine...

## R-bloggers

January 5, 2011
By

Just a quick FYI note in case you haven't seen this site.R-bloggers is an awesome site, bringing together more than 140 blogs (including mine) about R in a single location. See Tal Galili's motivation for creating the site, and his notes on the site here.

## New approach to analysis of phylogenetic community structure

January 5, 2011
By

Anthony Ives, of University of Wisconsin-Madison, and Matthew Helmus of the Xishuangbanna Tropical Botanical Garden, present a new statistical method for analyzing phylogenetic community structure in an early view paper in Ecological Monographs. See th...

January 5, 2011
By

Another exciting day at Adap’skiii!!! Yves Atchadé presented a very recent work on the fundamental issue of estimating the asymptotic variance estimation for adaptive MCMC algorithms, with an intriguing experimental observation that a non-converging bandwidth with rate 1/n was providing better coverage than the converging rate. (I always found the issue of estimating the asymptotic

## How many gifts did my true love give to me on all twelve nights…

January 5, 2011
By

How many gifts did my true love give to me on all twelve nights of Christmas? After seeing Information is Beautiful’s recent information animation, I decided that I’d make my own. I used R to generate a PDF and then did a screencast of the PDF with ffmpeg to make a video with the appropriate timings. The code for the PDF...

## Customizing the Theme of Your R HTML Help

January 4, 2011
By

R’s default theme of the HTML help pages is too plain for me to read, but we can easily modify the theme, which is essentially a CSS file. You can find the file under: file.path(R.home('doc'), 'html', 'R.css') Simply replace this file with my version: which looks like: Of course you can design your own R.css

## The R Journal: December 2010

January 4, 2011
By

Issue 2 of The R Journal (the peer-reviewed journal devoted to R) was published over the Christmas break. In addition to news about the latest release of R, it also includes contributed articles on using GPU processing to fit Bayesian models in R, processing text data in R, solving differential equations in R, and much more. Follow the link...

## Revised market prediction distributions

January 4, 2011
By

This provides revised plots of the prediction distributions published yesterday.  The previous plots of prediction distributions should be ignored — they are not doing as advertised. We show the prediction distribution of levels of several equity indices (plus oil price) at the end of 2011 assuming nothing happens.  That is, we’ve taken out market trends … Continue reading...

## Creating prediction distributions

January 4, 2011
By

Here we give details and code for the prediction distributions exhibited in yesterday’s blog post Tis the season to predict. Eight years of returns The equity indices use daily closing levels from the start of 2003.  This data comes from Yahoo. A roughly equivalent technique of selecting the last 2000 daily prices is used for … Continue reading...

## Customizing the Theme of Your R HTML Help

January 4, 2011
By

R's default theme of the HTML help pages is too plain for me to read, but we can easily modify the theme, which is essentially a CSS file. You can find the file under: file.path(R.home('doc'), 'html', 'R.css') Simply replac...

## How many of each gift did my true love give to me?

January 4, 2011
By

How many of each gift did my true love give to me?

## Gifts my true love gave to me on all twelve days of Christmas

January 4, 2011
By

Gifts my true love gave to me on all twelve days of Christmas

## First meeting of Toronto R User Group this Friday

January 3, 2011
By

The Greater Toronto Area R User Group is having its inaugural meeting this Friday, January 7. The meeting will feature Ben Bolker, professor and author of Ecological Models and Data in R, who will speak on Generalized Linear Models. It's also your opportunity to join the group and help define its future path, so if you're in the area...

## Example 8.19: Referencing lists of variables

January 3, 2011
By

In section 1.11.4 (p. 50), we discuss referring to lists of variables in a data set. In SAS, this can be done for variable stored in adjacent columns with the "var_x -- var_y" syntax and for variables with sequentially enumerated suffixes with the "va...

## Exporting KML from R

January 2, 2011
By

Google Earth has become a popular way of disseminating spatial data. KML is the data format required to do this. It is possible to load almost any type of spatial...

## LaTeX Typesetting – Tables

January 2, 2011
By

One area of LaTeX that is not always straightforward is producing tabular material. In this post we will look at some of the basics of creating tables in a LaTeX document. Fast Tube by Casper When creating tables there is a tabular environment that can be used and we need to specify the number of

## Animation from Kaggle traffic prediction competition

January 2, 2011
By

Below: two animations (same one twice, just as different speeds) I made from the highway data for the Kaggle traffic prediction competition.The data are travel times for 61 different sections of Sydney's M4 freeway, in both directions. The data comes i...

## Animation from Kaggle traffic prediction competition

January 2, 2011
By

Below: two animations (same one twice, just as different speeds) I made from the highway data for the Kaggle traffic prediction competition.The data are travel times for 61 different sections of Sydney's M4 freeway, in both directions. The data comes i...

## Le Monde puzzle [52|solution]

January 1, 2011
By
$Le Monde puzzle [52|solution]$

I have now received the first issue of Le Monde magazine, including the solution to puzzle #52 I solved just in time by simulated annealing! The trick is in using the following theorem: Iter(1,x,y) is divisible by 10x-1 if and only if y is divisible by 10x-1. Then the value to be found is divisible

## Ten News Stories of 2010 – and the Statistics that Made Them

January 1, 2011
By

Significance MagazineAccording to Significance Magazine (jointly published by Royal Statistical Society and the American Statistical Association) the following are the top ten stories of 2010.1. Progress in the preven...

## Ten News Stories of 2010 – and the Statistics that Made Them

January 1, 2011
By

Significance MagazineAccording to Significance Magazine (jointly published by Royal Statistical Society and the American Statistical Association) the following are the top ten stories of 2010.1. Progress in the preven...