setwd("H:/BigLake/Paper/Figure2")

Data <- read.csv("Figure2Data.csv", header = TRUE)

mean.pud = mean(Data$summerpud, na.rm = T)
var.pud = var(Data$summerpud, na.rm = T)
sd.pud = sd(Data$summerpud, na.rm = T)
min.pud = min(Data$summerpud, na.rm = T)
max.pud = max(Data$summerpud, na.rm = T)
num.obs = nrow(Data)
list(mean.pud,sd.pud,min.pud,max.pud,num.obs)
a <- table(Data$summerpud)
list(a)

DataPlus <- ifelse(Data$summerpud<10,Data$summerpud,10)

b <- table(DataPlus)
list(b)

barplot(b,ylim=c(0,50000),xlab="Summer Photo-User Days (2005 to 2014)",
        ylab="Frequency")