GOL.R

counts <- read.csv('counts.csv', header = TRUE, sep = '\t')
counts <- transform(counts, logCount = log(Count))

ggplot(counts, aes(x = Length, y = logCount)) +
  geom_point() +
  geom_smooth(method = 'lm') +
  xlab('Length of Entire Word') +
  ylab('Log Number of Occurrences of Spelling') +
  opts(title = 'How Many O\'s Does It Take to Make a GOL?')

# C = a * e ^ l
# => log(C) = log(a) + l

lm.fit <- lm(logCount ~ Length, data = counts)
summary(lm.fit)