xy.data<-read.csv("03_xy.csv")

summary(xy.data)

pdf("histograms.pdf", width=8, height=5)
par(mfrow=c(1,2), mar=c(5,4,2,0.5))
hist(xy.data$x, xlim=c(0, 60), ylim=c(0,30), breaks=seq(0,60, 10),
     main=NULL, las=1, col="green", xlab="X values")
par(mar=c(5,2,2,2.5))
hist(xy.data$y, xlim=c(0, 60), ylim=c(0,30),main=NULL,
     yaxt='n', col="red", ylab=NULL, xlab="Y values")
axis(2,labels=F)
dev.off()

par(mfrow=c(1,1), mar=c(5,4,2,2))

plot(y~x, data=xy.data)
?par