Benford's Law

gistfile1.txt

install.packages("devtools")  # so we can install from github
library("devtools")
install_github("ropensci/plotly")  # plotly is part of ropensci
library(plotly)

py <- plotly(username="r_user_guide", key="mw5isa4yqp")  # open plotly connection

library(ggplot2)
d <- 1:9
prob = log10(1 + (1/d))
benford_dist <- data.frame(type = "Benford", digit = 1:9, prob = prob)
q <- ggplot(benford_dist, aes(x = factor(digit), y = prob, fill = type))
q <- q + geom_bar(stat = "identity") + labs(title = "Benford's Law")

py$ggplotly()