AmeliaMN

# Data preparation

# Not included in this code snippet to respect the current proprietary nature of this research.
# Once research is to be reproduced, data preparation procedures will be made available.

#Plotting

par(mar=c(5,4,4,5)+.1)
plot(newyears, pathtangent, type="l", col="white",xlim=c(1890, 2030), ylim=c(30, 225), xaxt="n", yaxt="n", ylab="", xlab="")
axis(1, at=seq(from=1890, to=2030, by=10))

AmeliaMN

# Data preparation

# Not included in this code snippet to respect the current proprietary nature of this research.
# Once research is to be reproduced, data preparation procedures will be made available.

#Plotting

plot(ypole, xpole, xlim=c(0.35, -0.12), ylim=c(0.15, -0.15), type="l", col="white", ylab="", xlab="", yaxt="n", xaxt="n")
axis(1, at=seq(from=0.35, to=-0.12, by=-0.05), labels=c(0.35, 0.3, 0.25, 0.2, 0.15, 0.1, 0.05, 0, -0.05, -0.1))
axis(2, at=seq(from=0.15, to=-0.15, by=-0.05), labels=c(0.15, 0.1, 0.05, 0, -0.05, -0.1, -0.15))