Daily ride counts

library(dplyr)
library(ggplot2)

# Summarize the number of taxi rides by pickup day

daily_count <- tx %>%
group_by(pickup_date) %>%
summarise(n_rides = n())

# Create a line plot

ggplot(daily_count, aes(pickup_date, n_rides)) +
geom_line()

Distribution of cab fare amount

library(ggplot2)

# Create a histogram of total_amount

ggplot(tx, aes(total_amount)) +
geom_histogram() +
scale_x_log10()

Distribution of payment type

library(ggplot2)

# Create a bar chart of payment_type

ggplot(tx, aes(payment_type)) +
geom_bar()

Relationship between trip duration and total fare

library(ggplot2)

# Create a hexagon-binned plot of total_amount vs. trip_duration

ggplot(tx, aes(trip_duration, total_amount)) +
geom_hex(bins = 75) +
scale_x_log10() +
scale_y_log10()

Faceting daily rides

library(dplyr)
library(ggplot2)

# Summarize taxi rides count by payment type, pickup date, pickup day of week
daily_count <- tx %>%
filter(payment_type %in% c("Card", "Cash")) %>%
group_by(pickup_date, pickup_dow, payment_type) %>%
summarise(n_rides = n())

# Plot the data

ggplot(daily_count, aes(pickup_date, n_rides)) +
geom_point() +
facet_grid(payment_type ~ pickup_dow) +
coord_fixed(ratio = 0.4)

Tip amount distribution faceted by payment type

library(ggplot2)

# Histogram of the tip amount faceted on payment type

ggplot(tx, aes(tip_amount + 0.01)) +
geom_histogram() +
 scale_x_log10() +
facet_wrap(~ payment_type, ncol = 1, scales = "free_y")

Comparing fare distribution by payment type

library(ggplot2)
library(dplyr)
library(tidyr)

# Get data ready to plot

amount_compare <- tx_pop %>%
mutate(total_no_tip = total_amount - tip_amount) %>%
select(total_amount, total_no_tip, payment_type) %>%
gather(amount_type, amount, -payment_type)

# Quantile plot
ggplot(amount_compare, aes(sample = amount, color = payment_type)) +
geom_qq(distribution = stats::qunif, shape = 21) +
facet_wrap(~ amount_type) +
ylim(c(3, 20))

Trelliscope faceting gapminder by country

library(ggplot2)
library(trelliscopejs)

# Create the plot

ggplot(gapminder, aes(year, lifeExp)) +
geom_line() +
# Facet on country and continent
facet_trelliscope(~ country + continent)
///////////////////
library(ggplot2)
library(trelliscopejs)

ggplot(gapminder, aes(year, lifeExp)) +

geom_line() +
facet_trelliscope(~ country + continent,
name = "lifeExp_by_country",
desc = "Life expectancy vs. year per country",
nrow = 1, ncol = 2)

//3

Customizing the gapminder display

library(trelliscopejs)
library(ggplot2)

# Create the plot

ggplot(gapminder, aes(year, lifeExp)) +
geom_point() +
geom_smooth(method = "lm", se = FALSE) +
facet_trelliscope(~ country + continent,
name = "lifeExp_by_country",
desc = "Life expectancy vs. year for 142 countries.",
nrow = 1, ncol = 2,
# Set the scales
scales = "sliced",
# Specify automatic cognostics
auto_cog = TRUE)
//3

Adding custom cognostics

library(ggplot2)
library(dplyr)
library(gapminder)
library(trelliscopejs)
space_to_dash <- function(x) gsub(" ", "-", x)

# Group by country and create the two new variables

gap <- gapminder %>%
group_by(country) %>%
mutate(
delta_lifeExp = tail(lifeExp, 1) - head(lifeExp, 1),
ihme_link = paste0("http://www.healthdata.org/", space_to_dash(country)))

# Add the description

gap$delta_lifeExp <- cog(gap$delta_lifeExp, desc = "Overall change in life
expectancy")
# Specify the default label
gap$ihme_link <- cog(gap$ihme_link, default_label = TRUE)

ggplot(gap, aes(year, lifeExp)) +

geom_point() +
facet_trelliscope(~ country + continent,
name = "lifeExp_by_country",
desc = "Life expectancy vs. year.",
nrow = 1, ncol = 2,
scales = c("same", "sliced"))

//1

Grouping and nesting

library(dplyr)
library(tidyr)

# Nest stocks by symbol

by_symbol <- stocks %>%
group_by(symbol) %>%
nest()
///////////
//////////2
library(dplyr)
library(tidyr)
library(purrr)

# Nest stocks by symbol

by_symbol <- stocks %>%
group_by(symbol) %>%
nest()

min_volume_fn <- function(x) min(x$volume)

# Create new column

by_symbol_min <- by_symbol %>%
mutate(min_volume = map_dbl(data, min_volume_fn))

Stock price display

library(trelliscopejs)
 library(dplyr)
library(purrr)
library(plotly)

ohlc_plot <- function(d) {

plot_ly(d, x = ~date, type = "ohlc",
open = ~open, close = ~close,
high = ~high, low = ~low)
}

by_symbol_plot <- by_symbol %>%

mutate(panel = map_plot(data, ohlc_plot))
////////////////
//3
////////////////
library(trelliscopejs)
library(dplyr)
library(purrr)
library(plotly)

ohlc_plot <- function(d) {

plot_ly(d, x = ~date, type = "ohlc",
open = ~open, close = ~close,
high = ~high, low = ~low)
}

by_symbol_plot <- by_symbol %>%

mutate(panel = map_plot(data, ohlc_plot))

trelliscope(by_symbol_plot, name = "ohlc_top500")

//2

Adding cognostics
library(trelliscopejs)
library(dplyr)

# Create market_cap_log
by_symbol <- mutate(by_symbol,
market_cap_log = cog(
val = log10(market_cap),
desc = "log base 10 market capitalization"
)
)

Cognostics from nested data frames

library(trelliscopejs)
library(dplyr)
library(purrr)

annual_return <- function(x)

100 * (tail(x$close, 1) - head(x$open, 1)) / head(x$open, 1)

# Compute by_symbol_avg
by_symbol_avg <- mutate(by_symbol,
stats = map(data, function(x) {
data_frame(
mean_close = mean(x$close),
 mean_volume = mean(x$volume),
annual_return = annual_return(x)
)
}))

//4

Customizing the stock display

library(trelliscopejs)

# Create the trelliscope display

p <- trelliscope(
by_symbol,
name = "ohlc_top500",
desc = "a plot",
width = 600,
height = 300
)

# View the display

p

Visualizing Pokemon
library(dplyr)
library(trelliscopejs)

# Create a new data frame to use for plotting

pokemon2 <- pokemon %>%
# Reduce the variables in the dataset
select(pokemon, type_1, attack, generation_id, url_image) %>%
mutate(
# Respecify pokemon
pokemon = cog(val = pokemon, default_label = TRUE),
# Create panel variable
panel = img_panel(url_image)
)

# Create the display

trelliscope(pokemon2, name = "pokemon", nrow = 3, ncol = 6)

//1

Number of daily rides

library(dplyr)
library(ggplot2)

# Compute daily counts

daily <- bike %>%
group_by(start_day, weekday) %>%
summarise(n = n())

# Plot the result

ggplot(daily, aes(start_day, n, color = weekday)) +
geom_point()

Examining time-of-day
library(dplyr)
library(ggplot2)
 # Compute week_hod
week_hod <- bike %>%
group_by(start_wk, start_hod, weekday) %>%
summarise(n = n())

# Plot the result

ggplot(week_hod, aes(start_wk, n, color = weekday)) +
geom_point() +
facet_grid(~ start_hod) +
scale_y_sqrt()

Effect of membership and weekday

library(dplyr)
library(ggplot2)

# Compute wk_memb_hod
wk_memb_hod <- bike %>%
group_by(start_wk, start_hod, weekday, membership) %>%
summarise(n = n())

# Plot the result

ggplot(wk_memb_hod, aes(start_wk, n, color = weekday)) +
geom_point() +
facet_grid(membership ~ start_hod) +
scale_y_sqrt()

Daily plots
library(dplyr)
library(ggplot2)

# Compute daily_may
daily_may <- bike %>%
filter(start_mon == 5) %>%
group_by(start_day, start_hod, membership) %>%
summarise(n = n())

# Plot the result

ggplot(daily_may, aes(start_hod, n, color = membership)) +
geom_point() +
facet_wrap(~ start_day, ncol = 7)

Looking at all days

library(ggplot2)

ggplot(daily_may, aes(start_hod, n, color = membership)) +

geom_point() +
# Facet on start_day
facet_trelliscope(~ start_day, nrow = 3, ncol = 7)

Augmenting the Data: Route Summary Statistics

library(trelliscopejs)
library(ggplot2)
library(dplyr)

# Function to construct a Google maps URL with cycling directions

make_gmap_url <- function(start_lat, start_lon, end_lat, end_lon) {
paste0("https://www.google.com/maps/dir/?api=1",
"&origin=", start_lat, ",", start_lon,
 "&destination=", end_lat, ",", end_lon,
"&travelmode=bicycling")
}

# Compute tot_rides, weekday_diff, and map_url

route_hod_updated <- route_hod %>%
group_by(start_station_code, end_station_code) %>%
mutate(
tot_rides = sum(n),
weekday_diff = mean(n[weekday == "workweek"]) - mean(n[weekday ==
"weekend"]),
map_url = make_gmap_url(start_lat, start_lon, end_lat, end_lon))

Visualizing the Data: Counts by Hour-of-Day

library(trelliscopejs)
library(ggplot2)

# Create the plot

ggplot(route_hod, aes(start_hod, n, color = weekday)) +
geom_point(size = 3) +
facet_trelliscope(~ start_station_name + end_station_name, nrow = 2, ncol =
4) +
theme(legend.position = "none")

//3