data(ozone,package="plyr")
n <- length(ozone)
br <- 231.5+cumsum(seq(from=0,by=1,length=19))
k <- length(br)
aa <- hist(ozone,breaks=br)
ib <- aa$breaks[2:k]-aa$breaks[1:(k-1)]
# sum(aa$dens*ib)
#[1] 1
# $intensities is exactly aa$counts/n/ib
# aa$dens is equal to aa$count/n/ib
# plot histogram
pdf("histo1.pdf",width=7,height=6)
par(mar=c(4.5,4.5,1,1))
plot(0,0,xlim=c(aa$breaks[1],aa$breaks[k]),ylim=c(0,max(aa$counts/n)),type="n",
xlab="Ozon-Daten",ylab="Relative Häufigkeiten",cex.lab=1.3)
for (i in 1:length(ib)){
x1 <- c(aa$mids[i]-ib[i]/2,aa$mids[i]+ib[i]/2,aa$mids[i]+ib[i]/2,
aa$mids[i]-ib[i]/2,aa$mids[i]-ib[i]/2)
y1 <- c(0,0,aa$counts[i]/n,aa$counts[i]/n,0)
polygon(x1,y1)
}
dev.off()
data(ozone,package="plyr")
n <- length(ozone)
br <- 231.5+cumsum(seq(from=0,by=1,length=19))
k <- length(br)
aa <- hist(ozone,breaks=br)
ib <- aa$breaks[2:k]-aa$breaks[1:(k-1)]
# sum(aa$dens*ib)
#[1] 1
# $intensities is exactly aa$counts/n/ib
# aa$dens is equal to aa$count/n/ib
# plot histogram
pdf("histo2.pdf",width=7,height=6)
par(mar=c(4.5,4.5,1,1))
plot(0,0,xlim=c(aa$breaks[1],aa$breaks[k]),ylim=c(0,max(aa$dens)),type="n",
xlab="Ozon-Daten",ylab="Relative Häufigkeiten/Intervallbreite",cex.lab=1.3)
for (i in 1:length(ib)){
x1 <- c(aa$mids[i]-ib[i]/2,aa$mids[i]+ib[i]/2,aa$mids[i]+ib[i]/2,
aa$mids[i]-ib[i]/2,aa$mids[i]-ib[i]/2)
y1 <- c(0,0,aa$dens[i],aa$dens[i],0)
polygon(x1,y1)
}
dev.off()