Golden Cheetah

# R script will run on selection.
##
## GC.activity(compare=FALSE)
## GC.metrics(all=FALSE,
## compare=FALSE)
##
## Get the current ride or metrics
##

## get data
compares <- GC.metrics(compare=TRUE)

## ylim/xlim set
maxy <- 0
miny <- 0
maxx <- 0
minx <- 0


# loop through each looking for max min
for (compare in compares) {

# get the metrics
metrics <- compare$metrics

# Duration on Y axis
y <- max(metrics$NP)
if (y >maxy) maxy = y
y <- min(metrics$NP)
if (y < miny) miny = y

# Distance on X axis
x <- max(metrics$hr)
if (x >maxx) maxx = x
x <- min(metrics$hr)
if (x < minx) minx = x

}

# set first plot
plot(x=metrics$hr, y=metrics$NP,
col="black", pch=20, main="", xlab="",
ylim=c(miny,maxy), xlim=c(minx,maxx))

## now plot for real
for (compare in compares) {
metrics <- compare$metrics
points(x=metrics$hr, y=metrics$NP,
col=compare$color, pch=20)
}

## title
title(main="HR v Power", xlab="HR (Bpm)",
ylab="NP (Watt)")

compares <- GC.metrics(compare=TRUE)

for (compare in compares) {
metrics <- compare$metrics
x <- metrics$date
y <- metrics$HrNp_Ratio

# plot points
plot(x, y, col=compare$color, pch=20)
# linear model
model <- lm(y ~ x)
# plot line model
abline(model, col='green')
}

## title
title(main="HrNp_Ratio", xlab="Date", 
ylab="Ratio")

compares <- GC.metrics(compare=TRUE)

## ylim/xlim set
maxy <- 0
miny <- 0
maxx <- Sys.Date()
minx <- Sys.Date()

# loop through each looking for max min
for (compare in compares) {
# get the metrics
metrics <- compare$metrics

# HrNp_Ratio on Y axis
tmp <- max(metrics$HrNp_Ratio)
if (tmp > maxy) maxy = tmp
tmp <- min(metrics$HrNp_Ratio)
if (tmp < miny) miny = tmp

# Date on X axis
tmp <- max(metrics$date)
if (tmp > maxx) maxx = tmp
tmp <- min(metrics$date)
if (tmp < minx) minx = tmp
}

# set first plot
plot(x=metrics$date, y=metrics$HrNp_Ratio, col="black", pch=20, main="", xlab="", ylim=c(miny,maxy), xlim=c(minx,maxx))

for (compare in compares) {
metrics <- compare$metrics
data_x <- metrics$date
data_y <- metrics$HrNp_Ratio

# plot points
points(x = data_x, y = data_y, col=compare$color)
# linear model
model <- lm(data_y ~ data_x)
# plot line model
abline(model, col=compare$color)
}

## title
title(main="HrNp_Ratio", xlab="Date", 
ylab="Ratio")