Case Study: Running activities
Megan Beckett
16 August 2018
Introduction
This follows on from investigating the lubridate package within the tidyverse and putting some of the features into practice by looking at running activity data over time.
I was interested to investigate some ways of plotting my running activities over time. The data is a CSV export from Garmin Connect.
Load packages needed
library(lubridate)
library(dplyr)
library(ggplot2)
library(skimr)
library(summarytools)
library(RColorBrewer)Read in data
activities <- read.csv("../data/Activities.csv", stringsAsFactors = FALSE)Initial view of data
Notes:
- I normally use a combination of str(), summary(), glimpse() to do this and recently explored skim().
- For skim(), I like how it splits the summary statistics across variable type and shows a basic histogram for numerical types.
- Missing values are also shown and this, along with the histograms, gives some indication of what to explore further.
- Subsequent to this, I have now used the package summarytools and specifically, dfSummary(). The histogram plots are even better (although render quite small) and the developers are working to improve the HTML output to work with R Markdown.
- See this vignette: https://cran.r-project.org/web/packages/summarytools/vignettes/Recommendations-rmarkdown.html
dim(activities)[1] 660 42#str(activities)
#summary(activities)
glimpse(activities)Observations: 660
Variables: 42
$ Activity.Type <chr> "running", "running", "running", "running", "running", "running", "running", "running...
$ Date <chr> "2018-08-01 18:19:10", "2018-07-31 06:23:37", "2018-07-29 08:11:23", "2018-07-28 10:2...
$ Favorite <chr> "false", "false", "false", "false", "false", "false", "false", "false", "false", "fal...
$ Title <chr> "Cape Town Running", "Cape Town Running", "Cape Town Running", "Cape Town Running", "...
$ Distance <dbl> 10.27, 7.36, 24.92, 14.13, 7.21, 8.13, 10.82, 10.16, 1.68, 18.97, 10.23, 12.01, 9.07,...
$ Calories <chr> "351", "311", "1,731", "1,308", "243", "262", "472", "510", "125", "1,235", "610", "6...
$ Time <chr> "50:36", "31:40", "4:17:17", "2:49:24", "38:03", "40:48", "49:51", "51:09", "12:15", ...
$ Avg.HR <chr> "130", "159", "139", "150", "123", "127", "153", "159", "160", "156", "166", "164", "...
$ Max.HR <chr> "154", "180", "185", "188", "141", "157", "178", "174", "173", "179", "190", "187", "...
$ Aerobic.TE <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
$ Avg.Run.Cadence <chr> "168", "180", "140", "136", "172", "170", "172", "168", "152", "152", "170", "168", "...
$ Max.Run.Cadence <chr> "180", "194", "248", "248", "178", "178", "248", "176", "168", "246", "184", "182", "...
$ Avg.Pace <chr> "4:56", "4:18", "10:20", "11:59", "5:17", "5:01", "4:36", "5:02", "7:18", "7:49", "5:...
$ Best.Pace <chr> "1:19", "3:15", "4:11", "5:00", "4:16", "3:39", "3:03", "3:42", "5:00", "0:20", "3:49...
$ Elev.Gain <chr> "91", "7", "1,779", "1,190", "36", "42", "18", "91", "32", "1,288", "170", "206", "9"...
$ Elev.Loss <chr> "92", "6", "1,581", "1,130", "37", "42", "18", "89", "148", "1,179", "171", "207", "8...
$ Avg.Stride.Length <dbl> 1.21, 1.29, 0.68, 0.61, 1.10, 1.17, 1.25, 1.18, 0.90, 0.84, 1.13, 1.18, 1.33, 1.07, 1...
$ Avg.Vertical.Ratio <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
$ Avg.Vertical.Oscillation <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
$ Avg.Ground.Contact.Time <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Avg.GCT.Balance <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Avg.Run.Cadence.1 <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Max.Run.Cadence.1 <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Normalized.Power...NP.. <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ L.R.Balance <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Training.Stress.Score. <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
$ Max.Avg.Power..20.min. <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Power <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Max.Power <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Total.Strokes <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Avg..Swolf <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Avg.Stroke.Rate <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Max.Depth <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Bottom.Time <chr> "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00...
$ Min.Water.Temp <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
$ Gas.Type <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Surface.Interval <chr> "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00", "0:00...
$ Decompression <chr> "No", "No", "No", "No", "No", "No", "No", "No", "No", "No", "No", "No", "No", "No", "...
$ Weight <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Current <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Surface.Conditions <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...
$ Water.Type <chr> "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "--", "...#skim(activities)Tidy up data
Extract variables of interest
Many of the columns are empty, so let’s create a smaller subset of those we are interested in.
# Create subset of variables
running <- subset(activities, select = c(Activity.Type, Date, Distance, Calories, Time,
Avg.Pace, Elev.Gain))
# View the classes of the subset of data
lapply(running, class)$Activity.Type
[1] "character"
$Date
[1] "character"
$Distance
[1] "numeric"
$Calories
[1] "character"
$Time
[1] "character"
$Avg.Pace
[1] "character"
$Elev.Gain
[1] "character"Transform data types
read.csv imported some of the the numbers, dates and times as characters, so let’s fix this. Also, convert Activity.Type to factor to see if there is more than one level.
# Numeric values
# Substitute 1000 comma separator then convert to numeric
running$Calories <- as.numeric(gsub(",", "", running$Calories))
running$Elev.Gain <- as.numeric(gsub(",", "", running$Elev.Gain))
# Date-time columns (Date, Time, Avg.Pace, Best.Pace)
running$Date = ymd_hms(activities$Date)
# TODO Figure out how to convert different time durations to actual times
#running$Duration <- ms(running$Time)
#running$Avg.Pace <- ms(running$Avg.Pace)
#running$Best.Pace <- ms(running$Best.Pace)
# Convert activity type to factor
running$Activity.Type <- as.factor(running$Activity.Type)
levels(running$Activity.Type)[1] "running" "street_running" "track_running" "trail_running" View data again using summarytools package
view(dfSummary(running), method = 'render')Data Frame Summary
running
N: 660| No | Variable | Stats / Values | Freqs (% of Valid) | Graph | Valid | Missing |
|---|---|---|---|---|---|---|
| 1 | Activity.Type [factor] | 1. running 2. street_running 3. track_running 4. trail_running | 653 (98.9%) 2 (0.3%) 2 (0.3%) 3 (0.4%) |  |
660 (100%) | 0 (0%) |
| 2 | Date [POSIXct, POSIXt] | min : 2014-12-24 07:43:13 med : 2017-03-25 12:42:45 max : 2018-08-01 18:19:10 range : 3y 7m 8d 10H 35M 57S | 660 distinct val. |  |
660 (100%) | 0 (0%) |
| 3 | Distance [numeric] | mean (sd) : 10.43 (7.67) min < med < max : 0.23 < 8.14 < 65.6 IQR (CV) : 4.24 (0.74) | 465 distinct values |  |
660 (100%) | 0 (0%) |
| 4 | Calories [numeric] | mean (sd) : 627.37 (406.27) min < med < max : 14 < 516 < 3033 IQR (CV) : 281.5 (0.65) | 460 distinct values |  |
659 (99.85%) | 1 (0.15%) |
| 5 | Time [character] | 1. 30:21 2. 30:45 3. 39:04 4. 39:15 5. 40:52 6. 1:12:03 7. 1:13:37 8. 1:20:13 9. 1:26:24 10. 20:17 [ 577 others ] | 3 (0.4%) 3 (0.4%) 3 (0.4%) 3 (0.4%) 3 (0.4%) 2 (0.3%) 2 (0.3%) 2 (0.3%) 2 (0.3%) 2 (0.3%) 635 (95.2%) |  |
660 (100%) | 0 (0%) |
| 6 | Avg.Pace [character] | 1. 5:05 2. 5:14 3. 5:02 4. 4:56 5. 5:03 6. 5:08 7. 4:57 8. 5:07 9. 4:55 10. 5:01 [ 205 others ] | 19 (2.9%) 17 (2.6%) 16 (2.4%) 15 (2.3%) 15 (2.3%) 14 (2.1%) 13 (2.0%) 13 (2.0%) 12 (1.8%) 12 (1.8%) 514 (77.6%) |  |
660 (100%) | 0 (0%) |
| 7 | Elev.Gain [numeric] | mean (sd) : 194.71 (266.72) min < med < max : 1 < 87 < 1779 IQR (CV) : 202 (1.37) | 311 distinct values |  |
645 (97.73%) | 15 (2.27%) |
Generated by summarytools 0.8.7 (R version 3.4.4)
2018-09-15
Subset again to only look at last 3 full years of data
running <- subset(running, running$Date >= "2016-01-01")Create new variables
running$year <- year(running$Date)
running$month <- month(running$Date, label=TRUE)
running$week <- week(running$Date)
running$wday <- wday(running$Date, label=TRUE, abbr=FALSE)
running$hour <- hour(running$Date)
week <- c("Sunday", "Saturday", "Friday", "Thursday", "Wednesday", "Tuesday", "Monday")
running$wday <- factor(running$wday, levels = week)Exploratory visualisation
Use dplyr to group and summarise and create some plots.
# Create summary
year_month <- running %>%
group_by(year, month) %>%
summarise(total_runs = n(), total_distance = sum(Distance))# Plot data
# TODOHeatmaps
Weekly pattern - most popular times to run during the days of the week
# Create a complete grid to include all possible combinations of weekday and hour even if no data
hour = c(0:23)
grid = expand.grid(week, hour)
# Group by week day, then hour and summarise by counting number of activities per group
heatmap_runs <- running %>%
group_by(wday, hour) %>%
summarise(no_runs = n())
# Join grid and heatmap_runs
# If you don't do this, the heatmap won't show the full week, and have blank tiles for missing values instead of rather showing 0 runs/activity
heatmap_runs_full = heatmap_runs %>%
right_join(grid, by = c("wday"="Var1", "hour"="Var2"))Column `wday`/`Var1` has different attributes on LHS and RHS of join# Replace NA with 0
heatmap_runs_full[is.na(heatmap_runs_full)] <- 0
# Plot heatmap using ggplot
ggplot(heatmap_runs_full, aes(x=hour, y=wday)) +
geom_tile(aes(fill=no_runs), colour="white") +
scale_fill_gradient(name = "No. of runs", low = "white", high = "red") +
scale_x_continuous(breaks=seq(0,23,1)) + ggtitle("Heatmap showing weekly activity pattern")
# Using ColorBrewer to generate palette
display.brewer.all()
pal <- colorRampPalette(brewer.pal(9, "YlOrRd"))(100)
ggplot(heatmap_runs_full, aes(x=hour, y=wday)) +
geom_tile(aes(fill=no_runs), colour="white") +
scale_fill_gradientn(name = "No. of runs", colours=c("white", pal)) +
scale_x_continuous(breaks=seq(0,23,1)) + ggtitle("Heatmap showing weekly activity pattern")
# Save plot
ggsave("../results/heatmap_weekly_counts.png", width = 10, height = 7)Weekly pattern - average distance run at different times during the week
# Group by week day, then hour and summarise by calculating mean distance per group
heatmap_dist <- running %>%
group_by(wday, hour) %>%
summarise(mean_dist = mean(Distance))
# Join grid and heatmap_dist
heatmap_dist_full = heatmap_dist %>%
right_join(grid, by = c("wday"="Var1", "hour"="Var2"))Column `wday`/`Var1` has different attributes on LHS and RHS of join# Replace NA with 0
heatmap_dist_full[is.na(heatmap_dist_full)] <- 0
# Plot heatmap using ggplot
pal <- colorRampPalette(brewer.pal(9, "YlGnBu"))(100)
ggplot(heatmap_dist_full, aes(x=hour, y=wday)) +
geom_tile(aes(fill=mean_dist), colour="white") +
scale_fill_gradientn(name = "Mean distance", colours=c("white", pal)) +
scale_x_continuous(breaks=seq(0,23,1)) + ggtitle("Heatmap showing weekly activity pattern")
Mmmm…that doesn’t look right - what happened at 2pm on a Sunday?! Let’s investigate.
running_outliers = subset(running, Distance >=50)
running_outliersFirstly, look at the average pace of entries 345 and 352. This is rather quick for a run! As this is my own data, I know what these outliers are from! These were actually cycles where I didn’t classify the activities correctly. Secondly, the 65km entry on 13 August was when I used my Garmin watch to measure how far we travelled on a boat in Brazil :) So, let’s remove these three entries.
running <- subset(running, Distance <=60)
# And create the graphs again above.
# Repeated here just for ease of reference.
# Group by week day, then hour and summarise by calculating mean distance per group
heatmap_dist <- running %>%
group_by(wday, hour) %>%
summarise(mean_dist = mean(Distance))
# Join grid and heatmap_dist
heatmap_dist_full = heatmap_dist %>%
right_join(grid, by = c("wday"="Var1", "hour"="Var2"))Column `wday`/`Var1` has different attributes on LHS and RHS of join# Replace NA with 0
heatmap_dist_full[is.na(heatmap_dist_full)] <- 0
# Plot heatmap using ggplot
pal <- colorRampPalette(brewer.pal(9, "YlGnBu"))(100)
ggplot(heatmap_dist_full, aes(x=hour, y=wday)) +
geom_tile(aes(fill=mean_dist), colour="white") +
scale_fill_gradientn(name = "Mean distance", colours=c("white", pal)) +
scale_x_continuous(breaks=seq(0,23,1)) + ggtitle("Heatmap showing weekly activity pattern")
# Save plot
ggsave("../results/heatmap_weekly_dist.png", width = 10, height = 7)Calendar heat map
source ("calendarHeat.R")
calendarHeat(running$Date, running$Distance, varname = "running activities and distances", color = "y2b", ncolors=20)Loading required package: lattice
Loading required package: grid
Loading required package: chron
Attaching package: ‘chron’
The following objects are masked from ‘package:lubridate’:
days, hours, minutes, seconds, years
the condition has length > 1 and only the first element will be used
<<<<<<< HEAD
======= >>>>>>> b1a1cff6c76adb84db272af543df1ec3940595db
LS0tCnRpdGxlOiAiQ2FzZSBTdHVkeTogUnVubmluZyBhY3Rpdml0aWVzIgphdXRob3I6ICJNZWdhbiBCZWNrZXR0IgpkYXRlOiAiMTYgQXVndXN0IDIwMTgiCm91dHB1dDoKICBodG1sX25vdGVib29rOgogICAgZmlnX2hlaWdodDogNwogICAgZmlnX3dpZHRoOiAxMAogICAgdG9jOiB5ZXMKLS0tCiMgSW50cm9kdWN0aW9uCgpUaGlzIGZvbGxvd3Mgb24gZnJvbSBpbnZlc3RpZ2F0aW5nIHRoZSBsdWJyaWRhdGUgcGFja2FnZSB3aXRoaW4gdGhlIHRpZHl2ZXJzZSBhbmQgcHV0dGluZyBzb21lIG9mIHRoZSBmZWF0dXJlcyBpbnRvIHByYWN0aWNlIGJ5IGxvb2tpbmcgYXQgcnVubmluZyBhY3Rpdml0eSBkYXRhIG92ZXIgdGltZS4KCkkgd2FzIGludGVyZXN0ZWQgdG8gaW52ZXN0aWdhdGUgc29tZSB3YXlzIG9mIHBsb3R0aW5nIG15IHJ1bm5pbmcgYWN0aXZpdGllcyBvdmVyIHRpbWUuIFRoZSBkYXRhIGlzIGEgQ1NWIGV4cG9ydCBmcm9tIEdhcm1pbiBDb25uZWN0LgoKIyBMb2FkIHBhY2thZ2VzIG5lZWRlZApgYGB7ciBtZXNzYWdlPUZBTFNFLHdhcm5pbmc9RkFMU0V9CmxpYnJhcnkobHVicmlkYXRlKQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KGdncGxvdDIpCmxpYnJhcnkoc2tpbXIpCmxpYnJhcnkoc3VtbWFyeXRvb2xzKQpsaWJyYXJ5KFJDb2xvckJyZXdlcikKYGBgCgojIFJlYWQgaW4gZGF0YQpgYGB7cn0KYWN0aXZpdGllcyA8LSByZWFkLmNzdigiLi4vZGF0YS9BY3Rpdml0aWVzLmNzdiIsIHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkKYGBgCgojIEluaXRpYWwgdmlldyBvZiBkYXRhCk5vdGVzOgoKLSBJIG5vcm1hbGx5IHVzZSBhIGNvbWJpbmF0aW9uIG9mIHN0cigpLCBzdW1tYXJ5KCksIGdsaW1wc2UoKSB0byBkbyB0aGlzIGFuZCByZWNlbnRseSBleHBsb3JlZCBza2ltKCkuCi0gRm9yICoqc2tpbSgpKiosIEkgbGlrZSBob3cgaXQgc3BsaXRzIHRoZSBzdW1tYXJ5IHN0YXRpc3RpY3MgYWNyb3NzIHZhcmlhYmxlIHR5cGUgYW5kIHNob3dzIGEgYmFzaWMgaGlzdG9ncmFtIGZvciBudW1lcmljYWwgdHlwZXMuCi0gTWlzc2luZyB2YWx1ZXMgYXJlIGFsc28gc2hvd24gYW5kIHRoaXMsIGFsb25nIHdpdGggdGhlIGhpc3RvZ3JhbXMsIGdpdmVzIHNvbWUgaW5kaWNhdGlvbiBvZiB3aGF0IHRvIGV4cGxvcmUgZnVydGhlci4KLSBTdWJzZXF1ZW50IHRvIHRoaXMsIEkgaGF2ZSBub3cgdXNlZCB0aGUgcGFja2FnZSAqKnN1bW1hcnl0b29scyoqIGFuZCBzcGVjaWZpY2FsbHksIGRmU3VtbWFyeSgpLiBUaGUgaGlzdG9ncmFtIHBsb3RzIGFyZSBldmVuIGJldHRlciAoYWx0aG91Z2ggcmVuZGVyIHF1aXRlIHNtYWxsKSBhbmQgdGhlIGRldmVsb3BlcnMgYXJlIHdvcmtpbmcgdG8gaW1wcm92ZSB0aGUgSFRNTCBvdXRwdXQgdG8gd29yayB3aXRoIFIgTWFya2Rvd24uCi0gU2VlIHRoaXMgdmlnbmV0dGU6IGh0dHBzOi8vY3Jhbi5yLXByb2plY3Qub3JnL3dlYi9wYWNrYWdlcy9zdW1tYXJ5dG9vbHMvdmlnbmV0dGVzL1JlY29tbWVuZGF0aW9ucy1ybWFya2Rvd24uaHRtbAoKYGBge3J9CmRpbShhY3Rpdml0aWVzKQojc3RyKGFjdGl2aXRpZXMpCiNzdW1tYXJ5KGFjdGl2aXRpZXMpCmdsaW1wc2UoYWN0aXZpdGllcykKI3NraW0oYWN0aXZpdGllcykKYGBgCgojIFRpZHkgdXAgZGF0YQojIyBFeHRyYWN0IHZhcmlhYmxlcyBvZiBpbnRlcmVzdApNYW55IG9mIHRoZSBjb2x1bW5zIGFyZSBlbXB0eSwgc28gbGV0J3MgY3JlYXRlIGEgc21hbGxlciBzdWJzZXQgb2YgdGhvc2Ugd2UgYXJlIGludGVyZXN0ZWQgaW4uCgpgYGB7cn0KIyBDcmVhdGUgc3Vic2V0IG9mIHZhcmlhYmxlcwpydW5uaW5nIDwtIHN1YnNldChhY3Rpdml0aWVzLCBzZWxlY3QgPSBjKEFjdGl2aXR5LlR5cGUsIERhdGUsIERpc3RhbmNlLCBDYWxvcmllcywgVGltZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBBdmcuUGFjZSwgRWxldi5HYWluKSkKCiMgVmlldyB0aGUgY2xhc3NlcyBvZiB0aGUgc3Vic2V0IG9mIGRhdGEKbGFwcGx5KHJ1bm5pbmcsIGNsYXNzKQpgYGAKCiMjIFRyYW5zZm9ybSBkYXRhIHR5cGVzCnJlYWQuY3N2IGltcG9ydGVkIHNvbWUgb2YgdGhlIHRoZSBudW1iZXJzLCBkYXRlcyBhbmQgdGltZXMgYXMgY2hhcmFjdGVycywgc28gbGV0J3MgZml4IHRoaXMuIEFsc28sIGNvbnZlcnQgQWN0aXZpdHkuVHlwZSB0byBmYWN0b3IgdG8gc2VlIGlmIHRoZXJlIGlzIG1vcmUgdGhhbiBvbmUgbGV2ZWwuCgpgYGB7ciBtZXNzYWdlPUZBTFNFLHdhcm5pbmc9RkFMU0V9CiMgTnVtZXJpYyB2YWx1ZXMKIyBTdWJzdGl0dXRlIDEwMDAgY29tbWEgc2VwYXJhdG9yIHRoZW4gY29udmVydCB0byBudW1lcmljCnJ1bm5pbmckQ2Fsb3JpZXMgPC0gYXMubnVtZXJpYyhnc3ViKCIsIiwgIiIsIHJ1bm5pbmckQ2Fsb3JpZXMpKQpydW5uaW5nJEVsZXYuR2FpbiA8LSBhcy5udW1lcmljKGdzdWIoIiwiLCAiIiwgcnVubmluZyRFbGV2LkdhaW4pKQoKIyBEYXRlLXRpbWUgY29sdW1ucyAoRGF0ZSwgVGltZSwgQXZnLlBhY2UsIEJlc3QuUGFjZSkKcnVubmluZyREYXRlID0geW1kX2htcyhhY3Rpdml0aWVzJERhdGUpCgojIFRPRE8gRmlndXJlIG91dCBob3cgdG8gY29udmVydCBkaWZmZXJlbnQgdGltZSBkdXJhdGlvbnMgdG8gYWN0dWFsIHRpbWVzCiNydW5uaW5nJER1cmF0aW9uIDwtIG1zKHJ1bm5pbmckVGltZSkKI3J1bm5pbmckQXZnLlBhY2UgPC0gbXMocnVubmluZyRBdmcuUGFjZSkKI3J1bm5pbmckQmVzdC5QYWNlIDwtIG1zKHJ1bm5pbmckQmVzdC5QYWNlKQoKIyBDb252ZXJ0IGFjdGl2aXR5IHR5cGUgdG8gZmFjdG9yCnJ1bm5pbmckQWN0aXZpdHkuVHlwZSA8LSBhcy5mYWN0b3IocnVubmluZyRBY3Rpdml0eS5UeXBlKQpsZXZlbHMocnVubmluZyRBY3Rpdml0eS5UeXBlKQpgYGAKClZpZXcgZGF0YSBhZ2FpbiB1c2luZyBzdW1tYXJ5dG9vbHMgcGFja2FnZQpgYGB7cn0KdmlldyhkZlN1bW1hcnkocnVubmluZyksIG1ldGhvZCA9ICdyZW5kZXInKQpgYGAKCiMjIFN1YnNldCBhZ2FpbiB0byBvbmx5IGxvb2sgYXQgbGFzdCAzIGZ1bGwgeWVhcnMgb2YgZGF0YQpgYGB7cn0KcnVubmluZyA8LSBzdWJzZXQocnVubmluZywgcnVubmluZyREYXRlID49ICIyMDE2LTAxLTAxIikKYGBgCgojIENyZWF0ZSBuZXcgdmFyaWFibGVzCmBgYHtyfQpydW5uaW5nJHllYXIgPC0geWVhcihydW5uaW5nJERhdGUpCnJ1bm5pbmckbW9udGggPC0gbW9udGgocnVubmluZyREYXRlLCBsYWJlbD1UUlVFKQpydW5uaW5nJHdlZWsgPC0gd2VlayhydW5uaW5nJERhdGUpCnJ1bm5pbmckd2RheSA8LSB3ZGF5KHJ1bm5pbmckRGF0ZSwgbGFiZWw9VFJVRSwgYWJicj1GQUxTRSkKcnVubmluZyRob3VyIDwtIGhvdXIocnVubmluZyREYXRlKQoKd2VlayA8LSBjKCJTdW5kYXkiLCAiU2F0dXJkYXkiLCAiRnJpZGF5IiwgIlRodXJzZGF5IiwgIldlZG5lc2RheSIsICJUdWVzZGF5IiwgIk1vbmRheSIpCnJ1bm5pbmckd2RheSA8LSBmYWN0b3IocnVubmluZyR3ZGF5LCBsZXZlbHMgPSB3ZWVrKQpgYGAKCiMgRXhwbG9yYXRvcnkgdmlzdWFsaXNhdGlvbgpVc2UgZHBseXIgdG8gZ3JvdXAgYW5kIHN1bW1hcmlzZSBhbmQgY3JlYXRlIHNvbWUgcGxvdHMuCgpgYGB7cn0KIyBDcmVhdGUgc3VtbWFyeQp5ZWFyX21vbnRoIDwtIHJ1bm5pbmcgJT4lCiAgZ3JvdXBfYnkoeWVhciwgbW9udGgpICU+JQogIHN1bW1hcmlzZSh0b3RhbF9ydW5zID0gbigpLCB0b3RhbF9kaXN0YW5jZSA9IHN1bShEaXN0YW5jZSkpCmBgYAoKYGBge3J9CiMgUGxvdCBkYXRhCiMgVE9ETwpgYGAKCiMgSGVhdG1hcHMKIyMgV2Vla2x5IHBhdHRlcm4gLSBtb3N0IHBvcHVsYXIgdGltZXMgdG8gcnVuIGR1cmluZyB0aGUgZGF5cyBvZiB0aGUgd2VlawpgYGB7cn0KIyBDcmVhdGUgYSBjb21wbGV0ZSBncmlkIHRvIGluY2x1ZGUgYWxsIHBvc3NpYmxlIGNvbWJpbmF0aW9ucyBvZiB3ZWVrZGF5IGFuZCBob3VyIGV2ZW4gaWYgbm8gZGF0YQpob3VyID0gYygwOjIzKQpncmlkID0gZXhwYW5kLmdyaWQod2VlaywgaG91cikKCiMgR3JvdXAgYnkgd2VlayBkYXksIHRoZW4gaG91ciBhbmQgc3VtbWFyaXNlIGJ5IGNvdW50aW5nIG51bWJlciBvZiBhY3Rpdml0aWVzIHBlciBncm91cApoZWF0bWFwX3J1bnMgPC0gcnVubmluZyAlPiUKICBncm91cF9ieSh3ZGF5LCBob3VyKSAlPiUKICBzdW1tYXJpc2Uobm9fcnVucyA9IG4oKSkKCiMgSm9pbiBncmlkIGFuZCBoZWF0bWFwX3J1bnMKIyBJZiB5b3UgZG9uJ3QgZG8gdGhpcywgdGhlIGhlYXRtYXAgd29uJ3Qgc2hvdyB0aGUgZnVsbCB3ZWVrLCBhbmQgaGF2ZSBibGFuayB0aWxlcyBmb3IgbWlzc2luZyB2YWx1ZXMgaW5zdGVhZCBvZiByYXRoZXIgc2hvd2luZyAwIHJ1bnMvYWN0aXZpdHkKaGVhdG1hcF9ydW5zX2Z1bGwgPSBoZWF0bWFwX3J1bnMgJT4lCiAgcmlnaHRfam9pbihncmlkLCBieSA9IGMoIndkYXkiPSJWYXIxIiwgImhvdXIiPSJWYXIyIikpCgojIFJlcGxhY2UgTkEgd2l0aCAwCmhlYXRtYXBfcnVuc19mdWxsW2lzLm5hKGhlYXRtYXBfcnVuc19mdWxsKV0gPC0gMAoKIyBQbG90IGhlYXRtYXAgdXNpbmcgZ2dwbG90CmdncGxvdChoZWF0bWFwX3J1bnNfZnVsbCwgYWVzKHg9aG91ciwgeT13ZGF5KSkgKwogIGdlb21fdGlsZShhZXMoZmlsbD1ub19ydW5zKSwgY29sb3VyPSJ3aGl0ZSIpICsKICBzY2FsZV9maWxsX2dyYWRpZW50KG5hbWUgPSAiTm8uIG9mIHJ1bnMiLCBsb3cgPSAid2hpdGUiLCBoaWdoID0gICJyZWQiKSArCiAgc2NhbGVfeF9jb250aW51b3VzKGJyZWFrcz1zZXEoMCwyMywxKSkgKyBnZ3RpdGxlKCJIZWF0bWFwIHNob3dpbmcgd2Vla2x5IGFjdGl2aXR5IHBhdHRlcm4iKQpgYGAKCgpgYGB7cn0KIyBVc2luZyBDb2xvckJyZXdlciB0byBnZW5lcmF0ZSBwYWxldHRlCgpkaXNwbGF5LmJyZXdlci5hbGwoKQoKcGFsIDwtIGNvbG9yUmFtcFBhbGV0dGUoYnJld2VyLnBhbCg5LCAiWWxPclJkIikpKDEwMCkKCmdncGxvdChoZWF0bWFwX3J1bnNfZnVsbCwgYWVzKHg9aG91ciwgeT13ZGF5KSkgKwogIGdlb21fdGlsZShhZXMoZmlsbD1ub19ydW5zKSwgY29sb3VyPSJ3aGl0ZSIpICsKICBzY2FsZV9maWxsX2dyYWRpZW50bihuYW1lID0gIk5vLiBvZiBydW5zIiwgY29sb3Vycz1jKCJ3aGl0ZSIsIHBhbCkpICsKICBzY2FsZV94X2NvbnRpbnVvdXMoYnJlYWtzPXNlcSgwLDIzLDEpKSArIGdndGl0bGUoIkhlYXRtYXAgc2hvd2luZyB3ZWVrbHkgYWN0aXZpdHkgcGF0dGVybiIpCgpgYGAKCmBgYHtyfQojIFNhdmUgcGxvdApnZ3NhdmUoIi4uL3Jlc3VsdHMvaGVhdG1hcF93ZWVrbHlfY291bnRzLnBuZyIsIHdpZHRoID0gMTAsIGhlaWdodCA9IDcpCmBgYAojIyBXZWVrbHkgcGF0dGVybiAtIGF2ZXJhZ2UgZGlzdGFuY2UgcnVuIGF0IGRpZmZlcmVudCB0aW1lcyBkdXJpbmcgdGhlIHdlZWsKCmBgYHtyfQojIEdyb3VwIGJ5IHdlZWsgZGF5LCB0aGVuIGhvdXIgYW5kIHN1bW1hcmlzZSBieSBjYWxjdWxhdGluZyBtZWFuIGRpc3RhbmNlIHBlciBncm91cApoZWF0bWFwX2Rpc3QgPC0gcnVubmluZyAlPiUKICBncm91cF9ieSh3ZGF5LCBob3VyKSAlPiUKICBzdW1tYXJpc2UobWVhbl9kaXN0ID0gbWVhbihEaXN0YW5jZSkpCgojIEpvaW4gZ3JpZCBhbmQgaGVhdG1hcF9kaXN0CmhlYXRtYXBfZGlzdF9mdWxsID0gaGVhdG1hcF9kaXN0ICU+JQogIHJpZ2h0X2pvaW4oZ3JpZCwgYnkgPSBjKCJ3ZGF5Ij0iVmFyMSIsICJob3VyIj0iVmFyMiIpKQoKIyBSZXBsYWNlIE5BIHdpdGggMApoZWF0bWFwX2Rpc3RfZnVsbFtpcy5uYShoZWF0bWFwX2Rpc3RfZnVsbCldIDwtIDAKCiMgUGxvdCBoZWF0bWFwIHVzaW5nIGdncGxvdApwYWwgPC0gY29sb3JSYW1wUGFsZXR0ZShicmV3ZXIucGFsKDksICJZbEduQnUiKSkoMTAwKQpnZ3Bsb3QoaGVhdG1hcF9kaXN0X2Z1bGwsIGFlcyh4PWhvdXIsIHk9d2RheSkpICsKICBnZW9tX3RpbGUoYWVzKGZpbGw9bWVhbl9kaXN0KSwgY29sb3VyPSJ3aGl0ZSIpICsKICBzY2FsZV9maWxsX2dyYWRpZW50bihuYW1lID0gIk1lYW4gZGlzdGFuY2UiLCBjb2xvdXJzPWMoIndoaXRlIiwgcGFsKSkgKwogIHNjYWxlX3hfY29udGludW91cyhicmVha3M9c2VxKDAsMjMsMSkpICsgZ2d0aXRsZSgiSGVhdG1hcCBzaG93aW5nIHdlZWtseSBhY3Rpdml0eSBwYXR0ZXJuIikKYGBgCgpNbW1tLi4udGhhdCBkb2Vzbid0IGxvb2sgcmlnaHQgLSB3aGF0IGhhcHBlbmVkIGF0IDJwbSBvbiBhIFN1bmRheT8hIExldCdzIGludmVzdGlnYXRlLgoKYGBge3J9CnJ1bm5pbmdfb3V0bGllcnMgPSBzdWJzZXQocnVubmluZywgRGlzdGFuY2UgPj01MCkKcnVubmluZ19vdXRsaWVycwoKYGBgCgpGaXJzdGx5LCBsb29rIGF0IHRoZSBhdmVyYWdlIHBhY2Ugb2YgZW50cmllcyAzNDUgYW5kIDM1Mi4gVGhpcyBpcyByYXRoZXIgcXVpY2sgZm9yIGEgcnVuISBBcyB0aGlzIGlzIG15IG93biBkYXRhLCBJIGtub3cgd2hhdCB0aGVzZSBvdXRsaWVycyBhcmUgZnJvbSEgVGhlc2Ugd2VyZSBhY3R1YWxseSBjeWNsZXMgd2hlcmUgSSBkaWRuJ3QgY2xhc3NpZnkgdGhlIGFjdGl2aXRpZXMgY29ycmVjdGx5LiBTZWNvbmRseSwgdGhlIDY1a20gZW50cnkgb24gMTMgQXVndXN0IHdhcyB3aGVuIEkgdXNlZCBteSBHYXJtaW4gd2F0Y2ggdG8gbWVhc3VyZSBob3cgZmFyIHdlIHRyYXZlbGxlZCBvbiBhIGJvYXQgaW4gQnJhemlsIDopIFNvLCBsZXQncyByZW1vdmUgdGhlc2UgdGhyZWUgZW50cmllcy4KCmBgYHtyfQpydW5uaW5nIDwtIHN1YnNldChydW5uaW5nLCBEaXN0YW5jZSA8PTYwKQoKIyBBbmQgY3JlYXRlIHRoZSBncmFwaHMgYWdhaW4gYWJvdmUuCiMgUmVwZWF0ZWQgaGVyZSBqdXN0IGZvciBlYXNlIG9mIHJlZmVyZW5jZS4KCiMgR3JvdXAgYnkgd2VlayBkYXksIHRoZW4gaG91ciBhbmQgc3VtbWFyaXNlIGJ5IGNhbGN1bGF0aW5nIG1lYW4gZGlzdGFuY2UgcGVyIGdyb3VwCmhlYXRtYXBfZGlzdCA8LSBydW5uaW5nICU+JQogIGdyb3VwX2J5KHdkYXksIGhvdXIpICU+JQogIHN1bW1hcmlzZShtZWFuX2Rpc3QgPSBtZWFuKERpc3RhbmNlKSkKCiMgSm9pbiBncmlkIGFuZCBoZWF0bWFwX2Rpc3QKaGVhdG1hcF9kaXN0X2Z1bGwgPSBoZWF0bWFwX2Rpc3QgJT4lCiAgcmlnaHRfam9pbihncmlkLCBieSA9IGMoIndkYXkiPSJWYXIxIiwgImhvdXIiPSJWYXIyIikpCgojIFJlcGxhY2UgTkEgd2l0aCAwCmhlYXRtYXBfZGlzdF9mdWxsW2lzLm5hKGhlYXRtYXBfZGlzdF9mdWxsKV0gPC0gMAoKIyBQbG90IGhlYXRtYXAgdXNpbmcgZ2dwbG90CnBhbCA8LSBjb2xvclJhbXBQYWxldHRlKGJyZXdlci5wYWwoOSwgIllsR25CdSIpKSgxMDApCmdncGxvdChoZWF0bWFwX2Rpc3RfZnVsbCwgYWVzKHg9aG91ciwgeT13ZGF5KSkgKwogIGdlb21fdGlsZShhZXMoZmlsbD1tZWFuX2Rpc3QpLCBjb2xvdXI9IndoaXRlIikgKwogIHNjYWxlX2ZpbGxfZ3JhZGllbnRuKG5hbWUgPSAiTWVhbiBkaXN0YW5jZSIsIGNvbG91cnM9Yygid2hpdGUiLCBwYWwpKSArCiAgc2NhbGVfeF9jb250aW51b3VzKGJyZWFrcz1zZXEoMCwyMywxKSkgKyBnZ3RpdGxlKCJIZWF0bWFwIHNob3dpbmcgd2Vla2x5IGFjdGl2aXR5IHBhdHRlcm4iKQpgYGAKCmBgYHtyfQojIFNhdmUgcGxvdApnZ3NhdmUoIi4uL3Jlc3VsdHMvaGVhdG1hcF93ZWVrbHlfZGlzdC5wbmciLCB3aWR0aCA9IDEwLCBoZWlnaHQgPSA3KQpgYGAKCiMjIENhbGVuZGFyIGhlYXQgbWFwCgpgYGB7cn0Kc291cmNlICgiY2FsZW5kYXJIZWF0LlIiKQoKY2FsZW5kYXJIZWF0KHJ1bm5pbmckRGF0ZSwgcnVubmluZyREaXN0YW5jZSwgdmFybmFtZSA9ICJydW5uaW5nIGFjdGl2aXRpZXMgYW5kIGRpc3RhbmNlcyIsIGNvbG9yID0gInkyYiIsIG5jb2xvcnM9MjApCgpgYGAKPDw8PDw8PCBIRUFECgoKPT09PT09PQo+Pj4+Pj4+IGIxYTFjZmY2Yzc2YWRiODRkYjI3MmFmNTQzZGYxZWMzOTQwNTk1ZGIK