Eccentricity, Obliquity and Precession

The Earth climate is primarily influenced by the Sun,
thus changes in the orbital variations are generally
considered in paleoceanographic studies. Thus,
my first step was to download this data.
In this code, you download the data directly from the net,
and plot it in R. No extra packages required.

#  Orbital Variations and Insolation Database: Insolation.readme file
 orbit=read.table("ftp://ftp.ncdc.noaa.gov/pub/data/paleo/insolation/orbit91", 
      skip=3, col.names=c("age", "ECC","OMEGA","OBL","PREC","I65NJul","I65SJan","I15NJul","I15SJan"))
 #summary(orbit)
 orbit$age=abs(orbit$age)
 par(mfrow=c(2,2), mar=c(4,4,1,1), cex=.7)
 plot(ECC~age, orbit[orbit$age<450,], typ="l", col=4, main="Eccentricity")
 plot(OMEGA~age, orbit[orbit$age<450,], typ="l", col=2)
 plot(OBL~age, orbit[orbit$age<450,], typ="l", col=3, main="Obliquity")
 plot(PREC~age, orbit[orbit$age<450,], typ="l", col=5, main="Precession")

Created by Pretty R at inside-R.org

 And than you get something like this: