Radiation levels at Fukushima

[This article was first published on Back Side Smack » R Stuff, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

From BWR

The above graph is derived from data scraped from TEPCO press releases. Every hour or so for the first few days of the crisis, a TEPCO van would record radiation (probably Beta/Gamma, but the translation is unclear) at the front gate to the Fukushima complex. The series is less continuous than the graph would lead you to believe, but I don’t feel the actual series is so different as to warrant concern. I pulled the data down from here (some explanation of the original docs on a reddit thread) and converted the dates from a ragged MM/DD/YYYY HR:MN AM format to something a little more tractable using Google Refine. I used Google Refine because making such a change programmatically was much more challenging, with the irregular format trouncing my understanding of strptime(). As you will see from the code below I actually only used the dates for labeling. But the y axis is completely true to the original data.

rad<-read.csv(file="~/Downloads/radtext.csv",as.is=TRUE)
radplot<- ggplot(data=rad,aes(x=1:619,y=cumsum(microSv))) + geom_line(size=I(1),linetype=I(2)) + scale_y_log10(name="microSieverts") + scale_x_continuous(labels=c(rad$time[seq(from=100,to=600,by=100)]),breaks=c(seq(from=100,to=600,by=100)),name='')
radplot<- radplot + geom_line(aes(x=1:619,y=microSv),colour="blue") + geom_hline(aes(yintercept=50000),colour="red")
radplot<- radplot + annotate(geom="text",x=575,y=35000,label="Yearly Limit",colour="red") + annotate(geom="text",x=575,y=1000,label="Hourly Dose",colour="blue") + annotate(geom="text",x=575,y=200000,label="Cumulative")
view raw fukushima.R hosted with ❤ by GitHub

To leave a comment for the author, please follow the link and comment on their blog: Back Side Smack » R Stuff.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)