A script to process and chart step count data exported from Pedometer++ for iOS.

step-count-analysis.r

#! /usr/bin/env Rscript

library(zoo)

raw_data = read.csv('Export.csv')
raw_data$Date <- as.Date(raw_data$Date, "%m/%d/%y")
readings = data.frame(
    dates=raw_data$Date,
    steps=raw_data$Steps
)

readings = subset(readings, dates > "2012-01-01" )
readings.steps = read.zoo(readings, drop = FALSE)
readings$rolling_average <- rollmean(
    readings.steps$steps,
    30,
    partial=TRUE,
    fill=NA,
    align="left"
)

mean_line <- rep(mean(readings$steps), length(readings$steps))
is_streak <- function(x) {return(if (x > 6500) 10 else NA)}


# Output Setup

png(
  "step-counts.png",
  width     = 1600,
  height    = 900,
)
par(bg = 'white')
par(mar = c(5, 4, 4, 4) + 0.3)

# Plots

plot(
    mean_line,
    axes=FALSE,
    bty='n',
    xlab="",
    ylab="",
    type='l',
    col='orange3',
    lwd=3.0,
    ylim=c(0, max(readings$steps) + 20)
)
par(new = TRUE)
plot(
    readings$dates,
    readings$steps,
    xlab="",
    ylab="",
    col='blue3',
    type='l',
    lwd=2.0,
    ylim=c(0, max(readings$steps))
)
mtext("Steps (blue)", line=3, side=2)
mtext("Date", line=3, side=1)
par(new = TRUE)
plot(
    readings$dates,
    readings$rolling_average,
    axes=FALSE,
    bty='n',
    xlab="",
    ylab="",
    col='green4',
    lwd=2.0,
    cex=0.6,
    pch=21,
    bg="green",
    ylim=c(0, max(readings$steps) + 20)
)
par(new = TRUE)
plot(
    readings$dates,
    sapply(readings$steps, FUN=is_streak),
    axes=FALSE,
    bty='n',
    xlab="",
    ylab="",
    type='l',
    col='red3',
    lwd=3.0,
    ylim=c(0, max(readings$steps) + 20)
)

dev.off()