bxb_vo2_analysis

library(ggplot2);
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(pastecs);

## Loading required package: boot

## 
## Attaching package: 'pastecs'

## The following objects are masked from 'package:dplyr':
## 
##     first, last

library(psych);

## 
## Attaching package: 'psych'

## The following object is masked from 'package:boot':
## 
##     logit

## The following objects are masked from 'package:ggplot2':
## 
##     %+%, alpha

library(nlme);

## 
## Attaching package: 'nlme'

## The following object is masked from 'package:dplyr':
## 
##     collapse

setwd ("~/maxVO2_RSA/VMax_bxb_csv")     #Set the working directory to the file folder for this paper
id130 = read.csv("130.csv",header = TRUE)   #Import the data for analysis
id130 <- filter(id130,vo2_L < 7.5)

scatter2 <- ggplot(id130,aes(tod, vo2_L)) + labs(title = "VO2 during a max aerobic test", x = "Time",y="VO2 (L/min" )
scatter2 + geom_point()

scatter <- ggplot(id130,aes(vo2_L, vco2_L)) + labs(title = "VO2 and VCO2 during a max aerobic test", x = "VO2 (L/min)",y="VCO2 (L/min" )
scatter + geom_point()

scatter3 <- ggplot(id130,aes(hr, vo2_L)) + labs(title = "HR and VO2 during a max aerobic test", x = "HR",y="VO2 (L/min" )
scatter3 + geom_point()

## Warning: Removed 3 rows containing missing values (geom_point).