blog.R

library(ggplot2)

F=read.csv("PLOSONE2013Data.csv")
CY=(F$CYNewWorldOrder+F$CYMLK+F$CYMoon+F$CYJFK+F$CY911+F$CYDiana)/6
CLIM=(F$CNatFluct+F$CdueGHG+F$CseriousDamage+F$CO2causesCC+F$HumansInsign)/5
LGO13=data.frame(cbind(CLIM,CY))

F=read.csv("LskyetalPsychSciClimate.csv")
CY=(F$CYNewWorldOrder+F$CYSARS+F$CYPearlHarbor+F$CYMLK+F$CYMoon+F$CYJFK+F$CY911+F$CYDiana+F$CYOkla+F$CYCoke+F$CYArea51+F$CYRoswell)/12
CLIM=(F$CO2TempUp+F$CO2AtmosUp+F$CO2WillNegChange+F$CO2HasNegChange+F$CauseCO2)/5
LOG13=data.frame(cbind(CLIM,CY))

CLIM=CLIM[CY<=2.5]
CY=CY[CY<=2.5]
LOG13tr=data.frame(cbind(CLIM,CY))

ggplot(data=LGO13,aes(x=CY,y=CLIM)) + geom_point(size=0.8,position = position_jitter(w = 0.02, h = 0.02)) + geom_smooth(method="loess")
ggsave(file="blog1a.png",width=3,height=3,dpi=400)
ggplot(data=LGO13,aes(x=CLIM,y=CY)) + geom_point(size=0.8,position = position_jitter(w = 0.02, h = 0.02)) + geom_smooth(method="loess")
ggsave(file="blog1b.png",width=3,height=3,dpi=400)

ggplot(data=LOG13) + geom_histogram(aes(x=CLIM),binwidth=1/5,origin=1-1/10) +ggtitle("LOG13-blogs")
ggsave(file="blog2a.png",width=3,height=3,dpi=400)
ggplot(data=LGO13) + geom_histogram(aes(x=CLIM),binwidth=1/5,origin=1-1/10) +ggtitle("LGO13-panel")
ggsave(file="blog2b.png",width=3,height=3,dpi=400)

ggplot(data=LOG13,aes(x=CLIM,y=CY)) + geom_point(size=0.8,position = position_jitter(w = 0.02, h = 0.02)) + geom_smooth(method="loess")
ggsave(file="blog3.png",width=3,height=3,dpi=400)

ggplot(data=LOG13tr,aes(x=CLIM,y=CY)) + geom_point(size=0.8,position = position_jitter(w = 0.02, h = 0.02)) + geom_smooth(method="loess")
ggsave(file="blog4.png",width=3,height=3,dpi=400)

Back to list of R-files