Solutions to Selected Exercises for Braun and Murdoch's a First

Home , Timor

Solutions to Selected Exercises for Braun and Murdoch’s A First Course in Statistical Programming with R

Kristy Alexander, Yiwen Diao, Qiang Fu, and Yu Han W. John Braun and Duncan J. Murdoch

November 1, 2007 Chapter 3

Programming Statistical Graphics

3.1 High level plots

1. (a) > hist(islands)

Histogram of islands 40 30 20 Frequency 10 0

0 5000 10000 15000

islands

6 Sturges Scott 40 40 30 30 20 20 Frequency Frequency 10 10 0 0

0 5000 10000 15000 0 5000 10000 15000 20000

islands islands

> par(mfrow=c(1,2)) > hist(log(islands),breaks="Sturges",main="Sturges") > hist(log(islands),breaks="Scott",main="Scott")

Sturges Scott 30 15 25 20 10 15 Frequency Frequency 10 5 5 0 0

2 4 6 8 10 2 4 6 8 10

log(islands) log(islands)

(e) > par(mfrow=c(1,2)) > dotchart(islands,main="islands") > dotchart(log(islands),main="log(islands)")

7 islands log(islands)

Victoria ● Victoria ● Vancouver ● Vancouver ● Timor ● Timor ● Tierra del Fuego ● Tierra del Fuego ● Tasmania ● Tasmania ● Taiwan ● Taiwan ● Sumatra ● Sumatra ● Spitsbergen ● Spitsbergen ● Southampton ● Southampton ● South America ● South America ● Sakhalin ● Sakhalin ● Prince of Wales ● Prince of Wales ● Novaya Zemlya ● Novaya Zemlya ● North America ● North America ● Newfoundland ● Newfoundland ● New Zealand (S) ● New Zealand (S) ● New Zealand (N) ● New Zealand (N) ● New Guinea ● New Guinea ● New Britain ● New Britain ● Moluccas ● Moluccas ● Mindanao ● Mindanao ● Melville ● Melville ● Madagascar ● Madagascar ● Luzon ● Luzon ● Kyushu ● Kyushu ● Java ● Java ● Ireland ● Ireland ● Iceland ● Iceland ● Honshu ● Honshu ● Hokkaido ● Hokkaido ● Hispaniola ● Hispaniola ● Hainan ● Hainan ● Greenland ● Greenland ● Europe ● Europe ● Ellesmere ● Ellesmere ● Devon ● Devon ● Cuba ● Cuba ● Celon ● Celon ● Celebes ● Celebes ● Britain ● Britain ● Borneo ● Borneo ● Banks ● Banks ● Baffin ● Baffin ● Axel Heiberg ● Axel Heiberg ● Australia ● Australia ● Asia ● Asia ● Antarctica ● Antarctica ● Africa ● Africa ●

0 5000 10000 15000 4 6 8 10

3. (a) > plot(pressure$temperature,pressure$pressure,xlab="temperature",ylab="pressure")

● 800

600 ●

400 ● pressure

●

200 ●

● ● ● ● ● ● ● ● ● ● ● ● ● ● 0

0 50 100 150 200 250 300 350

temperature

8 ●

●

2.5 ●

●

● 2.0 ●

●

1.5 ● pressure ●

●

● 1.0 ● ● ● ● 0.5 ● ●

0 50 100 150 200 250 300 350

temperature

Chapter Exercises

1. (a) > hist(log(islands,10),breaks="Scott",axes=FALSE,xlab="area",main="Histogram of Island Areas") > axis(1,at=1:5,labels=10^(1:5)) > axis(2) > box()

Histogram of Island Areas 35 30 25 20 15 Frequency 10 5 0

10 100 1000 10000 1e+05

area

(c) > hist(round(log(islands,10)),breaks="Sturges",axes=FALSE,xlab="area", + main="Histogram of Island Areas") > axis(1,at=1:5,labels=10^(1:5)) > axis(2) > box()

9 Histogram of Island Areas 20 15 10 Frequency 5 0

10 100 1000 10000

area