If you’re not using version control, you should be. Learn git. If you’re not on github, you should be. That’s real open source. To help some colleagues get started with git and github, I wrote a minimal tutorial. There are lots of git and github resources available, but I thought I’d give just the bare […]
There was lots of discussion on twitter yesterday about Graham Coop’s paper with Peter Ralph (or vice versa), on The geography of recent genetic ancestry across Europe, particularly regarding the FAQ they’d created. I was eager to take a look, and, it’s slightly embarrassing to say, I first did a search to see if they’d […]
It’s great that, in a step towards improved reproducibility, the Nature journals are removing page limits on Methods sections: To allow authors to describe their experimental designs and methods in enough detail for others to interpret and replicate them, the participating journals are removing length restrictions on Methods sections. But couldn’t they include the Methods […]
When physicists do mathematics, they don’t say they’re doing “number science”. They’re doing math. If you’re analyzing data, you’re doing statistics. You can call it data science or informatics or analytics or whatever, but it’s still statistics. If you say that one kind of data analysis is statistics and another kind is not, you’re not […]
Another R tip: beware of as.character applied to a list. Really, beware of grep with a list: You might have thought that the result would be just 1, but grep expects a vector of character strings. If the input is not that, it uses as.character(). Since the result of that starts with "c(", grep finds […]
It’s widely understood that, in R programming, one should avoid for loops and always try to use apply-type functions. But this isn’t entirely true. It may have been true for Splus, back in the day: As I recall, that had to do with the entire environment from each iteration being retained in memory. Here’s a […]
![x[[c(5,3)]]](http://feeds.wordpress.com/1.0/comments/kbroman.wordpress.com/1600/)
An R tip: Did you know that x[[c(5,3)]] is the same as x[[5]][[3]]? I should make more thorough use of this. In the help file for [[: [[ can be applied recursively to lists, so that if the single index i is a vector of length p, alist[[i]] is equivalent to alist[[i1]]...[[ip]] providing all but […]
The Statistics Department at UW-Madison has a course on statistical consulting, offered each semester. I’m often asked to give a lecture, which I do in an informal way: summarizing my experiences and answering questions. I thought it might be useful write my thoughts on statistical consulting here: why, how, and difficulties. This will be a […]
I’ve created another D3 example, of QTL analysis for a phenotype measured over time. (Click on the image for the interactive version.) The code is on github. It took me about a day. The hardest part was figuring out the right data structures. A pixel here is linked to curves over there and over there […]
I attended the ENAR meeting in Orlando, Florida, last week. (ENAR = “Eastern North American Region” of the International Biometric Society.) I had a great time, but I did come to the strong realization that what I view as important is distinctly different from what the typical ENAR attendee views as important. (Rafa said, incredulously, […]
