#ozone data
#always look first at the raw data file

ozone.df_read.table("ozone.dat")

data.class(ozone.df)
dim(ozone.df)
names(ozone.df)
ozone.df[1:10,]

summary(ozone.df)

boxplot(ozone.df)
par(mfrow=c(5,2))
for(i in 1:10) boxplot(ozone.df[,i])
par(mfrow=c(1,1))
for(i in 1:10) boxplot(ozone.df[,i],xlab=names(ozone.df)[i],ask=T)

plot(ozone.df$doy, ozone.df$temp)
lines(ozone.df$doy,ozone.df$temp)

attach(ozone.df)
plot(doy,temp)
plot(doy,ozone)
pairs(ozone.df)
cor(ozone.df)