R posts in other languages
MLB Gameday stores its game data in XML format, with the players denoted in ID numbers. To find out who is who, the codes are stored in pitchers.xml or batters.xml of each game. My DownloadPitchFX.R script can download the ID numbers, but it doesn’t look to see who the ID is because of the extra 
Good long weekend... Ready to get back into R and data mining...After a personal appointment this morning, I'm now reading an early draft of Data Mining with R: Learning with Case Studies (Chapman & Hall/CRC Data Mining and Knowledge Discovery Seri...
Maybe the title should have been prepended with “Don’t…” The source code of R is littered with “attribute_hidden” declarations. These declarations attempt to ensure that the variable or function may only be accessed by code in the core R distribution, and not by R extension packages. Generally there is a good reason for this. For
![Vanilla Rao-Blackwellisation [re]revised](http://l.wordpress.com/latex.php?latex=%5Cleft%7C+%5Cbegin%7Bmatrix%7D+%5Ctau+%26%5Ctext%7Bmedian%7D+%26%5Ctext%7Bmean+%7D%26q_%7B.8%7D+%26q_%7B.9%7D+%26%5Ctext%7Btime%7D%5C%5C+0.25+%260.0+%268.85+%264.9+%2613+%264.2%5C%5C+0.50+%260.0+%266.76+%264+%2611+%262.25%5C%5C+1.00+%260.25+%266.15+%264+%2610+%262.5%5C%5C+2.00+%260.20+%265.90+%263.5+%268.5+%264.5%5C%5C%5Cend%7Bmatrix%7D+%5Cright%7C&bg=000000&fg=B0B0B0&s=0)
Although the revision is quite minor, it took us two months to complete from the time I received the news in the Atlanta airport lounge… The vanilla Rao-Blackwellisation paper with Randal Douc has thus been resubmitted to the Annals of Statistics. And rearXived. The only significant change is the inclusion of two tables detailing computing 
After downloading some pitch f/x data using my R script, we can finally have some fun. But because the pitch f/x data is very elaborate, R can easily get overwhelmed by copying the dataset back and forth in memory, as you manipulate the data. So the natural progression is to use relational database systems. Here, 
It's official: R 2.11.1 is out. Source code and binaries for Windows and MacOS are available at the master CRAN mirror, and will be available for download from your local mirror soon. As anticipated, this is an update release focussing mainly on bugfixes and with just one new feature. According to R core team member Peter Dalgaard one fix...
In light of my recent attempt at aRt, Tal from R bloggers suggested I submit a T-shirt design for this contest. That got me thinking that R needs a logo freshining in general, so I dusted off my technical pens and drafted something. I’ll explain why I think this makes a good logo for R,
I've released version 0.2-0 of highlight to CRAN This version brings some more additions to the sweave driver that uses highlight to produce nice looking vignettes with color coded R chunks The driver gains new arguments boxes, bg and border to c...
Memorial Day weekend is also time for the annual Bike The Drive in Chicago. This time only half the household got up bright and early and enjoyed Lakeshore Drive free of cars. A highly recommended event.
I was reading over at math-blog.com about a concept called numeri ritardatari. This sounds a lot like “retarded numbers” in Italian, but apparently “retarded” here is used in the sense of “late” or “behind” and not in the short bus sense. I barely scanned the page, but I think I got the gist of it:
Details matter (at least, they do for me), but we don't yet have a systematic way of going back and forth between the structure of a graph, its details, and the underlying questions that motivate our visualizations. (Cleveland, Wilkinson, and....
It's Memorial Day weekend so it was time for the Chicago's JP Morgan Chase Corporate Challenge on Thursday. The weather was glorious, the usual 20-some thousand runners participated and a good time was had. Work had arranged for a nice tent, food, mu...
In our previous example, we demonstrated how to calculate the Kaplan-Meier estimate of the survival function for time to event data. A related quantity is the Nelson-Aalen estimate of cumulative hazard. In addition to summarizing the hazard incurred ...
In the GCaP class earlier this month, we talked about the meaning of the load average (in Unix and Linux) and simulating a grocery store checkout lane, but I didn't actually do it. So, I decided to take a shot at constructing a discrete-event simulatio...
This is the talk I am giving at the workshop on model uncertainty organised by the Centre for Research in Statistical Methodology (CRiSM) at the University of Warwick, on May 30-June 1. Careful readers will notice there is not much difference with my previous talk on the topic, as I only included the Savage-Dickey slides 
(This is the second of a series of ongoing posts on using Graph Algebra in the Social Sciences.) First-order linear difference equations are powerful, yet simple modeling tools. They can provide access to useful substantive insights to real-world phenomena. They can have powerful predictive ability when used appropriately. Additionally, they may be classified in any number
Notice that even though output is in a log scale, output is shooting up in an exponential way. DATA from Brad DeLong
R's interactive programming style is similar to what I have seen in other environments (e.g. ruby's irb and Oracle's SQL*Plus, etc). There are a few commands that you need to be aware of to get up and running with developing R programs.To identify yo...
As usual click on the image for a full-size version. Code: par(bg="black") par(mar=c(0,0,0,0)) plot(c(0,1),c(0,1),col="white",pch=".",xlim=c(0,1),ylim=c(0,1)) iters = 500 for(i in 1:iters) { center = runif(2) size = rbeta(2,1,50) # Let's create random HTML-style colors color = sample(c(0:9,"A","B","C","D","E","F"),12,replace=T) fill = paste("#", paste(color[1:6],collapse=""),sep="") brdr = paste("#", paste(color[7:12],collapse=""),sep="") rect(center[1]-size[1], center[2]-size[2], center[1]+size[1], center[2]+size[2], col=fill, border=brdr, density=NA, lwd=1.5) }
The version 0.1-8 of highlight introduced a small bug in the latex renderer. This is now fixed in version 0.1-9 and the latex renderer also gains an argument "minipage" which wraps the latex code in a minipage environment. I've used this to make...
Motivation In the past few months I have been using DropBox for syncing my work files between my home and work computer. It has saved me from numerous mistakes and from sending the files to myself via e-mail. Recently I found this service highly useful for sharing files with 4 other people with whom I am working on a...
Western Fresno Soil Hierarchy: partial view of the hierarchy within the US Soil Taxonomic system Maybe this is just craziness, but wouldn't be neat to have an XML formatted version of the Keys to Soil Taxonomy? The format might look something like the ...
Coded as follows: x = seq(.001,50,.001) par(bg="black") par(mar=c(0,0,0,0)) plot(x,sin(1/x)*rpois(length(x),x),pch=20,col="blue")
COUNTERINTUITIVE PROBLEM, INTUITIVE REPRESENTATION Blog posts about counterintuitive probability problems generate lots of opinions with a high probability. Andrew Gelman and readers have been having a lot of fun with the following probability problem: I have two children. One is a boy born on a Tuesday. What is the probability I have two boys? The
(This is the first in a series on the use of Graph Algebraic models for Social Science.) Linear Difference models are a hugely important first step in learning Graph Algebraic modeling. That said, linear difference equations are a completely independent thing from Graph Algebra. I’ll get into the Graph algebra stuff in the next post or
Having worked with Unix (BSD, HPUX, IRIX, Linux and OSX), Windows (NT4, 2000, XP, Vista and 7) for quite a while I have seen a lot of different software tools. I would like to quickly exhibit my “must have” list. These are the packages that I find to be the single “must have offerings” in Related posts:
A 3d wireframe plot is a type of graph that is used to display a surface – geographic data is an example of where this type of graph would be used or it could be used to display a fitted model with more than one explanatory variable. These plots are related to contour plots which
According to internet lore, there's a mathematical equation that governs the lower bound for the socially acceptable age of a potential dating partner: half your age plus 7, or, in mathematical terms, if x is your age then the lower bound is f(x) = x/2 + 7. Seems simple, right? if you're 20, then the minimum socially acceptable age...
Zero Inflated Models and Generalized Linear Mixed Models with R.
Zuur, Saveliev, Ieno (2012).
