# inst/doc/eddington.R In eddington: Compute a Cyclist's Eddington Number

```## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk\$set(
collapse = TRUE,
comment = "#>",
fig.width = 7,
fig.height = 5,
fig.align = "center"
)

# Putting this here too so that we don't have pkg startup messages.
library(dplyr)

## ----setup--------------------------------------------------------------------
library(eddington)

## ----xform--------------------------------------------------------------------
library(dplyr)

days <- rides %>%
group_by(ride_date) %>%
summarize(n = n(), total = sum(ride_length))

## ----summary------------------------------------------------------------------
summary(days)

## ---- echo=FALSE--------------------------------------------------------------
hist(
as.integer(days\$total),
breaks = 30,
main = "Histogram of Daily Mileages",
xlab = "Miles"
)

abline(v = E_num(days\$total), col = "darkred")

legend(
"topright",
legend = "Eddington Number",
col = "darkred",
bty = "n",
lty = 1L
)

## ----enum---------------------------------------------------------------------
E_num(days\$total)

## ----ecum---------------------------------------------------------------------
days\$E <- E_cum(days\$total)

## ----needle, echo=FALSE-------------------------------------------------------
E <- E_num(days\$total)

E_contribs <- days[days\$total >= E,]

plot(
y = days\$total,
x = days\$ride_date,
type = "h",
main = "Ride Mileages in 2009",
xlab = "Ride Day",
ylab = "Total Miles",
bty = "n",
ylim = c(0, 90)
)

lines(
y = c(0, days\$E),
x = c(as.Date("2009-01-01"), days\$ride_date),
type = "s",
col = "darkred"
)

abline(h = E, lty = 2L, col = "darkred")

text(
E_contribs[,c("ride_date","total")],
labels = as.integer(E_contribs[["total"]]),
pos = 3,
cex = 0.7
)

legend(
"topleft",
title = "Eddington Number",
legend = c("Cumulative", "Summary"),
col = "darkred",
bty = "n",
lty = c(1L, 2L)
)

## ----enext--------------------------------------------------------------------
E_next(days\$total)

## ----ereq---------------------------------------------------------------------
E_req(days\$total, 50)

## ----esat---------------------------------------------------------------------
E_sat(days\$total, 30)
```

## Try the eddington package in your browser

Any scripts or data that you put into this service are public.

eddington documentation built on March 24, 2020, 5:08 p.m.