--- title: 'MATH 301: Activity 8' author: "Brad Thiessen" date: "October 26, 2014" output: html_document --- ***** ###### Load packages ```{r, message=FALSE} library(mosaic) library(dplyr) library(ggplot2) library(ggvis) library(parallel) ```
#### Chi-square goodness-of-fit: Benford's Law ```{r , message=FALSE} ## Load Benford's Law package and ## Sino Forest Data library(benford.analysis) sino.forest <- read.csv("http://www.bradthiessen.com/html5/data/sinoforest.csv") ## Extract first digit first_digit <- extract.digits(sino.forest$value, number.of.digits = 1, sign="positive", second.order = FALSE, discrete=TRUE, round=3) ## Get frequencies first_digit %>% group_by(data.digits) %>% summarize(n=n(), rel_freq=n()/772) # Run the chi-squared analysis benford <- benford(first_digit$data.digits, 1, sign="both") #generates benford object benford #prints plot(benford) #plots ## Roll a die 120 times die <- sample(1:6, 120, replace=T) table(die) ```