R code for Shifting Bell Curves

A commentator named Jonathan asked for the code by which I produced my bell curve graph way back in 2013. Here it is. Note: I'm using the Cowtan and Way 2.0 temperature reconstruction in this example rather than NASA GISS as in the original post as the Cowtan and Way data is more accessible for use in R.

>Cowtan<-read.table("http://www-users.york.ac.uk/~kdc3/papers/coverage2013/had4_krig_v2_0_0.txt", header=F)
> names(Cowtan)<-c("Year", "Temperature", "Uncertainty1", "Uncertainty2", "Uncertainty3") #Name the columns
> summary(Cowtan) #Check to see if the column names look right and the data imported correctly
> S1850s<-subset(Cowtan, Year<1860) #Get subsets of each decade
> S1950s<-subset(Cowtan, Year>=1950 & Year<1960)
> S2007<-subset(Cowtan, Year>=2007 & Year<2017)
> D1850s=density(S1850s$Temperature) #Get the density kernals
> D1950s<-density(S1950s$Temperature)
> D2007<-density(S2007$Temperature)
> plot(D1850s, main="Bell curve of global temperatures during selected decades", xlab="Temperature anomalies (ºC)",xlim=c(-1.1,1.3), ylim=c(0,3), lwd=1.5)
> points(D2007, type="l",col="red", lwd=1.5)
> points(D1950s, type="l",col="blue", lwd=1.5)
> legend("topleft", legend=c("1850s", "1950s", "2007-2016"), col=c("black", "blue", "red"), lwd=2)


Comments

Popular posts from this blog

James Taylor gets polar ice wrong—as usual

Tom Luongo's multiple lies about climate change