### aids.S ### Read in the  AIDS  data -- and write it to the
### ~~~~~~ ### 				  datasets.Data  collection

 attach("/u/sfs/S/datasets.Data", pos = 1)

 aids <- read.table("/u/sfs/ueb/datasets/aids.dat", header = T)

 aids $ Time <- aids[,"Jahr"] + .25* aids[,"Quartal"]

 attach("/u/maechler/S/goodies.Data") #-->  To get a  write.data.frame method

 aids.ok <- !is.na(aids[,"kumul.Faelle"])
 write(aids[aids.ok, c("Time","kumul.Faelle")],file="/u/sfs/ueb/datasets/aids2.dat")

 #-- Nun noch das File aids2.dat etwas besser formatiert (mehr "space") .. --
 
 #-- Statt der folgenden 2 Linien koennte man auch
 #--    attach(aids)	verwenden....
 aids2 <- read.table("/u/sfs/ueb/datasets/aids2.dat")
 attach(aids2)

 plot(Time, kumul.Faelle, type='b')
 detach()

## Plot aufgeteilt nach Geschlecht
par(pch=16,mkh=0.05)
attach(aids2)
plot(Time, log(aids[,5]), type='b',ylim=c(0,6),ylab="log.Anzahl")
lines(Time, log(aids[,6]), type='b',lty=2)
legend(1990,1.5,lty=1:2,c("Maenner","Frauen"),cex=1.3)
mtext("AIDS-Daten: log. Quartalszahlen",cex=1.3)
detach()