# this is a comment
getwd() ## shows working directory
dir() ## lists directory contents

dat <- read.delim("./data/milestone.csv", header = TRUE, sep = ',')

head(dat)
tail(dat)

plot(percent ~ year, dat, subset = milestone == "independent")
barplot(percent ~ year, dat, subset = milestone == "independent")

plot(percent ~ year, dat, subset = milestone == "independent")
mtext("Percent living independently decreases over time", side = 3, font = 2)

plot(percent ~ year, dat, subset = milestone == "independent", xlab = "Year", ylab = "Percent Independent", type = 'o', pch = 19, las = 1, xlim = c(1980, 2030), ylim = c(50, 100))
mtext("Percent living independently decreases over time", side = 3, font = 2)

plot(1:25, pch = 1:25, col = 1:25, cex = 1.5)