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

ozone_read.table("ozone.dat")

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

summary(ozone)

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

plot(ozone$temp)
lines(ozone$doy*(330/365),ozone$temp)

attach(ozone)
plot(doy,temp)
plot(doy,ozone$ozone)

pairs(ozone)
cor(ozone)

#=========================================================================
#lunatics data
#first check raw data file

lunatics_read.table("lunatics.dat",header=T,skip=8)

summary(lunatics)
attach(lunatics)
cor(lunatics)

motif()
pairs(lunatics)
plot(DIST,PHOME)
identify(DIST,PHOME)

lunatics2_lunatics[-13,]

detach()
attach(lunatics2)

plot(DIST,PHOME)

lun.lm.1_lm(PHOME~DIST)
summary(lun.lm.1)
summary(lun.lm.1,cor=F)
lines(DIST,fitted(lun.lm.1))

plot(log(DIST),PHOME)
lun.lm.2_lm(PHOME~log(DIST))
lines(log(DIST),fitted(lun.lm.2))