cpfiffer

# Misc functions for conservation biology stats in R.

shannon <- function (x) {
  proportion_x = x/sum(x)
  print(proportion_x)
  print(log(proportion_x))
  print(proportion_x * log(proportion_x))
  print(sum(proportion_x * log(proportion_x)))
  return(-sum(proportion_x*log(proportion_x)))
}

cpfiffer

# Cameron Pfiffer
# July 31st, 2017
# Julia Does Numbers 2: The Mathening

# See the link below for information about the challenge.
# https://fivethirtyeight.com/features/pick-a-number-any-number/

#=
From Itay Bavly, a chain-link number problem:

cpfiffer

#The riddle can be found at https://fivethirtyeight.com/features/how-long-will-it-take-to-blow-out-the-birthday-candles/

tests = 10000000
blows <- rep(0,tests)
i = 1

while (i <=tests) {
  candles = 30
  
  while (candles > 0) {