###########Figure 1.3 from Lavine's book#################

par(mfrow=c(2,1))

mydata=read.table("OceanTemps.txt",header=T)
 
lat=mydata$latitude
lon=mydata$longitude
temp=mydata$temperature

hist(temp,prob=T)
lines(density(temp))


############################################################
#Load the data:
data(discoveries)

#Plot histogram
hist (discoveries, xlab="discoveries", ylab="", main="",
prob=T, breaks=seq(-.5,12.5,by=1))

#Plot Poisson density function
lines(0:12, dpois(0:12, 3.1), type="b"