Friday, August 19, 2011

It warms my heart

It really warms my heart when things go right in the world.  Twenty-eight years and $1 million is about right for a judge preying on vulnerable children.  This part gets me:
The sentence was four times the 87 months sentence that Ciavarella and federal prosecutors had agreed to when he pleaded guilty to charges in 2009.
The prosecutors had a good case against this guy and they wanted to agree to _what_!?!?!? Many thanks to the judge who threw out that plea-deal.

Saturday, August 13, 2011

Improving statistical understanding in biology: teh ANOVA, it's robust to stuff

Recently there was a blog post at Expansed about how R might improve the practice of statistics.  The basic idea is that when you have to face the code which implements your statistical analysis you might think it through a little more than clicking.  I think that logic fail pretty quickly, faster than you can say cut-and-paste.  The practice of cut-and-pasting R code for a whole analysis (or worse, BUGS code!) is common even with practitioners who have worked with those tools for multiple years.  Sooner or later it catches up with you and you might learn something, but it's not inevitable.

For kicks, I coded up the following to express my dismay at the state of art:

files <- dir()
output <- list()
for ( f in files ) {
   dat <- read.csv(f)
   output[[f]] <- lm( y ~ ., data = dat )
   test <- summary(output[[f]])
   ps <- test$coefficients[,4]
   if (any(ps < 0.05)) {
      cat("Eureka! There's gold in ", f, ".", sep = "")

In a few weeks I'll try to code up something that uses stepwise significance testing to do variable selection.  Now how long before this makes it into somebody’s analysis...? Drop me a line if you use it!