> mdata <-read.table("beer.csv", header=T, sep=",")
> tmp <-mdata[,1]
> bh  <-mdata[,4]
> cbind(tmp,bh)
       tmp      bh
 [1,]  7.4 360.539
 [2,]  7.9 434.770
 [3,] 12.2 549.474
 [4,] 16.1 601.685
 [5,] 18.4 592.449
 [6,] 21.6 747.202
 [7,] 28.0 712.046
 [8,] 28.0 657.574
 [9,] 23.1 525.572
[10,] 19.0 541.173
[11,] 11.6 517.968
[12,]  7.2 689.865

> summary(cbind(tmp,bh))
      tmp              bh       
 Min.   : 7.20   Min.   :360.5  
 1st Qu.:10.68   1st Qu.:523.7  
 Median :17.25   Median :571.0  
 Mean   :16.71   Mean   :577.5  
 3rd Qu.:21.98   3rd Qu.:665.6  
 Max.   :28.00   Max.   :747.2  

> sd(tmp)
[1] 7.563844

> sd(bh)
[1] 113.8329

> cor(cbind(tmp,bh))
          tmp        bh
tmp 1.0000000 0.5794202
bh  0.5794202 1.0000000

> var(cbind(tmp,bh))
          tmp         bh
tmp  57.21174   498.8889
bh  498.88892 12957.9204

> plot(tmp,bh)