monsterLakeArea.R

la <- c(Okanagan=348, Seneca=173, Cayuga=172, Ballston=1.10884, Reindeer=5650, Somyn=NA, Thetis=8.336524, Manitoba=4624, Cold=373, Winnipegosis=5370, Crescent=18.401866, Simcoe=722, Muskrat=12.02, Superior=127700, Ontario=7340, Timiskaming=295, Memphremagog=110, Tahoe=490, Tarpon=10, Fulks=NA, Flathead=510, Hopatcong=10, Champlain=1331, Willoughby=1687, Erie=25667, Bear=280, Raystown=33.6, GreatSalt=4400, Arenal=85, NahuelHuapi=530, Dilolo=12.35, MamfePool=NA, DiloloSwamps=NA, Edward=2325, MaikaMarshes=NA, Victoria=68800, Makele=NA, Tangankyika=32900, Bangweulu=15100, KafueFlats=NA, Mweru=5120, IshibaNgandu=NA, Kariba=5580, Bassenthwaite=5.1, Ness=56, Morar=26.7, Maree=28.6, Lagarfljot=53, LoughLeane=19, Muckross=2.7, UpperLake=1.7, Como=146, Seljordsvatnet=16.7, Storsjon=464, Gryttjen=NA, Brosno=7.55, Labynkyr=61, LakeVan=3755, Kanas=45.73, TianChiHeaven=9.82, Ikeda=11, Kussharo=79.3, Pinatubo=1.83, Modewarre=NA, Utopia=14) # area in km^2; Dilolo I measured on google maps and I excluded the NE corner that looked to be marsh; measured Brosno on google maps, following the not-earth view depiction of water; 
ra <- c(WhiteRiver=1162, Payette=134, HudsonRiver=507) # length in km; Payette is listed as lake, but could only find river length (is probably reservoir)

pdf("~/Desktop/MonsterLakeArea.pdf")

#Fig 1
hist(log(la))

#Fig 2
hist(log(la), breaks=15)

#Fig 3
ML_area <- seq(from=log(1), to=log(1.5E5), length.out=100)
ML_area_ecdf <- ecdf(log(la)) 
plot(ML_area, ML_area_ecdf(ML_area), xlab="Monster Lake Area (log(km^2))", ylab="Empirical CDF")

dev.off()